USGS Science⌕ Search

SEARCH · USGS Science

Results for “Agricultural Systems”

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 1,135 records · Page 63Linked to original sources

Bathymetric survey and digital elevation model of Little Holland Tract, Sacramento-San Joaquin Delta, California

The U.S. Geological Survey conducted a bathymetric survey in Little Holland Tract, a flooded agricultural tract, in the northern Sacramento-San Joaquin Delta (the “Delta”) during the summer of 2015. The new bathymetric data were combined with existing data to generate a digital elevation model (DEM) at 1-meter resolution. Little Holland Tract (LHT) was historically diked off for agricultural uses and has been tidally inundated since an accidental levee breach in 1983. Shallow tidal regions such as LHT have the potential to improve habitat quality in the Delta. The DEM of LHT was developed to support ongoing studies of habitat quality in the area and to provide a baseline for evaluating future geomorphic change. The new data comprise 138,407 linear meters of real-time-kinematic (RTK) Global Positioning System (GPS) elevation data, including both bathymetric data collected from personal watercraft and topographic elevations collected on foot at low tide. A benchmark (LHT15_b1) was established for geodetic control of the survey. Data quality was evaluated both by comparing results among surveying platforms, which showed systematic offsets of 1.6 centimeters (cm) or less, and by error propagation, which yielded a mean vertical uncertainty of 6.7 cm. Based on the DEM and time-series measurements of water depth, the mean tidal prism of LHT was determined to be 2,826,000 cubic meters. The bathymetric data and DEM are available at http://dx.doi.org/10.5066/F7RX9954 .

California↗

Landsat benefiting society for fifty years

Since 1972, data acquired by the Landsat series of satellites have become integral to land management for both government and the private sector, providing scientists and decision makers with key information about agricultural productivity, ice sheet dynamics, urban growth, forest monitoring, natural resource management, water quality, and supporting disaster response. Landsat 9 continues the mission of unrivaled space-based Earth observation and will lead the Landsat program into its second half century of Earth imagery provided to users, worldwide, at no charge. More than 8 million Landsat scenes held in the USGS archive to date are used in conjunction with advanced geographic information systems, image processing software, and cloud computing capabilities to enable Landsat users to study changes on the Earth’s surface across continental regions and extended time periods. The Operational Land Imager 2 (OLI-2) and the Thermal Infrared Sensor 2 (TIRS-2) instruments onboard Landsat 9 —which replicate the technologically-advanced instruments introduced onboard Landsat 8—allow for the collection of continuous high-quality data required for advancing Earth applications, including our ability to map surface temperature and surface water quality. Landsat 9 will build on the Landsat legacy, achieving a half-century record of global Earth observations.

Report↗

Evaluating land cover influences on model uncertainties—A case study of cropland carbon dynamics in the Mid-Continent Intensive Campaign region

Quantifying spatial and temporal patterns of carbon sources and sinks and their uncertainties across agriculture-dominated areas remains challenging for understanding regional carbon cycles. Characteristics of local land cover inputs could impact the regional carbon estimates but the effect has not been fully evaluated in the past. Within the North American Carbon Program Mid-Continent Intensive (MCI) Campaign, three models were developed to estimate carbon fluxes on croplands: an inventory-based model, the Environmental Policy Integrated Climate (EPIC) model, and the General Ensemble biogeochemical Modeling System (GEMS) model. They all provided estimates of three major carbon fluxes on cropland: net primary production (NPP), net ecosystem production (NEP), and soil organic carbon (SOC) change. Using data mining and spatial statistics, we studied the spatial distribution of the carbon fluxes uncertainties and the relationships between the uncertainties and the land cover characteristics. Results indicated that uncertainties for all three carbon fluxes were not randomly distributed, but instead formed multiple clusters within the MCI region. We investigated the impacts of three land cover characteristics on the fluxes uncertainties: cropland percentage, cropland richness and cropland diversity. The results indicated that cropland percentage significantly influenced the uncertainties of NPP and NEP, but not on the uncertainties of SOC change. Greater uncertainties of NPP and NEP were found in counties with small cropland percentage than the counties with large cropland percentage. Cropland species richness and diversity also showed negative correlations with the model uncertainties. Our study demonstrated that the land cover characteristics contributed to the uncertainties of regional carbon fluxes estimates. The approaches we used in this study can be applied to other ecosystem models to identify the areas with high uncertainties and where models can be improved to reduce overall uncertainties for regional carbon flux estimates.

Mid-Continent Intensive Campaign region↗

Simulation of groundwater conditions and streamflow depletion to evaluate water availability in a Freeport, Maine, watershed

