USGS ScienceSearch

SEARCH · USGS Science

Results for “Science of the Total Envionrment”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,687 records · Page 32Linked to original sources

Effect of mineral deposit data on predictions from the three-part approach to quantitative mineral resource assessment—A study of 16 previous U.S. Geological Survey assessments

The three-part approach to quantitative mineral resource assessment requires information about the properties of undiscovered mineral deposits in an assessment area. These properties are unknown, so the properties of discovered mineral deposits of the same mineral deposit type are used instead. In the three-part approach, these discovered mineral deposits come from around the world, and their properties constitute the pooled data for that mineral deposit type. Alternatively, these discovered mineral deposits could come from the assessment area, and their properties constitute the tract data for that mineral deposit type. Tract data may be more representative of the undiscovered mineral deposits in the assessment area than the pooled data. The goal of this study was to determine whether resource predictions using pooled data are equivalent to resource predictions using tract data. To this end, 16 previous U.S Geological Survey assessments were studied. For each assessment, resources were predicted for one undiscovered mineral deposit in the assessment area. One set of predictions used pooled data, and another used tract data. The two sets of predictions were compared with an equivalence test, using the six assessment statistics that are commonly reported for mineral resource assessments. Practical equivalence is the condition that two corresponding assessment statistics are within a factor of 1.5 of one another. For each of 2 assessments, all 6 assessment statistics were practically equivalent. For both assessments, the assessment statistics from the pooled data, relative to the corresponding assessment statistics from the tract data, ranged from 1.30 times smaller to 1.03 times larger. For each of 14 assessments, 1 or more of the 6 assessment statistics were not practically equivalent. The assessment statistics from the pooled data, relative to the corresponding assessment statistics from the tract data, ranged from 26.6 times smaller to 5.53 times larger. The use of pooled data has been a standard procedure in the three-part approach since at least 1986. The 16 assessments in this study are not a representative sample of those prior assessments that used pooled data. So, it is inappropriate to use the study results to infer whether pooled data affected the resource predictions for those prior assessments.

Scientific Investigations Report

Review and synthesis of the applications of machine learning to coalbed methane recovery

Over the last 30 years, a substantial literature has evolved on the use of machine learning (ML) to assess, predict, and improve the efficiency of coalbed methane (CBM) recovery. In the United States, the production of CBM declined as shale gas production matured, but CBM continues to be an important energy resource in other parts of the world. ML applications that have the potential to improve CBM reservoir management and production forecasts, and to increase exploration and operational efficiency, are still of significant interest. The integration of geostatistical techniques into the CBM ML applications has been largely absent but represents an opportunity for improvement. The literature demonstrates the widespread interest in, and applicability of, ML algorithms applied to CBM problems, and that they continue to result in improvements in predictive performance. However, (1) much of the research is more academic than operational, (2) many results are based on simulations, or small or proprietary datasets, (3) ML performance information can be inconsistent and sometimes entirely omitted, (4) most methodologies are unique to the specific CBM situation and likely not generalizable, (5) no standard data repositories are available to directly compare the performance of competing algorithms, and (6) the spatial component is often omitted. Finally, relatively new ML protocols involving causality analysis and reinforced learning, as well as hybrid workflows combining both supervised and unsupervised learning, are anticipated to dominate the future investigations. Integration of geostatistical and geospatial analysis with ML should enhance performance.

Book chapter

Effect of deicing chemicals on the hydrologic environment in Massachusetts; evaluation of surface-water data collection

The objective of this study is to develop predictive relationships that can be used to describe the impact of highway deicing salts on the hydrologic environment. A method, presented in an earlier report, for estimating yearly mean chloride concentrations from estimated or actual runoff and salt-application data was applied to the period 1971 through1974 and the results compared with yearly mean chloride concentrations computed from records of specific conductance. Estimates were found to be from minus 63 percent to plus 56 percent in error when the error was expressed as the difference between estimated and computed values as a percentage of the computed values. Chloride concentrations used in this study are computed from specific conductance/chloride concentration relationships and records of specific conductance. Preliminary results from graphic analyses and least-squares regression analysis show that relationships between measured values of specific conductance and chloride concentration have correlation coefficients ranging from 0.30 to 0.97. Analysis of streamflow for all major dissolved constituents is recommended with the purpose of attempting to describe the variations in the specific conductance/chloride concentration relationships.

