Quantifying the effects of temperature and concentration on variable-density flow in numerical modeling of groundwater systems: Implications for predictive uncertainty and data collection

Groundwater systems of different densities are often mathematically modeled to understand and predict environmental behavior such as seawater intrusion or submarine groundwater discharge. Additional data collection may be justified if it will cost-effectively aid in reducing the uncertainty of a model's prediction. The collection of salinity, as well as, temperature data could aid in reducing predictive uncertainty in a variable-density model. However, before numerical models can be created, rigorous testing of the modeling code needs to be completed. This research documents the benchmark testing of a new modeling code, SEAWAT Version 4. The benchmark problems include various combinations of density-dependent flow resulting from variations in concentration and temperature. The verified code, SEAWAT, was then applied to two different hydrological analyses to explore the capacity of a variable-density model to guide data collection. ^ The first analysis tested a linear method to guide data collection by quantifying the contribution of different data types and locations toward reducing predictive uncertainty in a nonlinear variable-density flow and transport model. The relative contributions of temperature and concentration measurements, at different locations within a simulated carbonate platform, for predicting movement of the saltwater interface were assessed. Results from the method showed that concentration data had greater worth than temperature data in reducing predictive uncertainty in this case. Results also indicated that a linear method could be used to quantify data worth in a nonlinear model. ^ The second hydrological analysis utilized a model to identify the transient response of the salinity, temperature, age, and amount of submarine groundwater discharge to changes in tidal ocean stage, seasonal temperature variations, and different types of geology. The model was compared to multiple kinds of data to (1) calibrate and verify the model, and (2) explore the potential for the model to be used to guide the collection of data using techniques such as electromagnetic resistivity, thermal imagery, and seepage meters. Results indicated that the model can be used to give insight to submarine groundwater discharge and be used to guide data collection. ^

Coastal groundwater discharge – an additional source of phosphorus for the oligotrophic wetlands of the Everglades

In this manuscript we define a new term we call coastal groundwater discharge (CGD), which is related to submarine groundwater discharge (SGD), but occurs when seawater intrudes inland to force brackish groundwater to discharge to the coastal wetlands. A hydrologic and geochemical investigation of both the groundwater and surface water in the southern Everglades was conducted to investigate the occurrence of CGD associated with seawater intrusion. During the wet season, the surface water chemistry remained fresh. Enhanced chloride, sodium, and calcium concentrations, indicative of brackish groundwater discharge, were observed in the surface water during the dry season. Brackish groundwaters of the southern Everglades contain 1–2.3μM concentrations of total phosphorus (TP). These concentrations exceed the expected values predicted by conservative mixing of local fresh groundwater and intruding seawater, which both have TPμM. The additional source of TP may be from seawater sediments or from the aquifer matrix as a result of water–rock interactions (such as carbonate mineral dissolution and ion exchange reactions) induced by mixing fresh groundwater with intruding seawater. We hypothesize that CGD maybe an additional source of phosphorus (a limiting nutrient) to the coastal wetlands of the southern Everglades.

Sea level rise and South Florida coastal forests

Coastal ecosystems lie at the forefront of sea level rise. We posit that before the onset of actual inundation, sea level rise will influence the species composition of coastal hardwood hammocks and buttonwood (Conocarpus erectus L.) forests of the Everglades National Park based on tolerance to drought and salinity. Precipitation is the major water source in coastal hammocks and is stored in the soil vadose zone, but vadose water will diminish with the rising water table as a consequence of sea level rise, thereby subjecting plants to salt water stress. A model is used to demonstrate that the constraining effect of salinity on transpiration limits the distribution of freshwater-dependent communities. Field data collected in hardwood hammocks and coastal buttonwood forests over 11 years show that halophytes have replaced glycophytes. We establish that sea level rise threatens 21 rare coastal species in Everglades National Park and estimate the relative risk to each species using basic life history and population traits. We review salinity conditions in the estuarine region over 1999–2009 and associate wide variability in the extent of the annual seawater intrusion to variation in freshwater inflows and precipitation. We also examine species composition in coastal and inland hammocks in connection with distance from the coast, depth to water table, and groundwater salinity. Though this study focuses on coastal forests and rare species of South Florida, it has implications for coastal forests threatened by saltwater intrusion across the globe.