In order to evaluate water availability in the State of Maine, the U.S. Geological Survey (USGS) and the Maine Geological Survey began a cooperative investigation to provide the first rigorous evaluation of watersheds deemed "at risk" because of the combination of instream flow requirements and proportionally large water withdrawals. The study area for this investigation includes the Harvey and Merrill Brook watersheds and the Freeport aquifer in the towns of Freeport, Pownal, and Yarmouth, Maine. A numerical groundwater- flow model was used to evaluate groundwater withdrawals, groundwater-surface-water interactions, and the effect of water-management practices on streamflow. The water budget illustrates the effect that groundwater withdrawals have on streamflow and the movement of water within the system. Streamflow measurements were made following standard USGS techniques, from May through September 2009 at one site in the Merrill Brook watershed and four sites in the Harvey Brook watershed. A record-extension technique was applied to estimate long-term monthly streamflows at each of the five sites. The conceptual model of the groundwater system consists of a deep, confined aquifer (the Freeport aquifer) in a buried valley that trends through the middle of the study area, covered by a discontinuous confining unit, and topped by a thin upper saturated zone that is a mixture of sandy units, till, and weathered clay. Harvey and Merrill Brooks flow southward through the study area, and receive groundwater discharge from the upper saturated zone and from the deep aquifer through previously unknown discontinuities in the confining unit. The Freeport aquifer gets most of its recharge from local seepage around the edges of the confining unit, the remainder is received as inflow from the north within the buried valley. Groundwater withdrawals from the Freeport aquifer in the study area were obtained from the local water utility and estimated for other categories. Overall, the public-supply withdrawals (105.5 million gallons per year (Mgal/yr)) were much greater than those for any other category, being almost 7 times greater than all domestic well withdrawals (15.3 Mgal/yr). Industrial withdrawals in the study area (2.0 Mgal/yr) are mostly by a company that withdraws from an aquifer at the edge of the Merrill Brook watershed. Commercial withdrawals are very small (1.0 Mgal/yr), and no irrigation or other agricultural withdrawals were identified in this study area. A three-dimensional, steady-state groundwater-flow model was developed to evaluate stream-aquifer interactions and streamflow depletion from pumping, to help refine the conceptual model, and to predict changes in streamflow resulting from changes in pumping and recharge. Groundwater levels and flow in the Freeport aquifer study area were simulated with the three-dimensional, finite-difference groundwater-flow modeling code, MODFLOW-2005. Study area hydrology was simulated with a 3-layer model, under steady-state conditions. The groundwater model was used to evaluate changes that could occur in the water budgets of three parts of the local hydrologic system (the Harvey Brook watershed, the Merrill Brook watershed, and the buried aquifer from which pumping occurs) under several different climatic and pumping scenarios. The scenarios were (1) no pumping well withdrawals; (2) current (2009) pumping, but simulated drought conditions (20-percent reduction in recharge); (3) current (2009) recharge, but a 50-percent increase in pumping well withdrawals for public supply; and (4) drought conditions and increased pumping combined. In simulated drought situations, the overall recharge to the buried valley is about 15 percent less and the total amount of streamflow in the model area is reduced by about 19 percent. Without pumping, infiltration to the buried valley aquifer around the confining unit decreased by a small amount (0.05 million gallons per day (Mgal/d)), and discharge to the streams increased by about 8 percent (0.3 Mgal/d). A 50-percent increase in pumping resulted in a simulated decrease in streamflow discharge of about 4 percent (0.14 Mgal/d). Streamflow depletion in Harvey Brook was evaluated by use of the numerical groundwater-flow model and an analytical model. The analytical model estimated negligible depletion from Harvey Brook under current (2009) pumping conditions, whereas the numerical model estimated that flow to Harvey Brook decreased 0.38 cubic feet per second (ft 3 /s) because of the pumping well withdrawals. A sensitivity analysis of the analytical model method showed that conducting a cursory evaluation using an analytical model of streamflow depletion using available information may result in a very wide range in results, depending on how well the hydraulic conductivity variables and aquifer geometry of the system are known, and how well the aquifer fits the assumptions of the model. Using the analytical model to evaluate the streamflow depletion with an incomplete understanding of the hydrologic system gave results that seem unlikely to reflect actual streamflow depletion in the Freeport aquifer study area. In contrast, the groundwater-flow model was a more robust method of evaluating the amount of streamflow depletion that results from withdrawals in the Freeport aquifer, and could be used to evaluate streamflow depletion in both streams. Simulations of streamflow without pumping for each measurement site were compared to the calibratedmodel streamflow (with pumping), the difference in the total being streamflow depletion. Simulations without pumping resulted in a simulated increase in the steady-state flow rate of 0.38 ft 3 /s in Harvey Brook and 0.01 ft 3 /s in Merrill Brook. This translates into a streamflow-depletion amount equal to about 8.5 percent of the steady-state base flow in Harvey Brook, and an unmeasurable amount of depletion in Merrill Brook. If pumping was increased by 50 percent and recharge reduced by 20 percent, the amount of streamflow depletion in Harvey Brook could reach 1.41 ft 3 /s.