Massachusetts

Prediction of regional broadband strong ground motions using a teleseismic source model of the 18 April 2014 Mw 7.3 Papanoa, Mexico, earthquake

To estimate predicted ground motion from a teleseismic slip model, we use a low‐ and high‐frequency hybrid method to simulate the regional, strong ground motions observed following the 18 April 2014 moment magnitude ( ⁠M w ⁠ ) 7.3 Papanoa, Mexico, earthquake. To generate the regional ground motion at low frequencies (<1 Hz), a teleseismically derived, finite‐fault, kinematic model is used to define the earthquake source, taking into account slip‐model variations identified with a parameter sampling approach that considers possible errors in the fault geometry, the hypocenter depth, and the rupture velocity. A 3D crustal model is used to calculate the low‐frequency ground motions using a finite‐element calculation that includes topography and considers variations in the source model to estimate the uncertainty in the calculations. High frequencies (>1 Hz) are added using a 1D full‐wave propagation code that estimates uncertainties by considering multiple random distributions of slip with different spatial correlation lengths. The synthetic, broadband (0.05–10.0 Hz) ground motions are obtained by combining the low‐ and high‐frequency portions match filtered at 1 Hz. These synthetic ground motions are compared with the regional observations using velocity records, peak ground acceleration, and medians of the orientation‐independent response spectra of the horizontal components (RotD50) calculated at periods of 0.2, 0.3, 0.5, 1.0, 2.0, 3.0, 5.0, 7.5, and 10.0 s. The results indicate that ground motions estimated at these periods using our hybrid approach based primarily on a teleseismically derived source model are comparable to the values observed for the 2014 Papanoa earthquake at regional distances. The approach could be used to estimate strong‐motion spectral levels expected for regions with limited local and regional recordings and could also fill in magnitude or distance gaps in ground‐motion prediction relations utilized in the assessment of seismic hazard.

Papanoa

Active volcanism on the Arabian Shield—Geology, volcanology, and geophysics of northern Harrat Rahat and vicinity, Kingdom of Saudi Arabia

Active volcanic systems pose serious hazards to people and property including inundation and incineration by lava, blanketing by tephra (volcanic ash), exposure to noxious volcanic gases, and damage from shallow earthquakes triggered by ascending molten material (magma). To improve understanding of volcanism and associated seismicity on the western Arabia Plate, the Saudi Geological Survey and the U.S. Geological Survey conducted a multi-year investigation of the northern Harrat Rahat volcanic field adjacent to the city of Al Madīnah al Munawwarah, Kingdom of Saudi Arabia. Project components included creation of a high-resolution digital topographic base; interpretation of eruptive history supported by detailed geologic mapping, paleomagnetism, and abundant high-precision geochronology of volcanic deposits; assessments of eruptive styles and volcanic hazards by physical volcanology; investigation of the origins of magmas in the mantle and of their differentiation in the crust revealed by chemical and isotopic petrology; gravity and magnetotelluric surveys to reveal crustal structures and to search for magma reservoirs; and regional and local seismic tomography and analyses of seismic hazards. Project results are presented in this Professional Paper as chapters written for technical scientific audiences.

Professional Paper

The abandoned mine inventory of the United States—A brief summary

The Infrastructure Investment and Jobs Act of 2021 required the Secretary of the Interior to establish a program to inventory abandoned hard-rock mines in the United States. The Department of the Interior’s Office of Environmental Policy and Compliance asked the U.S. Geological Survey’s Mineral Deposit Database project (USMIN) to use existing data sources to build an inventory of all individual abandoned mine features in the United States. In addition to feature locations, this new database documents the surface land management agency, associated physical and environmental hazards, and any completed mitigation efforts. This information will improve risk assessment and support land management efforts, including hazard mitigation, ecosystem restoration planning, and reclamation. This fact sheet provides a brief overview of USMIN’s collaborative efforts with Federal, State, and Tribal agencies to build a comprehensive and authoritative national inventory of abandoned mine features.

Fact Sheet

A northeast-dipping zone of low frequency earthquakes at the southern edge of Cascadia subduction

