SEARCH · Geology Reports

Results for “Algorithms”

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

SourceResearch connectionFrom yearThrough year

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

At least 721 records · Page 40Linked to original sources

USGS2016

Spectrally based mapping of riverbed composition

Remote sensing methods provide an efficient means of characterizing fluvial systems. This study evaluated the potential to map riverbed composition based on in situ and/or remote measurements of reflectance. Field spectra and substrate photos from the Snake River, Wyoming, USA, were used to identify different sediment facies and degrees of algal development and to quantify their optical characteristics. We hypothesized that accounting for the effects of depth and water column attenuation to isolate the reflectance of the streambed would enhance distinctions among bottom types and facilitate substrate classification. A bottom reflectance retrieval algorithm adapted from coastal research yielded realistic spectra for the 450 to 700 nm range; but bottom reflectance-based substrate classifications, generated using a random forest technique, were no more accurate than classifications derived from above-water field spectra. Additional hypothesis testing indicated that a combination of reflectance magnitude (brightness) and indices of spectral shape provided the most accurate riverbed classifications. Convolving field spectra to the response functions of a multispectral satellite and a hyperspectral imaging system did not reduce classification accuracies, implying that high spectral resolution was not essential. Supervised classifications of algal density produced from hyperspectral data and an inferred bottom reflectance image were not highly accurate, but unsupervised classification of the bottom reflectance image revealed distinct spectrally based clusters, suggesting that such an image could provide additional river information. We attribute the failure of bottom reflectance retrieval to yield more reliable substrate maps to a latent correlation between depth and bottom type. Accounting for the effects of depth might have eliminated a key distinction among substrates and thus reduced discriminatory power. Although further, more systematic study across a broader range of fluvial environments is needed to substantiate our initial results, this case study suggests that bed composition in shallow, clear-flowing rivers potentially could be mapped remotely.

Geomorphology↗

USGS2025

Favorability mapping for hydrothermal power resource assessments of the Great Basin, USA

The U.S. Geological Survey (USGS) is updating the 2008 assessment of conventional hydrothermal resources for the Great Basin in the western United States. As part of this work, the workflow for hydrothermal resource favorability maps is being modified to integrate modern data-driven machine learning (ML) methods. Improvements include: [1] using new and refined evidence layers (features); [2] using an order of magnitude more training sites (labeled examples); [3] utilizing simple but non-linear supervised ML algorithms; [4] representing positive training sites (wells with measured heat flow) with their ordinal value proportional to the magnitude of convective upflow (i.e., low, high, or very high convective signals instead of past strategies using positive-negative labels); [5] supplementing training sites with additional sites with low convective signals to represent diverse under-sampled areas where hydrothermal systems are unlikely to exist; [6] comparing with competing approaches; and [7] utilizing Monte Carlo cross-validation to estimate and evaluate prediction uncertainty. For the new favorability map, over half of the power-producing systems (i.e., 15 of 28) are predicted in the 99th percentile of most favorable locations (i.e., the highest 1 % of favorability, corresponding to 1 % of the map area), exceeding the performance of past models that have explicitly used power plants as training sites. Previous favorability maps predicted approximately half of the power-producing hydrothermal systems above the 80th percentile (i.e., 20 % of the map area). For the new favorability map, 93 % of power-producing systems (i.e., 26 of 28) are above the 80th percentile. The power-producing systems for which the new model does not perform well are either comparatively small, low-temperature systems or systems also not predicted well by prior modeling approaches, suggesting that these few systems are unusual when compared with most power-producing systems. Focusing research on these known, seemingly different systems may yield new insights and subsequent discovery of new prospects.

California, Idaho, Nevada, Oregon, Utah↗

USGS2024

Hyperspectral (VNIR-SWIR) analysis of roll front uranium host rocks and industrial minerals from Karnes and Live Oak Counties, Texas Coastal Plain

