Comparison of EPIC-Simulated and MODIS-Derived Leaf Area Index (LAI) across Multiple Spatial Scales.

Modeled leaf area index (LAI) in conjunction with satellite-derived LAI data streams may be used to support various regional and local scale air quality models for retrospective and future meteorological assessments. The Environmental Policy Integrated Climate (EPIC) model holds promise for providing LAI within a dynamic range for input into climate and air quality models, improving on current LAI distribution assumptions typical within atmospheric modeling. To assess the potential use of EPIC LAI, we first evaluated the Moderate Resolution Imaging Spectroradiometer (MODIS) LAI product collections 5 and 6 (i.e., Mc5, Mc6) with in situ LAI estimates upscaled at four 1.0 km resolution research sites distributed over the Albemarle-Pamlico Basin in North Carolina and Virginia, USA. We then compared the EPIC modeled 12.0 km resolution LAI to aggregated MODIS LAI (Mc5, Mc6) over a 3 × 3 grid (or 36 km × 36 km) centered over the same four research sites. Upscaled in situ LAI comparison with MODIS LAI showed improvement with the newer collection where the Mc5 overestimate of +2.22 LAI was reduced to +0.97 LAI with the Mc6. On three of the four sites, the EPIC/MODIS LAI comparison at 12.0 km resolution grid showed similar weighted mean LAI differences (LAI 1.29-1.34), with both Mc5 and Mc6 exceeding EPIC LAI across most dates. For all four research sites, both MODIS collections showed a positive bias when compared to EPIC LAI, with Mc6 (LAI = 0.40) aligning closer to EPIC than the Mc5 (LAI = 0.61) counterpart. Despite modest differences between both MODIS collections and EPIC LAI, the overestimation trend suggests the potential for EPIC to be used for future meteorological alternative management applications on a regional or national scale.

Projections of Atmospheric Nitrogen Deposition to the Chesapeake Bay Watershed.

Atmospheric deposition is among the largest pathways of nitrogen loading to the Chesapeake Bay Watershed (CBW). The interplay between future climate and emission changes in and around the CBW will likely shift the future nutrient deposition abundance and chemical regime (e.g., oxidized vs. reduced nitrogen). In this work, a Representative Concentration Pathway (RCP) from the Community Earth System Model is dynamically downscaled using the Weather Research and Forecasting (WRF) and Community Multiscale Air Quality (CMAQ) model coupled to the agro-economic Environmental Policy Integrated Climate (EPIC) model. The relative impacts of emission and climate changes on atmospheric nutrient deposition are explored for a recent historical period and a period centered on 2050. The projected regional emissions in CMAQ reflect current federal and state regulations, which use baseline and projected emission years 2011 and 2040, respectively. The historical simulations of 2-m temperature and precipitation have cool and dry biases, and temperature and precipitation are projected to both increase. Ammonium wet deposition agrees well with observations, but nitrate wet deposition is underpredicted. Climate and deposition changes increase simulated future ammonium fertilizer application. In the CBW at 2050, these changes (along with widespread decreases in anthropogenic nitrogen oxide and sulfur oxide emissions, and relatively constant NH3 emissions) decrease total nitrogen deposition by 21%, decrease annual average oxidized nitrogen deposition by 44%, and increase reduced nitrogen deposition by 10%. These results emphasize the importance of decreased anthropogenic emissions on the control of future nitrogen loading to the Chesapeake Bay in a changing climate.

Modeling the combined effects of changing land cover, climate, and atmospheric deposition on nitrogen transport in the Neuse River Basin.

STUDY REGION: The SWAT model was used to estimate the combined effects of changing land cover, climate and Clean Air Act (CAAA)-related atmospheric nitrogen (N) deposition to watershed nitrogen fate and transport for two watersheds in North Carolina, USA. STUDY FOCUS: Two different model simulation scenarios were applied: one included CAAA-related atmospheric N deposition, climate and land cover (CAAD+C+L) and the other only included CAAA-related N deposition (CAAD) in simulation. NEW HYDROLOGICAL INSIGHTS FOR THE REGION: Results show both scenarios generated overall decreasing trends for nearly all N outputs between 2010 and 2070 which resulted primarily from CAAA-related reductions in oxidized N deposition. In both watersheds, including climate and land cover change in simulation resulted in a relative 30% higher NO3 load, 30% higher denitrification, 10% higher organic N load and a 20% smaller level of plant N uptake in year 2070 compared to not including climate and landcover changes in simulation. The increases in N transport for both watersheds indicates the combined impacts from climate and land cover change may offset benefits provided by the CAAA regulations; however, future NO3 loads for the Little River watershed were small relative to current N loading rates. Conversely, the increasing NO3 and organic N loads for the nearby Nahunta watershed were significant compared to current rates demonstrating that watershed nutrient responses to climate and land cover changes may vary significantly over relatively small spatial scales.