Tectonic tremor monitoring occasionally detects events in an anomalous zone in southern Cascadia, 50–100 km west of the main tremor band, near the expected southern edge of the subducting Gorda slab at the Mendocino triple junction. To investigate the geometry and temporal behavior of this tremor, we examine its constituent low-frequency earthquakes (LFEs) by developing 27 stacked LFE waveform templates that we use to detect events from 2018 to 2024. We then relocate LFE sources together with regional seismicity. We find that LFE hypocenters form a northeast-dipping alignment at 22–29 km depth, extending eastward from a zone of micro-earthquakes, ∼15 km south of the southern edge of Gorda slab seismicity. These LFE families exhibit small bursts of activity every few days. Considering the strong world-wide association of tremor and LFEs with high slip-rate, plate-bounding faults, we hypothesize these LFEs may demark the southern edge of Cascadia subduction.

California

Deep syntectonic burial of the Anthracite belt, Eastern Pennsylvania

Fluid inclusion microthermometry and Raman spectroscopy of fluid inclusions in quartz veins from the Pennsylvanian rocks of the Anthracite belt, eastern Pennsylvania support a deep burial model of coalification in favor of focused orogenic hot fluid flow. High-temperature (250 to 255 °C) trapping of CH 4 ± CO 2 saturated aqueous fluids and CH 4 ± CO 2 inclusions indicate fluid trapping at depths of 11.5 to 13.4 km under a cover of Pennsylvanian to Permian(?) syntectonic load. In the folded rocks to the south of the Anthracite belt, CH 4 ± CO 2 fluid inclusions indicate a sediment load that was up to 16.3 km thick. Re -equilibrated aqueous fluid inclusions from veins in Silurian through Devonian rocks give the same range of trapping conditions but a wide range of fluid salinities suggesting that folding, fracturing, and meteoric recharge resulted in the intermixing of fluids from throughout the stratigraphic succession.

Pennsylvania

Geologic map and structure sections along the southern part of the Bartlett Springs Fault Zone and adjacent areas from Cache Creek to Lake Berryessa, northern Coast Ranges, California

Introduction Located in the Coast Ranges of northern California, the Bartlett Springs Fault Zone is the easternmost fault in the San Andreas Fault system in northern California. The fault is a right-lateral, strike-slip structure considered capable of producing an earthquake of moment magnitude 7. The purpose of this mapping is to better characterize the geology and earthquake hazards associated with the southern part of the Bartlett Springs Fault Zone and to help identify any evidence of active uplift on the faults bounding the Coast Ranges. Although the area immediately surrounding the Bartlett Springs Fault Zone is sparsely populated, its southern segment presents a potential seismic hazard to northern California communities as far away as the San Francisco Bay region and Sacramento. There are also nearby water resources, mineral resources, and public lands used for public recreation. The Coast Ranges of northern California are a series of northwest-southeast-oriented mountain ranges and valleys located north of the San Francisco Bay region, between the Pacific Ocean to the west and the Sacramento Valley to the east. The region has rugged terrain, high mountain peaks that reach more than 2,400 meters above sea level, isolated and narrow valley bottoms on which most human settlements are located, and large drainage systems that tend to follow the northwest-southeast-oriented topographic grain. The physiographic character of the region is shaped by its bedrock geology, deformational history, and active faulting. The basement rocks of the northern Coast Ranges consist of the Franciscan Complex and the Great Valley complex, the latter of which consists of two informal units, the Coast Range ophiolite and the Great Valley sequence. The Franciscan Complex and the Great Valley complex are in structural contact along the Coast Range Fault, a regional-scale structure and fundamental crustal boundary. The Franciscan Complex and the Great Valley complex are superposed by active, northwest-southeast-striking strike-slip faults that are associated with seismicity swarms. These active strike-slip faults can produce moderate to large earthquakes that have moment magnitudes of 7–8. In places, these active structures bound large ranges and valleys, suggesting that much of the modern topographic expression is the result of active deformation processes. This report contains new 1:24,000-scale geologic mapping along the southern part of the Bartlett Springs Fault Zone between Clear Lake and Lake Berryessa. The map area spans 738 square kilometers in northern Napa County, southern Lake County, and parts of Yolo and Colusa Counties. The south and east borders of the map are 90 kilometers north of San Francisco and 70 kilometers west of Sacramento, respectively. The map area is within the Knoxville mining district, which has a history of mercury and gold mining dating back to the mid-19th century. The two main towns in the region, Lower Lake and Clearlake, California, are west-northwest of the map area. Approximately 71,000 people live in the cities and rural communities located within a 40-kilometer radius of the center of the map area. The bedrock geology, cross sections, and structural data presented herein are critical for evaluating the long-term evolution of the Bartlett Springs Fault Zone. This work will supplement studies on local seismic hazards, liquefaction potential, landslide hazards, earthquake geology, natural resources, groundwater resources, engineering geology, and tectonic history by providing the background information for site-specific investigations on these subjects.