VNIR-SWIR (400–2500 nm) reflectance measurements were made on the surfaces of various cores, cuttings and sample splits of sedimentary rocks from the Tertiary Jackson Group, and Catahoula, Oakville and Goliad Formations. These rocks vary in composition and texture from mudstone and claystone to sandstone and are known host rocks for roll front uranium occurrences in Karnes and Live Oak Counties, Texas. Spectral reflectance profiles, 569 in total, were reduced to 125 representative spectral signatures, which were analyzed using the U.S. Geological Survey's (USGS) Material Identification and Characterization Algorithm (MICA). MICA uses an automated continuum-removal procedure together with a least-squares linear regression to determine the fit of observed sample spectral absorption features to those of reference mineral standards in a spectral library. The reference minerals include various clay, mica, carbonate, ferric and ferrous iron minerals and their mixtures. In addition, absorption feature band-depth analysis was done to identify rock surfaces exhibiting absorption features related to uranium and zeolite minerals, which were not included in the command files used to execute MICA. Rocks from each of the four geologic units produced broadly similar spectral signatures as a result of comparable mineral compositions, but there were some notable differences. For example, Ca- and Na-montmorillonite was matched most frequently to the spectral absorption features in 2-μm (∼2000–2500 nm) wavelengths, while goethite occurred often at 1-μm (∼400–1000 nm) wavelengths. The latter is related to limonitic iron-staining in and around oxidized zones of the uranium roll front as described in previous papers. Rocks of the Jackson Group differed from those of the Catahoula, Oakville and Goliad units in that the former exhibited spectral features we interpret as being due to the presence of lignite-bearing mudstone layers. Goliad rocks exhibit spectral features related to dolomite, gypsum, anhydrite, and an unidentified green clay mineral that is possibly glauconite . Jackson Group rocks also exhibit weak but well-resolved absorption features at 964 and 1157 nm related to either or both zeolite minerals clinoptilolite and heulandite. These zeolite minerals and a few spectra exhibiting hydrous silica absorption features are indicative of alteration of volcanic glass in tuffaceous mudstone and claystone layers. A few sample spectra exhibited strong absorption features at around 1135 nm related to the uranium mineral coffinite. Both the 1135 nm coffinite and 1157 nm zeolite absorption features overlap somewhat, potentially making them difficult to distinguish without additional hyperspectral field, laboratory or remote sensing data. The results of this study were compared to mixtures of minerals described for ore, gangue and alteration minerals in deposit models for sandstone-hosted uranium, sedimentary bentonite and sedimentary zeolite. Use of these spectra can help facilitate mapping of both waste materials from the legacy mining of the above commodities, as well as future exploration and resource assessment activities.

Texas↗

USGS2024

Use of Doppler velocity radars to monitor and predict debris and flood wave velocities and travel times in post-wildfire basins

The magnitude and timing of extreme events such as debris and floodflows (collectively referred to as floodflows) in post-wildfire basins are difficult to measure and are even more difficult to predict. To address this challenge, a sensor ensemble consisting of noncontact, ground-based (near-field), Doppler velocity (velocity) and pulsed (stage or gage height) radars, rain gages, and a redundant radio communication network was leveraged to monitor flood wave velocities, to validate travel times, and to compliment observations from NEXRAD weather radar. The sensor ensemble (DEbris and Floodflow Early warNing System, DEFENS) was deployed in Waldo Canyon, Pike National Forest, Colorado, USA, which was burned entirely (100 percent burned) by the Waldo Canyon fire during the summer of 2012 ( MTBS, 2020 ). Surface velocity, stage, and precipitation time series collected during the DEFENS deployment on 10 August 2015 were used to monitor and predict flood wave velocities and travel times as a function of stream discharge (discharge; streamflow). The 10 August 2015 event exhibited spatial and temporal variations in rainfall intensity and duration that resulted in a discharge equal to 5.01 cubic meters per second (m 3 /s). Discharge was estimated post-event using a slope-conveyance indirect discharge method and was verified using velocity radars and the probability concept algorithm. Mean flood wave velocities – represented by the kinematic celerity ( c k = 2.619 m e t e r s p e r s e c o n d , m / s ± 0.556 p e r c e n t "> 𝑐 𝑘 = 2.619 𝑚𝑒𝑡𝑒𝑟𝑠 𝑝𝑒𝑟 𝑠𝑒𝑐𝑜𝑛𝑑, m/s ± 0.556 𝑝𝑒𝑟𝑐𝑒𝑛𝑡) and dynamic celerity ( c d = 3.533 m / s ± 0.181 p e r c e n t a n d t h e i r u n c e r t a i n t i e s "> 𝑐 𝑑 = 3.533 m/s ± 0.181 𝑝𝑒𝑟𝑐𝑒𝑛𝑡) and their uncertainties were computed. L-moments were computed to establish probability density functions (PDFs) and associated statistics for each of the at-a-section hydraulic parameters to serve as a workflow for implementing alert networks in hydrologically similar basins that lack data. Measured flood wave velocities and travel times agreed well with predicted values. Absolute percent differences between predicted and measured flood wave velocities ranged from 1.6 percent to 49 percent and varied with water slope, hydraulic radius, and depth. The kinematic celerity was a better predictor for steep slopes and wide flood plains associated with the Upper Waldo and Middle Waldo radar streamgages; whereas, the dynamic celerity was a better surrogate for shallow slopes and incised channels such as the Lower Waldo radar streamgage. The method demonstrates the potential extensibility of a post-wildfire warning system by (1) leveraging multiple systems (i.e., weather radar, near-field velocity and stage radars, and rain gages) for accurate and timely warnings of debris and floodflows, (2) establishing an order of operations to site, install, and operate near-field radars and conventional rain gages to record floodflows, forecast travel times, and document geomorphic change in this basin and hydrologically similar basins that lack data, and (3) communicating data operationally with the Colorado Department of Transportation engineering staff, National Weather Service forecasters, and emergency managers.