Maine↗

Macropolygon morphology, development, and classification on North Panamint and Eureka playas, Death Valley National Park CA

Panamint and Eureka playas, both located within Death Valley National Park, exhibit a host of surficial features including fissures, pits, mounds, and plant-covered ridges, representing topographic highs and lows that vary up to 2 m of relief from the playa surface. Aerial photographs reveal that these linear strands often converge to form polygons, ranging in length from several meters to nearly a kilometer. These features stand out in generally dark contrast to the brighter intervening expanse of flat, plant-free, desiccated mud of the typical playa surface. Ground-truth mapping of playa features with differential GPS (Global Positioning System) was conducted in 1999 (North Panamint Valley) and 2002 (Eureka Valley). High-resolution digital maps reveal that both playas possess macropolygons of similar scale and geometry, and that fissures may be categorized into one of two genetic groups: (1) shore-parallel or playa-interior desiccation and shrinkage; and (2) tectonic-induced cracks. Early investigations of these features in Eureka Valley concluded that their origin may have been related to agricultural activity by paleo-Indian communities. Although human artifacts are abundant at each locale, there is no evidence to support the inference that surface features reported on Eureka Playa are anthropogenic in origin. Our assumptions into the genesis of polygons on playas is based on our fortuitous experience of witnessing a fissure in the process of formation on Panamint Playa after a flash flood (May 1999); our observations revealed a paradox that saturation of the upper playa crusts contributes to the establishment of some desiccation features. Follow-up visits to the same feature over 2 yrs' time are a foundation for insight into the evolution and possible longevity of these features. ?? 2005 Elsevier B.V. All rights reserved.

Earth-Science Reviews↗

Death Valley regional ground-water flow system, Nevada and California -- hydrogeologic framework and transient ground-water flow model

A numerical three-dimensional (3D) transient ground-water flow model of the Death Valley region was developed by the U.S. Geological Survey for the U.S. Department of Energy programs at the Nevada Test Site and at Yucca Mountain, Nevada. Decades of study of aspects of the ground-water flow system and previous less extensive ground-water flow models were incorporated and reevaluated together with new data to provide greater detail for the complex, digital model. A 3D digital hydrogeologic framework model (HFM) was developed from digital elevation models, geologic maps, borehole information, geologic and hydrogeologic cross sections, and other 3D models to represent the geometry of the hydrogeologic units (HGUs). Structural features, such as faults and fractures, that affect ground-water flow also were added. The HFM represents Precambrian and Paleozoic crystalline and sedimentary rocks, Mesozoic sedimentary rocks, Mesozoic to Cenozoic intrusive rocks, Cenozoic volcanic tuffs and lavas, and late Cenozoic sedimentary deposits of the Death Valley Regional Ground-Water Flow System (DVRFS) region in 27 HGUs. Information from a series of investigations was compiled to conceptualize and quantify hydrologic components of the ground-water flow system within the DVRFS model domain and to provide hydraulic-property and head-observation data used in the calibration of the transient-flow model. These studies reevaluated natural ground-water discharge occurring through evapotranspiration and spring flow; the history of ground-water pumping from 1913 through 1998; ground-water recharge simulated as net infiltration; model boundary inflows and outflows based on regional hydraulic gradients and water budgets of surrounding areas; hydraulic conductivity and its relation to depth; and water levels appropriate for regional simulation of prepumped and pumped conditions within the DVRFS model domain. Simulation results appropriate for the regional extent and scale of the model were provided by acquiring additional data, by reevaluating existing data using current technology and concepts, and by refining earlier interpretations to reflect the current understanding of the regional ground-water flow system. Ground-water flow in the Death Valley region is composed of several interconnected, complex ground-water flow systems. Ground-water flow occurs in three subregions in relatively shallow and localized flow paths that are superimposed on deeper, regional flow paths. Regional ground-water flow is predominantly through a thick Paleozoic carbonate rock sequence affected by complex geologic structures from regional faulting and fracturing that can enhance or impede flow. Spring flow and evapotranspiration (ET) are the dominant natural ground-water discharge processes. Ground water also is withdrawn for agricultural, commercial, and domestic uses. Ground-water flow in the DVRFS was simulated using MODFLOW-2000, a 3D finite-difference modular ground-water flow modeling code that incorporates a nonlinear least-squares regression technique to estimate aquifer parameters. The DVRFS model has 16 layers of defined thickness, a finite-difference grid consisting of 194 rows and 160 columns, and uniform cells 1,500 m on each side. Prepumping conditions (before 1913) were used as the initial conditions for the transient-state calibration. The model uses annual stress periods with discrete recharge and discharge components. Recharge occurs mostly from infiltration of precipitation and runoff on high mountain ranges and from a small amount of underflow from adjacent basins. Discharge occurs primarily through ET and spring discharge (both simulated as drains) and water withdrawal by pumping and, to a lesser amount, by underflow to adjacent basins, also simulated by drains. All parameter values estimated by the regression are reasonable and within the range of expected values. The simulated hydraulic heads of the final calibrated transient model gener

