Geology ReportsSearch

SEARCH · Geology Reports

Results for “Energy”

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.

At least 433 records · Page 24Linked to original sources

Soil disturbance as a driver of increased stream salinity in a semiarid watershed undergoing energy development

Salinization is a global threat to the quality of streams and rivers, but it can have many causes. Oil and gas development were investigated as one of several potential causes of changes in the salinity of Muddy Creek, which drains 2470 km2 of mostly public land in Wyoming, U.S.A. Stream discharge and salinity vary with seasonal snowmelt and define a primary salinity-discharge relationship. Salinity, measured by specific conductance, increased substantially in 2009 and was 53-71% higher at low discharge and 33-34% higher at high discharge for the years 2009-2012 compared to 2005-2008. Short-term processes (e.g., flushing of efflorescent salts) cause within-year deviations from the primary relation but do not obscure the overall increase in salinity. Dissolved elements associated with increased salinity include calcium, magnesium, and sulfate, a composition that points to native soil salts derived from marine shales as a likely source. Potential causes of the salinity increase were evaluated for consistency by using measured patterns in stream chemistry, slope of the salinity-discharge relationship, and inter-annual timing of the salinity increase. Potential causes that were inconsistent with one or more of those criteria included effects from precipitation, evapotranspiration, reservoirs, grazing, irrigation return flow, groundwater discharge, discharge of energy co-produced waters, and stream habitat restoration. In contrast, surface disturbance of naturally salt-rich soil by oil and gas development activities, such as pipeline, road, and well pad construction, is a reasonable candidate for explaining the salinity increase. As development continues to expand in semiarid lands worldwide, the potential for soil disturbance to increase stream salinity should be considered, particularly where soils host substantial quantities of native salts.

Journal of Hydrology: Regional Studies

Evaluation of streamflow predictions from LSTM models in water- and energy-limited regions in the United States

The application of Long Short-Term Memory (LSTM) models for streamflow predictions has been an area of rapid development, supported by advancements in computing technology, increasing availability of spatiotemporal data, and availability of historical data that allows for training data-driven LSTM models. Several studies have focused on improving the performance of LSTM models; however, few studies have assessed the applicability of these LSTM models across different hydroclimate regions. This study investigated the single-basin trained local (one model for each basin), multi-basin trained regional (one model for one region), and grand (one model for several regions) models for predicting daily streamflow in water-limited Great Basin (18 basins) and energy-limited New England (27 basins) regions in the United States using the CAMELS (Catchment Attributes and Meteorology for Large-sample Studies) data set. The results show a general pattern of higher accuracy in daily streamflow predictions from the regional model when compared to local or grand models for most basins in the New England region. For the Great Basin region, local models provided smaller errors for most basins and substantially lower for those basins with relatively larger errors from the regional and grand models. The evaluation of one-layer and three-layer LSTM network architectures trained with 1-day lag information indicates that the addition of model complexity by increasing the number of layers may not necessarily increase the model skill for improving streamflow predictions. Findings from our study highlight the strengths and limitations of LSTM models across contrasting hydroclimate regions in the United States, which could be useful for local and regional scale decisions using standalone or potential integration of data-driven LSTM models with physics-based hydrological models.

Machine Learning with Applications

Estimating wildfire fuel consumption with multitemporal airborne laser scanning data and demonstrating linkage with MODIS-derived fire radiative energy

