SEARCH · Search USGS

Results for “Agricultural Systems Journal”

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.

137 records · Page 8Linked to original sources

USGS2006

Nest survival of clay-colored and vesper sparrows in relation to woodland edge in mixed-grass prairies

The quantity and quality of northern mixed-grass prairie continues to decline because of conversion to agriculture, invasion of woody and exotic plants, and disruption of important ecological processes that shape grasslands. Declines in grassland bird populations in North Dakota, USA, have coincided with these largely anthropogenic alterations to prairie habitat. In grasslands of north-central and northwestern North Dakota, woody plants have increased due primarily to fire suppression, extirpation of bison (Bos bison), and widescale planting of tree shelter belts. In northern grasslands, effects of woody vegetation on survival of grassland birds are poorly understood, and conclusions are based mainly on studies conducted outside the region. We examined nest survival of clay-colored sparrows (Spizella pallida) and vesper sparrows (Pooecetes gramineus) relative to the distance nests were located from aspen (Populus tremuloides,) woodland edges and relative to other habitat features near the nest. Clay-colored and vesper sparrow nest survival was higher for nests located near woodland edges, nests with greater cover of Kentucky bluegrass (Poa pratensis), and nests more concealed by vegetation. Vesper sparrow nest survival increased as the percent cover of tall shrubs near the nest increased. Based on video-camera data, the 13-lined ground squirrel (Spermophilus tridecemlineatus,) was the most common predator of sparrow eggs and young. Thirteen-lined ground squirrels were more common far from woodland edges than near, and this pattern may, in part, explain clay-colored and vesper sparrow nest survival in relation to woodland edges. In contrast to our results, studies conducted in other grassland systems generally report lower nest survival for grassland birds nesting near trees and shrubs. This disparity in results demonstrates the need to identify specific nest predators and their distributions with respect to important habitat features because these data can be important in explaining-and perhaps predicting-patterns of nest predation.

Journal of Wildlife Management↗

USGS2014

Analysis of projected water availability with current basin management plan, Pajaro Valley, California

The projection and analysis of the Pajaro Valley Hydrologic Model (PVHM) 34 years into the future using MODFLOW with the Farm Process (MF-FMP) facilitates assessment of potential future water availability. The projection is facilitated by the integrated hydrologic model, MF-FMP that fully couples the simulation of the use and movement of water from precipitation, streamflow, runoff, groundwater flow, and consumption by natural and agricultural vegetation throughout the hydrologic system at all times. MF-FMP allows for more complete analysis of conjunctive-use water-resource systems than previously possible with MODFLOW by combining relevant aspects of the landscape with the groundwater and surface-water components. This analysis is accomplished using distributed cell-by-cell supply-constrained and demand-driven components across the landscape within “water-balance subregions” (WBS) comprised of one or more model cells that can represent a single farm, a group of farms, watersheds, or other hydrologic or geopolitical entities. Analysis of conjunctive use would be difficult without embedding the fully coupled supply-and-demand into a fully coupled simulation, and are difficult to estimate a priori. The analysis of projected supply and demand for the Pajaro Valley indicate that the current water supply facilities constructed to provide alternative local sources of supplemental water to replace coastal groundwater pumpage, but may not completely eliminate additional overdraft. The simulation of the coastal distribution system (CDS) replicates: 20 miles of conveyance pipeline, managed aquifer recharge and recovery (MARR) system that captures local runoff, and recycled-water treatment facility (RWF) from urban wastewater, along with the use of other blend water supplies, provide partial relief and substitution for coastal pumpage (aka in-lieu recharge). The effects of these Basin Management Plan (BMP) projects were analyzed subject to historical climate variations and assumptions of 2009 urban water demand and land use. Water supplied directly from precipitation, and indirectly from reuse, captured local runoff, and groundwater is necessary but inadequate to satisfy agricultural demand without coastal and regional storage depletion that facilitates seawater intrusion. These facilities reduce potential seawater intrusion by about 45% with groundwater levels in the four regions served by the CDS projected to recover to levels a few feet above sea level. The projected recoveries are not high enough to prevent additional seawater intrusion during dry-year periods or in the deeper aquifers where pumpage is greater. While these facilities could reduce coastal pumpage by about 55% of the historical 2000–2009 pumpage for these regions, and some of the water is delivered in excess of demand, other coastal regions continue to create demands on coastal pumpage that will need to be replaced to reduce seawater intrusion. In addition, inland urban and agricultural demands continue to sustain water levels below sea level causing regional landward gradients that also drive seawater intrusion. Seawater intrusion is reduced by about 45% but it supplies about 55% of the recovery of groundwater levels in the coastal regions served by the CDS. If economically feasible, water from summer agricultural runoff and tile-drain returnflows could be another potential local source of water that, if captured and reused, could offset the imbalance between supply and demand as well as reducing discharge of agricultural runoff into the National Marine Sanctuary of Monterey Bay. A BMP update (2012) identifies projects and programs that will fund a conservation program and will provide additional, alternative water sources to reduce or replace coastal and inland pumpage, and to replenish the aquifers with managed aquifer recharge in an inland portion of the Pajaro Valley.