Colorado↗

USGS2007

Saturn's icy satellites investigated by Cassini-VIMS. I. Full-disk properties: 350-5100 nm reflectance spectra and phase curves

Saturn's icy satellites are among the main scientific objectives of the Cassini-VIMS (Visual and Infrared Mapping Spectrometer) experiment. This paper contains a first systematic and comparative analysis of the full-disk spectral properties of Dione, Enceladus, Epimetheus, Hyperion, Iapetus, Mimas, Phoebe, Rhea and Tethys as observed by VIMS from July 2004 to June 2005. The disk integrated properties (350-5100 nm reflectance spectra and phase curves at 550-2232 nm) and images of satellites are reported and discussed in detail together with the observed geometry. In general, the spectra in the visible spectral range are almost featureless and can be classified according to the spectral slopes: from the bluish Enceladus and Phoebe to the redder Iapetus, Hyperion and Epimetheus. In the 1000-1300 nm range the spectra of Enceladus, Tethys, Mimas and Rhea are characterized by a negative slope, consistent with a surface largely dominated by water ice, while the spectra of Iapetus, Hyperion and Phoebe show a considerable reddening pointing out the relevant role played by darkening materials present on the surface. In between these two classes are Dione and Epimetheus, which have a flat spectrum in this range. The main absorption bands identified in the infrared are the 1520, 2020, 3000 nm H2O/OH bands (for all satellites), although Iapetus dark terrains show mostly a deep 3000 nm band while the 1520 and 2020 nm bands are very faint. In this spectral range, the Iapetus spectrum is characterized by a strong reddening. The CO2 band at 4260 nm and the Fresnel ice peak around 3100 nm are evident only on Hyperion, Phoebe and Iapetus. The phase curves at 550 and at 2232 nm are reported for all the available observations in the 0??-144?? range; Rhea shows an opposition surge at visible wavelengths in the 0.5??-1.17?? interval. The improvement on the retrieval of the full-disk reflectance spectra can be appreciated by a direct comparison with ground-based telescopic data available from literature. Finally, data processing strategies and recent upgrades introduced in the VIMS-V calibration pipeline (flat-field and destriping-despiking algorithm) are discussed in appendices. ?? 2006 Elsevier Inc. All rights reserved.

Icarus↗

USGS2015

Enumeration of Mars years and seasons since the beginning of telescopic exploration

A clarification for the enumeration of Mars Years prior to 1955 is presented, along with a table providing the Julian dates associated with Ls = 0° for Mars Years -183 (beginning of the telescopic study of Mars) to 100. A practical algorithm for computing Ls as a function of the Julian Date is provided. No new science results are presented