Characterizing pre- and post-fire fuels remains a key challenge for estimating biomass consumption and carbon emissions from wildfires. Airborne laser scanning (ALS) data have demonstrated effectiveness for estimating canopy, and to a lesser degree, surface fuel components at fine-scale (i.e., 30 m) across landscapes. Using pre- and post-fire ALS data and corresponding field data, this study estimated consumption of canopy fuel (ΔCF), understory fuel (ΔUF), total fuel (ΔTF), and canopy bulk density (ΔCBD) for the 2012 Pole Creek fire in Oregon, USA (10,760 ha), and portions of the 2011 Las Conchas fire in New Mexico, USA (4,934 ha). Additionally, the feasibility of predicting fuel consumption was tested using separate pre- and post-fire models (PrePost), models combining all pre- and post-fire data (Pooled), and models using all data from both fires (Global). Estimates of ΔTF were then compared to fire radiative energy (FRE, units: MJ) derived from Fire Radiative Power (FRP, units: MW) observations from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor onboard NASA Terra and Aqua satellites to mechanistically derive a biomass combustion coefficient (BCC, units: kg MJ −1 ). The PrePost and Pooled approaches yielded similar results at Las Conchas, but at Pole Creek insufficient pre-fire field data resulted in erroneous fuel consumption estimates outside the fire perimeter using the PrePost models. These results demonstrated that pre-fire field data were less important for these models than having field data which represent the full range of fuel conditions likely to exist across the landscape. Estimated total biomass consumed for the PrePost, Pooled, and Global models were 226 Gg, 224 Gg, and 224 Gg at Las Conchas, and 581 Gg, 713 Gg, and 552 Gg at Pole Creek. Comparisons between estimated ΔTF and FRE yielded an average BCC for both fires of 0.367 (s.d. ± 0.049) kg MJ −1 based on pixels with at least five MODIS observations. Both higher MODIS observations per pixel and accounting for canopy occlusion of FRE improved the relationship between ΔTF and MODIS-FRE. This study suggested a practical modelling approach for future efforts using only post-fire field observations and quantified a landscape-scale relationship between MODIS-derived FRE and fine-scale fuel consumption consistent with prior experiments.

New Mexico, Oregon

Risk assessment of brine contamination to aquatic resources from energy development in glacial drift deposits: Williston Basin, USA

Contamination to aquatic resources from co-produced water (brine) associated with energy development has been documented in the northeastern portion of the Williston Basin; an area mantled by glacial drift. The presence and magnitude of brine contamination can be determined using the contamination index (CI) value from water samples. Recently, the U.S. Geological Survey published a section (~ 2.59 km 2 ) level risk assessment of brine contamination to aquatic resources for Sheridan County, Montana, using oilfield and hydrogeological parameters. Our goal was to improve the Sheridan County assessment (SCA) and evaluate the use of this new Williston Basin assessment (WBA) across 31 counties mantled by glacial drift in the Williston Basin. To determine if the WBA model improved the SCA model, results from both assessments were compared to CI values from 37 surface and groundwater samples collected to evaluate the SCA. The WBA (R 2 = 0.65) outperformed the SCA (R 2 = 0.52) indicating improved model performance. Applicability across the Williston Basin was evaluated by comparing WBA results to CI values from 123 surface water samples collected from 97 sections. Based on the WBA, the majority (83.5%) of sections lacked an oil well and had minimal risk. Sections with one or more oil wells comprised low (8.4%), moderate (6.5%), or high (1.7%) risk areas. The percentage of contaminated water samples, percentage of sections with at least one contaminated sample, and the average CI value of contaminated samples increased from low to high risk indicating applicability across the Williston Basin. Furthermore, the WBA performed better compared to only the contaminated samples (R 2 = 0.62) versus all samples (R 2 = 0.38). This demonstrates that the WBA was successful at identifying sections, but not individual aquatic resources, with an increased risk of contamination; therefore, WBA results can prioritize future sampling within areas of increased risk.

Montana, North Dakota

Corrigendum to “A comprehensive analysis of interseasonal and interannual energy and water balance dynamics in semiarid shrubland and forest ecosystems” [Sci. Total Environ. 651 (2019) 381–398]

The authors regret the omission of an author, funding sources, and key support staff. The omitted author and their affiliation is: Matthew J Germino US Geological Survey, Forest and Rangeland Ecosystem Science Center, Boise ID 83706. The correct citation should be: Valayamkunnath, P., Sridhar, V., Zhao, W. Allen, R.G., and Germino, M.J., 2019. A comprehensive analysis of interseasonal and interannual energy and water balance dynamics in semiarid shrubland and forest ecosystems. Science of The Total Environment, 651, 381–398. Additions to the Acknowledgements include: Funding was provided by the U.S. Bureau of Land Management and US Department of Agriculture National Research Initiative. Jeremy Greth, Clarence Robinson, Ricardo Trezza, Jeppe Kjaersgaard, and many other field technicians from Idaho State University and Boise State University assisted in setting up the experiment. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.

Science of the Total Environment

Radionuclides in fly ash and bottom ash: Improved characterization based on radiography and low energy gamma-ray spectrometry