California

Magnitude, depth and methodological variations of spectral stress drop within the SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest Earthquake Sequence

We present the first ensemble analysis of the 56 different sets of results submitted to the ongoing Community Stress Drop Validation Study using the 2019 Ridgecrest, California, earthquake sequence. Different assumptions and methods result in different estimation of the source contribution to recorded seismograms, and hence to the source parameters (principally corner frequency, f c ⁠ , spectral stress drop, Δσ, and seismic moment, M 0 ⁠ ) obtained from modeling calculated source spectra. For earthquakes smaller than magnitude (M) 2.5 there is negligible correlation between the f c values obtained by different studies, implying that no present method is reliable using available data. For larger magnitude events, correlation between f c measurements of different studies, within even a small M range is always higher than spectral ⁠Δσ , because the f c measurements simply reflect the underlying physical decrease in f c with increasing M. We model the observed trends of submitted f c with both magnitude and depth. Most methods report an increase in spectral Δσ with M, although a magnitude‐invariant spectral Δσ is within the confidence limits. The depth dependence is smaller and depends on whether a study allows attenuation to vary with source depth; a combination of depth‐dependent attenuation correction, and depth‐dependent shear‐wave velocity can compensate for reported depth trends. We model the submitted values to remove differing M and depth variation to investigate the relative interevent variability. We find consistent relative variation between individual events, and also lower relative spectral Δσ in the northwest of the aftershock sequence, and higher on the cross fault and in the region of main fault intersection. This large‐scale comparison implies that absolute spectral Δσ estimates are dependent on the methods used; studies of different regions or using different methods should not be directly compared and improved constraints on path and site corrections are needed to resolve these absolute spectral Δσ differences.

California

Peri-Gondwanan sediment in the Arkoma Basin derived from the north: The detrital zircon record of a uniquely concentrated non-Laurentian source signal in the late Paleozoic

During the assembly of Pangea, peri-Gondwanan terranes collided with the eastern and southern margins of Laurentia and brought with them unique detrital zircon U-Pb signatures. Discriminating between individual peri-Gondwanan terranes in the detrital record is difficult due to their similar geologic histories. However, characterization of this provenance is critical for understanding late Paleozoic sediment routing during development of Pangea. Along southeastern Laurentia, in the Arkoma Basin (present-day Arkansas and eastern Oklahoma, southeastern United States), we identified Middle Pennsylvanian (Desmoinesian) strata that exhibit a concentrated peri-Gondwanan detrital zircon signature (e.g., ca. 800–550 Ma). Although several southern peri-Gondwanan terranes (e.g., Maya, Suwannee) are closer to the Arkoma Basin, geologic data, such as predominantly north-to-south paleocurrents and proximal-to-distal facies relationships in these Desmoinesian strata, support a northern source (e.g., Ganderia, Avalonia, Meguma). Further evidence of a northern source comes from detrital zircon source mapping, which reveals the persistence of this peri-Gondwanan signal in depocenters to the north of the basin after the signal had diminished in the Arkoma Basin. To this end, bottom-up detrital zircon source modeling, source mapping, regional stratigraphy, paleocurrent data, and sandstone petrography allow us to reconstruct the evolution of this Middle Pennsylvanian (Desmoinesian) sediment pathway in the context of intraplate and plate-margin tectonic activity. This reconstruction documents processes affecting Earth’s surface (e.g., tectonics, climate) during the assembly of Pangea and describes in detail part of a dynamic continental-scale drainage system.