Critical Review of Eutrophication Models for Life Cycle Assessment.

This paper evaluates the current state of life cycle impact assessment (LCIA) methods used to estimate potential eutrophication impacts in freshwater and marine ecosystems and presents a critical review of the underlying surface water quality, watershed, marine, and air fate and transport (F&T) models. Using a criteria rubric, we assess the potential of each method and model to contribute to further refinements of life cycle assessment (LCA) eutrophication mechanisms and nutrient transformation processes as well as model structure, availability, geographic scope, and spatial and temporal resolution. We describe recent advances in LCIA modeling and provide guidance on the best available sources of fate and exposure factors, with a focus on midpoint indicators. The critical review identifies gaps in LCIA characterization modeling regarding the availability and spatial resolution of fate factors in the soil compartment and identifies strategies to characterize emissions from soil. Additional opportunities are identified to leverage detailed F&T models that strengthen existing approaches to LCIA or that have the potential to link LCIA modeling more closely with the spatial and temporal realities of the effects of eutrophication.

A Comparison of Simulated and Field-Derived Leaf Area Index (LAI) and Canopy Height Values from Four Forest Complexes in the Southeastern USA.

Vegetative leaf area is a critical input to models that simulate human and ecosystem exposure to atmospheric pollutants. Leaf area index (LAI) can be measured in the field or numerically simulated, but all contain some inherent uncertainty that is passed to the exposure assessments that use them. LAI estimates for minimally managed or natural forest stands can be particularly difficult to develop as a result of interspecies competition, age and spatial distribution. Satellite-based LAI estimates hold promise for retrospective analyses, but we must continue to rely on numerical models for alternative management analysis. Our objective for this study is to calculate and validate LAI estimates generated from the USDA Environmental Policy Impact Climate (EPIC) model (a widely used, field-scale, biogeochemical model) on four forest complexes spanning three physiographic provinces in Virginia and North Carolina. Measurements of forest composition (species and number), LAI, tree diameter, basal area, and canopy height were recorded at each site during the 2002 field season. Calibrated EPIC results show stand-level temporally resolved LAI estimates with R2 values ranging from 0.69 to 0.96, and stand maximum height estimates within 20% of observation. This relatively high level of performance is attributable to EPIC's approach to the characterization of forest stand biogeochemical budgets, stand history, interspecies competition and species-specific response to local weather conditions. We close by illustrating the extension of this site-level approach to scales that could support regional air quality model simulations.

Integrating multimedia models to assess nitrogen losses from the Mississippi River basin to the Gulf of Mexico.

This study describes and implements an integrated, multimedia, process-based system-level approach to estimating nitrogen (N) fate and transport in large river basins. The modeling system includes the following components: (1) Community Multiscale Air Quality (CMAQ),(2) Weather Research and Forecasting Model (WRF), (3) Environmental Policy Integrated Climate (EPIC), and (4) Soil and Water Assessment Tool (SWAT). The previously developed Fertilizer Emission Scenario Tool for CMAQ (FEST-C), an advanced user interface, integrated EPIC with the WRF model and CMAQ. The FEST-C system, driven by process-based WRF weather simulations, includes atmospheric N additions to agricultural cropland and agricultural cropland contributions to ammonia emissions. This study focuses on integrating the watershed hydrology and water quality model with FEST-C system so that a full multimedia assessment on water quality in large river basins to address impacts of fertilization, meteorology, and atmospheric N deposition on water quality can be achieved. Objectives of this paper are to describe how to expand the previous effort by integrating the SWAT model with the FEST-C (CMAQ/WRF/EPIC) modeling system, as well as to demonstrate application of the Integrated Modeling System (IMS) to the Mississippi River basin (MRB) to simulate streamflow and dissolved N loadings to the Gulf of Mexico (GOM). IMS simulation results generally agree with US Geological Survey (USGS) observations/estimations; the annual simulated streamflow is 218.9 mm and USGS observation is 211.1 mm and the annual simulated dissolved N is 2.1 kg ha-1 and the USGS estimation is 2.8 kg ha-1. Integrating SWAT with the CMAQ/WRF/EPIC modeling system allows for its use within large river basins without losing EPIC's more detailed biogeochemistry processes, which will strengthen the assessment of impacts of future climate scenarios, regulatory and voluntary programs for N oxide air emissions, and land use and land management on N transport and transformation in large river basins.

