SEARCH · Geology Reports
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.
A U.S. Geological Survey program to monitor and assess channel reconfiguration activity is described. A data base available on the world wide web will enable land-management agencies and other interested parties to evaluate the long-term success of various channel reconfiguration projects. A demonstration project on the Lake Fork of the Gunnison River, Colorado, illustrates the program objectives and approach.
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.
As part of U.S. Geological Survey's (USGS) efforts to identify and assess geothermal energy resources of the US, a three-dimensional (3D) geologic and thermal model has been constructed for the Williston Basin, USA. The geologic model consists of all sedimentary units above the Proterozoic and Archean crystalline rock (called basement herein), with a total sedimentary thickness of up to 5 km near the basin center. Twenty-nine geologic units were mapped from interpreted formation tops from 16,465 wells. A 3D temperature model was constructed to a depth of 7 km by constructing a 3D heat flow model for the sedimentary units, followed by estimating underlying temperature using a one-dimensional (1D) analytic solution for heat flow within the underlying crystalline basement. Using the sedimentary basin model, heat flow was simulated in 3D and was calibrated using three temperature datasets: 1) 24 high-confidence static temperature logs (equilibrium thermal profiles), 2) more than15,000 drill stem test (DST) measurements from >7,000 wells, and 3) more than 45,000 bottomhole temperature (BHT) measurements from >14,000 wells. The DST and BHT datasets provide broad spatial coverage, but are lower confidence, primarily because measurements were made prior to attaining thermal equilibrium. DST and BHT measurements were binned regionally to develop representative thermal profiles that generally agree with these lower quality data (hereafter called pseudowell temperature profiles). Layer properties (primarily thermal conductivity and compaction curves) were set to best estimate values, then the heat flow model was calibrated to fit pseudowell and static temperature logs primarily by adjusting basal heat flow to approximate the overall temperature profile. Minor adjustments to thermal conductivity allowed adjusting changes in slope at lithologic contacts. Resulting maps include 3D temperature and basal (bottom of sedimentary units) heat flow estimates, which are used as input for the temperature model of the basement. The crystalline basement temperature model uses an analytic 1D solution to the heat flow equation that requires estimates of heat flow and temperature at the upper boundary (i.e., the sediment/basement contact), radiogenic heat production within the crystalline basement, and reference thermal conductivity (i.e., uncorrected for temperature). Two regions of high heat flow are identified: 1) in western North Dakota along the North American Central Plains Conductivity Anomaly and 2) in eastern Montana near the Poplar dome. Within the sedimentary column in the center of the basin of the basin, an area of approximately 100,000 km2 is predicted to have moderate- to high-temperature geothermal resources (>90 °C) under the thickest sequences of sediments. Where thick insulation and high heat flow coincide, electric-grade resources can be less than 4 km deep. Assuming a maximum feasible drilling depth of 7 km, temperatures are predicted to be as high as 175 °C. The geologic model may be used to identify strata at sufficient temperatures that may have natural permeability or that may have conditions that favor development of enhanced/engineered geothermal systems resources.
Underserved communities, especially those in coastal areas in Puerto Rico, face significant threats from natural hazards such as hurricanes and rising sea levels. Limited funding hinders the investment in costly mitigation measures, increasing exposure to natural disasters. Providing coastal resources and data products through effective communication mechanisms is fundamental to improving the well-being of these underserved coastal communities. The overall objectives of the pilot effort to engage and connect with underserved coastal communities in Puerto Rico were the following: (1) compile a comprehensive database of the projects and resources relevant to natural hazards in Puerto Rico; (2) foster connections with Puerto Rican interested parties to better understand their priorities regarding coastal hazards and provide them with pertinent U.S. Geological Survey (USGS) resources; and (3) identify knowledge gaps to guide future USGS projects in Puerto Rico. Here we outline our participatory engagement framework and process, along with two specific resources developed with the information collected from this effort. These resources are available in English and Spanish and consist of user-friendly, non-technical information products. Among them are: (1) a website where users can learn about USGS research on landslides, hurricanes, earthquakes, water resources, coastal hazards, tsunamis, and ecosystem hazards and environmental contaminants, and (2) a geonarrative highlighting shoreline changes in Puerto Rico with sections on historical shoreline trends, hurricane impacts, and potential solutions that could help protect communities and mitigate coastal hazards. Continuing participatory engagement in future projects could enhance the accessibility and usability of natural hazards resources within the community.
The United States Geological Survey (USGS) provides scientific information for the Department of Interior and the nation, consistent with its original mission expressed in the Organic Act of 1879 (43 U.S.C. 31): “the classification of the public lands and examination of the geological structure, mineral resources, and products within and outside the national domain.” Legislation such as the Energy Act of 2020 and the 2022 Infrastructure Investment and Jobs Act (43 USC 31l) and recent executive actions (Executive Orders 14154 , 14153, 14241 Secretary’s Orders 3417, 3418, 3422, 3436) underscore the importance of mineral resources and focus USGS activities on mapping and assessing mineral resources, with a particular focus on those presently identified as critical, both in ground and above ground in mine wastes. This article reviews selected activities and accomplishments by the USGS Mineral Resources Program related to critical minerals in 2025. Highlights include a new List of Critical Minerals, a first-ever national mine waste inventory, international minerals partnerships, mineral resource assessment advancements, and national data collection activities and outcomes of the Earth Mapping Resources Initiative (Earth MRI). The selected contributions are not comprehensive but are intended to demonstrate USGS leadership in critical mineral mapping and assessment, the importance of domestic and global partnerships, and the breadth of research activities that are responsive to national needs and priorities.
Techniques for management of drainage basins that use water budgets to balance available water resources with actual or anticipated water use require accurate and precise estimates of basin withdrawals, interbasin transfers of freshwater, unaccounted-for use, water use, consumptive use, inflow and infiltration, basin return flow, and interbasin transfers of wastewater. Frequently, interbasin transfers of freshwater and wastewater are not included in basin water budgets because they occur within public water-delivery and wastewater-collection systems. A new 10-step method was developed to improve estimates of inflow and infiltration and interbasin transfers using readily available statewide data. The accuracy and precision of water-use estimates determined by this method are improved through careful application of coefficients for small users and the use of metered values for large users. The method was developed and tested with data for the Ten Mile River Basin in southeastern Massachusetts. This report uses examples from the basin to illustrate each step of the method.
Geochemical data from sedimentary rocks are the primary source of information regarding Earth's surface evolution through time, including its air and water envelopes and interactions with life and deep Earth processes. The Sedimentary Geochemistry and Paleoenvironments Project (SGP) is a scientific consortium centered around open data and community-driven development of cyberinfrastructure tools and resources for sedimentary geochemistry and Earth history. Here we describe the SGP Phase 2 data release, which focused on incorporating Paleoproterozoic and Mesoproterozoic (2500–1000 million years ago) data and better accommodating carbonate data. This data release was built through the involvement of >200 researchers worldwide in academia, government, and industry, and provides the largest available public data resource for our user community in the academic fields of geochemistry, sedimentology, tectonics, paleontology, Earth history, and paleoclimate, as well as the petroleum and minerals industries. The dataset now encompasses 126,006 samples and 4,132,371 geochemical analyses. In addition to direct entry by SGP Team Members, we have ingested and incorporated datasets from the Geoscience Australia OZCHEM database, the Alberta Geological Survey, and the Deep-Time Marine Sedimentary Element Database (DM-SED) compilation. This paper details sampling in the Phase 2 dataset with respect to age, geography, lithology, and other geological characteristics, documents access via our search website and API, discusses possible issues and/or biases in the dataset that could impact analyses, describes plans for governance and stewardship of data from Indigenous lands, and serves as the citable reference paper for the data release.
Winter snow accumulation and summer snow, firn, and ice melt were measured at South Cascade Glacier, Washington, to determine the winter and net balances for the 1999 balance year. The 1999 winter snow balance, averaged over the glacier, was 3.59 meters, and the net balance was 1.02 meters. Since the winter balance record began in 1959, only three winters have had a higher winter balance. Since the net balance record began in 1953, only 2 years have had a greater positive net balance than 1999. Runoff was measured from the glacier and an adjacent non-glacierized basin. Air temperature, precipitation, and humidity were measured nearby, and ice speed was measured. This report makes these data available to the glaciological and climatological community.
The U.S. Geological Survey (USGS), in cooperation with the Monterey County Water Resources Agency, conducted studies to help evaluate the surface-water and groundwater resources of the Salinas Valley study area, consisting of the entire Salinas River watershed and several smaller, adjacent coastal watersheds draining into Monterey Bay. The Salinas Valley study area is a highly productive agricultural region that depends on the coordinated use of surface water and groundwater to meet demand for irrigation and public water supply. To continue to meet these demands, a better understanding of the historical water balance and the effects of water-resource development on the long-term sustainability of water resources in the Salinas Valley study area is needed.
A digital numerical model of the ground-water-flow system in the Virginia Coastal Plain was constructed as part of the Regional Aquifer System Analysis program of the U.S. Geological Survey and subsequently revised based on reinterpretations of the hydrogeologic framework. The revised model has been incorporated by the Virginia Department of Environmental Quality as a means to evaluate the effects of existing and proposed ground-water withdrawals. The revised model, which is documented here, includes many design features of the original model. Principal modifications are (1) reassignment of some model layers to represent different aquifers, (2) replacement of applied recharge and the uppermost surface-water layer by constant-head cells to represent the unconfined aquifer, and (3) revision of the lateral extents and transmissivity and vertical leakance values assigned to model layers. Ground-water levels simulated by the revised model under steady-state conditions, using withdrawal rates representing the period 1978-80, generally are similar to the water levels simulated by the original model under transient conditions using the same withdrawal rates. Differences probably result largely from changes in transmissivity and vertical leakance values. T he revised model enables the continued evaluation of development of the large and complex aquifer system. Improved information on parts of some aquifers has been incorporated. Only large, regional-scale trends in ground-water water levels, flow directions, and flow rates are simulated, however, and not local trends and (or) short-term changes in water level and flow. Changing hydraulic stresses, knowledge of hydrogeologic conditions, and resource-management needs require periodic revision of the model to maintain its usefulness for the management of ground-water resources in the Virginia Coastal Plain.
A comprehensive scientific research roadmap is essential to bridge knowledge gaps and deepen the understanding of key geological, geochemical, and geophysical aspects of natural hydrogen (H 2 ) as a potential new energy resource. This paper reviews major scientific uncertainties on natural H 2 , suggesting research priorities, as a guide for defining exploration strategies, techniques, and data interpretation. The uncertainties concern all phases of the natural H 2 cycle, from generation (source rocks) through migration (advection and diffusion) and accumulation (reservoir and cap rocks) to the application and interpretation of subsurface and surface geochemical and geophysical exploration techniques. Understanding H 2 sources and generation rates (the amount of H 2 generated by a given volume of rock over time) is crucial for determining whether a geological H 2 system operates as a short-term dynamic system with rapid H 2 production and release, or as a conventional gas system with long-term accumulations, analogous to petroleum reservoirs. Preliminary estimates for serpentinisation, radiolysis, and organic matter degradation suggest that H 2 generation is not inherently fast, especially for non-hydrothermal continental systems (crystalline basement of shields, ophiolites, peridotite massifs, sedimentary basins), and long-term accumulations, like those of fossil natural gas systems, represent the most likely scenario. The mechanisms of H 2 migration through geological formations require application of fundamental principles of fluid-flow physics, distinguishing advection and diffusion, as well as their forms (from gas-phase, bubble flows to aqueous solutions). Additional studies of H 2 accumulation and retention in subsurface reservoirs could improve understanding of mechanisms of H 2 migration by focusing on the rock fluid-bearing properties and the factors affecting H 2 preservation, such as the presence of cap rocks impermeable to H 2 , pressure conditions, residence times, and microbial or abiotic consumption. Advanced techniques, including reservoir modelling, flow simulations, 3D imaging (micro-CT) of H 2 -bearing rocks, and extraction and analysis of gas occluded in rocks, can provide insights into the stability and potential recoverability of H 2 accumulations. The interpretation of surface exploration techniques, including gas geochemistry, geophysics, and remote sensing, long employed in mineral and energy resource exploration, is now being adapted for natural H 2 studies, but challenges remain in the data interpretation. Distinguishing H 2 seepage due to geological degassing from H 2 produced near the surface by modern microbial processes or artificial sources, such as hammering or drilling for soil-gas sampling, drilling into aquifers, and corrosion in boreholes, is an essential step in exploration. The simple detection of H 2 in soils, even in morphological structures like sub-circular depressions or “fairy circles”, cannot be cursorily interpreted as a signal of natural H 2 seepage from a deep source. A holistic geochemical approach, including isotopic analyses of gases associated with H 2 , is recommended to distinguish among the variety of possible H 2 origins. Observations of H 2 in wells should be interrogated to rule out possible artifacts such as corrosion and drill bit metamorphism. The integration of multiple geophysical methods, including seismic, gravimetric, magnetic, and electro-magnetic surveys, is recommended to mitigate interpretation ambiguities regarding the structure of a subsurface H 2 system (source and reservoir rocks, including fluid and gas storage), due to the non-uniqueness of rock-specific physical properties.
Mellette and Todd Counties are located in south-central South Dakota and have a combined area of 2,694 square miles. The White River and its tributaries, which include the Little White River, drain Mellette County and about one-half of Todd County. Tributaries to the Niobrara River, which include the Keya Paha River, drain the other one-half of Todd County. The average discharge of the Little White River is about 56 cubic feet per second as the river enters Todd County and is about 131 cubic feet per second as it discharges to the White River in northern Mellette County. The average discharge of the Keya Paha River just outside Todd County is about 39 cubic feet per second. The average annual runoff for Mellette and Todd Counties ranges from 0.94 to 2.36 inches based on records from nine streamflow-gaging stations in and near the counties. The average annual runoff is 1.62 inches, which compares with the average annual precipitation of about 19 inches. In Todd County, shallow wells completed in the alluvial, Ogallala, Arikaree, and White River aquifers generally can supply water that has low concentrations of dissolved solids, is fresh, and is soft to moderately hard. Ground water from shallow aquifers is limited in Mellette County; therefore, deep wells, often greater than 1,000 feet, are sometimes installed. The Pierre Shale often is used to supply rural domestic and stock wells in Mellette County even though well yields are low and the water has high dissolved solids, is moderately saline, and is very hard. Alluvial aquifers are present in both counties and store an estimated 1.6 million acre-feet of water. The water quality of the alluvial aquifers is dependent on the underlying deposits, and generally the water has low concentrations of dissolved solids, is fresh, and is soft to moderately hard where underlain by the Ogallala and Arikaree Formations; has moderate concentrations of dissolved solids, is slightly saline, and is hard where underlain by the White River Group; and has high concentrations of dissolved solids, is saline, and is very hard where underlain by the Pierre Shale. Also, yields often are lower where the alluvial aquifers are underlain by the Pierre Shale. The Ogallala aquifer is present in only Todd County, and the Arikaree aquifer is present throughout most of Todd County and southwestern and south-central Mellette County. The Ogallala aquifer contains an estimated 17 million acre-feet of water in storage, and the Arikaree aquifer contains an estimated 50 million acre-feet of water in storage. Both aquifers generally are suitable for irrigation, and yields from these aquifers are sometimes greater than 1,000 gallons per minute. Nitrate concentrations in 13 out of 92 water samples collected from the Ogallala aquifer exceeded the Primary Drinking Water Maximum Contaminant Level (MCL) of 10 milligrams per liter. In 11 out of 46 samples collected from the Arikaree aquifer, arsenic concentrations exceeded the MCL of 50 micrograms per liter. The White River aquifer, where present, is usually the shallowest source of ground water in Mellette County. The White River aquifer also is used in northern Todd County where the Ogallala and Arikaree aquifers are not present. The White River aquifer contains an estimated 50 million acre-feet of water in storage. Reported yields from the aquifer range from 1 to 30 gallons per minute, which generally is insufficient to support irrigation in most areas. However, yields are sufficient for livestock-watering and rural-domestic purposes. In both counties, the Pierre Shale is the shallowest bedrock aquifer and is exposed at the land surface throughout most of Mellette County. This aquifer is used primarily in Mellette County. Although the aquifer contains an estimated maximum of 1.5 million acre-feet of water in storage, it is not a viable source of ground water because the aquifer is relatively impermeable, yields are low, and water usually can be obtained from shallower sources, especially in Todd County. Reported yields from the Pierre Shale aquifer range from 1 to 8 gallons per minute. Because few test holes and wells penetrate below the Pierre Shale, little is known about the extent of the deeper bedrock aquifers. All wells completed in the Dakota Sandstone, Inyan Kara, and Minnelusa and Madison aquifers in the counties are used for stock-watering purposes. High concentrations of dissolved solids and hard water are characteristic of the water quality in the bedrock aquifers. Depths to the top of the deeper bedrock aquifers range from 1,270 feet to greater than 2,000 feet below land surface.
Snowmobile use in Yellowstone National Park has increased substantially in the past three decades. In areas of greatest snowmobile use, elevated levels of by-products of gasoline combustion such as ammonium and benzene have been detected in snowpack samples. Annual snowpacks and snow-covered roadways trap deposition from local and regional atmospheric emissions. Snowpack samples representing most of the winter precipitation were collected at about the time of maximum annual snow accumulation at a variety of locations in the park to observe effects of a range of snowmobile traffic levels. Concentrations of organic and inorganic compounds in snow samples from pairs of sites located directly in and off snow-packed roadways used by snowmobiles were compared to concentrations in samples collected at nearby off-road sites. Concentrations of ammonium were 2 to 5 times higher for the in-road snow compared to off-road snow for each pair of sites. Thus, concentrations decreased rapidly with distance from roadways. In addition, concentrations of ammonium, nitrate, sulfate, benzene, and toluene in snow were positively correlated with snowmobile use.
Many individual measurement points are required to characterize groundwater velocity within an aquifer. Groundwater velocity is most commonly measured using a network of >5 cm diameter monitoring wells, which, if not already present at a site, are expensive and labor-intensive to install. Drive-point piezometers—simple, cost-effective wells that can be installed by hand—are a common tool for sampling groundwater in shallow, alluvial aquifers, but most groundwater velocity measurement techniques require equipment that is too large for these narrow (usually <2 cm inside diameter) piezometers. In this technical note, we introduce a low-cost sensor and well packer system (<$90 USD) for performing point dilution tests in narrow piezometers. Field data show that the magnitude of groundwater velocity measured with this technique agrees with velocities computed from natural gradient tracer tests. Additionally, with proper calibration, these sensors can be used to continuously monitor in-well specific conductance, either during inter-well tracer tests with saline tracers or for water quality monitoring. This system is a viable tool for rapid assessment of the magnitude of groundwater velocity in shallow aquifers.
Flood-peak discharges and frequencies are presented for 57 gaged sites in Puerto Rico for recurrence intervals ranging from 2 to 500 years. The log-Pearson Type III distribution, the methodology recommended by the United States Interagency Committee on Water Data, was used to determine the magnitude and frequency of floods at the gaged sites having 10 to 43 years of record. A technique is presented for estimating flood-peak discharges at recurrence intervals ranging from 2 to 500 years for unregulated streams in Puerto Rico with contributing drainage areas ranging from 0.83 to 208 square miles. Loglinear multiple regression analyses, using climatic and basin characteristics and peak-discharge data from the 57 gaged sites, were used to construct regression equations to transfer the magnitude and frequency information from gaged to ungaged sites. The equations have contributing drainage area, depth-to-rock, and mean annual rainfall as the basin and climatic characteristics in estimating flood peak discharges. Examples are given to show a step-by-step procedure in calculating a 100-year flood at a gaged site, an ungaged site, a site near a gaged location, and a site between two gaged sites.
Habitat characterization is part of a multidisciplinary approach to water-quality assessment implemented by the National Water-Quality Assessment Program. Habitat data were collected in the Georgia-Florida Coastal Plain study unit at 24 sites during 1993-95. Data were collected for habitat characteristics at three spatial scales: basin, segment, and reach. Basin data include physiography, land resource provinces, and land use, providing a description of the environmental setting at each site. Segment data include length, gradient, and sinuosity. A Kendall correlation analysis performed on segment characteristics and the log-of-basin area showed a correlation between segment gradient and the log-of-basin area and a correlation between sinuosity and segment length. Reach data consist of field-collected measurements of both instream and riparian habitats. Sand and detritus were the most common channel-bed substrates among the sampled sites. Measurements of channel width, water depth, and bank width and height were used to create cross-sectional profiles of each sampled area. Elevations of selected durations plotted on cross sections illustrated the percentage of time that the banks were inundated at each site. Sites were divided into two groups based on duration of bank inundation (less than or equal to 1 percent and greater than 1 percent). Bank woody vegetation was also sampled and a clustering algorithm known as Two-Way INdicator SPecies ANalysis (TWINSPAN) was used to analyze these data. TWINSPAN divided the sites into two groups based on their vegetation composition. A statistical comparison of the two types of site groups (duration of bank inundation and vegetation) was performed. The significant association between these groups was consistent with the hypothesis that inundation frequency affected riparian vegetation.
This report documents methods and guidelines for model calibration using inverse modeling. The inverse modeling and statistical methods discussed are broadly applicable, but are presented as implemented in the computer programs UCODE, a universal inverse code that can be used with any application model, and MODFLOWP, an inverse code limited to one application model. UCODE and MODFLOWP perform inverse modeling, posed as a parameter-estimation problem, by calculating parameter values that minimize a weighted least-squares objective function using nonlinear regression. Minimization is accomplished using a modified Gauss-Newton method, and prior, or direct, information on estimated parameters can be included in the regression. Inverse modeling in many fields is plagued by problems of instability and nonuniqueness, and obtaining useful results depends on (1) defining a tractable inverse problem using simplifications appropriate to the system under investigation and (2) wise use of statistics generated using calculated sensitivities and the match between observed and simulated values, and associated graphical analyses. Fourteen guidelines presented in this work suggest ways of constructing and calibrating models of complex systems such that the resulting model is as accurate and useful as possible.
Evaporation plays an essential role in the water cycle, influencing local and regional climates while directly impacting water availability in lakes. However, directly measuring evaporation over water bodies remains challenging due to the high costs of installing and maintaining the required in situ instrumentation. Although several remote sensing algorithms have been providing evaporation estimates, the lack of a global validation hinders our understanding of their relative uncertainties and performances across different regions. Here, we analyze the performance of a suite of models that leverage satellite data and meteorological reanalysis to estimate evaporation over lakes worldwide. We compare 3 remote sensing-based models, 1 reanalysis-driven model and 1 ensemble approach, using in situ observations from 27 lakes representing a diverse range of geographic and climatic regions. Our results demonstrate that, overall, the ensemble outperformed any individual model in terms of accuracy, with a RMSE and a bias of 1.3 and 0.3 mm day −1 , respectively. These findings highlight the benefits of using an ensemble approach to estimate open water evaporation with satellite-based models at the global scale, leveraging the unique strengths of each model. For the individual models, differences in the representation of heat storage changes and advection effects led to lower values of RMSE and bias, depending on the location and depth of the lakes. This study sets the path for future improvement of open water evaporation algorithms globally, while remote sensing techniques are proven satisfactory to monitoring of water loss in lakes globally, an essential step toward effective large-scale water resources management.