California↗

USGS2026

What is the (real) rate of soil health practice adoption? Making sense of three data sources

Conservation stakeholders looking to quantify the impact of their investments to increase soil health practice adoption over time often face challenges in interpreting practice adoption data due to discrepancies in language and results among data sources. Similarly, efforts to estimate environmental outcomes of practice adoption, such as water quality and greenhouse gas emissions, can vary depending on different practice adoption input data. To help make sense of different adoption data sources, we compared county-level adoption data for winter cover crops (WCC), no-till (NT), and reduced tillage (RT) in three areas of the United States with contrasting climates and production systems: central Illinois (CIL), southern Illinois (SIL), and western New York (WNY). We analyzed data available during 2015 through 2022 from the Operational Tillage Information System (OpTIS, remote sensing), US Census of Agriculture (AgCensus, a farmer survey), and, specifically in Illinois, the Illinois Soil Conservation Transect Survey (Transect, a roadside survey). The magnitude of differences between the datasets depended on the practice and geographic location. For example, OpTIS and AgCensus tillage data were much more similar in Illinois (average difference of less than 4 percentage points) compared to New York (average differences of 20 percentage points). Similarly, there was less variability and smaller differences between OpTIS and AgCensus WCC data in Illinois compared to WNY. AgCensus tended to report lower WCC adoption for Illinois and greater adoption in WNY compared to OpTIS. All data sources agreed that the rate of change in tillage practices is slow (mainly –1% to 1%) and that adoption of WCC is low (assuming linear growth, it could take nearly a century to reach 50% WCC adoption in CIL). Differences among the datasets were attributed to definitional inconsistencies for RT and NT and how WCC data were acquired. For example, the AgCensus asks if a WCC was planted, whereas OpTIS and Transect evaluate the presence of a standing WCC. Data sources also reflect different time periods (calendar years or crop years) and types of cropland assessed (corn [ Zea mays L.], soybean [ Glycine max {L.} Merr.], or all cropland). We propose two recommendations to improve interpretation and consistency: (1) a working group to harmonize definitions and protocols and develop educational materials for data users, and (2) a research effort that integrates different adoption data types and produces publicly available adoption data at HUC-10 and county scales. Such activities could help improve data access and utility for evidence-based conservation decision-making and enhance the accuracy of environmental models that rely on adoption data as input.

Journal of Soil and Water Conservation↗

USGS2023

Marshbird response to herbicide control of cattail in northwestern Minnesota

Wetlands provide essential habitat for a wide variety of wildlife species. In the once wetland-rich Prairie Pothole Region and adjacent areas of central North America, many wetlands have been converted to agricultural production. Many remaining wetlands experience ecological change via the invasion and spread of non-native plant species, such as non-native narrowleaf ( Typha angustifolia ) and hybrid cattail ( Typha x glauca ), which spread aggressively and displace native vegetation, especially in large, impounded wetlands. Management of wetlands in these landscapes often includes broad-scale herbicide application intended to break up mats of cattail and restore areas to more wildlife-friendly conditions. Although restoration of wildlife habitat is a common goal of such management, marshbird response to invasive cattail control is poorly understood. To evaluate the effects of cattail management on wetland wildlife, we conducted standardized call-broadcast surveys for 5 species of marshbirds at 9 study sites that included survey locations associated with areas treated with herbicide and paired areas not treated with herbicide in wetland impoundments in northwestern Minnesota, USA, using a before-after, control-impact study design. We surveyed American bitterns ( Botaurus lentiginosus ), least bitterns ( Ixobrychus exilis ), pied-billed grebes ( Podilymbus podiceps ), soras ( Porzana carolina ), and Virginia rails ( Rallus limicola ) during the breeding season prior to herbicide application (late summer and early autumn of 2015) and during the 3 breeding seasons after herbicide application (2016–2018). We modeled species counts using a generalized linear mixed model with year-by-treatment interactions as fixed effects and site as a random effect. Before herbicide application, expected mean counts did not differ between treatment and control survey locations. Three years post-treatment, we detected significant increases in expected mean counts at treatment compared to control survey locations for soras ( t 193 = −3.373, P = 0.020) and Virginia rails ( t 193 = −3.167, P = 0.037), and point estimates for all species except least bittern were higher at treatment survey locations. Overall, our results suggest that these marshbird species responded positively to herbicide control of invasive cattail and that breeding marshbirds in these and similar wetland systems may experience positive population response over a period of at least 3 years following treatment.