Icarus↗

USGS2018

A remote sensing-based model of tidal marsh aboveground carbon stocks for the conterminous United States

Remote sensing based maps of tidal marshes, both of their extents and carbon stocks, have the potential to play a key role in conducting greenhouse gas inventories and implementing climate mitigation policies. Our objective was to generate a single remote sensing model of tidal marsh aboveground biomass and carbon that represents nationally diverse tidal marshes within the conterminous United States (CONUS). We developed the first calibration-grade, national-scale dataset of aboveground tidal marsh biomass, species composition, and aboveground plant carbon content (%C) from six CONUS regions: Cape Cod, MA, Chesapeake Bay, MD, Everglades, FL, Mississippi Delta, LA, San Francisco Bay, CA, and Puget Sound, WA. Using the random forest machine learning algorithm, we tested whether imagery from multiple sensors, Sentinel-1 C-band synthetic aperture radar, Landsat, and the National Agriculture Imagery Program (NAIP), can improve model performance. The final model, driven by six Landsat vegetation indices and with the soil adjusted vegetation index as the most important (n = 409, RMSE = 310 g/m 2 , 10.3% normalized RMSE), successfully predicted biomass for a range of marsh plant functional types defined by height, leaf angle and growth form. Model results were improved by scaling field-measured biomass calibration data by NAIP-derived 30 m fraction green vegetation. With a mean plant carbon content of 44.1% (n = 1384, 95% C.I. = 43.99%–44.37%), we generated regional 30 m aboveground carbon density maps for estuarine and palustrine emergent tidal marshes as indicated by a modified NOAA Coastal Change Analysis Program map. We applied a multivariate delta method to calculate uncertainties in regional carbon densities and stocks that considered standard error in map area, mean biomass and mean %C. Louisiana palustrine emergent marshes had the highest C density (2.67 ± 0.004 Mg/ha) of all regions, while San Francisco Bay brackish/saline marshes had the highest C density of all estuarine emergent marshes (2.03 ± 0.004 Mg/ha). Estimated C stocks for predefined jurisdictional areas ranged from 1023 ± 39 Mg in the Nisqually National Wildlife Refuge in Washington to 507,761 ± 14,822 Mg in the Terrebonne and St. Mary Parishes in Louisiana. This modeling and data synthesis effort will allow for aboveground C stocks in tidal marshes to be included in the coastal wetland section of the U.S. National Greenhouse Gas Inventory. With the increased availability of free post-processed satellite data, we provide a tractable means of modeling tidal marsh aboveground biomass and carbon at the global extent as well.

ISPRS Journal of Photogrammetry and Remote Sensing↗

USGS2011

Monitoring of urban subsidence with SAR interferometric point target analysis: A case study in Suzhou, China

Ground subsidence, mainly caused by over exploitation of groundwater and other underground resources, such as oil, gas and coal, occurs in many cities in China. The annual direct loss associated with subsidence across the country is estimated to exceed 100 million US dollar. Interferometric SAR (InSAR) is a powerful tool to map ground deformation at an unprecedented level of spatial detail. It has been widely used to investigate the deformation resulting from earthquakes, volcanoes and subsidence. Repeat-pass InSAR, however, may fail due to impacts of spatial decorrelation, temporal decorrelation and heterogeneous refractivity of atmosphere. In urban areas, a large amount of natural stable radar reflectors exists, such as buildings and engineering structures, at which radar signals can remain coherent during a long time interval. Interferometric point target analysis (IPTA) technique, also known as persistent scatterers (PS) InSAR is based on these reflectors. It overcomes the shortfalls in conventional InSAR. This paper presents a procedure for urban subsidence monitoring with IPTA. Calculation of linear deformation rate and height residual, and the non-linear deformation estimate, respectively, are discussed in detail. Especially, the former is highlighted by a novel and easily implemented 2-dimensional spatial search algorithm. Practically useful solutions that can significantly improve the robustness of IPTA, are recommended. Finally, the proposed procedure is applied to mapping the ground subsidence in Suzhou city, Jiangsu province, China. Thirty-four ERS-1/2 SAR scenes are analyzed, and the deformation information over 38,881 point targets between 1992 and 2000 are generated. The IPTA-derived deformation estimates correspond well with leveling measurements, demonstrating the potential of the proposed subsidence monitoring procedure based on IPTA technique. Two shortcomings of the IPTA-based procedure, e.g., the requirement of large number of SAR images and assumed linear plus non-linear deformation model, are discussed as the topics of further research.