The Rise of Native Voices Against Natural Resource Extraction in Latin America [Student's Paper Series]

Latin America, a region rich in both energy resources and native heritage, faces a rising politico-social confrontation that has been growing for over two decades. While resources like oil and gas are exploited to enhance the state’s economic growth, indigenous groups feel threatened because the operations related to this exploitation are infringing on their homelands. Furthermore, they believe that the potential resource wealth found in these environmentally-sensitive regions is provoking an “intrusion” in their ancestral territory of either government agencies or corporations allowed by governmental decree. Indigenous groups, which have achieved greater political voice over the past decade, are protesting against government violations. These protests have reached the media and received international attention, leading the discourse on topics such as civil and human rights violations. When this happens, the State finds itself “between a rock and a hard place”: In a debate between indigenous groups’ rights and economic sustainability.

Online solid phase extraction liquid chromatography tandem mass spectrometry (SPE-LC-MS/MS) method for the determination of sucralose in reclaimed and drinking waters and its photo degradation in natural waters from South Florida

Background Sucralose has gained popularity as a low calorie artificial sweetener worldwide. Due to its high stability and persistence, sucralose has shown widespread occurrence in environmental waters, at concentrations that could reach up to several μg/L. Previous studies have used time consuming sample preparation methods (offline solid phase extraction/derivatization) or methods with rather high detection limits (direct injection) for sucralose analysis. This study described a faster and sensitive analytical method for the determination of sucralose in environmental samples. Results An online SPE-LC–MS/MS method was developed, being capable to quantify sucralose in 12 minutes using only 10 mL of sample, with method detection limits (MDLs) of 4.5 ng/L, 8.5 ng/L and 45 ng/L for deionized water, drinking and reclaimed waters (1:10 diluted with deionized water), respectively. Sucralose was detected in 82% of the reclaimed water samples at concentrations reaching up to 18 μg/L. The monthly average for a period of one year was 9.1 ± 2.9 μg/L. The calculated mass loads per capita of sucralose discharged through WWTP effluents based on the concentrations detected in wastewaters in the U. S. is 5.0 mg/day/person. As expected, the concentrations observed in drinking water were much lower but still relevant reaching as high as 465 ng/L. In order to evaluate the stability of sucralose, photodegradation experiments were performed in natural waters. Significant photodegradation of sucralose was observed only in freshwater at 254 nm. Minimal degradation (<20%) was observed for all matrices under more natural conditions (350 nm or solar simulator). The only photolysis product of sucralose identified by high resolution mass spectrometry was a de-chlorinated molecule at m/z 362.0535, with molecular formula C12H20Cl2O8. Conclusions Online SPE LC-APCI/MS/MS developed in the study was applied to more than 100 environmental samples. Sucralose was frequently detected (>80%) indicating that the conventional treatment process employed in the sewage treatment plants is not efficient for its removal. Detection of sucralose in drinking waters suggests potential contamination of surface and ground waters sources with anthropogenic wastewater streams. Its high resistance to photodegradation, minimal sorption and high solubility indicate that sucralose could be a good tracer of anthropogenic wastewater intrusion into the environment.

Development of effective approaches to the large-scale aerodynamic testing of low-rise buildings

Low-rise buildings are often subjected to high wind loads during hurricanes that lead to severe damage and cause water intrusion. It is therefore important to estimate accurate wind pressures for design purposes to reduce losses. Wind loads on low-rise buildings can differ significantly depending upon the laboratory in which they were measured. The differences are due in large part to inadequate simulations of the low-frequency content of atmospheric velocity fluctuations in the laboratory and to the small scale of the models used for the measurements. A new partial turbulence simulation methodology was developed for simulating the effect of low-frequency flow fluctuations on low-rise buildings more effectively from the point of view of testing accuracy and repeatability than is currently the case. The methodology was validated by comparing aerodynamic pressure data for building models obtained in the open-jet 12-Fan Wall of Wind (WOW) facility against their counterparts in a boundary-layer wind tunnel. Field measurements of pressures on Texas Tech University building and Silsoe building were also used for validation purposes. The tests in partial simulation are freed of integral length scale constraints, meaning that model length scales in such testing are only limited by blockage considerations. Thus the partial simulation methodology can be used to produce aerodynamic data for low-rise buildings by using large-scale models in wind tunnels and WOW-like facilities. This is a major advantage, because large-scale models allow for accurate modeling of architectural details, testing at higher Reynolds number, using greater spatial resolution of the pressure taps in high pressure zones, and assessing the performance of aerodynamic devices to reduce wind effects. The technique eliminates a major cause of discrepancies among measurements conducted in different laboratories and can help to standardize flow simulations for testing residential homes as well as significantly improving testing accuracy and repeatability. Partial turbulence simulation was used in the WOW to determine the performance of discontinuous perforated parapets in mitigating roof pressures. The comparisons of pressures with and without parapets showed significant reductions in pressure coefficients in the zones with high suctions. This demonstrated the potential of such aerodynamic add-on devices to reduce uplift forces.

Spatial and Temporal Changes in Groundwater Salinity in South Florida

A combination of statistical and interpolation methods and Geographic Information System (GIS) spatial analysis was used to evaluate the spatial and temporal changes in groundwater Cl− concentrations in Collier and Lee Counties (southwestern Florida), and Miami-Dade and Broward Counties (southeastern Florida), since 1985. In southwestern Florida, the average Cl− concentrations in the shallow wells (0–43 m) in Collier and Lee Counties increased from 132 mg L−1 in 1985 to 230 mg L−1 in 2000. The average Cl− concentrations in the deep wells (>43 m) of southwestern Florida increased from 392 mg L−1 in 1985 to 447 mg L−1 in 2000. Results also indicated a positive correlation between the mean sea level and Cl− concentrations and between the mean sea level and groundwater levels for the shallow wells. Concentrations in the Biscayne Aquifer (southeastern Florida) were significantly higher than those of southwestern Florida. The average Cl− concentrations increased from 159 mg L−1 in 1985 to 470 mg L−1 in 2010 for the shallow wells (<33 m) and from 1360 mg L−1 in 1985 to 2050 mg L−1 in 2010 for the deep wells (>33 m). In the Biscayne Aquifer, wells showed a positive or negative correlation between mean sea level and Cl− concentrations according to their location with respect to the saltwater intrusion line. Wells located inland behind canal control structures and west of the saltwater intrusion line showed negative correlation values, whereas wells located east of the saltwater intrusion line showed positive values. Overall, the results indicated that since 1985, there was a potential decline in the available freshwater resources estimated at about 12–17% of the available drinking-quality groundwater of the southeastern study area located in the Biscayne Aquifer.

Biogeochemical Effects of Simulated Sea Level Rise on Carbon Loss in an Everglades Mangrove Peat Soil

Saltwater intrusion and inundation can affect soil microbial activity, which regulates the carbon (C) balance in mangroves and helps to determine if these coastal forests can keep pace with sea level rise (SLR). This study evaluated the effects of increased salinity (+15 ppt), increased inundation (−8 cm), and their combination, on soil organic C loss from a mangrove peat soil (Everglades, Florida, USA) under simulated tides. Soil respiration (CO2 flux), methane (CH4) flux, dissolved organic carbon (DOC) production, and porewater nutrient concentrations were quantified. Soil respiration was the major pathway of soil organic C loss (94–98%) and was approximately 90% higher in the control water level than the inundated treatment under elevated salinity. Respiration rate increased with water temperature, but depended upon salinity and tidal range. CH4 flux was minimal, while porewater DOC increased with a concomitant, significant decline in soil bulk density under increased inundation. Porewater ammonium increased (73%) with inundation and soluble reactive phosphorus increased (32%) with salinity. Overall, the decline in soil organic C mineralization from combined saltwater intrusion and prolonged inundation was not significant, but results suggest SLR could increase this soil’s susceptibility to peat collapse and accelerate nutrient and DOC export to adjacent Florida Bay.

Compositional Effects of Sea-Level Rise in a Patchy Landscape: The Dynamics of Tree Islands in the Southeastern Coastal Everglades

The landscape structure of emergent wetlands in undeveloped portions of the southeastern coastal Everglades is comprised of two distinct components: scattered forest fragments, or tree islands, surrounded by a low matrix of marsh or shrub-dominated vegetation. Changes in the matrix, including the inland transgression of salt-tolerant mangroves and the recession of sawgrass marshes, have been attributed to the combination of sea level rise and reductions in fresh water supply. In this study we examined concurrent changes in the composition of the region’s tree islands over a period of almost three decades. No trend in species composition toward more salt-tolerant trees was observed anywhere, but species characteristic of freshwater swamps increased in forests in which fresh water supply was augmented. Tree islands in the coastal Everglades appear to be buffered from some of the short term effects of salt water intrusion, due to their ability to build soils above the surface of the surrounding wetlands, thus maintaining mesophytic conditions. However, the apparent resistance of tree islands to changes associated with sea level rise is likely to be a temporary stage, as continued salt water intrusion will eventually overwhelm the forests’ capacity to maintain fresh water in the rooting zone.

Ecohydrology, evapotranspiration and hydrogeochemistry of carbonate mangrove wetlands

Coastal environments can be highly susceptible to environmental changes caused by anthropogenic pressures and natural events. Both anthropogenic and natural perturbations may directly affect the amount and the quality of water flowing through the ecosystem, both in the surface and subsurface and can subsequently, alter ecological communities and functions. The Florida Everglades and the Sian Ka'an Biosphere Reserve (Mexico) are two large ecosystems with an extensive coastal mangrove ecotone that represent a historically altered and pristine environment, respectively. Rising sea levels, climate change, increased water demand, and salt water intrusion are growing concerns in these regions and underlies the need for a better understanding of the present conditions. The goal of my research was to better understand various ecohydrological, environmental, and hydrogeochemical interactions and relationships in carbonate mangrove wetlands. A combination of aqueous geochemical analyses and visible and near-infrared reflectance data were employed to explore relationships between surface and subsurface water chemistry and spectral biophysical stress in mangroves. Optical satellite imagery and field collected meteorological data were used to estimate surface energy and evapotranspiration and measure variability associated with hurricanes and restoration efforts. Furthermore, major ionic and nutrient concentrations, and stable isotopes of hydrogen and oxygen were used to distinguish water sources and infer coastal groundwater discharge by applying the data to a combined principal component analysis-end member mixing model. Spectral reflectance measured at the field and satellite scales were successfully used to estimate surface and subsurface water chemistry and model chloride concentrations along the southern Everglades. Satellite imagery indicated that mangrove sites that have less tidal flushing and hydrogeomorphic heterogeneity tend to have more variable evapotranspiration and soil heat flux in response to storms and restoration. Lastly, water chemistry and multivariate analyses indicated two distinct fresh groundwater sources that discharge to the phosphorus-limited estuaries and bays of the Sian Ka'an Biopshere Reserve; and that coastal groundwater discharge was an important source for phosphorus. The results of the study give us a better understanding of the ecohydrological and hydrogeological processes in carbonate mangrove environments that can be then be extrapolated to similar coastal ecosystems in the Caribbean.

Uniquely Identifiable Tamper-Evident Device Using Coupling between Subwavelength Gratings

Reliability and sensitive information protection are critical aspects of integrated circuits. A novel technique using near-field evanescent wave coupling from two subwavelength gratings (SWGs), with the input laser source delivered through an optical fiber is presented for tamper evidence of electronic components. The first grating of the pair of coupled subwavelength gratings (CSWGs) was milled directly on the output facet of the silica fiber using focused ion beam (FIB) etching. The second grating was patterned using e-beam lithography and etched into a glass substrate using reactive ion etching (RIE). The slightest intrusion attempt would separate the CSWGs and eliminate near-field coupling between the gratings. Tampering, therefore, would become evident. Computer simulations guided the design for optimal operation of the security solution. The physical dimensions of the SWGs, i.e. period and thickness, were optimized, for a 650 nm illuminating wavelength. The optimal dimensions resulted in a 560 nm grating period for the first grating etched in the silica optical fiber and 420 nm for the second grating etched in borosilicate glass. The incident light beam had a half-width at half-maximum (HWHM) of at least 7 µm to allow discernible higher transmission orders, and a HWHM of 28 µm for minimum noise. The minimum number of individual grating lines present on the optical fiber facet was identified as 15 lines. Grating rotation due to the cylindrical geometry of the fiber resulted in a rotation of the far-field pattern, corresponding to the rotation angle of moiré fringes. With the goal of later adding authentication to tamper evidence, the concept of CSWGs signature was also modeled by introducing random and planned variations in the glass grating. The fiber was placed on a stage supported by a nanomanipulator, which permitted three-dimensional displacement while maintaining the fiber tip normal to the surface of the glass substrate. A 650 nm diode laser was fixed to a translation mount that transmitted the light source through the optical fiber, and the output intensity was measured using a silicon photodiode. The evanescent wave coupling output results for the CSWGs were measured and compared to the simulation results.

Pirometamorfismo em calcários da formação Jandaíra, bacia potiguar, nordeste do Brasil

The present study aims the characterization of thermally affected carbonate rocks from Jandaíra Formation in contact with Paleogene and Neogene basic intrusions in the region of the Pedro Avelino and Jandaíra municipalities (RN), northeastern Brazil. For this study, field, petrographic, x-ray diffraction, electron microprobe, and whole rock litogeochemistry data of carbonates were undertaken. The thermally unaffected limestones are classified like wackstones, grainstones and packstones. They may constitute carbonates grains of benthic foraminifera, echinoderm spines, ostracods, algae, corals, bivalves, gastropods, peloids and intraclasts. The porosities are classified like vug, intraparticle, interparticle, intercrystal and moldic types. The major minerals are calcite, ankerite and dolomite; the detrital are montmorillonite, pyrite, limonite, quartz and microcline. The thermally affected limestones are very coarse to very fine-grained and light to dark gray color. The fossiliferous components totally disappear, and the porosity tends to disappear. With the data obtained, it can be inferred that the carbonate protoliths would be calciferous to dolomitic limestones, both with small amount of clay minerals. Crystalline carbonates from dolomitic protolith have rhombohedral calcite and iron oxides / hydroxides, making the rocks much darker. The carbonates from calciferous protolith have a wide variation of grain size according to the recrystallization degree, increasing toward contact with the basic bodies. In this group, it was identified the minerals lizardite and spinel in weakly to moderately affected samples, and spinel and spurrite in strongly affected rocks, as well as calcite, that occur everywhere. The geological context (shallow level diabase intrusions), the crystallization of the pyrometamorphic minerals spurrite and olivine, and comparison with diagrams from the literature allow estimating temperatures and pressures around 1050-1200 °C and 0.5-1.0 kbar, respectively, for PTOTAL=PCO2. The post-intrusion cooling would have afforded the releasing of metasomatic / hydrothermal fluids, allowing the opening of the metamorphic system, with possible contribution of chemical elements from host units (sandstones, shales) and from basic intrusions. This would induce hydration of previous phases, allowing the formation of serpentine, chlorite and brucite. The results discussed here reveal the strong influence of the heat from basic intrusions within the sedimentary pile. Whereas in the offshore portion of the basin occur sills with up to 1000 m thickness, the understanding of pyrometamorphism might be useful for understanding and measuring the thermally affected rocks.

Pirometamorfismo em calcários da formação Jandaíra, bacia potiguar, nordeste do Brasil

The present study aims the characterization of thermally affected carbonate rocks from Jandaíra Formation in contact with Paleogene and Neogene basic intrusions in the region of the Pedro Avelino and Jandaíra municipalities (RN), northeastern Brazil. For this study, field, petrographic, x-ray diffraction, electron microprobe, and whole rock litogeochemistry data of carbonates were undertaken. The thermally unaffected limestones are classified like wackstones, grainstones and packstones. They may constitute carbonates grains of benthic foraminifera, echinoderm spines, ostracods, algae, corals, bivalves, gastropods, peloids and intraclasts. The porosities are classified like vug, intraparticle, interparticle, intercrystal and moldic types. The major minerals are calcite, ankerite and dolomite; the detrital are montmorillonite, pyrite, limonite, quartz and microcline. The thermally affected limestones are very coarse to very fine-grained and light to dark gray color. The fossiliferous components totally disappear, and the porosity tends to disappear. With the data obtained, it can be inferred that the carbonate protoliths would be calciferous to dolomitic limestones, both with small amount of clay minerals. Crystalline carbonates from dolomitic protolith have rhombohedral calcite and iron oxides / hydroxides, making the rocks much darker. The carbonates from calciferous protolith have a wide variation of grain size according to the recrystallization degree, increasing toward contact with the basic bodies. In this group, it was identified the minerals lizardite and spinel in weakly to moderately affected samples, and spinel and spurrite in strongly affected rocks, as well as calcite, that occur everywhere. The geological context (shallow level diabase intrusions), the crystallization of the pyrometamorphic minerals spurrite and olivine, and comparison with diagrams from the literature allow estimating temperatures and pressures around 1050-1200 °C and 0.5-1.0 kbar, respectively, for PTOTAL=PCO2. The post-intrusion cooling would have afforded the releasing of metasomatic / hydrothermal fluids, allowing the opening of the metamorphic system, with possible contribution of chemical elements from host units (sandstones, shales) and from basic intrusions. This would induce hydration of previous phases, allowing the formation of serpentine, chlorite and brucite. The results discussed here reveal the strong influence of the heat from basic intrusions within the sedimentary pile. Whereas in the offshore portion of the basin occur sills with up to 1000 m thickness, the understanding of pyrometamorphism might be useful for understanding and measuring the thermally affected rocks.

Petrologia do Plutão Bom Jardim de Goiás (PBJG): implicação na evolução neoproterozoica da Província Tocantins

The Bom Jardim de Goiás Pluton (PBJG) is a semi-circular body, located in the central portion of the Tocantins Province, intrusive into orthogneisses and metassupracrustals of the Arenópolis Magmatic Arc. These metasupracrustals present a low to moderate dipping banding or schistosity, have a low to moderate angle of banding / foliation, defined by mica, andalusite and sillimanite and cordierite, which characterize an amphibolite facies metamorphism. This structure is crosscut by the emplacement of the PBJG rocks. The abrupt nature of the contacts and the absence of ductile structures indicate that the intrusion took place in a relatively cold crust. Under petrographic grounds, the pluton consists mainly of monzodiorites, tonalite and granodiorite, following the low to medium-K calk-alkaline alkaline trend. Rocks of the PBJG have hornblende and biotite as the main mafic phases, besides subordinate clinopyroxene, titanite, epidote and opaque. Late dikes of leucogranite contain only mineral biotite as relevant accessory mineral. One U-Pb zircon dating of a monzodiorite yielded an age of 550 ± 12 Ma (MSWD = 1.06). Whole-rock and mineral chemistry suggest that the studied rocks are calc-alkaline, having evolved by fractional crystallization of Ca- and Fe-Mg minerals under high oxygen fugacity. Using the amphibole-plagioclase geothermometer and the Al-in amphibole geobarometer, we calculate temperatures and pressures of, respectively, 692-791 °C e 2.4-5.0 kbar for the intrusion of the PBJG, which is corroborated by previous metamorphic assemblages in the country rocks. The geological, geochemical and geochronological features of PBJG demonstrate their post-tectonic or post-collisional nature, with emplacement into an already uplifted and relatively cool crust at the end of brasiliano orogeny in this portion of the Tocantins Province.

Petrologia do Plutão Bom Jardim de Goiás (PBJG): implicação na evolução neoproterozoica da Província Tocantins

The Bom Jardim de Goiás Pluton (PBJG) is a semi-circular body, located in the central portion of the Tocantins Province, intrusive into orthogneisses and metassupracrustals of the Arenópolis Magmatic Arc. These metasupracrustals present a low to moderate dipping banding or schistosity, have a low to moderate angle of banding / foliation, defined by mica, andalusite and sillimanite and cordierite, which characterize an amphibolite facies metamorphism. This structure is crosscut by the emplacement of the PBJG rocks. The abrupt nature of the contacts and the absence of ductile structures indicate that the intrusion took place in a relatively cold crust. Under petrographic grounds, the pluton consists mainly of monzodiorites, tonalite and granodiorite, following the low to medium-K calk-alkaline alkaline trend. Rocks of the PBJG have hornblende and biotite as the main mafic phases, besides subordinate clinopyroxene, titanite, epidote and opaque. Late dikes of leucogranite contain only mineral biotite as relevant accessory mineral. One U-Pb zircon dating of a monzodiorite yielded an age of 550 ± 12 Ma (MSWD = 1.06). Whole-rock and mineral chemistry suggest that the studied rocks are calc-alkaline, having evolved by fractional crystallization of Ca- and Fe-Mg minerals under high oxygen fugacity. Using the amphibole-plagioclase geothermometer and the Al-in amphibole geobarometer, we calculate temperatures and pressures of, respectively, 692-791 °C e 2.4-5.0 kbar for the intrusion of the PBJG, which is corroborated by previous metamorphic assemblages in the country rocks. The geological, geochemical and geochronological features of PBJG demonstrate their post-tectonic or post-collisional nature, with emplacement into an already uplifted and relatively cool crust at the end of brasiliano orogeny in this portion of the Tocantins Province.