Search USGS⌕ Search

Geology topics

Research about Brazil

Source-linked reports with geographic coverage including Brazil.

At least 19 recordsLinked to original sources
USGS2026

From transcriptomic dysregulation to behavioral impairment: Mechanistic effects of weathered Brazilian oil on early life stage fish

The 2019–2020 oil spill along the northeastern coast of Brazil affected coastal reef ecosystems and the organisms that inhabit these environments. To characterize the toxicity of persistent weathered oil residues to early life stage fish, embryonic zebrafish ( Danio rerio ) were exposed to environmentally relevant concentrations of a water-accommodated fraction (WAF) prepared from weathered oil recovered following the spill. Embryos were exposed from 24 to 168 h postfertilization (hpf) to low (2.14 μg ∑PAH 60 L –1 ) and high (17.1 μg ∑PAH 60 L –1 ) concentrations representative of chronic postspill exposure conditions. Transcriptomic profiling revealed early alterations in pathways associated with connective tissue, cardiovascular, and respiratory development at 48 hpf, which corresponded with a significant reduction in heart rate at 72 hpf. By 96 and 168 hpf, molecular responses shifted toward visual system dysfunction, including pathways related to photoreceptor development and phototransduction. Histological evaluation confirmed significant thinning of retinal layers, consistent with impaired retinal development or photoreceptor degeneration. These structural changes translated to functional impairment, as larvae exhibited a 40% reduction in optokinetic response. These integrated analyses identify mechanistic targets of weathered oil toxicity and establish molecular to behavioral relationships that improve understanding of chronic oil spill impacts on early life stage fish.

Pernambuco↗
USGS2026

Assessment of undiscovered conventional oil and gas resources in the Santos Basin, Campos Basin, and Espírito Santo Basin provinces of Brazil, 2024

Using a geology-based assessment methodology, the U.S. Geological Survey estimated undiscovered, technically recoverable mean conventional resources of 10.4 billion barrels of oil and 53.3 trillion cubic feet of gas in the Santos Basin, Campos Basin, and Espírito Santo Basin provinces of Brazil.

Campos Basin, Espírito Santo Basin, Santos Basin↗
USGS2025

Experimental drought suppresses amphibian pathogen yet intensifies transmission and disrupts protective skin microbiome

Shifting precipitation regimes driven by global climate change can alter vertebrate behavior and host-symbiont relationships, potentially compromising host resistance to pathogen invasion. In Brazil's Atlantic Forest, a biodiversity hotspot, prior research identified drought as a key factor disrupting the skin microbiome, contributing to a die-off of pumpkin toadlets due to the invasive waterborne fungal pathogen Batrachochytrium dendrobatidis (Bd). However, observational studies cannot disentangle the direct effect of moisture on Bd growth from increased amphibian activity during wet breeding seasons. Using field enclosures, we experimentally tested the influence of drought conditions on host microhabitat use, Bd disease dynamics, and the composition and predicted Bd-inhibitory function of cutaneous bacterial communities. Each enclosure housed ecologically realistic densities of Brachycephalus pitanga , a micro-endemic pumpkin toadlet. We simulated a short-term drought in half of the enclosures using translucent tarp coverings. To track individual toadlets, we identified their unique markings and collected skin swabs biweekly over 3 months. We then implemented molecular techniques to quantify Bd loads and characterize skin bacterial diversity and composition over time. Our findings indicate that while drought may reduce overall Bd loads on hosts, this effect is partially offset by an increase in the use of water-filled areas of the enclosures and by a disruption of the protective host skin microbiome. This study provides valuable insights into the cascading impacts of climate change on animal behavior, host-symbiont interactions, and disease dynamics.

Atlantic Forest↗
USGS2025

Reproductive habitat mismatch influences chytrid infection dynamics in a tropical amphibian community