Jaingsu Province↗

USGS2013

Hierarchical Bayesian spatial models for predicting multiple forest variables using waveform LiDAR, hyperspectral imagery, and large inventory datasets

In this paper we detail a multivariate spatial regression model that couples LiDAR, hyperspectral and forest inventory data to predict forest outcome variables at a high spatial resolution. The proposed model is used to analyze forest inventory data collected on the US Forest Service Penobscot Experimental Forest (PEF), ME, USA. In addition to helping meet the regression model's assumptions, results from the PEF analysis suggest that the addition of multivariate spatial random effects improves model fit and predictive ability, compared with two commonly applied modeling approaches. This improvement results from explicitly modeling the covariation among forest outcome variables and spatial dependence among observations through the random effects. Direct application of such multivariate models to even moderately large datasets is often computationally infeasible because of cubic order matrix algorithms involved in estimation. We apply a spatial dimension reduction technique to help overcome this computational hurdle without sacrificing richness in modeling.

Maine↗

USGS2014

Performance and effects of land cover type on synthetic surface reflectance data and NDVI estimates for assessment and monitoring of semi-arid rangeland

Federal land management agencies provide stewardship over much of the rangelands in the arid andsemi-arid western United States, but they often lack data of the proper spatiotemporal resolution andextent needed to assess range conditions and monitor trends. Recent advances in the blending of com-plementary, remotely sensed data could provide public lands managers with the needed information.We applied the Spatial and Temporal Adaptive Reflectance Fusion Model (STARFM) to five Landsat TMand concurrent Terra MODIS scenes, and used pixel-based regression and difference image analyses toevaluate the quality of synthetic reflectance and NDVI products associated with semi-arid rangeland. Pre-dicted red reflectance data consistently demonstrated higher accuracy, less bias, and stronger correlationwith observed data than did analogous near-infrared (NIR) data. The accuracy of both bands tended todecline as the lag between base and prediction dates increased; however, mean absolute errors (MAE)were typically ≤10%. The quality of area-wide NDVI estimates was less consistent than either spectra lband, although the MAE of estimates predicted using early season base pairs were ≤10% throughout the growing season. Correlation between known and predicted NDVI values and agreement with the 1:1regression line tended to decline as the prediction lag increased. Further analyses of NDVI predictions,based on a 22 June base pair and stratified by land cover/land use (LCLU), revealed accurate estimates through the growing season; however, inter-class performance varied. This work demonstrates the successful application of the STARFM algorithm to semi-arid rangeland; however, we encourage evaluation of STARFM’s performance on a per product basis, stratified by LCLU, with attention given to the influence of base pair selection and the impact of the time lag.

Idaho, Utah, Wyoming↗

USGS2019

Mapping cropland extent of Southeast and Northeast Asia using multi-year time-series Landsat 30-m data using Random Forest classifier on Google Earth Engine