Examining the impacts of increased corn production on groundwater quality using a coupled modeling system.

This study demonstrates the value of a coupled chemical transport modeling system for investigating groundwater nitrate contamination responses associated with nitrogen (N) fertilizer application and increased corn production. The coupled Community Multiscale Air Quality Bidirectional and Environmental Policy Integrated Climate modeling system incorporates agricultural management practices and N exchange processes between the soil and atmosphere to estimate levels of N that may volatilize into the atmosphere, re-deposit, and seep or flow into surface and groundwater. Simulated values from this modeling system were used in a land-use regression model to examine associations between groundwater nitrate-N measurements and a suite of factors related to N fertilizer and groundwater nitrate contamination. Multi-variable modeling analysis revealed that the N-fertilizer rate (versus total) applied to irrigated (versus rainfed) grain corn (versus other crops) was the strongest N-related predictor variable of groundwater nitrate-N concentrations. Application of this multi-variable model considered groundwater nitrate-N concentration responses under two corn production scenarios. Findings suggest that increased corn production between 2002 and 2022 could result in 56% to 79% increase in areas vulnerable to groundwater nitrate-N concentrations ≥5mg/L. These above-threshold areas occur on soils with a hydraulic conductivity 13% higher than the rest of the domain. Additionally, the average number of animal feeding operations (AFOs) for these areas was nearly 5 times higher, and the mean N-fertilizer rate was 4 times higher. Finally, we found that areas prone to high groundwater nitrate-N concentrations attributable to the expansion scenario did not occur in new grid cells of irrigated grain-corn croplands, but were clustered around areas of existing corn crops. This application demonstrates the value of the coupled modeling system in developing spatially refined multi-variable models to provide information for geographic locations lacking complete observational data; and in projecting possible groundwater nitrate-N concentration outcomes under alternative future crop production scenarios.

A Nitrogen Physical Input-Output Model for Illinois.

Nitrogen (N) presents an important challenge for sustainability. Human intervention in the global nitrogen cycle has been pivotal in in providing goods and services to society. However, release of N beyond its intended societal use has many negative health and environmental consequences. Several systems modeling approaches have been developed to understand the trade-offs between the beneficial and harmful effects of N. These efforts include life cycle modeling, integrated management practices and sustainability metrics for individuals and communities. However, these approaches do not connect economic and ecological N flows in physical units throughout the system, which could better represent these trade-offs for decision-makers. Physical Input-Output Table (PIOT) based models present a viable complementary solution to overcome this limitation. We developed a N-PIOT for Illinois representing the interdependence of sectors in 2002, using N mass units. This allows studying the total N flow required to produce a certain amount of N in the final product. An Environmentally Extended Input Output (EEIO) based approach was used to connect the physical economic production to environmental losses; allowing quantification of total environmental impact to support agricultural production in Illinois. A bottom up approach was used to develop the N-PIOT using Material Flow Analysis (MFA) tracking N flows associated with top 3 commodities (Corn, Soybean and Wheat). These three commodities cover 99% of N fertilizer use in Illinois. The PIOT shows that of all the N inputs to corn production the state exported 68% of N embedded in useful products, 9% went to animal feed manufacturing and only 0.03% was consumed directly within the state. Approximately 35% of N input to soybean farming ended up in animal feed. Release of N to the environment was highest from corn farming, at about 21.8% of total N fertilizer inputs, followed by soybean (9.2%) and wheat farming (4.2%). The model also allowed the calculation of life cycle N use efficiency for N based on physical flows in the economy. Hence, PIOTs prove to be a viable tool for developing a holistic approach to manage disrupted biogeochemical cycles, since these provide a detailed insight into physical flows in economic systems and allow physical coupling with ecological N flows.

Alternative futures of dissolved inorganic nitrogen export from the Mississippi River Basin: influence of crop management, atmospheric deposition, and population growth.