Minnesota↗

USGS2016

The origin of shallow lakes in the Khorezm Province, Uzbekistan, and the history of pesticide use around these lakes

The economy of the Khorezm Province in Uzbekistan relies on the large-scale agricultural production of cotton. To sustain their staple crop, water from the Amu Darya is diverted for irrigation through canal systems constructed during the early to mid-twentieth century when this region was part of the Soviet Union. These diversions severely reduce river flow to the Aral Sea. The Province has >400 small shallow (<3 m deep) lakes that may have originated because of this intensive irrigation. Sediment cores were collected from 12 lakes to elucidate their origin because this knowledge is critical to understanding water use in Khorezm. Core chronological data indicate that the majority of the lakes investigated are less than 150 years old, which supports a recent origin of the lakes. The thickness of lacustrine sediments in the cores analyzed ranged from 20 to 60 cm in all but two of the lakes, indicating a relatively slow sedimentation rate and a relatively short-term history for the lakes. Hydrologic changes in the lakes are evident from loss on ignition and pollen analyses of a subset of the lake cores. The data indicate that the lakes have transitioned from a dry, saline, arid landscape during pre-lake conditions (low organic carbon content) and low pollen concentrations (in the basal sediments) to the current freshwater lakes (high organic content), with abundant freshwater pollen taxa over the last 50–70 years. Sediments at the base of the cores contain pollen taxa dominated by Chenopodiaceae and Tamarix , indicating that the vegetation growing nearby was tolerant to arid saline conditions. The near surface sediments of the cores are dominated by Typha/Sparganium , which indicate freshwater conditions. Increases in pollen of weeds and crop plants indicate an intensification of agricultural activities since the 1950s in the watersheds of the lakes analyzed. Pesticide profiles of DDT (dichlorodiphenyltrichloroethane) and its degradates and γ-HCH (gamma-hexachlorocyclohexane), which were used during the Soviet era, show peak concentrations in the top 10 cm of some of the cores, where estimated ages of the sediments (1950–1990) are associated with peak pesticide use during the Soviet era. These data indicate that the lakes are relatively young (mostly <150 years old) and that without irrigation and canal inputs from the Amu Darya, the lakes would not exist as freshwater lakes.

Khorezm Province↗

USGS2013

Appraising options to reduce shallow groundwater tables and enhance flow conditions over regional scales in an irrigated alluvial aquifer system

Some of the world’s key agricultural production systems face big challenges to both water quantity and quality due to shallow groundwater that results from long-term intensive irrigation, namely waterlogging and salinity, water losses, and environmental problems. This paper focuses on water quantity issues, presenting finite-difference groundwater models developed to describe shallow water table levels, non-beneficial groundwater consumptive use, and return flows to streams across two regions within an irrigated alluvial river valley in southeastern Colorado, USA. The models are calibrated and applied to simulate current baseline conditions in the alluvial aquifer system and to examine actions for potentially improving these conditions. The models provide a detailed description of regional-scale subsurface unsaturated and saturated flow processes, thereby enabling detailed spatiotemporal description of groundwater levels, recharge to infiltration ratios, partitioning of ET originating from the unsaturated and saturated zones, and groundwater flows, among other variables. Hybrid automated and manual calibration of the models is achieved using extensive observations of groundwater hydraulic head, groundwater return flow to streams, aquifer stratigraphy, canal seepage, total evapotranspiration, the portion of evapotranspiration supplied by upflux from the shallow water table, and irrigation flows. Baseline results from the two regional-scale models are compared to model predictions under variations of four alternative management schemes: (1) reduced seepage from earthen canals, (2) reduced irrigation applications, (3) rotational lease fallowing (irrigation water leased to municipalities, resulting in temporary dry-up of fields), and (4) combinations of these. The potential for increasing the average water table depth by up to 1.1 and 0.7 m in the two respective modeled regions, thereby reducing the threat of waterlogging and lowering non-beneficial consumptive use from adjacent fallow and naturally-vegetated lands, is demonstrated for the alternative management intervention scenarios considered. Net annual average savings of up to about 9.9 million m 3 (8000 ac ft) and 2.3 million m 3 (1900 ac ft) of non-beneficial groundwater consumptive use is demonstrated for the study periods in each of the two respective study regions. Alternative water management interventions achieve varying degrees of benefits in each of the two regions, suggesting a need to adopt region-specific interventions and avoid a ‘one-size-fits-all’ approach. Impacts of the considered interventions on return flows to the river were predicted to be significant, highlighting the need for flow augmentation to comply with an interstate river compact and portending beneficial impacts on solute loading.