Cropland extent maps are useful components for assessing food security. Ideally, such products are a useful addition to countrywide agricultural statistics since they are not politically biased and can be used to calculate cropland area for any spatial unit from an individual farm to various administrative unites (e.g., state, county, district) within and across nations, which in turn can be used to estimate agricultural productivity as well as degree of disturbance on food security from natural disasters and political conflict. However, existing cropland extent maps over large areas (e.g., Country, region, continent, world) are derived from coarse resolution imagery (250 m to 1 km pixels) and have many limitations such as missing fragmented and\or small farms with mixed signatures from different crop types and\or farming practices that can be, confused with other land cover. As a result, the coarse resolution maps have limited useflness in areas where fields are small (<1 ha), such as in Southeast Asia. Furthermore, coarse resolution cropland maps have known uncertainties in both geo-precision of cropland location as well as accuracies of the product. To overcome these limitations, this research was conducted using multi-date, multi-year 30-m Landsat time-series data for 3 years chosen from 2013 to 2016 for all Southeast and Northeast Asian Countries (SNACs), which included 7 refined agro-ecological zones (RAEZ) and 12 countries (Indonesia, Thailand, Myanmar, Vietnam, Malaysia, Philippines, Cambodia, Japan, North Korea, Laos, South Korea, and Brunei). The 30-m (1 pixel = 0.09 ha) data from Landsat 8 Operational Land Imager (OLI) and Landsat 7 Enhanced Thematic Mapper (ETM+) were used in the study. Ten Landsat bands were used in the analysis (blue, green, red, NIR, SWIR1, SWIR2, Thermal, NDVI, NDWI, LSWI) along with additional layers of standard deviation of these 10 bands across 1 year, and global digital elevation model (GDEM)-derived slope and elevation bands. To reduce the impact of clouds, the Landsat imagery was time-composited over four time-periods (Period 1: January- April, Period 2: May-August, and Period 3: September-December) over 3-years. Period 4 was the standard deviation of all 10 bands taken over all images acquired during the 2015 calendar year. These four period composites, totaling 42 band data-cube, were generated for each of the 7 RAEZs. The reference training data (N = 7849) generated for the 7 RAEZ using sub-meter to 5-m very high spatial resolution imagery (VHRI) helped generate the knowledge-base to separate croplands from non-croplands. This knowledge-base was used to code and run a pixel-based random forest (RF) supervised machine learning algorithm on the Google Earth Engine (GEE) cloud computing environment to separate croplands from non-croplands. The resulting cropland extent products were evaluated using an independent reference validation dataset (N = 1750) in each of the 7 RAEZs as well as for the entire SNAC area. For the entire SNAC area, the overall accuracy was 88.1% with a producer’s accuracy of 81.6% (errors of omissions = 18.4%) and user’s accuracy of 76.7% (errors of commissions = 23.3%). For each of the 7 RAEZs overall accuracies varied from 83.2 to 96.4%. Cropland areas calculated for the 12 countries were compared with country areas reported by the United Nations Food and Agriculture Organization and other national cropland statistics resulting in an R 2 value of 0.93. The cropland areas of provinces were compared with the province statistics that showed an R 2 = 0.95 for South Korea and R 2 = 0.94 for Thailand. The cropland products are made available on an interactive viewer at www.croplands.org and for download at National Aeronautics and Space Administration’s (NASA) Land Processes Distributed Active Archive Center (LP DAAC): https://lpdaac.usgs.gov/node/1281 .

International Journal of Applied Earth Observation↗

USGS2020

Surface soil temperature seasonal variation estimation in a forested area using combined satellite observations and in-situ measurements

Surface soil temperature is the soil temperature from the surface to 10 cm in depth. Surface soil temperature plays a significant role in agricultural drought monitoring, ecosystem energy transfer modeling, and global carbon cycle evaluation. Studies have been proposed to estimate surface soil temperature, but surface soil temperature monitoring within forested areas still poses a significant challenge. In this study, we proposed a surface soil temperature retrieval method using combined satellite observations and in-situ measurements for the Great Dismal Swamp (GDS). The GDS is a U.S. protected area managed and protected by the U.S. Fish and Wildlife Service. It is located along the boundary of Virginia and North Carolina, with maple gum, Atlantic white cedar, and pine pocosin as the main forest cover types. Ground-based surface soil temperature measurements were collected for these forest types from May 2015 to April 2017. Both the Land Remote Sensing Satellite (Landsat) Thermal Infrared Sensor (TIRS) and the Moderate Resolution Imaging Spectroradiometer (MODIS) carry two thermal infrared (TIR) channels. The TIR channels with similar corresponding wavelengths were first fused using an improved fusing model to generate high resolution TIR measurements. Then the enterprise algorithm was applied to calculate land surface temperature (LST) from the fused TIR bands. An improved soil temperature retrieval method was applied to generate surface soil temperature based on LST and vegetation index (VI) within the study area for the three forest types. In-situ measurements were used to build the surface soil temperature retrieval method, and results were then validated. The normalized difference vegetation index (NDVI) and enhanced vegetation index (EVI) were integrated separately as VIs in the model to monitor surface soil temperature. The R 2 for retrieved surface soil temperature through satellite observations was 0.76, and the RMSE was 1.96 ℃ ">℃ ℃ when NDVI was integrated in the model; the R 2 was 0.78, and the RMSE was 1.85 ℃ ">℃ ℃ when EVI was used.