Scientific Investigations Report↗

Hydrogeology and water quality of the Leetown area, West Virginia

The U.S. Geological Survey’s Leetown Science Center and the co-located U.S. Department of Agriculture’s National Center for Cool and Cold Water Aquaculture both depend on large volumes of cold clean ground water to support research operations at their facilities. Currently, ground-water demands are provided by three springs and two standby production wells used to augment supplies during periods of low spring flow. Future expansion of research operations at the Leetown Science Center is dependent on assessing the availability and quality of water to the facilities and in locating prospective sites for additional wells to augment existing water supplies. The hydrogeology of the Leetown area, West Virginia, is a structurally complex karst aquifer. Although the aquifer is a karst system, it is not typical of most highly cavernous karst systems, but is dominated by broad areas of fractured rock drained by a relatively small number of solution conduits. Characterization of the aquifer by use of fluorometric tracer tests, a common approach in most karst terranes, therefore only partly defines the hydrogeologic setting of the area. In order to fully assess the hydrogeology and water quality in the vicinity of Leetown, a multi-disciplinary approach that included both fractured rock and karst research components was needed. The U.S. Geological Survey developed this multi-disciplinary research effort to include geologic, hydrologic, geophysical, geographic, water-quality, and microbiological investigations in order to fully characterize the hydrogeology and water quality of the Leetown area, West Virginia. Detailed geologic and karst mapping provided the framework on which hydrologic investigations were based. Fracture trace and lineament analysis helped locate potential water-bearing fractures and guided installation of monitoring wells. Monitoring wells were drilled for borehole geophysical surveys, water-quality sampling, water-level measurements, and aquifer tests to characterize the quality of water and the hydraulic properties of the aquifer. Surface geophysical surveys provided a 3-dimensional view of bedrock resistivity in order to assess geologic and lithologic controls on ground-water flow. Borehole geophysical surveys were conducted in monitoring wells to assess the storage and movement of water in subsurface fractures. Numerous single-well, multi-well, and straddle packer aquifer tests and step-drawdown tests were conducted to define the hydraulic properties of the aquifer and to assess the role of bedrock fractures and solution conduits in the flow of ground water. Water samples collected from wells and springs were analyzed to assess the current quality of ground water and provide a baseline for future assessment. Microbiological sampling of wells for indicator bacteria and human and animal DNA provided an analysis of agricultural and suburban development impacts on ground-water quality. Light detection and ranging (LiDAR) data were analyzed to develop digital elevation models (DEMs) for assessing sinkhole distribution, to provide elevation data for development of a ground-water flow model, and to assess the distribution of major fractures and faults in the Leetown area. The flow of ground water in the study area is controlled by lithology and geologic structure. Bedrock, especially low permeability units such as the shale Martinsburg Formation and the Conococheague Limestone, act as barriers to water flowing down gradient and across bedding. This retardation of cross-strike flow is especially pronounced in the Leetown area, where bedding typically dips at steep angles. Highly permeable fault and fracture zones that disrupt the rocks in cross-strike directions provide avenues through which ground water can flow laterally across or through strata of low primary permeability. Significant strike parallel thrust faults and cross-strike faults typically coincide with larger solution conduits and act as drains for the more pervasive network of interconnected diffuse fractures. Results of borehole geophysical surveys indicate that although numerous fractures may intersect a borehole, only one or two of the fractures typically transmit most of the water to a well. The diffuse-flow dominated network of fractures that provides the majority of storage occupies only a small proportion of the total aquifer volume but constitutes the majority of porosity within the aquifer. Solution conduits, while occupying a relatively small volume of the overall aquifer, are especially important because they serve as primary drains for the ground-water flow system. Surface resistivity maps and cross-sectionsshow anomalous areas of low resistivities coincident with the prevailing geologic strike at N. 20º E., with major cross-strike faults, and with major springs in the region. Transmissivity derived from straddle packer tests was highly variable, and ranged over three orders of magnitude (1.8 x 10 -6 to 5.9 x 10 -3 ft 2 /d) in diffuse-flow fractures. A similar large variability in transmissivity was documented by single- and multi-well aquifer tests conducted in conduit-flow dominated portions of the aquifer (2.0 x 10 3 to 1.4 x 10 4 ft 2 /d) in lowland areas immediately adjacent to the Leetown Science Center. A stream-gaging station installed on Hopewell Run near the point where the stream exits the Leetown watershed indicates average daily streamflow for the Hopewell Run of approximately 11.2 ft 3 /s, and ranged from a minimum of 1.80 ft 3 /s on September 28, 2005, to a maximum of 73.0 ft 3 /s on December 11, 2003. Base-flow (ground-water) discharge surveys identified numerous small seeps adjacent to streams in the area. Hydrographs of the stage of Balch Spring show rapid response to individual storms. Strong correlation of the flow of Hopewell Run and Balch Spring indicates the nearby losing stream reach is partly responsible for higher fluctuations in the stage of Balch Spring. A water budget for the study period (2003-2005), based on measured precipitation and hydrograph analyses, is expressed as Precipitation (38.60 in/yr) = Surface Runoff (1.36 in/yr) + Ground-Water Discharge (17.73 in/yr) + Evapotranspiration (24.23 in/yr) – Change in storage (4.72 in/yr). Flow of ground water through the epikarst, a shallow zone of intensely weathered rock and regolith, can be rapid (on the order of days or weeks) as flow is concentrated in solution conduits. Flow within the intermediate and deeper zones is typically much slower. Eight dye-tracer tests conducted in the Leetown area found ground-water flow patterns to be divergent, with velocities ranging from about 12.5 to 610 ft/day and a median velocity of 50 ft/day. Estimates of ground-water age in carbonate rocks in the region are on the order of 15 years in the shallower portions of the aquifer to 50 years or older for deeper portions of the aquifer. Shallow springs can have a significant component of fairly young water (< 5 years in age). Ground-water samples collected from 16 sites (12 wells and 4 springs) in the Leetown area were analyzed for more than 340 constituents. Only turbidity, indicator bacteria, and radon were typically present in concentrations exceeding U.S. Environmental Protection Agency (USEPA) drinking-water or aquatic life standards.