Geosphere

Chloride concentrations in groundwater from the western part of the Southern Hills regional aquifer system, Louisiana, 2021–22

Groundwater is heavily used for public supply and industrial uses in the Baton Rouge, Louisiana, area. Lowered water levels resulting from groundwater withdrawals have induced the movement of saltwater towards wells in East Baton Rouge and West Baton Rouge Parishes. Saltwater intrusion has the potential to affect water supply infrastructure, reduce water availability for some uses, and increase treatment costs. To document current conditions, samples were collected from 161 wells screened in 10 aquifers of the Southern Hills regional aquifer system during November 2021 through February 2022. The results were compared with historical data to identify where chloride concentrations are increasing, which could indicate that saltwater intrusion is occurring. Saltwater intrusion, to varying degrees and areal extents, was observed in most of the 10 aquifers. The limited availability of monitoring wells near or within some of the known saltwater plume areas restricts tracking of the movement or delineation of the plumes’ current extents.

Louisiana

Geologic models underpinning the 2024 U.S. Geological Survey assessment of undiscovered oil and gas resources in the Hosston and Travis Peak Formations of the onshore Gulf Coast region, U.S.A.

The Early Cretaceous (Berriasian–Hauterivian) Hosston Formation in Louisiana and eastward is time correlative to the Travis Peak Formation of Texas and southern Arkansas. The formation is a first-order clastic sequence with a regional carbonate transgressive surface as an upper contact. The Hosston and Travis Peak formations contain conventional natural gas and oil accumulations that have been produced for nearly a century. These mature reservoirs contain terrigenous fluvial-deltaic, shore-zone, and paralic deposits across the productive trend; organic-lean mudstone and siltstone lithologies are found outboard of the Lower Cretaceous shelf margin. Producing reservoirs exhibit normal pressure gradients from 0.43 to 0.55 psi/ft (9.7 to 12.4 kpa/m), depths from 4000 to over 20,000 ft (1220 to 6100 m), and temperatures from 150 to 385°F (65 to 196°C). Wells are primarily vertical completions. The number of new field wildcats has been declining since the late 1990s. This paper presents comprehensive geologic models, which include lithofacies maps, structure and isopach maps, burial history models, regional seismic interpretations, and events charts that underpin the recently completed U.S. Geological Survey assessment of undiscovered, technically recoverable hydrocarbons within the Hosston and Travis Peak formations. This study also provides geographic and stratigraphic distributions of Hosston–Travis Peak reservoir properties, including geopressure, reservoir temperature, porosity, permeability, API gravity, and gas-oil ratios. Results indicate estimated undiscovered, technically recoverable mean resources of 28 million barrels of oil and 35.8 trillion cubic ft of gas in conventional and continuous accumulations within the Lower Cretaceous Hosston and Travis Peak formations of the onshore U.S. Gulf Coast region. Quantitative assessment results are detailed in U.S. Geological Survey Fact Sheet 2025–3021 and associated Data Release.

GCAGS Journal

Geochemical monitoring of volcanic fluids in the twenty-first century

This chapter reviews the state-of-the-art of geochemical volcano monitoring techniques. We cover in-situ monitoring technologies that rely on sampling fluids (direct sampling) and on instrumental analysis of the composition of such fluids in real-time (remote sensing of volcanic fluids is covered elsewhere in this book). We first review key concepts and principles in the field, and then review the results of some selected case studies and applications. We cover the large variety of fluid categories emitted by volcanoes, in both the near-(crater fumaroles and lakes, and plumes) and far-(degassing soils, groundwaters) fields. Our aim is to demonstrate the utility of measuring the chemistry of fluids released by volcanoes, and how these can help characterize volcano unrest, and eventually the increased likelihood of eruption. We conclude with a brief discussion of current challenges and knowledge gaps, and on future directions in geochemical monitoring.

Book chapter

The 3D Elevation Program—Supporting the economy of Massachusetts