Two radiation-based techniques for determining the distribution and relative abundance of radionuclides are described, and applied to a suite of fly ash and bottom ash samples from a Kentucky power plant. The technique of fission-track radiography provides new observations of the variety of uranium hosts and of uranium distribution in individual particles of fly ash, and thus aids prediction of the leachability of uranium during long-term disposal or utilization of fly ash. Uranium is largely dispersed within glassy components of fly ash particles and shows little evidence for obvious surface enrichment that could be attributed to secondary adsorption. The technique of low energy gamma-ray spectrometry provides simultaneous, non-destructive determination of the relative abundance of 238 U, 226 Ra, 228 Ra and 210 Pb in representative 150–250 g samples. The measurements provide a means for screening samples to determine if the combustion process causes significant preferential redistribution of radionuclides that could affect their subsequent mobility. Results indicate that radium isotopes are not significantly (within 10–15%) fractionated from parent 238 U and 232 Th during coal combustion. In contrast, 210 Pb appears to be preferentially enriched in some samples of fly ash, and depleted in bottom ash relative to 238 U and 226 Ra. In this application 210 Pb acts as a tracer for elemental lead, and confirms the expected greater volatility of Pb compared to more refractory elements during coal combustion.

Fuel

Proximate composition and energy density of some North Pacific forage fishes

Mature pelagic forage fish species (capelin, sand lance, squid) had greater lipid concentrations than juvenile age-classes of large demersal and pelagic fish species (walleye pollock, Pacific cod, Atka mackerel, greenling, prowfish, rockfish, sablefish). Myctophids preyed on by puffins have at least twice as much lipid per gram compared to mature capelin, sand lance and squid, and an order of magnitude greater lipid concentrations than juvenile forage fish. Energy density of forage fishes was positively correlated with lipid content, and negatively correlated with water, ash-free lean dry mass (mostly protein), and ash contents.

Comparative Biochemistry and Physiology, Part A: M

Methane emissions from artificial waterbodies dominate the carbon footprint of irrigation: A study of transitions in the food-energy-water-climate nexus (Spain, 1900-2014)

Irrigation in the Mediterranean region has been used for millennia and has greatly expanded with industrialization. Irrigation is critical for climate change adaptation, but it is also an important source of greenhouse gas emissions. This study analyzes the carbon (C) footprint of irrigation in Spain, covering the complete historical process of mechanization. A 21-fold total, 6-fold area-based, and 4-fold product-based increase in the carbon footprint was observed during the 20th century, despite an increase in water use efficiency. CH 4 emissions from waterbodies, which had not previously been considered in the C footprint of irrigation systems, dominated the emission budget during most of the analyzed period. Technologies to save water and tap new water resources greatly increased energy and infrastructure demand, while improvements in power generation efficiency had a limited influence on irrigation emissions. Electricity production from irrigation dams may contribute to climate change mitigation, but the amount produced in relation to that consumed in irrigation has greatly declined. High uncertainty in CH 4 emission estimates from waterbodies stresses a need for more spatially resolved data and an improved empirical knowledge of the links between water quality, water level fluctuations, and emissions at the regional scale.

Environmental Science & Technology

Linear solvation energy relationships: "rule of thumb" for estimation of variable values

For the linear solvation energy relationship (LSER), values are listed for each of the variables (V i /100, π*, &beta m , α m ) for fundamental organic structures and functional groups. We give the guidelines to estimate LSER variable values quickly for a vast array of possible organic compounds such as those found in the environment. The difficulty in generating these variables has greatly discouraged the application of this quantitative structure-activity relationship (QSAR) method. This paper present the first compilation of molecular functional group values together with a utilitarian set of the LSER variable estimation rules. The availability of these variable values and rules should facilitate widespread application of LSER for hazard evaluation of environmental contaminants.

Environmental Science & Technology

Method development for analysis of urban dust using scanning electron microscopy with energy dispersive x-ray spectrometry to detect the possible presence of world trade center dust constituents