West Virginia↗

Dynamic modeling of Tampa Bay urban development using parallel computing

Urban land use and land cover has changed significantly in the environs of Tampa Bay, Florida, over the past 50 years. Extensive urbanization has created substantial change to the region's landscape and ecosystems. This paper uses a dynamic urban-growth model, SLEUTH, which applies six geospatial data themes (slope, land use, exclusion, urban extent, transportation, hillside), to study the process of urbanization and associated land use and land cover change in the Tampa Bay area. To reduce processing time and complete the modeling process within an acceptable period, the model is recoded and ported to a Beowulf cluster. The parallel-processing computer system accomplishes the massive amount of computation the modeling simulation requires. SLEUTH calibration process for the Tampa Bay urban growth simulation spends only 10 h CPU time. The model predicts future land use/cover change trends for Tampa Bay from 1992 to 2025. Urban extent is predicted to double in the Tampa Bay watershed between 1992 and 2025. Results show an upward trend of urbanization at the expense of a decline of 58% and 80% in agriculture and forested lands, respectively.

Computers & Geosciences↗

POLARIS: A 30-meter probabilistic soil series map of the contiguous United States

A new complete map of soil series probabilities has been produced for the contiguous United States at a 30 m spatial resolution. This innovative database, named POLARIS, is constructed using available high-resolution geospatial environmental data and a state-of-the-art machine learning algorithm (DSMART-HPC) to remap the Soil Survey Geographic (SSURGO) database. This 9 billion grid cell database is possible using available high performance computing resources. POLARIS provides a spatially continuous, internally consistent, quantitative prediction of soil series. It offers potential solutions to the primary weaknesses in SSURGO: 1) unmapped areas are gap-filled using survey data from the surrounding regions, 2) the artificial discontinuities at political boundaries are removed, and 3) the use of high resolution environmental covariate data leads to a spatial disaggregation of the coarse polygons. The geospatial environmental covariates that have the largest role in assembling POLARIS over the contiguous United States (CONUS) are fine-scale (30 m) elevation data and coarse-scale (~ 2 km) estimates of the geographic distribution of uranium, thorium, and potassium. A preliminary validation of POLARIS using the NRCS National Soil Information System (NASIS) database shows variable performance over CONUS. In general, the best performance is obtained at grid cells where DSMART-HPC is most able to reduce the chance of misclassification. The important role of environmental covariates in limiting prediction uncertainty suggests including additional covariates is pivotal to improving POLARIS' accuracy. This database has the potential to improve the modeling of biogeochemical, water, and energy cycles in environmental models; enhance availability of data for precision agriculture; and assist hydrologic monitoring and forecasting to ensure food and water security.

