Insights into tectonic reactivation and landscape development processes at the Paraná Basin border using integrated geomorphometric and radiometric analyses

In the Bom Jardim de Goiás region, the Transbrasiliano Lineament shows evidence of a fault reactivation along NE-SW structures (main direction) and secondary directions related to subsequent events. To study the structural framework and reactivation inputs toward tectonic events, distinctive features of the northern Paraná Basin border were analyzed by integrating gamma spectrometry and radar data, both supported by field research. The remote sensing database obtained by the Advanced Land Observing Satellite (ALOS) Phased Array type L-band Synthetic Aperture Radar (PALSAR) sensor allowed detecting lineaments and morphometric derivatives. Besides, the developed morphogenetic process images attest to the importance of structural control on relief evolution while concentrating lineaments as a function of lithological and geomorphological domains. Airborne gamma-ray spectrometry data processing focused on Kd (Anomalous Potassium) and F-Parameter indices. In this context, the fuzzy logic and Index Overlay Method allowed multivalued integration that indicated hydrothermal products and weathering processes correlated with Kd and F-Parameter anomalies. Thus, the joint interpretation of lineaments, geomorphological domains, and Kd/F-Parameter anomalies allowed differentiating primary geological elements and products from surface dynamics, directly or indirectly related to tectonic reactivation processes.