The collapse of the World Trade Center Towers on September 11, 2001, sent dust and debris across much of Manhattan and in the surrounding areas. Indoor and outdoor dust samples were collected and characterized by U.S. Geological Survey (USGS) scientists using scanning electron microscopy with energy-dispersive spectrometry (SEM/EDS). From this characterization, the U.S. Environmental Protection Agency and USGS developed a particulate screening method to determine the presence of residual World Trade Center dust in the indoor environment using slag wool as a primary "signature". The method describes a procedure that includes splitting, ashing, and sieving of collected dust. From one split, a 10 mg/mL dust/ isopropanol suspension was prepared and 10-30 ??L aliquots of the suspension placed on an SEM substrate. Analyses were performed using SEM/EDS manual point counting for slag wool fibers. Poisson regression was used to identify some of the sources of uncertainty, which are directly related to the small number of fibers present on each sample stub. Preliminary results indicate that the procedure is promising for screening urban background dust for the presence of WTC dust. Consistent sample preparation of reference materials and samples must be performed by each laboratory wishing to use this method to obtain meaningful and accurate results. ?? 2009 American Chemical Society.

Environmental Science & Technology

Linear solvation energy relationship of the limiting partition coefficient of organic solutes between water and activated carbon

A linear solvation energy relationship has been found for 353 values of the limiting adsorption coefficients of diverse chemicals: log K = −0.37 + 0.0341 V i − 1.07β + D + 0.65 P with R = 0.951, s = 0.51, n = 353, and F = 818.0, where V i is the intrinsic molar volume; β is a measure of the hydrogen bond acceptor strength of the solute; D is an index parameter for the research group which includes the effects of the different types of carbon used, the temperature, and the length of time allowed for the adsorption equilibrium to be established; and P is an index parameter for the flatness of the molecule. P is defined to be unity if there is an aromatic system in the molecule or if there is a double bond or series of conjugated double bonds with no more that one non-hydrogen atom beyond the double bond and zero otherwise. A slightly better fit is obtained if the two-thirds power of V i is used as a measure of the surface area in place of the volume term: log K = −1.75 + 0.227 V 2/3 − 1.10β + D + 0.60 P with R = 0.954, s = 0.49, n = 353, and F = 895.39. This is the first quantitative measure of the effect of the shape of the molecule on its tendency to be adsorbed on activated carbon.

Environmental Science & Technology

Evaluating the effects of historical land cover change on summertime weather and climate in New Jersey: Land cover and surface energy budget changes

The 19th-century agrarian landscape of New Jersey (NJ) and the surrounding region has been extensively transformed to the present-day land cover by urbanization, reforestation, and localized areas of deforestation. This study used a mesoscale atmospheric numerical model to investigate the sensitivity of the warm season climate of NJ to these land cover changes. Reconstructed 1880s-era and present-day land cover data sets were used as surface boundary conditions for a set of simulations performed with the Regional Atmospheric Modeling System (RAMS). Three-member ensembles with historical and present-day land cover were compared to examine the sensitivity of surface air and dew point temperatures, rainfall, and the individual components of the surface energy budget to these land cover changes. Mean temperatures for the present-day landscape were 0.3-0.6??C warmer than for the historical landscape over a considerable portion of NJ and the surrounding region, with daily maximum temperatures at least 1.0??C warmer over some of the highly urbanized locations. Reforested regions, however, were slightly cooler. Dew point temperatures decreased by 0.3-0.6??C, suggesting drier, less humid near-surface air for the present-day landscape. Surface warming was generally associated with repartitioning of net radiation from latent to sensible heat flux, and conversely for cooling. While urbanization was accompanied by strong surface albedo decreases and increases in net shortwave radiation, reforestation and potential changes in forest composition have generally increased albedos and also enhanced landscape heterogeneity. The increased deciduousness of forests may have further reduced net downward longwave radiation. Copyright 2008 by the American Geophysical Union.

New Jersey

Distributed energy balance modeling of South Cascade Glacier, Washington and assessment of model uncertainty

We have developed a physically based, distributed surface energy balance model to simulate glacier mass balance under meteorological and climatological forcing. Here we apply the model to estimate summer ablation on South Cascade Glacier, Washington, for the 2004 and 2005 mass balance seasons. To arrive at optimal mass balance simulations, we investigate and quantify model uncertainty associated with selecting from a range of physical parameter values that are not commonly measured in glaciological mass balance field studies. We optimize the performance of the model by varying values for atmospheric transmissivity, the albedo of surrounding topography, precipitation-elevation lapse rate, surface roughness for turbulent exchange of momentum, and snow albedo aging coefficient. Of these the snow aging parameter and precipitation lapse rates have the greatest influence on the modeled ablation. We examined model sensitivity to varying parameters by performing an additional 103 realizations with parameters randomly chosen over a ??5% range centered about the optimum values. The best fit suite of model parameters yielded a net balance of -1.69??0.38 m water equivalent (WE) for the 2004 water year and -2.10??0.30 m WE up to 11 September 2005. The 2004 result is within 3% of the measured value. These simulations account for 91% and 93% of the variance in measured ablation for the respective years. Copyright 2008 by the American Geophysical Union.