Geoderma↗

Soil phosphorus cycling in tropical soils: An ultisol and oxisol perspective

Phosphorus (P) is essential for life. It is the backbone of our DNA, provides energy for biological reactions, and is an integral component of cell membranes. As such, it is no surprise that P availability plays a strong role in regulating ecosystem structure and function (Wassen et al. 2005, Elser et al. 2007, Condit et al. 2013), and in determining our capacity to grow food for a burgeoning human population (Sharpley et al. 1997, Sims and Sharpley 2005, Lal 2009). Concerns that P supplies are insufficient to meet our species’ growing demands are on the rise (Richardson and Simpson 2011) and scientific and media outlets increasingly discuss P as an element worthy of our attention and concern (e.g., Cordell et al. 2009, Lougheed 2011, Edixhoven et al. 2013, Ulrich et al. 2013). Indeed, a number of groups are calling for the explicit stewardship of our planet’s P stocks (Schipper 2014, Withers et al. 2015). Yet a focus on P as a vital and limited resource is not new in the tropics, where an abundance of soils characterized by low P has resulted in a substantial, longstanding reliance on P inputs for tropical ecosystem function in both unmanaged and agriculture settings (Table 1, Figure 2; Sanchez 1976, Swap et al. 1992, Chadwick et al. 1999, Okin et al. 2004, Lal 2009). Indeed, there is a long history of cultivation in the tropics, where for thousands of years land management practices have included methods that effectively modify P availability for plant growth (e.g., Giardina et al. 2000, Lawrence and Schlesinger 2001, Vitousek et al. 2004, Lewis et al. 2015). Nevertheless, low soil fertility in tropical systems where fertilizer is scarce has enduringly been recognized as a major source of hunger and starvation (Sanchez and Buol 1975, Sanchez 2002, Sanchez and Swaminathan 2005).

Book chapter↗

Inventory of glaciers and perennial snowfields of the conterminous USA

This report summarizes an updated inventory of glaciers and perennial snowfields of the conterminous United States. The inventory is based on interpretation of mostly aerial imagery provided by the National Agricultural I magery Program, US Department of Agriculture, with some satellite imagery in places where aerial imagery was not suitable. The inventory includes all perennial snow and ice features ≥ 0.01 km 2 . Due to aerial survey schedules and seasonal snow cover, imageries acquired over a number of years were required. The earliest date is 2013 and the latest is 2020, but more than 73 % of the outlines were acquired from 2015 imagery. The inventory is compiled as shapefiles within a geographic information system that includes feature classification, area, and location. The inventory identified 1331 (366.52 ± 14.34 km 2 ) glaciers, 1176 (31.01 ± 9.30 km 2 ) perennial snowfields, and 35 (3.57 km 2 ± no uncertainty) buried-ice features. The data including both the shapefiles and tabulated results are publicly available at https://doi.org/10.15760/geology-data.03 (Fountain and Glenn, 2022).

Earth System Science Data↗

The statistical power to detect cross-scale interactions at macroscales

Macroscale studies of ecological phenomena are increasingly common because stressors such as climate and land-use change operate at large spatial and temporal scales. Cross-scale interactions (CSIs), where ecological processes operating at one spatial or temporal scale interact with processes operating at another scale, have been documented in a variety of ecosystems and contribute to complex system dynamics. However, studies investigating CSIs are often dependent on compiling multiple data sets from different sources to create multithematic, multiscaled data sets, which results in structurally complex, and sometimes incomplete data sets. The statistical power to detect CSIs needs to be evaluated because of their importance and the challenge of quantifying CSIs using data sets with complex structures and missing observations. We studied this problem using a spatially hierarchical model that measures CSIs between regional agriculture and its effects on the relationship between lake nutrients and lake productivity. We used an existing large multithematic, multiscaled database, LAke multiscaled GeOSpatial, and temporal database (LAGOS), to parameterize the power analysis simulations. We found that the power to detect CSIs was more strongly related to the number of regions in the study rather than the number of lakes nested within each region. CSI power analyses will not only help ecologists design large-scale studies aimed at detecting CSIs, but will also focus attention on CSI effect sizes and the degree to which they are ecologically relevant and detectable with large data sets.

Ecosphere↗

Development of a benchmark eddy flux evapotranspiration dataset for evaluation of satellite-driven evapotranspiration models over the CONUS