Colorado↗

USGS2022

Upper Rio Grande Basin water-resource status and trends: Focus area study review and synthesis

The Upper Rio Grande Basin (URGB) is a critical international water resource under pressure from a myriad of climatic, ecological, infrastructural, water-use, and legal constraints. The objective of this study is to provide a comprehensive assessment of the spatial distribution and temporal trends of selected water-budget components (snow processes, evapotranspiration (ET), streamflow processes, and groundwater storage) using integrated analyses, such as watershed modeling and water availability and use data in the URGB over the past three decades. A spatially distributed snow evolution modeling system simulated snowpack processes over 34 years (1984–2017). It highlighted snow water equivalent declines from -35 to -77 mm/decade with widespread variability across elevation zones and land cover types. Gridded actual ET data from the SSEBop model were developed and tested for the URGB and demonstrated that all land-cover types had significant decreasing trends (1986-2015) ranging from -14 to -80 mm/decade. Conductivity-mass-balance (CMB) hydrograph separation results found that baseflow forms a large component of total streamflow, ranging from 29 to 69% (49% average) of total streamflow at 17 URGB sites upstream of Albuquerque, NM. Three of 4 graphical hydrograph separation methods in the U.S. Geological Survey Groundwater Toolbox were found to be inappropriate for estimating baseflow in the URGB; the most promising method, baseflow index (BFI) Standard, was optimized using CMB data and tested at three URGB sites, with resulting overestimation of 0 to 47%. Simulated changes in groundwater storage were extracted from historical and recent groundwater-flow models of select alluvial basins (San Luis, Española, Middle Rio Grande, and Tularosa-Hueco). In general, decreases in groundwater storage were observed from 1903 to 2013 except for the San Luis alluvial basin (Colorado), where periods of recovery are observed. The PRMS hydrologic model was successfully calibrated for 9 near-native subbasins (Nash-Sutcliffe efficiency 0.47 to 0.85) and parameters translated to the remaining subbasins; compared to simulated near-native flows (with minimal influence of reservoirs or diversions), observed Rio Grande streamgage flows demonstrated reductions of 40% or more for New Mexico and Texas areas of the basin. Significant decreasing trends (1980-2015) in precipitation, snowmelt rate, streamflow, and baseflow were observed at many of the 12 streamgage basins studied, which suggests that the decreasing trends for actual ET may be related to overall decreasing water availability in the basin, with negative implications for agricultural production and groundwater abstraction. Water security concerns arise from our findings of higher fraction precipitation as rain, slower snowmelt rates leading to decreasing streamflow production, and an increasing fraction of baseflow, all of which will affect the timing and magnitude of water available for human needs in the basin.

Upper Rio Grande Basin↗

USGS2018

Demographic consequences of conservation reserve program grasslands for lesser prairie‐chickens