North Carolina, Virginia↗

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↗

USGS2022

Bathymetry retrieval from CubeSat image sequences with short time lags

The rapid expansion of CubeSat constellations could revolutionize the way inland and nearshore coastal waters are monitored from space. This potential stems from the ability of CubeSats to provide daily imagery with global coverage at meter-scale spatial resolution. In this study, we explore the unique opportunity to improve the retrieval of bathymetry offered by CubeSats, specifically those of the PlanetScope constellation. The orbital design of the PlanetScope constellation enables the acquisition of image sequences with short time lags (from seconds to hours). This characteristic allows multiple images to be captured during a short period of steady bathymetric conditions, especially in dynamic environments like rivers. We hypothesize that taking the ensemble mean of a CubeSat image sequence can enhance bathymetry retrieval compared to standard single-image analysis. Along with the existing optimal band ratio analysis (OBRA) algorithm, we also use a new neural network-based depth retrieval (NNDR) technique to infer bathymetry from both individual and time-averaged images. The two methodologies are evaluated using field data from five different river reaches with depths up to 15 m and both top-of-atmosphere (TOA) radiance and bottom-of-atmosphere (BOA) surface reflectance PlanetScope data products. Despite low spectral resolution and concerns about the radiometric quality of CubeSat imagery, accuracy assessment based on in-situ comparisons indicates the potential (0.52 < R 2 < 0.7 for the NNDR method) of PlanetScope imagery to retrieve depths up to ∼ 10 m in clear water conditions. The proposed image averaging consistently improves bathymetry retrieval over single image analysis. The NNDR technique was found to outperform OBRA, illustrating the importance of leveraging all spectral bands through machine learning approaches. TOA data provided more robust bathymetry results than BOA data for the OBRA technique, but the NNDR technique was minimally impacted by the type of data product.

International Journal of Applied Earth Observation↗

USGS2013

Reevaluation of lake trout and lake whitefish bioenergetics models

Using a corrected algorithm for balancing the energy budget, we reevaluated the Wisconsin bioenergetics model for lake trout (Salvelinus namaycush) in the laboratory and for lake whitefish (Coregonus clupeaformis) in the laboratory and in the field. For lake trout, results showed that the bioenergetics model slightly overestimated food consumption by the lake trout when they were fed low and intermediate rations, whereas the model predicted food consumption by lake trout fed ad libitum without any detectable bias. The slight bias in model predictions for lake trout on restricted rations may have been an artifact of the feeding schedule for these fish, and we would therefore recommend application of the Wisconsin lake trout bioenergetics model to lake trout populations in the field without any revisions to the model. Use of the Wisconsin bioenergetics model for coregonids resulted in overestimation of food consumption by lake whitefish both in the laboratory and in the field by between 20 and 30%, on average. This overestimation of food consumption was most likely due to overestimation of respiration rate. We therefore adjusted the respiration component of the bioenergetics model to obtain a good fit to the observed consumption in our laboratory tanks. The adjusted model predicted the consumption in the laboratory and the field without any detectable bias. Until a detailed lake whitefish respiration study can be conducted, we recommend application of our adjusted version of the Wisconsin generalized coregonid bioenergetics model to lake whitefish populations in the field.

Journal of Great Lakes Research↗

USGS2019

A comparison of chlorophyll a values obtained from an autonomous underwater vehicle to satellite-based measures for Lake Michigan