Batrachochytrium dendrobatidis ( Bd ) has been decimating amphibian populations globally; previous work indicates that infection risk increases with moisture and thermal mismatch from a host’s optimum. We hypothesized that, in addition to these abiotic influences, mismatch of hosts from their reproductive habitat heightens infection risk via exposure and/or susceptibility mechanisms. We evaluated this “reproductive habitat mismatch hypothesis” by quantifying the interplay of host breeding mode, habitat, and rainfall on Bd infection dynamics using two years of frog survey data—including swab data for 3427 captures representing 44 species—from Brazil’s Atlantic Forest. We modeled infection prevalence, infection intensity, and the number of frogs captured as a function of rainfall, reproductive mode (aquatic or terrestrial), and habitat (aquatic or terrestrial) using hierarchical models. High rainfall was associated with increases in infection prevalence and infection intensity; however, these increases were particularly apparent for species in habitats that were mismatched from the species’ reproductive habitat. Tropical regions experiencing increases in precipitation will likely see higher Bd risk, and our results indicate that such increases in rainfall will be particularly problematic for species that are forced to move from their reproductive habitats by factors such as habitat loss or thermal stress.

Parque Estadual da Serra do Mar–Núcleo Santa Virgí↗
USGS2025

Validation of the U37K' paleotemperature proxy in the South Brazilian Bight from core-top sediments