Introduction Massachusetts extends from the mountains of the Appalachian system in the west of the State to the sandy beaches and rocky shorelines of the Atlantic coast in the east. Inland topographic data support a wide range of important activities, including geologic mapping, transportation planning, forest and wildlife management, quantifying ecological services, water supply protection, commonwealth-wide infrastructure planning, local site planning, and flood-plain management. Nearshore bathymetry can be used to support coastal portions of the Commonwealth by addressing the combined threats of ocean warming, strong storm surge, and rising sea levels. The maintenance and (or) expansion of Massachusetts ports (for instance, Boston, New Bedford) and Cape Cod sediment management depends upon the accurate mapping of bathymetry and the frequent influx of sediment and redeposition. Critical applications that address the broad range of requirements depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features. The 3D Elevation Program (3DEP) is managed by the U.S. Geological Survey (USGS) in partnership with Federal, State, Tribal, U.S. territorial, and local agencies to acquire consistent lidar coverage at quality level 2 or better to meet the many needs of the Nation and Massachusetts. The status of available and in-progress 3DEP baseline lidar data in Massachusetts is shown in figure 1. 3DEP baseline lidar data include quality level 2 or better, 1-meter or better digital elevation models, and lidar point clouds, and must meet the Lidar Base Specification version 1.2 ( https://www.usgs.gov/3dep/lidarspec ) or newer requirements. The National Enhanced Elevation Assessment identified user requirements and conservatively estimated that availability of lidar data would result in at least $1.23 million in new benefits annually to Massachusetts. The top 10 Massachusetts business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Massachusetts

Shallow geologic framework of the Mississippi Sound and the potential for sediment resources

The Mississippi Sound, an estuarine environment located between the mainland and barrier islands bordering the northern Gulf of America (formerly the Gulf of Mexico), serves as a vital ecosystem for the States of Mississippi and Alabama. Spanning approximately 100 kilometers from east to west and covering 1,400 square kilometers, the sound is home to marine industry and ports, and its shallow and brackish waters sustain a diverse array of marine life. Barrier islands along the southern edge of the sound separate the microtidal estuary from the Gulf of America. This protection from gulf wave action mediates current flow within the sound, resulting in predominantly fine-grained sediment deposition along the seafloor. This study, conducted by the U.S. Geological Survey in cooperation with the U.S. Army Corps of Engineers, provides insight on fluvial and tidal processes spanning the past 5,000 years. The report synthesizes existing research to provide a comprehensive overview of the sound geology, from Pleistocene origins to present-day morphology, and utilizes high-resolution single channel seismic profiles and sediment data to identify and map sedimentary deposits and morphologic features at and below the seafloor. Despite its ecological significance, the Mississippi Sound faces environmental challenges, including water-quality issues, habitat degradation, storm-induced erosion, and the ongoing threats of sea-level rise and environmental changes. This study uses the present-day understanding of the sound's geology to inform coastal management decisions, hazard assessment, and potential mineral resources.

Louisiana, Mississippi

Horizontal transport of Picture Gorge Basalt magma through the Monument Dike Swarm determined by magnetic fabric

Flood basalts of the mid-Miocene Columbia River Basalt Group (CRBG) cover 210,000 km 2 of Washington, Oregon, and Idaho. The source of CRBG melt is debated; widely spaced feeder dike swarms can be projected toward hypothetical sources near the Oregon-Idaho border. In this study, we use anisotropy of magnetic susceptibility (AMS) to track magma flow in the Monument dike swarm (MDS), the feeder dikes of the Picture Gorge Basalt (PGB). This small formation of the main-phase CRBG eruptions allows us to explore in detail the localized dynamics of a large igneous province feeder system, with implications for the larger CRBG picture. We measured the magnetic fabric of 205 oriented paleomagnetic specimens subsampled from 97 samples collected from 15 dikes of the MDS. Thermal demagnetization and hysteresis loops show that the magnetic minerals are a mixture of single domain and multidomain sized titanomagnetites. At three dikes, the paleodepth of sampling was determined to be shallow (<350 m). Magma flowing through dikes has been shown—in most cases— to acquire an anisotropic magnetic fabric with an AMS ellipsoid minimum axis perpendicular to the wall and maximum axis aligned in the direction of flow. Of 15 dikes, 12 show horizontal flow directions in the plane of the dike. Only one dike displayed imbricated fabrics, showing westward flow away from the Oregon-Idaho border. We conclude that magma flow in the MDS was sub-horizontal from a distal source.

Oregon

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