Accurate methods to track changes in lake productivity through time and space are critical to fisheries management. Chlorophyll a is the most widely studied proxy for ecosystem primary production, and has been the topic of many studies. The main sources of chlorophyll a measurements are ship-based measures or multi-spectral satellite data. Autonomous underwater vehicles can survey large spatial extents approaching the scale of satellite data, but with the accuracy of ship-based water sampling methods. We use several statistical measures to compare measures of chlorophyll a collected in Lake Michigan with spatiotemporally matched satellite-derived measures of chlorophyll a from the MODIS Aqua multi-spectral sensor using NASA’s OC3 and the Great Lakes Fit algorithms. Our findings show a near one to one relationship between AUV data and both satellite-derived data sets when the AUV data are coarsened to the resolution of the satellite data. A comparison of satellite-based chlorophyll a to AUV-derived chlorophyll summarized in discrete water depth bins suggested that, based on decreasing coefficients of determination, satellite estimates of chlorophyll accounted for the most variability in chlorophyll a concentrations in the upper 10 m of the water column, even though satellite sensors may detect past this depth.

Lake Michigan↗

USGS2016

Uncertainty analysis of the Operational Simplified Surface Energy Balance (SSEBop) model at multiple flux tower sites

Evapotranspiration (ET) is an important component of the water cycle – ET from the land surface returns approximately 60% of the global precipitation back to the atmosphere. ET also plays an important role in energy transport among the biosphere, atmosphere, and hydrosphere. Current regional to global and daily to annual ET estimation relies mainly on surface energy balance (SEB) ET models or statistical and empirical methods driven by remote sensing data and various climatological databases. These models have uncertainties due to inevitable input errors, poorly defined parameters, and inadequate model structures. The eddy covariance measurements on water, energy, and carbon fluxes at the AmeriFlux tower sites provide an opportunity to assess the ET modeling uncertainties. In this study, we focused on uncertainty analysis of the Operational Simplified Surface Energy Balance (SSEBop) model for ET estimation at multiple AmeriFlux tower sites with diverse land cover characteristics and climatic conditions. The 8-day composite 1-km MODerate resolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) was used as input land surface temperature for the SSEBop algorithms. The other input data were taken from the AmeriFlux database. Results of statistical analysis indicated that the SSEBop model performed well in estimating ET with an R2 of 0.86 between estimated ET and eddy covariance measurements at 42 AmeriFlux tower sites during 2001–2007. It was encouraging to see that the best performance was observed for croplands, where R2 was 0.92 with a root mean square error of 13 mm/month. The uncertainties or random errors from input variables and parameters of the SSEBop model led to monthly ET estimates with relative errors less than 20% across multiple flux tower sites distributed across different biomes. This uncertainty of the SSEBop model lies within the error range of other SEB models, suggesting systematic error or bias of the SSEBop model is within the normal range. This finding implies that the simplified parameterization of the SSEBop model did not significantly affect the accuracy of the ET estimate while increasing the ease of model setup for operational applications. The sensitivity analysis indicated that the SSEBop model is most sensitive to input variables, land surface temperature (LST) and reference ET ( ET o ); and parameters, differential temperature ( dT ), and maximum ET scalar ( K max ), particularly during the non-growing season and in dry areas. In summary, the uncertainty assessment verifies that the SSEBop model is a reliable and robust method for large-area ET estimation. The SSEBop model estimates can be further improved by reducing errors in two input variables ( ET o and LST) and two key parameters ( K max and dT ).

Journal of Hydrology↗

USGS2016

An overview of current applications, challenges, and future trends in distributed process-based models in hydrology

Process-based hydrological models have a long history dating back to the 1960s. Criticized by some as over-parameterized, overly complex, and difficult to use, a more nuanced view is that these tools are necessary in many situations and, in a certain class of problems, they are the most appropriate type of hydrological model. This is especially the case in situations where knowledge of flow paths or distributed state variables and/or preservation of physical constraints is important. Examples of this include: spatiotemporal variability of soil moisture, groundwater flow and runoff generation, sediment and contaminant transport, or when feedbacks among various Earth’s system processes or understanding the impacts of climate non-stationarity are of primary concern. These are situations where process-based models excel and other models are unverifiable. This article presents this pragmatic view in the context of existing literature to justify the approach where applicable and necessary. We review how improvements in data availability, computational resources and algorithms have made detailed hydrological simulations a reality. Avenues for the future of process-based hydrological models are presented suggesting their use as virtual laboratories, for design purposes, and with a powerful treatment of uncertainty.

Journal of Hydrology↗