Knowledge of landscape and regional circumstances where conservation programs are successful on working lands in agricultural production are needed. Converting marginal croplands to grasslands using conservation programs such as the United States Department of Agriculture Conservation Reserve Program (CRP) should be beneficial for many grassland‐obligate wildlife species; however, addition of CRP grasslands may result in different population effects based on regional climate, characteristics of the surrounding landscape, or species planted or established. Within landscapes occupied by lesser prairie‐chickens (Tympanuchus pallidicinctus), CRP may provide habitat only for specific life stages and habitat selection for CRP may vary between wet and dry years. Among all study sites, we captured and fitted 280 female lesser prairie‐chickens with very high frequency (VHF)‐ and global positioning system (GPS) transmitters during the spring lekking seasons of 2013–2015 to monitor habitat selection for CRP in regions of varying climate. We also estimated vital rates and habitat selection for 148 individuals, using sites in northwest Kansas, USA. The greatest ecological services of CRP became apparent when examining habitat selection and densities. Nest densities were approximately 3 times greater in CRP grasslands than native working grasslands (i.e., grazed), demonstrating a population‐level benefit (CRP = 6.0 nests/10 km2 ± 1.29 [SE], native working grassland = 1.7 nests/10 km2 ± 0.62). However, CRP supporting high nest density did not provide brood habitat; 85% of females with broods surviving to 7 days moved their young to other cover types. Regression analyses indicated lesser prairie‐chickens were approximately 8 times more likely to use CRP when 5,000‐ha landscapes were 70% rather than 20% grassland, indicating variation in the level of ecological services provided by CRP was dependent upon composition of the larger landscape. Further, CRP grasslands were 1.7 times more likely to be used by lesser prairie‐chickens in regions receiving 40 cm compared to 70 cm of average annual precipitation and during years of greater drought intensity. Demographic and resource selection analyses revealed that establishing CRP grasslands in northwest Kansas can increase the amount nesting habitat in a region where it may have previously been limited, thereby providing refugia to sustain populations through periods of extreme drought. Nest survival, adult survival during breeding, and nonbreeding season survival did not vary between lesser prairie‐chickens that used and did not use CRP grasslands. The finite rate of population growth was also similar for birds using CRP and using only native working grasslands, suggesting that CRP provides habitat similar to that of native working grassland in this region. Overall, lesser prairie‐chickens may thrive in landscapes that are a mosaic of native working grassland, CRP grassland, with a minimal amount of cropland, particularly when nesting and brood habitat are in close proximity.

Journal of Wildlife Management↗

USGS2025

Component identification of solid biomass fuels using reflected light microscopy: Interlaboratory study 2

As nations transition toward sustainable energy systems , biomass has become a vital component of global energy portfolios. Derived from organic materials such as wood, agricultural residues, forestry byproducts, and organic waste, biomass is a renewable energy source with significant environmental and economic benefits. Responsible biomass energy production can improve waste management, reduce emissions of greenhouse gases, and mitigate environmental pollution. However, as the diversity of biomass-derived fuels increases, robust quality assessment methods are essential to ensure their efficiency, safety, and minimal environmental impact. Reflected light microscopy (RLM) is one such technique with the potential to complement conventional physico-chemical analyses by enabling a rapid identification of material constituents and impurities. To refine this methodology and evaluate the reproducibility of solid biomass component identification using RLM, an interlaboratory study (ILS) was conducted. The study involved the recognition of 58 components across 45 photomicrographs, with the participation of 65 scientists and students from 25 countries. The participants faced high difficulty identifying some of the marked components, and as a result, the percentage of correct answers ranged from 19.0 % to 98.3 %, with an average correct identification rate of 62.7 %. The most challenging aspects of the identification process included distinguishing between woody and non-woody (agro) biomass, accurately identifying petroleum-derived materials, and differentiating agro biomass from inorganic matter. The results suggest that while RLM is an important tool for characterizing solid biomass, further development of methodology guidelines and training are necessary to enhance its effectiveness. Future research should prioritize preparing detailed, image-rich, microscopic morphological descriptions of biomass fuel components, which could improve the accuracy and reliability of using RLM in biomass fuel characterization.

International Journal of Coal Geology↗

USGS2025

Multispectral red-edge indices accurately estimate nitrogen content in winter cereal cover crops