The paleothermometer based on the alkenone unsaturation index ( U 37 K ′ "> U37K′ ) is often used to reconstruct past sea surface temperatures (SST). In the SW Atlantic Ocean, however, a limited understanding of the seasonal and depth distribution of coccolithophores, which generates the U 37 K ′ "> U37K′ signal preserved in ocean sediments, hinders accurate estimates of past regional SSTs. We analyzed 45 core-top sediment samples from the continental shelf of the South Brazilian Bight (SBB, 23 °S to 28 °S) to assess regional spatial U 37 K ′ "> U37K′ -SST gradients and improve SST estimates. The U 37 K ′ "> U37K′ data were converted to SST using six published paleotemperature equations and compared to modern observational SST data from the World Ocean Atlas (2018) data set. Data indicate that the U 37 K ′ "> U37K′ signal is produced during the austral summer and autumn when regional episodic upwelling events occur over the inner-shelf and the nutricline shoals at the slope. Our analysis of core-top data shows that SST estimates derived from most calibration equations closely align with modern observed temperatures but are skewed toward warmer months associated with upwelling-derived nutrients. Our findings underscore the importance of considering regional and seasonal biases to improve the accuracy of paleotemperature reconstructions. Understanding the factors influencing the SBB U 37 K ′ "> U37K′ signal enables a more meaningful comparison between regional paleoceanographic studies, improving our understanding of past changes in the SW Atlantic Ocean and our ability to predict regional SST response to ongoing and future warming.

South Brazilian Bight↗
USGS2024

Recognition of artificial gases formed during drill-bit metamorphism using advanced mud gas

Drill-bit metamorphism (DBM) is the process of thermal degradation of drilling fluid at the interface of the bit and rock due to the overheating of the bit. The heat generated by the drill when drilling into a rock formation promotes the generation of artificial hydrocarbon and non-hydrocarbon gas, changing the composition of the gas. The objective of this work is to recognize and evaluate artificial gases originating from DBM in wells targeting oil accumulations in pre-salt carbonates in the Santos Basin, Brazil. For the evaluation, chromatographic data from advanced mud gas equipment, drilling parameters, drill type, and lithology were used. The molar concentrations of gases and gas ratios (especially ethene/ethene+ethane and dryness) were analyzed, which identified the occurrence of DBM. DBM is most severe when wells penetrate igneous and carbonate rocks with diamond-impregnated drill bits. The rate of penetration, weight on bit, and rotation per minute were evaluated together with gas data but did not present good correlations to assist in identifying DBM. The depth intervals over which artificial gases formed during DBM are recognized should not be used to infer pay zones or predict the composition and properties of reservoir fluids because the gas composition is completely changed.

Santos Basin↗
USGS2023

Intraspecific trait variability facilitates tree species persistence along riparian forest edges in Southern Amazonia

Tropical forest fragmentation from agricultural expansion alters the microclimatic conditions of the remaining forests, with effects on vegetation structure and function. However, little is known about how the functional trait variability within and among tree species in fragmented landscapes influence and facilitate species’ persistence in these new environmental conditions. Here, we assessed potential changes in tree species’ functional traits in riparian forests within six riparian forests in cropland catchments (Cropland) and four riparian forests in forested catchments (Forest) in southern Amazonia. We sampled 12 common functional traits of 123 species across all sites: 64 common to both croplands and forests, 33 restricted to croplands, and 26 restricted to forests. We found that forest-restricted species had leaves that were thinner, larger, and with higher phosphorus (P) content, compared to cropland-restricted ones. Tree species common to both environments showed higher intraspecific variability in functional traits, with leaf thickness and leaf P concentration varying the most. Species turnover contributed more to differences between forest and cropland environments only for the stem-specific density trait. We conclude that the intraspecific variability of functional traits (leaf thickness, leaf P, and specific leaf area) facilitates species persistence in riparian forests occurring within catchments cleared for agricultural expansion in Amazonia.

Scientific Reports↗
USGS2023

Trends and population estimate of the threatened Buff-breasted Sandpiper Calidris subruficollis wintering in coastal grasslands of southern Brazil

Information about population sizes, trends, and habitat use is key for species conservation and management. The Buff-breasted Sandpiper Calidris subruficollis (BBSA) is a long-distance migratory shorebird that breeds in the Arctic and migrates to south-eastern South America, wintering in the grasslands of southern Brazil, Uruguay, and Argentina. Most studies of Nearctic migratory species occur in the Northern Hemisphere, but monitoring these species at non-breeding areas is crucial for conservation during this phase of the annual cycle. Our first objective was to estimate trends of BBSA at four key areas in southern Brazil during the non-breeding season. We surveyed for BBSA and measured vegetation height in most years from 2008/09 to 2019/20. We used hierarchical distance sampling models in which BBSA abundance and density were modelled as a function of vegetation height and corrected for detectability. Next, we used on-the-ground surveys combined with satellite imagery and habitat classification models to estimate BBSA population size in 2019/20 at two major non-breeding areas. We found that abundance and density were negatively affected by increasing vegetation height. Abundance fluctuated five- to eight-fold over the study period, with peaks in the middle of the study (2014/15). We estimated the BBSA wintering population size as 1,201 (95% credible interval [CI]: 637–1,946) birds in Torotama Island and 2,232 (95% CI: 1,199–3,584) in Lagoa do Peixe National Park during the 2019/20 austral summer. Although no pronounced trend was detected, BBSA abundance fluctuated greatly from year to year. Our results demonstrate that only two of the four key areas hold high densities of BBSA and highlight the positive effect of short grass on BBSA numbers. Short-grass coastal habitats used by BBSA are strongly influenced by livestock grazing and climate, and are expected to shrink in size with future development and climatic changes.

Bird Conservation International↗
USGS2023

Estimating geomagnetically induced currents in southern Brazil using 3-D Earth resistivity model

Geomagnetically induced currents (GICs) result from the interaction of the time variation of ground magnetic field during a geomagnetic disturbance with the Earth's deep electrical resistivity structure. In this study, we simulate induced GICs in a hypothetical representation of a low-latitude power transmission network located mainly over the large Paleozoic Paraná basin (PB) in southern Brazil. Two intense geomagnetic storms in June and December 2015 are chosen and geoelectric fields are calculated by convolving a three-dimensional (3-D) Earth resistivity model with recorded geomagnetic variations. The dB / dt proxy often used to characterize GIC activity fails during the June storm mainly due to the relationship of the instantaneous geoelectric field to previous magnetic field values. Precise resistances of network components are unknown, so assumptions are made for calculating GIC flows from the derived geoelectric field. The largest GICs are modeled in regions of low conductance in the 3-D resistivity model, concentrated in an isolated substation at the northern edge of the network and in a cluster of substations in its central part where the east-west (E-W) oriented transmission lines coincide with the orientation of the instantaneous geoelectric field. The maximum magnitude of the modeled GIC was obtained during the main phase of the June storm, modeled at a northern substation, while the lowest magnitudes were found over prominent crustal anomalies along the PB axis and bordering the continental margin. The simulation results will be used to prospect the optimal substations for installation of GIC monitoring equipment.

Space Weather↗
USGS2023

Using simulated GEDI waveforms to evaluate the effects of beam sensitivity and terrain slope on GEDI L2A relative height metrics over the Brazilian Amazon Forest

The vertical structure of forests provides important parameters for estimating aboveground biomass (AGB) and it can be measured by LiDAR sensors. The Global Ecosystem Dynamics Investigation (GEDI) full-waveform LiDAR sensor collects data systematically over the Earth’s surface from the International Space Station. Since GEDI became operational, it has collected billions of ~25 m diameter footprints. This massive dataset has been used to create higher level gridded and non-grided products. However, GEDI’s ~25 m footprints can be subject to errors associated with effects of geolocation, terrain slope, and beam sensitivity, among others, which are likely transferred to the downstream products. This study aims to (1) evaluate the effect of beam sensitivity and terrain slope on the accuracy of relative heights (RH) of GEDI product L2A version 2 through comparison with discrete-return airborne LiDAR data collected over transects in the Brazilian Amazon Forest biome, (2) assess GEDI’s geolocation uncertainty and investigate its combined effects with beam sensitivity and terrain slope, and (3) re-evaluate beam sensitivity and terrain slope effects on the GEDI L2A RHs using footprints that were geolocation-adjusted through a simple novel approach. The analysis separates GEDI footprints by acquisition time, i.e., daytime, nighttime, and combined (all-data). The discrete-return airborne LiDAR point clouds are used to derive terrain slope within the GEDI footprints, and to simulate GEDI waveforms and derive RHs comparable to the GEDI L2A product. Results indicate that terrain slope only causes significant effects on GEDI data collected during daytime because solar radiation affects waveform signal-to-noise ratio. Beam sensitivity causes significant effects on nighttime and all-data GEDI L2A RHs. If geolocation uncertainty is considered, the effects of beam sensitivity and terrain slope have just minor changes. Geolocation-adjusted data continue showing significant effects on nighttime and all-data RH differences caused by beam sensitivity but produce unbiased results. This study improves the understanding of how beam sensitivity, terrain slope and their combined effect with geolocation uncertainty may affect the GEDI L2A RH collected over the Brazilian Amazon forest during daytime and nighttime.

Brazilian Amazon Forest↗
USGS2022

Floodplain ecology: A novel wetland community of the Amazon

An expedition to the upper estuarine reaches of the Amazon River reveals intriguing overlap of tropical mangrove wetlands with riverine floodplain forests. This newly discovered type of forested wetland assemblage may provide a uniquely process-rich carbon hotspot.

Amazon River↗
USGS2022

Geochemical characterization of natural gases in the pre-salt section of the Santos Basin (Brazil) focused on hydrocarbons and volatile organic sulfur compounds

The objective of this work is to characterize the geochemistry of a suite of natural gas samples from five fields in order to improve the understanding of the lacustrine petroleum system of the pre-salt section from the Santos Basin (Brazil). Additionally, the distribution of volatile organic sulfur compounds (VOSC) in petroleum reservoirs was examined to investigate possible applications to petroleum systems assessments. The hydrocarbon gases were generated by thermogenic processes associated with the oil window. The 13 C-enriched values for C 1 (>−40‰) were interpreted as an organic source signature rather than an indication of thermal maturity, except for the oil occurrence (Field B), where a different fluid charge mainly composed of methane and CO 2 from a minor external kitchen area in the Santos Basin was identified. The molecular composition and the carbon and hydrogen isotopic data of the hydrocarbon gases, when combined with the VOSC molecular compositional data, allow the identification of four gas families associated with different kitchens and/or migrations pathways. The total VOSC concentrations range from 0.7 to 23.9 ppm by volume (ppmV). The organic sulfides are mainly composed of carbonyl sulfide (COS) and diethyl sulfide. The main thiol compound is ethanethiol. The cyclic VOSC are primarily composed of thiophene, with a negligible amount of branched thiophenes (<0.1 ppmV). H2S showed a strong positive Pearson's correlation with COS and methanethiol (MeSH) concentrations (r = 0.943 and 0.807, respectively). This suggests that COS and MeSH formation was linked to H 2 S generation and/or post-catagenetic interactions between hydrocarbons and H 2 S, mainly related to thermochemical sulfate reduction (TSR). In contrast, the distribution of higher molecular weight VOSC seems to be controlled by source rock facies, rather than H 2 S concentration. Principal component analysis of the VOSC compositional data identified some subgroups within the gas families mainly associated with TSR. The results presented in this work reveal that VOSC can be an important auxiliary tool in petroleum system studies.

Atlantic Ocean, Santos Basin↗
USGS2022

Early Neoproterozoic gold deposits of the Alto Guaporé province, southwestern Amazon craton, western Brazil

The Alto Guaporé gold province, southwestern Amazon craton, contains gold deposits that have been mined since the beginning of the 18th century and these deposits, together, have modern-day, pre-mining gold resources of at least 1.8 Moz. The ore is associated with quartz vein systems along the southeastern part of the Aguapei belt, a ~35-km-wide and ~500-km-long, NNW-trending shear zone formed due to reactivation of a terrane-bounding suture. The Aguapei belt evolved by ca. 1150 to 1100 Ma rifting and deposition of siliciclastic sediments in an aulacogen basin, followed by deformation and low-grade metamorphism of the sedimentary sequences during 1100 to 900 Ma terrane collision along the craton margin. The deformation was characterized by a compressional regime until ca. 950 Ma and transition to a transpressional setting during the final 50 m.y. The gold deposits are hosted in a variety of structures that are second-order to the main Aguapei shear zone. The Ernesto and Pau-a-Pique deposits are located ~40 km apart and at jogs along the Aguapei belt. They are marginal to pre-ore igneous rocks, with Ernesto hosted in the basal part of the metasedimentary Fortuna Formation that overlies tonalite and Pau-a-Pique at the contact between metasedimentary rocks and diorite. Three deformational phases comprise the compressional (D 1 to D 2 ) to transpressional (D 3 ) tectonic events. In the Pau-a-Pique deposit and the deeper level of the Ernesto deposit, the ore-bearing veins are bedding parallel and follow D 2 strike-slip and reverse fault zones, respectively. However, the veins formed during D 3 reactivation of the older structures by an array of oblique accommodation faults. In contrast, ores at shallower levels of Ernesto, both in discordant and bedding-parallel veins, are hosted within a ~20-m-thick rigid metaconglomerate with associated dilation due to the structural complexity as sedimentary rocks of the Aguapei Group were folded around the dome-shaped roof of the pre-ore tonalite. The ores in both deposits, as well as in many other deposits of the province, are characterized by disseminated and vein-hosted pyrite. Gold occurs mainly as inclusions in the pyrite, with other hydrothermal phases comprising muscovite, Fe-Ti oxides, and minor apatite, chalcopyrite, and galena. Fluid inclusion data, coupled with stable isotope geochemistry and geothermometry, indicate that gold precipitated from a low-salinity, CO 2 -rich fluid at ~300°C and ~2.5 kbar. The source for the fluid and gold was the interbedded pelites during devolatilization of the Aguapei Group sequence. The aqueous-carbonic fluid inclusions and the narrow range of δ 18 O values of quartz (12 ± 1 ‰ ) from many auriferous veins from the central part of the province represent a regional ore-forming fluid. The broad range of δ D for hydrous minerals (–116 to –55 ‰ ) reflects influx of small amounts of meteoric water into the steeply dipping shear zones during postgold exhumation. The 40 Ar/ 39 Ar geochronology from hydrothermal muscovite indicates a widespread hydrothermal event along the belt between 928 and 920 Ma. Collectively, the geological, geochronological, and geochemical data suggest that metamorphic fluids migrated laterally into and then upward along the Aguapei belt and deposited gold in lower-order structures where strain gradients existed between lithounits. The province has many characteristics of large orogenic gold provinces worldwide and represents a highly prospective and underexplored target region for early Neoproterozoic gold, a time period that generally is not well endowed in gold ores.

Alto Guaporé gold province, Amazon craton↗
USGS2021

Mapping of suspended sediment transport using acoustic methods in a Pantanal tributary

Generally, fluvial systems are used for different objectives including energy production, water supply, recreation, and navigation. Thus, many impacts must be considered with their use. An understanding of sediment dynamics in fluvial systems is often of value for a variety of objectives, given that erosion and depositional processes can change the fluvial system morphology and can substantially alter the fluvial environment. In this sense, sediment monitoring is important because it helps to explain and quantify sediment dynamics in the environment. Hence, this study presents an innovative sediment monitoring technique: the use of the acoustic Doppler current profiler, commonly used to obtain discharge measurements, to obtain suspended sediment concentration (SSC). This paper aims to describe the application of additional corrections to the ADP-M9 signal to obtain SSC from measurement campaigns that used the ADP only for discharge measurements. The analyses were based on traditional sediment sampling methods and discharge measurements, with the ADP-M9, from 7 field campaigns at the Taquari River, a major tributary from the Alto Paraguay Basin, in the Pantanal Biome, known as the largest freshwater wetland system in the world. The correlation was assessed considering the following: (a) the equipment frequency operation mode (Smart Pulse or Fixed Frequency) and (b) by checking the influence of the sediment attenuation coefficient. Furthermore, extrapolation was conducted in filtered and unmeasured areas of the ADP to map the suspended sediment concentration over the entire cross section. Results indicate that ADP correlations can be an effective tool for estimating SSC in the Taquari River when samples cannot be collected. Correlations could be applied to past and future ADP measurements made at the location where the correlation was created, as long as similar environmental conditions are present as when the correlation was developed. The described technique can expand the amount of sediment data available at a monitoring site even with reduced traditional sampling and by leveraging instruments used for other monitoring purposes.

Taquari River Basin↗
USGS2021

Land use change influences ecosystem function in headwater streams of the Lowland Amazon Basin

Intensive agriculture alters headwater streams, but our understanding of its effects is limited in tropical regions where rates of agricultural expansion and intensification are currently greatest. Riparian forest protections are an important conservation tool, but whether they provide adequate protection of stream function in these areas of rapid tropical agricultural development has not been well studied. To address these gaps, we conducted a study in the lowland Brazilian Amazon, an area undergoing rapid cropland expansion, to assess the effects of land use change on organic matter dynamics (OM), ecosystem metabolism, and nutrient concentrations and uptake (nitrate and phosphate) in 11 first order streams draining forested ( n = 4) or cropland ( n = 7) watersheds with intact riparian forests. We found that streams had similar terrestrial litter inputs, but OM biomass was lower in cropland streams. Gross primary productivity was low and not different between land uses, but ecosystem respiration and net ecosystem production showed greater seasonality in cropland streams. Although we found no difference in stream concentrations of dissolved nutrients, phosphate uptake exceeded nitrate uptake in all streams and was higher in cropland than forested streams. This indicates that streams will be more retentive of phosphorus than nitrogen and that if fertilizer nitrogen reaches streams, it will be exported in stream networks. Overall, we found relatively subtle differences in stream function, indicating that riparian buffers have thus far provided protection against major functional shifts seen in other systems. However, the changes we did observe were linked to watershed scale shifts in hydrology, water temperature, and light availability resulting from watershed deforestation. This has implications for the conservation of tens of thousands of stream kilometers across the expanding Amazon cropland region.

Tanguro Ranch↗
USGS2021

Riparian forests buffer the negative effects of cropland on macroinvertebrate diversity in lowland Amazonian streams

Riparian forests regulate stream ecosystems and biodiversity. Therefore, changes to riparian structure may threaten stream ecosystem function by triggering taxonomic and functional changes to aquatic communities. Because macroinvertebrate assemblages are sensitive to environmental changes, they can be effective indicators of stream integrity in disturbed landscapes. To assess the role of riparian forests in maintaining tropical stream communities in areas experiencing large-scale watershed disturbance, we quantified the taxonomic and functional response of stream macroinvertebrate communities to forest clearing in the southeastern Amazon’s agricultural frontier, a region experiencing widespread deforestation. Our results show that watershed deforestation can lead to significant changes in macroinvertebrate richness and community composition. We found a predominance of shredders in forested watersheds; scrapers in cropland watersheds with riparian forests; and collector-filterers in cropland watersheds without riparian forest buffers. Taxonomic composition was controlled by available organic material in forested watersheds and by periphyton in cropland sites regardless of whether they had a riparian buffer. Our results show that the clearing of riparian forests alters food sources supporting aquatic food webs, leading to ecosystem-level shifts through changes in light and temperature dynamics that affect aquatic communities in areas with intense land-use change such as the southeastern Amazon.

Tanguro Ranch↗
USGS2021

Remote sensing-based actual evapotranspiration assessment in a data-scarce area of Brazil: A case study of the Urucuia Aquifer System

The large groundwater reserves of the Urucuia Aquifer System (UAS) enabled agricultural development and economic growth in the western Bahia State, in northeastern Brazil. Over the last several years, concern has grown around the aquifer’s diminishing water levels, and water balance (WB) studies are in demand. Considering the lack of measured actual evapotranspiration (ET a ), a major component of the water cycle, this work uses the Operational Simplified Surface Energy Balance (SSEBop) model to estimate ET a , and compares it to basin-scale estimates from the Soil Moisture Accounting Procedure (SMAP) monthly model and from an annual WB closure method, based on gridded meteorological data and the Gravity Recovery and Climate Experiment (GRACE) product. Additionally, a comparative assessment of different versions of the SSEBop parameterization was performed. Moderate Resolution Imaging Spectroradiometer (MODIS) imagery was used to implement eight different versions of the SSEBop algorithm over the UAS between 2000 and 2013. SSEBop and SMAP ET a yielded similar seasonal patterns, with correlation coefficient (r) up to 0.65, mean difference (MD) of 0.8 mm/month and mean absolute difference (MAD) of 18.5 mm/month. Comparison of SSEBop annual ET a estimates to annual SMAP and WB closure estimates yielded low MD (12.1 and −7.3 mm/year, respectively) and MAD (82.5 and 82.8 mm/year, respectively), but also low r values (0.00 and 0.37, respectively). The comparison of the different SSEBop versions indicated the need to incorporate a calibration step of the aerodynamic heat resistance (r ah ) parameter. SSEBop results were also used for land cover and drought monitoring. Analysis indicates that agriculture, associated with an increasing trend of atmospheric evaporative demand, is responsible for the decrease in groundwater levels and streamflow in the studied time period.

Urucuia Aquifer System↗