INTRODUCTION
Unique structural and tectonic features have been identified in the northern region of the Paraná Basin, Brazil, by previous geophysical studies and other regional surveys (Soares et al. 1982, Seer 1985, Machado and Souza 1990, Zalán et al. 1991, Marques et al. 1993, Milani et al. 2007, Pinheiro et al. 2019. These features are characterized by the impressive magnetic signature of the Transbrasiliano Lineament (TBL) (Vidotti et al. 2011, Curto et al. 2014, 2015, Pinto and Vidotti 2019, which limits major tectonic domains: Paraná Basin, Arenópolis Magmatic Arc (Brasília Belt), and Paraguay Belt (Cordani et al. 2010). The TBL reactivation from Ordovician to Cretaceous, mainly during the opening of the Atlantic (Alvarenga and Trompette 1993), and associated structures influenced the development of the Paraná Basin. In this context, the increasing availability of open remote sensing and high-resolution geophysical data represent the potential for applying methodologies for processing and integrating data to structural mapping in this area, which still lacks research studies on scales larger than 1:50.000.
The quantitative and qualitative analyses of the dataset collected at different scales (from the synthetic scale of remotely sensed data to field outcrop surveys) in a Geographic Information System (GIS) environment allowed establishing relationships among radiometric signatures, morphogenetic compartments, and field features of the study area. Characterizing geomorphometric features (Pike et al. 2009) provides a better understanding of tectonic events from Paleozoic-Mesozoic to Pleistocene in the geological evolution of the northern Paraná Basin, besides correlating the surface dynamics with landscape development.
The main objective of this paper is to improve the understanding of the geomorphometry and structural settings of the northern Paraná Basin border by integrating airborne gamma-ray spectrometry, geological, and remote sensing data. To this end, we used lineament analysis, F-Parameter, and Anomalous Potassium (K d ) methods. The two integration techniques used were (i) the Index Overlay Method (IOM) to assign different weights to each input map according to their relevance in the analysis, and (ii) the fuzzy modeling approach to build reliable integrated data.

GEOLOGICAL SETTING
The study area -located at the border region of the states of Goiás and Mato Grosso, Brazil -includes part of the northern Paraná Basin and its basement, bounded by the Brasília and Paraguay Neoproterozoic belts. The area covers the basin edge, at the boundary between sedimentary and basement rocks (Fig. 1).

Arenópolis Magmatic Arc
The Arenópolis Magmatic Arc (AMA), extending for several thousand kilometers in the southern Goiás Magmatic Arc, is one of the most important components of the Neoproterozoic Brasília Belt in central Brazil (Pimentel et al. 2000b). In the study area, AMA, to the east of the Serra Negra Fault (Seer 1985), is the main surface expression of the TBL system. Metamorphism reflecting the final closure of the Goiás Ocean and the continental collision took place between 650-630 Ma (Cordani et al. 2003). The AMA consists of orthogneiss units related to a (meta) volcanic-sedimentary sequence, suggesting a tectonic origin similar to modern volcanic arcs in a subduction zone environment (Pimentel and Fuck 1987, Pimentel and Fuck 1992a, Pimentel 2016.
The Bom Jardim de Goiás Group, one of the complete volcano-sedimentary calcium-alkaline sequences of the AMA, is represented by mafic to felsic volcanic rocks close to the Serra Negra Fault domain in the study area. Quartz veins with tourmaline in contact with volcanic rocks indicate hydrothermalism associated with a fault region.

Paraguay Belt and Cuiabá Group
The Neoproterozoic Paraguay Belt, a Pan-African-Brasiliano orogeny along the Amazônia southeastern border, is characterized by low-grade metamorphic sedimentary rocks subjected to polyphase deformation Trompette 1993, Alvarenga et al. 2007). The succession recognized in the study area lies on the west of the TBL, represented mainly by the Cuiabá Group. These rocks were deposited over a marginal basin installed in a passive limit SE of the Amazonian Craton during the Cryogenian and later deformed. The lithostratigraphic division of the Cuiabá Group into Campina de Pedras, Acorizal, and Coxipó formations, including meta-rhythmites, metaconglomerates, sandstones, pelites, and massive subordinate diamictite horizons, has been described by Tokashiki and Saes (2008).
The pelitic metasedimentary rocks of the Cuiabá Group, in contact with AMA lithotypes, present crenulation overprinting an old foliation, shear zones, and fractures. Drag folds associated with fault zones and gouge, as well as cataclasites, evidence recent tectonics. Normal faults placed part of the Cuiabá Group at a relatively higher structural level compared to the Furnas Formation, illustrating a compositional geomorphological control.

Neoproterozoic post-tectonic granites
The study area has four known granite bodies -Serra Negra, Macacos, Bom Jardim, and Serra Verde. The first intrudes the AMA, the second intrudes both the AMA and the Cuiabá Group, and the last ones intrude the Cuiabá Group. Regional granitic intrusions create large calcium-alkaline bodies with K-rich compounds, comprising predominantly isotropic biotite and equigranular granites (Pimentel et al. 1999, Pimentel et al. 2000b).  Pimentel et al. (2000b) suggested two episodes of post-tectonic granite magmatism in Arenópolis -from 590 to 560 Ma and from 508 to 485 Ma -and classified the younger Serra Negra and Iporá granites as Type-A granites.
The location of the granitic intrusions is part of large calcium-alkaline bodies rich in potassium. Pegmatite veins show centimetric K-feldspar crystals. High magmatic concentration of potassium influences the geological context, even in a background naturally rich in potassium.

Transbrasiliano lineament
Magmatic arc rocks are affected by the TBL fault system, with an average trend of N45°E. (Seer 1985). This strike-slip fault system (Schobbenhaus et al. 1975) relates to the final stage of the Pan-African-Brasiliano orogeny (Marini et al. 1984, Cordani et al. 2010. The TBL predominant kinematics is dextral (Pimentel and Fuck 1992b); however, Seer (1985) inferred sinistral kinematics for this zone based on indicators and relative displacement of Macacos Granite. Extensional faults trending N60° and N70°E are superimposed on the older tectonic lineament in the Phanerozoic (Curto 2015).

Paraná Basin
This large intracratonic basin in South America developed entirely on the continental crust, filled with sedimentary and volcanic rocks from Silurian to Cretaceous. It covers more than 1,000,000 km 2 of the Brazilian territory while reaching depths of 7000 m (Zalán et al. 1987, Milani et al. 2007). The stratigraphic record reflects the high basin latitude throughout much of the Palaeozoic, consisting of five major depositional sequences (Silurian, Devonian, Permo-Carboniferous, Triassic, Juro-Cretaceous). The first four depositional sequences are predominantly siliciclastic in nature, while the fifth contains one of the largest continental flood basalts on Earth: the Cretaceous Serra Geral tholeiitic magmatism (Zalán et al. 1991). The study site included parts of the following basin formations -Furnas, Aquidauana, Vila Maria, and Ponta Grossa.

LOCAL STRUCTURAL SETTING, FLUID MOBILIZATION, AND REACTIVATION PROCESSES
According to Cordani et al. (2010), the structural and tectonic regimes affecting the study area involve several tectonic reactivations on the South American Platform, including the Cretaceous (Waldensian) event (Almeida 1967).
In the study area, Seer (1985) defines two transcurrent structures in the TBL system -the Serra Negra and Aldeia faults. The first, predominantly N25E, cuts through Phanerozoic sedimentary rock domains and is characterized by discontinuous gravitational and brittle faults. Using integrated geophysical and fieldwork data, Curto et al. (2014) indicate that N45-60E and N25-30E trending structures are deep crustal components of the TBL, from surface to lower crust levels.
Brittle faults and joints affect the basement N40-60W and N40-70E trends. Such structures overprint pre-existing structures as a result of reactivations that occurred until the Mesozoic Era (Seer 1985). This scenario is confirmed by the Carboniferous and Devonian sequences of the Paraná Basin (Curto 2015). The N75W to E-W azimuthal family is strongly developed on vertical outcrops, and the N10-15W joints have both vertical and sub-vertical dips filled with quartz. The joint filling is also directly related to reactivations in the Aldeia Fault zone, where reactivations are associated with metamorphic rock foliation in the N10W direction.
Hydrothermal fluids may have played a role in the reactivation processes (e.g., Seer 1985, Harcouet-Menou et al. 2009, Curto et al. 2014, Eleraki et al. 2017, Taillefer et al. 2017), causing the enrichment or depletion of eTh/K/eU-bearing minerals. In the NE-SW Meatiq Shear Zone regional system (Egypt), Eleraki et al. (2017) found that most hydrothermal alteration zones are present at the sheared tectonic contact between different lithological units. The authors concluded that increasing Th could be interpreted as evidence of alteration through fractures and faults, while increasing K is related to weathering products. Both eTh and K anomalous values are correlated with tectonic reactivation processes.
Fault zone intersections improve fracture connections and favor the formation of enhanced directions of fluid circulation (Lucianetti et al. 2017, Pischiutta et al. 2013. Also, Harcouet-Menou et al. (2009) show that permeable layers and magmatic intrusions play an important role in controlling fluid drainage. According to Taillefer et al. (2017), the distribution of actual hydrothermal activity in the fault scarp relief of the Eastern Pyrenees (France) suggests that footwall topography is an essential parameter, allowing recharge and hydrothermal fluid circulation path. The authors argue that fault-related topography is crucial to establishing fluid advection and controlling hydrothermal activity intensity.

REACTIVATION EVENTS AND BASIN FORMATION
At least three events are indicative of TBL reactivation in the Phanerozoic -the formation of the Jaibaras Basin (Oliveira 2001) in the Cambrian-Ordovician; the formation of the Água Bonita Graben in the Silurian-Devonian; and the activation of a series of extensional structures during the opening of the Atlantic Ocean in the Late Cretaceous. The first sediments of the Paraná Basin in the Ordovician were deposited along the NE-SW depocenters formed by TBL reactivated structures (Zalán et al. 1991, Curto et al. 2014. Besides, Moura (2017) pointed out reactivations along the TBL between 450 and 429 Ma.
In the Late Cretaceous, a new sedimentation cycle overlapped more than 50% of the Paraná Basin area, forming the Bauru Basin, which exhibits a major elongated shape, trending depocenter in the north-central portion of the Paraná Basin, with NNW-SSE and N-S preferred directions (Batezelli 2006). In the study area, small elongated basins related to the Bauru Basin have settled along the TBL reactivated structures.
The Pantanal Basin -from the Quaternary (Pleistocene) -is a neotectonic reactivation product along the TBL while its elongated elliptical depocenter is controlled by fault activity (Assine 2004).
The Araguaia Formation of the Bananal Basin (Middle Pleistocene) is also controlled by the TBL (Sousa 2017) so that the Araguaia River is embedded in the northern extension of the Serra Negra Fault.
In addition to the works mentioned above, Brito Neves and Fuck (2013) and Pinheiro et al. (2019) recorded other Phanerozoic activities in the Paraná and Parnaíba basins. However, the known facts about reactivations are still based on scarce local data, while several knowledge gaps regarding tectonic evolution need to be filled.

METHODS: DATABASE AND PROCESSING
The remote sensing database was obtained from the elevation data of the ALOS PALSAR sensor. It allowed extracting lineaments and morphometric derivatives, as well as outlining geomorphological units and morphogenetic processes. Additionally, gamma-ray spectrometry data provided K d and F-Parameter indices. For data integration, the IOM assigned different weights to each input map according to their relevance in the analysis. Later, the fuzzy modeling approach was used to build reliable integrated data.

Data sources
On January 24, 2006, ALOS was launched from the Tanegashima Space Center by the Japan Aerospace Exploration Agency (JAXA) carrying three instruments -the Panchromatic Remote-Sensing Instrument for Stereo Mapping (PRISM), the Advanced Visible and Near Infrared Radiometer type 2 (AVNIR-2), and the PALSAR. (JAXA, 2006). PALSAR data are used to produce a digital elevation model (DEM) with a 12.5 m spatial resolution. The pre-processed data was provided by JAXA.
The DEM generated using the ALOS PALSAR ( JAXA, 2006) data was assessed based on geomorphometric and structural lineament analyses. Two methods were used for remote sensing processing. In the first, morphometric parameters were extracted and combined in a color composite RGB image, visually analyzed on-screen to identify geomorphological units and morphogenetic processes. In the second, the lineaments at the basin edge were extracted, followed by the statistical analysis of the identified lineaments. We used the software ENVI (L3HARRIS) 5.2 and ArcMap 10.3 (ESRI) for data processing.
The study area was included in the systematic airborne geophysical survey conducted by the Geological Survey of Brazil (Serviço Geológico do Brasil -CPRM) and the Government of the State of Goiás. The survey acquired magnetic and gamma-ray spectrometry data along the parallel N-S flight lines spaced approximately 500 m, with 100 m ground clearance, and E-W control lines spaced 5,000 m.

Surface morphometric parameters
The altitude, slope, and curvature data derived from the surface topography were extracted and, subsequently, combined in the RGB color composition as follows: altitude in the red channel, slope in the green channel, and curvature in the blue channel. A similar approach was previously used by Castro (2010) and Hermuche (2002) to identify geomorphological units. Our analysis was based on color and textures, identifying similar features by visual inspection and interpretation.

Lineament analysis and directional statistics
Lineaments are alignments of tonal variation on the observed images (Cianfarra and Salvini 2014). The shaded relief technique consists of a 3D terrain surface model that uses simulated sunlight position. This method enhances the texture of linear features so that refined lineaments can be selected based on images that simulate the sunlight position.
The low vertical exaggeration (less than two times) has the purpose of avoiding displacement between lineaments extracted from different lighting angles.
Shaded images were generated according to radar-like synthetic lighting conditions, perpendicular to the N40-60W and N40-70E directions. The methodology developed sought to highlight two main directions, obtaining as many lineaments as possible. According to Seer (1985), these directions mark the principal lineaments known in the study area. This multiple combination of lightning conditions and morphological contrast is useful in shaded images to suppress apparent false rotation (Cianfarra and Salvini 2014). The lineaments were individually extracted using the ArcMap 10.3 vector tool.
We always identified lineaments by negative breaks, that is, in less illuminated features of the shaded image, in order to avoid duplication. Afterward, lineaments provided directional statistics in a wind rose diagram. These data got buffered into 50 m range polygons to increase the extent of the analysis displaying how tectonics acted.

Pre-processing of radiometric data
Spurious data and inconsistent values related to acquisition errors were removed from the database. The correction consisted of comparing the overall mean of each isotope with the mean of each isotope channel. After correction, a constant value corresponding to the difference between the overall and database means was added to each channel (IAEA 2002). Negative values were removed, followed by an interpolation of the dataset using the minimum curvature algorithm to fill the gaps left by these negative values. The software used was Oasis Montaj 9.4 and ArcMap 10.2.

Anomalous potassium and F-Parameter
The K d methodology was developed by Saunders et al. (1994) and later used by Pires (1995) to highlight K d accumulations related to the dynamics of source and depositional areas.
Considering thorium as a resistant accessory mineral, its assembly values can be used to predict standard potassium values for a given lithotype. This correlation (K-Th) is represented by a linear function, obtained from the corrected data whose values describe the ideal potassium equation (Ki) (Eq. 1). The difference between corrected and ideal potassium values (obtained from the linear function) is known as the K d (Eq. 2).
To determine K d , a regression analysis of individual geological or geomorphological units in the study area is recommended.
This methodology adopted to emphasize local anomalies aimed to suppress the regional background that may interfere with the analyses. The heterogeneous geological context of the study area allows adjusting the linear coefficient due to its compositional similarity with average crustal values. The adjustments correspond to overall lithological means used by the IAEA (2003): (1) In which: a = the linear coefficient; b = the angular coefficient. (2) In which: K d = the observed deviation of potassium; Ks = the original value, measured in the survey; Ki = the ideal value.
The process flow below synthesizes the steps followed (Fig. 2): The F-Parameter (Gnojek and Prichystal 1985, Barbuena et al. 2013, Pereira and Ferreira 2018 derives from the contribution of potassium and uranium/thorium ratio. This procedure allows identifying different geological scenarios since relatively higher values correspond to areas with chemically altered rocks, whereas relatively lower values indicate areas with weathering of K-bearing minerals. Equation 3 defines this relationship: (3) Gamma-ray spectrometry and remote sensing data were used for qualitative and visual analyses (including spatial data overlap). The process applies the IOM and fuzzy logic mathematical analysis.

Index Overlay Method
IOM was used to integrate different data layers of interest. Each class has a different weight within the layers, depending on their relative importance to the context (Bonham-Carter 1994) (Fig. 3). For this reason, K d values were first separated K i : ideal potassium; K d : anomalous potassium; a: linear coefficient.  into two layers -geological and morphogenetic processes. The following steps were taken for each layer: • K d values were reclassified quantitatively so that all units were in the same range. This step is required because the basement has a heterogeneous composition; • Jenks optimization method classified features by natural breaks in data values. The Jenks optimization method is also known as the goodness of variance fit (GVF) (ESRI 2019); • After the quantitative classification, the values were ranked according to three qualitative classes (low, medium, and high K d ); Finally, the entire input was integrated by applying the IOM to all qualitative classes exclusively for K d values along the interpreted lineaments to understand the relationship between tectonics and radiometric composition.

Fuzzy logic approach
While standard logic applies only to completely true or completely false concepts, fuzzy logic generalizes standard logic, providing a degree of truth anywhere between 0.0 and 1.0 and allowing a gradual transition from false to true (Alamaniots et al. 2013). This value is chosen experimentally.
• After the K d values were separated according to geological and morphogenetic process units, they were rescaled to 0-1 values. Therefore, an MS Large function (ESRI ® ) was chosen to highlight the highest values. Next, a mosaic was made to gather the obtained data; • A second rescale was required to standardize data in the software model. This step, called "fuzzification", is necessary to obtain the fuzzy membership. In the fuzzy membership processes, the MS Large function was used again; • Subsequently, in the merging process, we used the "AND" operator, considering the highest coincident values in both data (from geological and morphogenetic processes) (Fig. 4).

Geomorphometric and morphogenetic processes
Geomorphometric and morphogenetic processes were analyzed according to economic zoning (MacroZEE 2014) to provide regional geomorphological information about Brazil. In the study area, geomorphological data are classified into Guimarães Plateau and Araguaia Surface units, corresponding to structural and homogeneous tabular dissections.
Geomorphometry, the science of quantitative land-surface analysis, is a modern, analytical-cartographic approach used to represent bare-earth topography by manipulating terrain height digitally (Pike et al. 2009). The RGB (R: altimetry, G: slope, and B: minimum curvature) ternary composition shows distinct zones that combine colors and textures (Fig. 5) ranked according to their importance and each morphometric parameter contribution to geomorphological units.
The classification followed the geomorphological units adopted by the Technical Manual for Geomorphology (IBGE 2009), a general guide commonly used in Brazil. Additionally, fieldwork data supported the ranking process and allowed identifying seven geomorphological units (Fig. 6) Morphogenetic processes represent the landscape movements responsible for sculpting the relief patterns. The seven units can be grouped into three relief sculpting agents classified as slow movements, deposition, and erosion, which controlled the development of geomorphological units (Fig. 6B). Slow movements sculpted transitional relief areas, erosion sculpted those affected by an intense erosive activity like weather and mass movements, and deposition processes sculpted areas that received residual sediments or were formed by resistant materials such as plateaus or residual hills in the basin.
The basin border/edge indicates a strong relationship between lithological and geomorphological units, e.g., the plains and their scarps form a system in which erosive processes act mostly on the basin boundary. Indeed, different relief types are found within the same geological domain, but geology is not the only factor influencing the landscape. Lineaments, as the record of tectonic events, are especially correlated with K d accumulation. Table 1 summarizes the data from Figure 6.

Anomalous potassium trend analysis
Lineaments were interpreted from the shaded relief image. Only negative breaks were demarcated to identify better the structures while different azimuth directions highlighted the data. Directions perpendicular to N40-70E marked lines in the NE direction, which showed the main reactivations in the area. Lighting at 45º, perpendicular to N40-60W (N40-60E), showed lines associated with events after the TBL. Our lineament analysis results are supported by a previous finding of Seer (1985). Finally, azimuths perpendicular to the TBL (N25E) helped identify the oldest features (Fig. 7): The rosette diagram shows that lineaments often follow the NE-SW trend, which is the predominant trend observed in the TBL and during the Brazilian Orogeny. At the same time, NW-SE trending lineaments are widely distributed in the area and also important. From this point, K d was integrated into lineament buffers, considering a 125 m margin (size of a pixel).
K d was calculated for the different lithological and morphogenetic units (Fig. 8). The basin border and basement environments are lithologically complex and heterogeneous, showing very different radiometric backgrounds. Each unit has distinct K d values since rocks with high K and low K background concentrations do not contain the same anomalous range of K d values. Figure 8A shows the K d classes of each geological unit, allowing us to analyze other possible factors that influenced the K d concentration, beyond the primary mineral contribution. Likewise, K d was classified according to the three basic morphogenetic processes (erosion, deposition, and slow movements) to assist in understanding the K mobilization dynamics. The structural control of morphological processes ends up crossing both pieces of information while examining different parameters.
To simplify the K d data analyses, the results were regrouped into three qualitative classes, low, medium, and high K d values (Fig. 8). In both maps, most of the medium to high K values are located along basement lineaments, many near the basin boundary and along the Serra Negra Fault, sites with known erosion and reactivation processes.

Index Overlay Method
Geological and morphogenetic process data were integrated using the IOM, a flexible technique able to indicate priorities in spatial units under the same influence. Greater weight was assigned to morphogenetic process data, whereas terrain was strongly correlated with structural control and K d anomalies, especially for local processes. Furthermore, a morphogenetic process unit can include several lithological units. A new map was produced by intersecting the IOM resulting image and lineament buffers (Fig. 9).
The very high class corresponds to the highest potassium enrichment due to the structures. Positive anomalies are located mostly along N60E, N45E, and between N30W and N45W structures.

Fuzzy logic approach
Fuzzy logic requires rescaling data since anomalous values are different for each geological or morphogenetic unit. Afterward, K d fuzzy membership images were intersected by 125 m lineament buffers, using the MS Large function in both products (Fig. 10).
As expected, the fuzzy logic final product is more restrictive than the IOM result. Despite the lower number of results, fuzzy logic probably shows the most important lineaments with positive anomalies in the area. Nonetheless, the ambiguity related to the origin of K accumulations was solved by using the F-Parameter to indicate areas dominated by either weathering or hydrothermal processes.

F-Parameter analysis
According to Gnojek and Prichystal (1985), low F-Parameter values indicate slightly altered rocks, whereas high values (more than 2) suggest a high concentration of potassium and uranium compared to thorium. Hydrothermalism is expected in such areas (Fig. 11A). In contrast, very low values (less than 1.2 and 1.4), as observed in the Paraná Basin (Fig. 11B), allow assuming that physical weathering and erosive processes are preeminent since these rocks are less chemically altered than others in their surroundings.

DISCUSSION
Technically, both IOM and fuzzy logic were efficient and could be used as complementary techniques. The first method requires fewer processing steps since weights do not need to be assigned in a predetermined way. On the other hand, the second method allows several types of observations, enabling a restrictive analysis. Furthermore, as it uses integrated geological data as a simple validation, fuzzy logic presented more accurate results compared to IOM. At the regional scale, the N25-30E direction is the most important TBL expression, being less expressed in the study area for its inherited crustal tectonic features.
We were able to correlate these lineaments to the structural framework generated by tectonic reactivations of the Aldeia and Serra Negra faults in the N-S and N60E directions, respectively. This structural set shows K, eTh, and eU enrichment, mostly related to geological units of the Cuiabá Group and post-tectonic granite intrusions.
Fuzzy logic integration results have considered gamma-ray spectrometric responses associated with geological and structural data, highlighting potential areas under the influence of tectonism in landscape configuration and isotopic remobilization by lineaments. High K d values related to high F-Parameter values determined for structural features allowed associating these data with those of previous studies by Seer (1985), evidencing tectonic reactivations in the study area and their connection with current relief settings. Hydrothermal activity distribution in the fault scarp relief suggests that footwall topography is an essential parameter for allowing recharge and hydrothermal fluid circulation path (Taillefer et al. 2017).
Locally, relief aspects have a more direct relationship with K d accumulation than pure lithology. Different morphogenetic processes are observed (the reverse is not always true) in the same lithological domain, as relief also results from morphostructural and tectonic control.
Geological and morphogenetic processes bind these tectonic dynamics since stable areas correspond to basin units where the environment is quite enduring. Similarly, deposition areas correspond to basement units where tectonics have operated energetically before basin formation and during reactivation in later periods.
Erosion units correspond to transitional activities, with intense erosive activity in slow movement units and material relocation in deposition unit areas.
In the primary tectonic analysis, morphogenesis also included events that formed the Goiás Magmatic Arc and the TBL, establishing the current basement configuration of the morphogenetic deposition unit basin. Occurrences following the basin formation are related to possible reactivations of old basement structures, subsidence in old lineament gutters, and reactivations of shear zones that eventually "tear" the edge of the Paraná Basin, as seen in Seer (1985).
Gamma-ray spectrometry is an excellent hydrothermalism indicator through F-Parameter analysis, which can reflect the abundance ratios of K and Th/U or U and Th/K. The integrated K d data and F-Parameter extracted from the lineaments     evidenced the hydrothermalism associated with possible tectonic reactivations, highlighted by the combined high K d values and F-Parameter. This result is corroborated by Eleraki et al. (2017), who correlated increased Th and K concentrations with tectonic reactivation processes.
Field data and isotope analysis of the lineaments show high concentrations of K, Th, and U, especially in basement areas. In this context, chemical availability is more diverse due to rock composition. Additionally, field data indicated clay mineral formations and quartz intrusions with associated tourmaline of hydrothermal origin, corroborating the F-Parameter analysis.
The direct relationship between morphogenesis and gamma-ray spectrometry is demonstrated by the dynamics and surface analysis of channels, derived from a shallow method. Morphogenetic processes are often responsible for the isotopic mobilization of elements and their concentration in structural features, thus allowing analyzing the origin of each radio element and inter-relationships. In the study area, this fact is evidenced by the K concentration in depositional units, especially in areas of lineament concentrations: the channels served as channels for K concentration, probably from the reactivation or subsidence of these structures. Since K concentrations are predominantly located in Cuiabá Group domains, the percolation of K-rich fluids in a sedimentary unit supports the reactivation hypothesis.
Also, anomalous values indicate hydrothermalism is concentrated in the basement, with few lineaments belonging to the basin. The highest hydrothermal activity is found in rocks of the Cuiabá and Bom Jardim groups and granitic intrusions of the region. The distribution pattern of the lineaments does not favor any particular lithology, indicating homogeneous hydrothermal activity throughout the basement, a product of reactivations, as seen in Seer (1985).
Moreover, Barbuena et al. (2013) noted that, according to K d data, well-defined lineaments with anomalous values indicate K-rich hydrothermal alteration aligned along faults and lineaments. Similarly, Th and U accumulations are also often related to post-magmatic evidence and hydrothermal processes associated with metamorphism and geomorphological processes (Ulbrich et al. 2009, Ribeiro et al. 2013). The abundance of K and the eU/eTh ratio depend on geomorphological circumstances that favor certain geochemical behavior of elements. These hydrothermal processes can be easily identified (Ostrovskiy 1975, Portnov 1987) by establishing relationships between the F-Parameter, (Efimov 1978 apud Gnojek andPrichystal 1985) and K d (Saunders et al. 1994). Figures 12A and 12B show that most lineaments with high values of fuzzy membership of K d integration also present high F-Parameter.
The results indicate that data produced from the integrated analysis of different datasets generated from various processes allow a more cohesive interpretation of established correlations. Pereira and Ferreira (2018) also discussed the correlation between the K d and F-Parameter and identified areas favorable to hydrothermal mineralization. The analysis confirmed high favorability in areas where the K d and F-Parameter correlation is high.
Furthermore, drag folds, clay-filled fractures, normal faults, and breccia with obliterated quartz and kaolin, as well as the displacement of the Furnas Formation in relation to the Cuiabá Group due to normal faulting, were indicators of fault reactivation (Assumpção and Sacek 2013) found in the field (Fig. 1, point 0BJ-14B).
Tectonics, temporal agents, and fluid surface dynamics govern the landscape evolution by directly affecting the isotopic composition measured in the gamma-ray spectrometry data. The tectonic processes forming the described structures also play a role in the chemical composition of the circulating fluids, through hydrothermal alterations and the weathering generated by erosive morphogenetic relief processes (Pinheiro et al. 2019).
The changing surface dynamics also reflects the formation of lateritic crusts in the study area. X-ray analyses show that lateritic crusts in the Aquidauana Formation contain illite, a secondary mineral originated from the weathering of muscovite-rich material. These processes also remobilize K, which may result in anomalous concentration. Sandstones of the Vila Maria Formation also present surface accumulation of muscovite.

Conclusions
The characterized geomorphometric features allow a deeper understanding of tectonic events while elucidating the consequences on surface dynamics and landscape shaping during the geological evolution of the northern Paraná Basin from the Paleozoic-Mesozoic to the Pleistocene. The correlation between anomalous targets and the dynamic context of landscape modeling indicates how tectonics may have promoted the remobilization of chemical elements. This work demonstrates that many lineaments are located in erosive units, confirming that tectonic reactivation processes can enhance the anomalous concentration of potassium. Although the predominant NE-SW trend of the Serra Negra Fault is broadly observed, NW-SE lineaments are also noticed, suggesting that structures can act together in terrain dynamics and obliterated sediment deposits. The multidisciplinary approach adopted has successfully proven the role of tectonic reactivations in the study area, thus representing a useful tool that can be applied to other regions with similarly enigmatic geotectonic settings.