A large sample of ground-based evapotranspiration (ET) measurements made in the United States, primarily from eddy covariance systems, were post-processed to produce a benchmark ET dataset. The dataset was produced primarily to support the intercomparison and evaluation of the OpenET satellite-based remote sensing ET (RSET) models and could also be used to evaluate ET data from other models and approaches. OpenET is a web-based service that makes field-delineated and pixel-level ET estimates from well-established RSET models readily available to water managers, agricultural producers, and the public. The benchmark dataset is composed of flux and meteorological data from a variety of providers covering native vegetation and agricultural settings. Flux footprint predictions were developed for each station and included static flux footprints developed based on average wind direction and speed, as well as dynamic hourly footprints that were generated with a physically based model of upwind source area. The two footprint prediction methods were rigorously compared to evaluate their relative spatial coverage. Data from all sources were post-processed in a consistent and reproducible manner including data handling, gap-filling, temporal aggregation, and energy balance closure correction. The resulting dataset included 243,048 daily and 5,284 monthly ET values from 194 stations, with all data falling between 1995 and 2021. We assessed average daily energy imbalance using 172 flux sites with a total of 193,021 days of data, finding that overall turbulent fluxes were understated by about 12% on average relative to available energy. Multiple linear regression analyses indicated that daily average latent energy flux may be typically understated slightly more than sensible heat flux. This dataset was developed to provide a consistent reference to support evaluation of RSET data being developed for a wide range of applications related to water accounting and water resources management at field to watershed scales.

continental United States↗

Sediment loads in the Moore Creek drainage‐basin, Idaho 1939–40

The Boise River Project in southwestern Idaho comprises an area of about 333,000 acres of highly developed agricultural land. Precipitation in the irrigated valley averages about ten inches a year which is too low to support any but desert vegetation. Water for irrigation during the growing season is obtained from Arrowrock Reservoir located on the Boise River four miles above Moore Creek, and from Deer Flat Reservoir near Nampa. The deposition of sediment in the river‐channel, in irrigation‐canals, and in farm‐laterals has interfered with the normal operation of the irrigation‐system and has occasionally required the expenditure of considerable time and money for removal of the deposited material.

Idaho↗

Evaluation of chemical, biological, and physical conditions in the Winter Haven chain of lakes, Florida, March-June 1976

Reconnaissance of water-quality conditions of 14 interconnected navigable lakes, in and around Winter Haven, Fla., revealed that in March and May, 1976 most were eutrophic, on the basis of high nutrient (nitrogen and phosphorus) concentrations. Lakes Lulu and Shipp were the most enriched as a result of surface runoff from residential, agricultural, and highly urbanized areas, and many years of municipal and industrial waste effluent input. Phytoplankton counts were greater than a million cells per milliliter in some lakes sampled; algal blooms have ocurred, and water clarity was low. The level of Lake Howard fell to the lowest stage recorded in 31 years during May 1976. The record low was likely due to rainfall deficiency. Leakage of water through the lake beds to the ground-water system is also possible, but determination of the escaping water volume would require additional study. (Woodard-USGS)

Florida↗

Groundwater flow model for the Little Plover River basin in Wisconsin’s Central Sands