Nitrogen (N) export from the Mississippi River Basin contributes to seasonal hypoxia in the Gulf of Mexico (GOM). We explored monthly dissolved inorganic N (DIN) export to the GOM for a historical year (2002) and two future scenarios (year 2022) by linking macroeonomic energy, agriculture market, air quality, and agriculture land management models to a DIN export model. Future scenarios considered policies aimed at encouraging bioenergy crop production and reducing atmospheric N-emissions, as well as the effect of population growth and the states' infrastructure plans on sewage fluxes. Model-derived DIN export decreased by about 9% (from 279 to 254 kg N km-2 year-1) between 2002 and 2022 due to a 28% increase in area planted with corn, 24% improvement in crop N-recovery efficiency (NRE, to 0.52), 22% reduction in atmospheric N deposition, and 23% increase in sewage inputs. Changes in atmospheric and sewage inputs had a relatively small effect on DIN export and the effect of bioenergy crop production depended on nutrient management practices. Without improved NRE, increased production of corn would have increased DIN export by about 14% (to 289 kg N km-2 year-1) between 2002 and 2022. Model results suggest that meeting future crop demand while reducing the areal extent of hypoxia could require aggressive actions, such improving basin-level crop NRE to 0.62 or upgrading N-removal capabilities in waste water treatment plants beyond current plans. Tile-drained cropland could contribute up to half of DIN export; thus, practices that reduce N losses from tile drains could also have substantial benefit.

The role of the atmosphere in the provision of ecosystem services.

Solving the environmental problems that we are facing today requires holistic approaches to analysis and decision making that include social and economic aspects. The concept of ecosystem services, defined as the benefits people obtain from ecosystems, is one potential tool to perform such assessments. The objective of this paper is to demonstrate the need for an integrated approach that explicitly includes the contribution of atmospheric processes and functions to the quantification of air-ecosystem services. First, final and intermediate air-ecosystem services are defined. Next, an ecological production function for clean and clear air is described, and its numerical counterpart (the Community Multiscale Air Quality model) is introduced. An illustrative numerical example is developed that simulates potential changes in air-ecosystem services associated with the conversion of evergreen forest land in Mississippi, Alabama and Georgia to commercial crop land. This one-atmosphere approach captures a broad range of service increases and decreases. Results for the forest to cropland conversion scenario suggest that although such change could lead to increased biomass (food) production services, there could also be coincident, seasonally variable decreases in clean and clear air-ecosystem services (i.e., increased levels of ozone and particulate matter) associated with increased fertilizer application. Metrics that support the quantification of these regional air-ecosystem changes require regional ecosystem production functions that fully integrate biotic as well as abiotic components of terrestrial ecosystems, and do so on finer temporal scales than are used for the assessment of most ecosystem services.

A regional atmospheric fate and transport model for atrazine. 2. Evaluation.

The Community Multiscale Air Quality (CMAQ) modeling system has been adapted to simulate the regional fate and transport of atrazine. Model modifications and simulations spanning April to mid-July 1995 are described in a previous paper. CMAQ results for atrazine concentrations in air and rainfall are evaluated against field observations taken along the Mississippi River and the shores of Lake Michigan in 1995. CMAQ results agree within 10% of published annual wet deposition load estimates for Lake Michigan and predicted annual dry deposition lies within published error bounds. Comparisons of weekly observed and predicted air and rainfall concentrations along the Mississippi River yield order-of-magnitude differences. Precipitation weighting of concentrations in rainfall good agreement for seasonal time frames. Weekly ambient gas form concentrations tend to be overpredicted by the CMAQ and semivolatile particulate fractons are underpredicted. Uncertainty in CMAQ predictions of air and rainfall concentrations for atrazine appear to derive primarily from uncertainty in emissions estimates, simulated precipitation, and spatial scale.

A regional atmospheric fate and transport model for atrazine. 1. Development and implementation.

The Community Multiscale Air Quality (CMAQ) modeling system is adapted to simulate the regional transport and fate of atrazine, one of the most widely used herbicides in the United States. Model chemistry and deposition are modified, and a gas-to-particle partitioning algorithm is added to accommodate semivolatile behavior. The partitioning algorithm depends on humidity, temperature, and particulate matter concentration and composition. Results indicate that gaseous atrazine will usually dominate warm season atmospheric concentrations, but particulate form can surpass gas forms when atmospheric humidity is high (> 70%) and less-acidic (pH > 2.5) aqueous aerosol component is present. Implementation of the modified CMAQ for atrazine is illustrated, and, within the limits of our current understanding, preliminary transport and fate patterns appear to be reasonable. This research represents one of the first attempts to include a gas-to-particulate matter partitioning mechanism in an Eulerian grid-model.