Journal of Geophysical Research F: Earth Surface

Climatic effects of 30 years of landscape change over the Greater Phoenix, Arizona, region: 1. Surface energy budget changes

This paper is part 1 of a two-part study that evaluates the climatic effects of recent landscape change for one of the nation's most rapidly expanding metropolitan complexes, the Greater Phoenix, Arizona, region. The region's landscape evolution over an approximate 30-year period since the early 1970s is documented on the basis of analyses of Landsat images and land use/land cover (LULC) data sets derived from aerial photography (1973) and Landsat (1992 and 2001). High-resolution, Regional Atmospheric Modeling System (RAMS), simulations (2-km grid spacing) are used in conjunction with consistently defined land cover data sets and associated biophysical parameters for the circa 1973, circa 1992, and circa 2001 time periods to quantify the impacts of intensive land use changes on the July surface temperatures and the surface radiation and energy budgets for the Greater Phoenix region. The main findings are as follows: since the early 1970s the region's landscape has been altered by a significant increase in urban/suburban land area, primarily at the expense of decreasing plots of irrigated agriculture and secondarily by the conversion of seminatural shrubland. Mean regional temperatures for the circa 2001 landscape were 0.12??C warmer than the circa 1973 landscape, with maximum temperature differences, located over regions of greatest urbanization, in excess of 1??C. The significant reduction in irrigated agriculture, for the circa 2001 relative to the circa 1973 landscape, resulted in dew point temperature decreases in excess of 1??C. The effect of distinct land use conversion themes (e.g., conversion from irrigated agriculture to urban land) was also examined to evaluate how the most important conversion themes have each contributed to the region's changing climate. The two urbanization themes studied (from an initial landscape of irrigated agriculture and seminatural shrubland) have the greatest positive effect on near-surface temperature, increasing maximum daily temperatures by 1??C. Overall, sensible heat flux differences between the circa 2001 and circa 1973 landscapes result in a 1 W m -2 increase in domain-wide sensible heating, and a similar order of magnitude decrease in latent heating, highlighting the importance of surface repartitioning in establishing near-surface temperature trends. In part 2 of this study, we address the role of the surface budget changes on the mesoscale dynamics/thermodynamics, in context of the large-scale environment. Copyright 2009 by the American Geophysical Union.

Journal of Geophysical Research D: Atmospheres

Carbon, water, and energy fluxes in a semiarid cold desert grassland during and following multiyear drought

The net exchanges of carbon dioxide, water vapor, and energy were examined in a perennial Colorado Plateau grassland for 5 years. The study began within a multiyear drought and continued as the drought ended. The grassland is located near the northern boundary of the influence of the North American monsoon, a major climatic feature bringing summer rain. Following rain, evapotranspiration peaked above 8 mm d -1 but was usually much smaller (2-4 mm d -1 ). Net productivity of the grassland was low compared to other ecosystems, with peak hourly net CO 2 uptake in the spring of 4 (mu or u)mol m -2 s -1 and springtime carbon gain in the range of 42 + or - 11 g C m -2 (based on fluxes) to 72 + or - 55 g C m -2 (based on carbon stocks; annual carbon gain was not quantified). Drought decreased gross ecosystem productivity (GEP) and total ecosystem respiration, with a much larger GEP decrease. Monsoon rains led to respiratory pulses, lasting a few days at most, and only rarely resulted in net CO 2 gain, despite the fact that C 4 grasses dominated plant cover. Minor CO 2 uptake was observed in fall following rain. Spring CO 2 uptake was regulated by deep soil moisture, which depended on precipitation in the prior fall and winter. The lack of CO 2 uptake during the monsoon and the dependence of GEP on deep soil moisture are in contrast with arid grasslands of the warm deserts. Cold desert grasslands are most likely to be impacted by future changes in winter and not summer precipitation.