The Little Plover River is a groundwater-fed stream in the sand plains region of central Wisconsin. In this region, sandy sediment deposited during or soon after the last glaciation forms an important unconfined sand and gravel aquifer. This aquifer supplies water for numerous high-capacity irrigation, municipal, and industrial wells that support a thriving agricultural industry. In recent years the addition of many new wells, combined with observed diminished flows in the Little Plover and other nearby rivers, has raised concerns about the impacts of the wells on groundwater levels and on water levels and flows in nearby lakes, streams, and wetlands. Diverse stakeholder groups, including well operators, Growers, environmentalists, local land owners, and regulatory and government officials have sought a better understanding of the local groundwater-surface water system and have a shared desire to balance the water needs of the he liagricultural, industrial, and urban users with the maintenance and protection of groundwater-dependent natural resources. To help address these issues, the Wisconsin Department of Natural Resources requested that the Wisconsin Geological and Natural History Survey and U.S. Geological Survey cooperatively develop a groundwater flow model that could be used to demonstrate the relationships among groundwater, surface water, and well withdrawals and also be a tool for testing and evaluating alternative water management strategies for the central sands region. Because of an abundance of previous studies, data availability, local interest, and existing regulatory constraints the model focuses on the Little Plover River watershed, but the modeling methodology developed during this study can apply to much of the larger central sands of Wisconsin. The Little Plover River groundwater flow model simulates three-dimensional groundwater movement in and around the Little Plover River basin under steady-state and transient conditions. This model explicitly includes all high-capacity wells in the model domain and simulates seasonal variations in recharge and well pumping. The model represents the Little Plover River, and other significant streams and drainage ditches in the model domain, as fully connected to the groundwater system, computes stream base flow resulting from groundwater discharge, and routes the flow along the stream channel. A separate soil-water-balance (SWB) model was used to develop groundwater recharge arrays as input for the groundwater flow model. The SWB model uses topography, soils, land use, and climatic data to estimate recharge as deep drainage from the soil zone. The SWB model explicitly includes recharge originating as irrigation water, and computes irrigation using techniques similar to those used by local irrigation operators. The groundwater flow model uses the U.S. Geological Survey’s MODFLOW modeling code which is freely available, widely accepted, and commonly used by the groundwater community. The groundwater flow model and the SWB model use identical high-resolution numerical grids having model cells 100 feet on a side, with physical properties assigned to each grid cell. This grid allows accurate geographic placement of wells, streams, and other model features. The 3-dimensional grid has three layers; layers 1 and 2 represent the sand and gravel aquifer and layer 3 represents the underlying sandstone. The distribution of material properties in the model (hydraulic conductivity, aquifer thickness, etc.) comes from previous published geologic studies of the region, updated by calibration to recent streamflow and groundwater level data. The SWB model operates on a daily time step. The groundwater flow model was calibrated to monthly stress periods with time steps ranging from 1 to 16 days. More detailed time discretization is possible. The groundwater model was calibrated to water-level and streamflow data collected during 2013 and 2014 by adjusting model parameters (primarily hydraulic conductivity, storage, and recharge) until the model produced a conditionally optimal fit between field observations and model output, subject to consistency with previously published geologic studies. Calibration was performed under both steady and transient conditions, and used a sophisticated parameter-estimation procedure (PEST) for the calibration process and to identify important model parameters. For the Little Plover River, the two most important parameters are the global recharge multiplier and the hydraulic conductivity of the stream bed. The calibrated model produces water-level and mass-balance results that are consistent with field observations and previous studies of the area. The completed model is a powerful tool for testing and demonstrating alternative water-management scenarios. Example model applications described in this report include simulating how the cumulative impacts of pumping and land-use change have affected average baseflow in the Little Plover River. Depletion-potential mapping represents a method for predicting which wells and well locations have the greatest impact on nearby surface-water resources. The completed model is publicly available, along with a companion user’s guide to assist with its operation, at http://wgnhs.org/littleplover- river-groundwater-model.

Wisconsin↗

Fungicides: An overlooked pesticide class?

Fungicides are indispensable to global food security and their use is forecasted to intensify. Fungicides can reach aquatic ecosystems and occur in surface water bodies in agricultural catchments throughout the whole growing season due to their frequent, prophylactic application. However, in comparison to herbicides and insecticides, the exposure to and effects of fungicides have received less attention. We provide an overview of the risk of fungicides to aquatic ecosystems covering fungicide exposure (i.e., environmental fate, exposure modelling, and mitigation measures) as well as direct and indirect effects of fungicides on microorganisms, macrophytes, invertebrates, and vertebrates. We show that fungicides occur widely in aquatic systems, that the accuracy of predicted environmental concentrations is debatable, and that fungicide exposure can be effectively mitigated. We additionally demonstrate that fungicides can be highly toxic to a broad range of organisms and can pose a risk to aquatic biota. Finally, we outline central research gaps that currently challenge our ability to predict fungicide exposure and effects, promising research avenues, and shortcomings of the current environmental risk assessment for fungicides.

Environmental Science & Technology↗

Beyond the edge: Linking agricultural landscapes, stream networks, and best management practices

Despite much research and investment into understanding and managing nutrients across agricultural landscapes, nutrient runoff to freshwater ecosystems is still a major concern. We argue there is currently a disconnect between the management of watershed surfaces (agricultural landscape) and river networks (riverine landscape). These landscapes are commonly managed separately, but there is limited cohesiveness between agricultural landscape-focused research and river science, despite similar end goals. Interdisciplinary research into stream networks that drain agricultural landscapes is expanding but is fraught with problems. Conceptual frameworks are useful tools to order phenomena, reveal patterns and processes, and in interdisciplinary river science, enable the joining of multiple areas of understanding into a single conceptual–empirical structure. We present a framework for the interdisciplinary study and management of agricultural and riverine landscapes. The framework includes components of an ecosystems approach to the study of catchment–stream networks, resilience thinking, and strategic adaptive management. Application of the framework is illustrated through a study of the Fox Basin in Wisconsin, USA. To fully realize the goal of nutrient reduction in the basin, we suggest that greater emphasis is needed on where best management practices (BMPs) are used within the spatial context of the combined watershed–stream network system, including BMPs within the river channel. Targeted placement of BMPs throughout the riverine landscape would increase the overall buffering capacity of the system to nutrient runoff and thus its resilience to current and future disturbances.

Journal of Environmental Quality↗