Winter cover crops reduce erosion and nutrient runoff from agricultural systems. Although cereal cover crops can decrease field nitrate leaching by 50%–95%, the magnitude of this reduction varies within and between fields, making it challenging to monitor the impact of cover crops on nitrate leaching at large spatial extents. Satellite remote sensing using red-edge bands has been shown to effectively estimate crop nitrogen (N) content (kg ha −1 ) in later growth-stage crops with a closed canopy. In this study, we evaluated 15 spectral indices derived from Sentinel-2 imagery to estimate N concentration (%) and content (kg ha −1 ) of cereal cover crops, using 1627 destructive samples collected from 2018 to 2023 in Maryland. Observed N content ranged from 0.1 to 214.7 kg ha −1 , while N concentration ranged from 0.6% to 5.5%. The 15 indices considered were poor predictors of N concentration (adj. R 2 = 0.089, root mean squared error [RMSE] = 0.802%), but were more successful at measuring N content (biomass × N concentration). Delta red-edge (ΔRE) was the best predictor of N content (adj. R 2 = 0.748, RMSE = 13.10 kg ha −1 from cross-validation with 80% train and 20% test splits iterated 100 times) using samples with imagery collected within ±4 days of destructive sampling ( n = 1110). Our findings indicate that longer red-edge wavelengths (783 and 740 nm) are more suited for estimating N content in cereal cover crops compared to shorter red-edge wavelengths, which have been shown to be more sensitive to biomass. Leave-one-year-out cross-validation demonstrated that the relationship between ΔRE and N content was robust across all four cover crop sampling years included in the study (adj. R 2 = 0.700–0.769, RMSE = 10.70–15.40 kg ha −1 ). Regression model performance improved with the addition of multiple predictors, including biomass (estimated from Normalized Difference Vegetation Index), weather variables (adj. R 2 = 0.765, RMSE = 12.37 kg ha −1 ), management variables (species, season, adj. R 2 = 0.772, and RMSE = 12.13 kg ha −1 ), and biophysical variables (height, fractional ground cover, adj. R 2 = 0.818, and RMSE = 10.29 kg ha −1 ). These findings demonstrate the feasibility of quantifying N content in cereal cover crops using a red-edge-based spectral index across large geographic extents and indicate the inclusion of additional predictors, such as weather and management data, improves model accuracy. This work has implications for quantifying reductions in N leaching associated with cover crops, aiding in policymaking and evaluation of conservation programs that impact water bodies such as Chesapeake Bay.

Maryland↗

USGS2018

Metamodeling and mapping of nitrate flux in the unsaturated zone and groundwater, Wisconsin, USA

Nitrate contamination of groundwater in agricultural areas poses a major challenge to the sustainability of water resources. Aquifer vulnerability models are useful tools that can help resource managers identify areas of concern, but quantifying nitrogen (N) inputs in such models is challenging, especially at large spatial scales. We sought to improve regional nitrate (NO 3 − ) input functions by characterizing unsaturated zone NO 3 − transport to groundwater through use of surrogate, machine-learning metamodels of a process-based N flux model. The metamodels used boosted regression trees (BRTs) to relate mappable landscape variables to parameters and outputs of a previous “vertical flux method” (VFM) applied at sampled wells in the Fox, Wolf, and Peshtigo (FWP) river basins in northeastern Wisconsin. In this context, the metamodels upscaled the VFM results throughout the region, and the VFM parameters and outputs are the metamodel response variables. The study area encompassed the domain of a detailed numerical model that provided additional predictor variables, including groundwater recharge, to the metamodels. We used a statistical learning framework to test a range of model complexities to identify suitable hyperparameters of the six BRT metamodels corresponding to each response variable of interest: NO 3 − source concentration factor (which determines the local NO 3 − input concentration); unsaturated zone travel time; NO 3 − concentration at the water table in 1980, 2000, and 2020 (three separate metamodels); and NO 3 − “extinction depth”, the eventual steady state depth of the NO 3 − front. The final metamodels were trained to 129 wells within the active numerical flow model area, and considered 58 mappable predictor variables compiled in a geographic information system (GIS). These metamodels had training and cross-validation testing R 2 values of 0.52 – 0.86 and 0.22 – 0.38, respectively, and predictions were compiled as maps of the above response variables. Testing performance was reasonable, considering that we limited the metamodel predictor variables to mappable factors as opposed to using all available VFM input variables. Relationships between metamodel predictor variables and mapped outputs were generally consistent with expectations, e.g. with greater source concentrations and NO 3 − at the groundwater table in areas of intensive crop use and well drained soils. Shorter unsaturated zone travel times in poorly drained areas likely indicated preferential flow through clay soils, and a tendency for fine grained deposits to collocate with areas of shallower water table. Numerical estimates of groundwater recharge were important in the metamodels and may have been a proxy for N input and redox conditions in the northern FWP, which had shallow predicted NO 3 − extinction depth. The metamodel results provide proof-of-concept for regional characterization of unsaturated zone NO 3 − transport processes in a statistical framework based on readily mappable GIS input variables.

Wisconsin↗