Colorado Plateau

Moments, magnitudes, and radiated energies of non-volcanic tremor near Cholame, CA, from ground motion spectra at UPSAR

By averaging the spectra of events within two episodes of tremor (on Jan. 21 and 24, 2005) across the 12 stations of UPSAR, we improved the S/N sufficiently to define source spectra. Analysis of eleven impulsive events revealed attenuation-corrected spectra of displacement similar to those of earthquakes, with a low-frequency plateau, a corner frequency, and a high frequency decay proportional to f −2 . Seismic moments, M 0 , estimated from these spectra range from about 3 to 10 × 10 11 N-m or moment magnitudes in the range 1.6 to 1.9. The corner frequencies range from 2.6 to 7.2 Hz and, if interpreted in the same way as for earthquakes, indicate low stress drops that vary from 0.001 to 0.04 MPa. Seismic energies, estimated from the ground motion spectra, vary from 0.2 × 10 5 to 4.4 × 10 5 J, or apparent stresses in the range 0.002 to 0.02 MPa. The low stress parameters are consistent with a weak fault zone in the lower crust at the depth of tremor. In contrast, the same analysis on a micro-earthquake, located near Cholame (depth = 10.3 km), revealed a stress drop of 0.5 MPa and an apparent stress of 0.02 MPa. Residual spectra from ω −2 model fits to the displacement spectra of the non-volcanic tremor events show peaks near 4 Hz that are not apparent in the spectra for the microearthquake nor for the spectrum of earth noise. These spectral peaks may indicate that tremor entails more than shear failure reminiscent of mechanisms, possibly entailing fluid flow, associated with volcanic tremor or deep volcanic earthquakes.

Geophysical Research Letters

Resolving small-scale forest snow patterns using an energy-balance snow model with a 1-layer canopy

Modelling spatiotemporal dynamics of snow in forests is challenging, as involved processes are strongly dependent on small-scale canopy properties. In this study, we explore how local canopy structure information can be integrated in a medium-complexity energy-balance snow model to replicate observed snow patterns at very high spatial resolutions. Snow depth distributions simulated with the Flexible Snow Model (FSM2) were tested against extensive experimental data acquired in discontinuous subalpine forest stands in Eastern Switzerland over three winters. While the default canopy implementation in FSM2 fails to capture the observed snow depth variability, performance is considerably improved when local canopy cover fraction and hemispherical sky view fraction are additionally accounted for (30% reduction in RMSE). However, realistic snow depth distribution patterns throughout the season are only achieved if effective temperatures of near and distant canopy elements are discerned, and if a mechanism to mimic preferential deposition of snow in canopy gaps is included. We demonstrate that by diversifying the canopy structure input in order to reflect respective portions of the canopy relevant to different processes, even a simple model based on widely used process parametrizations and canopy metrics can be applied for high-resolution simulations of the sub-canopy snow cover with just a few modifications. The presented approaches could be implemented in commonly used land surface models, allowing upscaling experiments and development of sub-grid parametrizations without necessitating complex high-resolution models.

Water Resources Research

Evaluation of the energy budget method of determining evaporation at Williams Lake, Minnesota, using alternative instrumentation and study approaches

Best estimates of evaporation at Williams Lake, north central Minnesota, were determined by the energy budget method using optimum sensors and optimum placement of sensors. These best estimates are compared with estimates derived from using substitute data to determine the effect of using less accurate sensors, simpler methods, or remotely measured data. Calculations were made for approximately biweekly periods during five open water seasons. For most of the data substitutions that affected the Bowen ratio, new values of evaporation differed little from best estimates. The three data substitution methods that caused the largest deviations from the best evaporation estimates were (1) using changes in the daily average surface water temperature as an indicator of the lake heat storage term, (2) using shortwave radiation, air temperature, and atmospheric vapor pressure data from a site 110 km away, and (3) using an analog surface water temperature probe. Recalculations based on these data substitutions resulted in differences from the best estimates as much as 89%, 21%, and 10%, respectively. The data substitution method that provided evaporation values that most closely matched the best estimates was measurement of the lake heat storage term at one location in the lake, rather than at 16 locations. Evaporation values resulting from this substitution method usually were within 2% of the best estimates.

Minnesota