Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mean”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,063 records · Page 59

SSP‐Based Land‐Use Change Scenarios: A Critical Uncertainty in Future Regional Climate Change Projections

Abstract To better understand the role projected land‐use changes (LUCs) may play in future regional climate projections, we assess the combined effects of greenhouse‐gas (GHG)‐forced climate change and LUCs in regional climate model (RCM) simulations. To do so, we produced RCM simulations that are complementary to the North‐American Coordinated Regional Downscaling Experiment (NA‐CORDEX) simulations, but with future LUCs that are consistent with particular Shared Socioeconomic Pathways (SSPs) and related to a specific Representative Concentration Pathway (RCP). We examine the state of the climate at the end of the 21st century with and without two urban and agricultural LUC scenarios that follow SSP3 and SSP5 using the Weather Research and Forecasting (WRF) model forced by one global climate model, the MPI‐ESM, under the RCP8.5 scenario. We find that LUCs following different societal trends under the SSPs can significantly affect climate projections in different ways. In regions of significant cropland expansion over previously forested area, projected annual mean temperature increases are diminished by around 0.5°C–1.0°C. Across all seasons, where urbanization is high, projected temperature increases are magnified. In particular, summer mean temperature projections are up to 4°C–5°C greater and minimum and maximum temperature projections are increased by 2.5°C–6°C, amounts that are on par with the warming due to GHG‐forced climate change. Warming is also enhanced in the urban surroundings. Future urbanization also has a large influence on precipitation projections during summer, increasing storm intensity, event length, and the overall amount over urbanized areas, and decreasing precipitation in surrounding areas.

Bukovsky, Melissa S.↗

A Retrospective and Prospective Examination of the 1960s U.S. Northeast Drought

As the most severe drought over the Northeastern United States (NEUS) in the past century, the 1960s drought had pronounced socioeconomic impacts. Although it was followed by a persisting wet period, the conditions leading to the 1960s extreme drought could return in the future, along with its challenges to water management. To project the potential consequences of such a future drought, pseudo-global warming simulations using the Weather Research and Forecasting Model are performed to simulate the dynamical conditions of the historical 1960s drought, but with modified thermodynamic conditions under the shared socioeconomic pathway SSP585 scenario in the early (2021–2027), middle (2041–2047), and late (2091–2097) 21st century. Our analysis focuses on essential hydroclimatic variables including temperature, precipitation, evapotranspiration, soil moisture, snowpack, and surface runoff. In contrast to the historical 1960s drought, similar dynamical conditions will generally produce more precipitation, increased soil moisture and evapotranspiration, and reduced snowpack. However, we also find that although wet months get much wetter, dry months may become drier, meaning that wetting trends are most significant in wet months but are essentially negligible for extremely dry months with negative monthly mean net precipitation. For these months, the trend toward wetting conditions provides little relief from the effects of extreme dry months. These conditions may even aggravate water shortages due to an increasingly rapid transition from wet to dry conditions. Other challenges emerge for residents and stakeholders in this region, including more extreme hot days, record-low snow pack, frozen ground degradation and subsequent decreases in surface runoff.

54 ENVIRONMENTAL SCIENCES↗

CMIP5 Climate Models Overestimate Cooling by Volcanic Aerosols

We compare projections of the observed hemispherical mean surface temperature (HadCRUT4.6.0.0) and the ensemble mean of CMIP5 climate models' simulations on a set of standard regression model forcing variables. We find that the volcanic aerosol regression coefficients of the CMIP5 simulations are consistently significantly larger (by 40–49%) than the volcanic aerosol coefficients of the observed temperature. The probability that the observed differences are caused just by chance is much less than 0.01. The overestimate is due to the climate models' response to volcanic aerosol radiative forcing. The largest overestimate occurs in the winter season of each hemisphere. We hypothesize that the models' parameterization of aerosol-cloud interactions within ice and mixed phase clouds is a likely source of this discrepancy. Furthermore, the models significantly underestimate the effect of solar variability on temperature for both hemispheres.

54 ENVIRONMENTAL SCIENCES↗

Weakening of the Extratropical Storm Tracks in Solar Geoengineering Scenarios

Solar geoengineering that aims to offset global warming could nonetheless alter atmospheric temperature gradients and humidity and thus affect the extratropical storm tracks. Here, we first analyze climate model simulations from experiment G1 of the Geoengineering Model Intercomparison Project, in which a reduction in incoming solar radiation balances a quadrupling of CO 2 . The Northern Hemisphere extratropical storm track weakens by a comparable amount in G1 as it does for increased CO2 only. The Southern Hemisphere storm track also weakens in G1, in contrast to a strengthening and poleward shift for increased CO 2 . Using mean available potential energy, we show that the changes in zonal-mean temperature and humidity are sufficient to explain the different responses of storm-track intensity. We also demonstrate similar weakening in a more complex geoengineering scenario. Here, our results offer insight into how geoengineering affects storm tracks, highlighting the potential for geoengineering to induce novel climate changes.

54 ENVIRONMENTAL SCIENCES↗

Analysis of Electric and Magnetic Lightning-Generated Wave Amplitudes Measured by the Van Allen Probes

Herein we provide a statistical analysis of both electric and magnetic field wave amplitudes of very low frequency lightning-generated waves (LGWs) based on the equivalent of 11.5 years of observations made by the Van Allen Probes encompassing ~24.6 × 10 6 survey mode measurements. We complement this analysis with data from the ground-based World Wide Lightning Location Network to explore differences between satellite and ground-based measurements. LGW mean amplitudes are found to be low compared with other whistler mode waves (1 ± 1.6 pT and 19 ± 59 μV/m). Extreme events (1/5,000) can reach 100 pT and contributes strongly to the mean power below $L$ = 2. We find excellent correlations between World Wide Lightning Location Network-based power and wave amplitudes in space at various longitudes. We reveal strong dayside ionospheric damping of the LGW electric field. LGW amplitudes drop for $L$ < 2, contrary to the Earth's intense equatorial lightning activity. We conclude that it is difficult for equatorial LGW to propagate and remain at $L$ < 2.

Electric wave power↗

Seasonally Modulated Stratospheric Aerosol Geoengineering Alters the Climate Outcomes

Abstract By reflecting some incoming solar radiation, stratospheric aerosol intervention using SO 2 would reduce global mean temperature. Previous research has shown that multiple injection latitudes can be used to maintain not only global mean temperature, but also interhemispheric and equator‐to‐pole temperature gradients. However, the regional climate response depends not only on where the SO 2 is injected, but also on when. We show here that even with these same objectives and same choices of latitudes, injecting in only one season instead of continuously throughout the year results in significant differences in regional climate, for instance in the magnitude of precipitation changes over India. The differential outcomes highlight the potential for underlying trade‐offs, with different choices regarding deployment leading to a different distribution of benefits or harms. This aspect of climate engineering should be considered in developing governance and emphasizes the need for a deeper understanding of the mechanisms underlying the regional responses.

Visioni, Daniele↗

Quantifying the Occurrence of Record Hot Years Through Normalized Warming Trends

Surface air temperature trends and extreme events are of global concern and they are related. Here we show that the occurrence of record hot years over different latitudes from 1960-2019 are more strongly correlated with the observational annual mean temperature trends normalized by internal variability. Compared with the raw trends showing Arctic amplification, the normalized trends show a tropical amplification over land. Two hot spots with more frequent occurrence of record hot years are identified: northern hemisphere ocean (versus land) and southern hemisphere tropical land (versus mid- and high-latitude lands). Ensemble mean results from 32 Earth system models agree with observations better than individual models, but they do not reproduce observed large differences in correlations across latitudes between normalized trends and record-breaking events over land versus ocean. Furthermore, our results enable the quantification of record hot year occurrence through normalized warming trends and provide new metrics for model evaluation and improvement.

58 GEOSCIENCES↗

Notable Contributions of Aerosols to the Predictability of Hail Precipitation

There is an increasing concern of the uncertainty produced by aerosols in forecasting precipitation including hail precipitation. This study provides an assessment of the uncertainties in hail and total precipitation by varying initial cloud condensation nuclei (CCN) number concentration (CCNC) and meteorological conditions based on 1200 cloud-resolving simulations of an idealized hailstorm. Although the meteorological perturbations produce large uncertainties in hail precipitation (including rate and maximum hail size) as well as total precipitation, varying CCNC by an order of magnitude can cause even larger uncertainties, especially pairing with the thermodynamics perturbation (i.e., potential temperature and water vapor). Changing CCNC modifies the predictability of hail precipitation, with a higher predictability in moderate polluted environments compared with the very clean and polluted environments. Increasing CCNC consistently leads a non-monotonic response of ensemble mean with an optimal CCNC for hail precipitation but a monotonic decreasing response of total precipitation with the various meteorological perturbations, meaning the initial meteorological perturbations does not qualitatively change the aerosol effects. Investigation with 10-fold reduced initial perturbation further supports the large CCN effects are not dependent of metrological perturbations. The findings suggest the importance of considering CCN effects in severe weather simulations and forecasting.

54 ENVIRONMENTAL SCIENCES↗

Groundwater regulates interannual variations in evapotranspiration in a riparian semiarid ecosystem

Dryland ecosystems are an important component of Earth’s carbon, water, and energy budgets. Changes in hydroclimatic conditions can significantly alter land surface fluxes in these ecosystems. In regions with groundwater-river water exchange, land surface fluxes are strongly modulated by groundwater table fluctuations in response to river stage variations; however, interactions between fluxes and groundwater remain poorly understood in drylands. Here, we use eddy covariance data collected in paired upland and riparian semiarid ecosystems during two years with contrasting river flow conditions to examine interactions among river stage conditions, water table dynamics, meteorological drivers, and ecosystem fluxes. Elevated groundwater availability resulting from groundwater-surface water exchange at the riparian site enhanced the latent heat flux (LE, the energy form of evapotranspiration), which was sustained longer in the growing season during both years. At the upland site, mean LE during the dry season differed by less than 2% between the two years. In contrast, at the riparian site contrasting river flow conditions resulted in a 41% difference in peak monthly LE between the two years. The two-year mean LE in the summer at the riparian site was twice that at the upland site. Carbon uptake (NEP) and LE were well correlated at the upland site during both years, but at the riparian site enhanced summertime LE during the high river flow year did not correspond t enhanced NEP. Our results demonstrate that enhanced groundwater availability modulated by river stage variability is a critical driver of land surface fluxes in semiarid riparian ecosystems.

Missik, Justine E.↗

Saharan Dust Transport Predictability Utilizing a Subseasonal Experiment (SubX) Model

Prediction of Saharan dust customarily requires complex aerosols models and observations. A previous study of the Miami, Florida dust record in conjunction to reanalysis data discovered a possible source of subseasonal predictability using a dust-transport-efficiency (DTE) index. Development of the Subseasonal Forecast Experiment (SubX) has expanded global forecast products; producing multi-model ensemble forecasts out to 45 days. Retrospective forecast data from the Community Climate System Model version 4.0 (CCSM4) is used in direct comparison to National Centers for Environmental Prediction (NCEP) reanalysis to evaluate the CCSM4 subseasonal forecast and DTE index prediction skill of weekly dust variability. Successful prediction of weekly dust transport using the DTE index is variable year-to-year. The DTE most successfully predicts dust when there is high variability in the tropical winds, likely due to a fluctuating subtropical high, and is not dependent on the overall mean flow or total dust mass transported. The CCSM4 SubX retrospective forecast well represents North Atlantic meteorology out to week-3 in both mean flow and variability. Dust transport can be predicted using the DTE index and CCSM4 SubX retrospective forecasts at week-1 leads to the same success as contemporaneous NCEP reanalysis.

54 ENVIRONMENTAL SCIENCES↗

Consideration of Initial Cloud Droplet Size Distribution Shapes in Quantifying Different Entrainment-Mixing Mechanisms

Entrainment-mixing process significantly affects cloud micro/macrophysics and the evaluation of aerosol indirect effects. However, it remains an open question as to how to quantify entrainment-mixing mechanism for broad cloud droplet size distributions (CDSDs). Here, CDSDs with different spectral widths are used to initialize the Explicit Mixing Parcel Model to address this problem, and microphysical properties are compared for different CDSD widths. For relatively broad CDSDs, as the number concentration and liquid water content decrease, the volume-mean radius and mean radius increase, which is opposite to the scenario for relatively narrow CDSDs and the homogeneous/inhomogeneous conceptual model. For the relatively broad CDSDs, such relationships could be mistaken as inhomogeneous mixing with subsequent ascent in analysis of the conventional microphysical mixing diagram. This then causes traditional homogenous mixing degrees to be unrealistically negative and not applicable for relatively broad CDSDs. New measures are introduced to explicitly account for the effect of CDSD spectral shapes on quantifying entrainment-mixing mechanism. The new measures yield reasonable ranges of values that conform to the dynamical measures such as the Damköhler and transitional scale numbers, and their temporal variation can be explained physically. This study extends the measures of homogeneous mixing degrees from narrow to broad CDSDs, and sheds new light on the consideration of CDSD shapes in parameterizations of entrainment-mixing mechanism.

54 ENVIRONMENTAL SCIENCES↗

Collar Properties and Measurement Time Confer Minimal Bias Overall on Annual Soil Respiration Estimates in a Global Database

Abstract Measuring the soil‐to‐atmosphere carbon dioxide (CO 2 ) flux (soil respiration, R S ) is important to understanding terrestrial carbon balance and to forecasting climate change. Such measurements are frequently made using measurement collars permanently inserted into the soil surface. However, differences in measurement duration and frequency, as well as collar properties, may lead to biases in the estimation of annual R S . Using a newly updated global R S database (SRDB‐V5), we investigated the annual R S bias associated with five methodological factors: collar height, collar coverage area, collar insertion depth, measurement duration, and measurement frequency. We found that annual R S was negatively correlated with collar insertion depth, consistent with the idea that collar insertion cuts roots and thus reduces R S . Annual R S was also negatively related with collar height and collar coverage area, perhaps because uniform head‐space mixing is difficult to achieve in larger volume chambers; however, these effects were quantitatively small (bias of ~2% to 10% of mean R S ). We found no correlation of measurement duration or measurement frequency with annual R S . These findings suggest that variation in R S methodology generally introduces minimal bias overall. Therefore, compilations of minimally adjusted annual R S measurements provide a reliable resource for synthesis studies, global annual R S modeling, and investigation of how soil carbon responds to climate change.

Jian, Jinshi↗

Projecting Permafrost Thaw of Sub-Arctic Tundra With a Thermodynamic Model Calibrated to Site Measurements

Northern circumpolar permafrost thaw affects global carbon cycling, as large amounts of stored soil carbon becomes accessible to microbial breakdown under a warming climate. The magnitude of carbon release is linked to the extent of permafrost thaw, which is locally variable and controlled by soil thermodynamics. Soil thermodynamic properties, such as thermal diffusivity, govern the reactivity of the soil-atmosphere thermal gradient, and are controlled by soil composition and drainage. In order to project permafrost thaw for an Alaskan tundra experimental site, we used seven years of site data to calibrate a soil thermodynamic model using a data assimilation technique. The model reproduced seasonal and interannual temperature dynamics for shallow (5–40 cm) and deep soil layers (2–4 m), and simulations of seasonal thaw depth closely matched observed data. The model was then used to project permafrost thaw at the site to the year 2100 using climate forcing data for three future climate scenarios (RCP 4.5, 6.0, and 8.5). Minimal permafrost thawing occurred until mean annual air temperatures rose above the freezing point, after which we measured over a 1 m increase in thaw depth for every 1 °C rise in mean annual air temperature. Under no projected warming scenario was permafrost remaining in the upper 3 m of soil by 2100. We demonstrated an effective data assimilation method that optimizes parameterization of a soil thermodynamic model. The sensitivity of local permafrost to climate warming illustrates the vulnerability of sub-Arctic tundra ecosystems to significant and rapid soil thawing.

54 ENVIRONMENTAL SCIENCES↗

Distortions of the Rain Distribution With Warming, With and Without Self-Aggregation

We investigate how mesoscale circulations associated with convective aggregation can modulate the sensitivity of the hydrologic cycle to warming. We quantify changes in the full distribution of rain across radiative-convective equilibrium states in a cloud-resolving model. For a given Sea Surface Temperature (SST), the shift in mean rainfall between disorganized and organized states is associated with a shift in atmospheric radiative cooling, and is roughly analogous to the effect of a 4K SST increase. With rising temperatures, the increase in mean rain rate is insensitive to the presence of organization, while extremes can intensify faster in the aggregated state, leading to a faster amplification in the sporadic nature of rain. When convection aggregates, heavy rain is enhanced by 20%–30% and nonlinear behaviors are observed as a function of SST and strength of aggregation feedbacks. First, radiative- and surface-flux aggregation feedbacks have multiplicative effects on extremes, illustrating a non-trivial sensitivity to the degree of organization. Second, alternating Clausius-Clapeyron and super-Clausius-Clapeyron regimes in extreme rainfall are found as a function of SST, corresponding to varying thermodynamic and dynamic contributions, and a large sensitivity to precipitation efficiency variations in some SST ranges. The potential for mesoscale circulations in amplifying the hydrologic cycle is established. However, these nonlinear distortions question the quantitative relevance of idealized self-aggregation. This calls for a deeper investigation of relationships which capture the coupling between global energetics, aggregation feedbacks and local convection, and for systematic tests of their sensitivity to domain configurations, surface boundary conditions, microphysics, and turbulence schemes.

54 ENVIRONMENTAL SCIENCES↗

Dynamics of Fungal and Bacterial Biomass Carbon in Natural Ecosystems: Site-level Applications of the CLM-Microbe Model

Explicitly representing microbial processes has been recognized as a key improvement to Earth system models for the realistic projections of soil carbon (C) and climate dynamics. The CLM-Microbe model builds upon the CLM4.5 and explicitly represents two major soil microbial groups, fungi and bacteria. Based on the compiled time-series data of fungal (FBC) and bacterial (BBC) biomass C from nine biomes, we parameterized and validated the CLM-Microbe model, and further conducted sensitivity analysis and uncertainty analysis for simulating C cycling. The model performance was evaluated with mean absolute error (MAE), root mean square error (RMSE), and coefficient of determination (R2) for relative change in FBC and BBC. The CLM-Microbe model is able to reasonably capture the seasonal dynamics of FBC and BBC across biomes, particularly for tropical/subtropical forest, temperate broadleaf forest, and grassland, with MAE < 0.49 for FBC and <0.36 for BBC and RMSE <0.52 FBC and <0.39 for BBC, while R2 values are relatively smaller in some biomes (e.g., shrub) due to small sample sizes. We found good consistencies between simulated and observed FBC (R 2 =0.70, P<0.001) and BBC (R 2 =0.26, P<0.05) on average across biomes, but the model is not able to fully capture the large variation in observed FBC and BBC. Sensitivity analysis shows the most critical parameters are turnover rate, carbon-to-nitrogen ratio of fungi and bacteria, and microbial assimilation efficiency. This study confirms that the explicit representation of soil microbial mechanisms enhances model performance in simulating C variables such as heterotrophic respiration and soil organic C density. The further application of the CLM-Microbe model would deepen our understanding of microbial contributions to the global C cycle.

54 ENVIRONMENTAL SCIENCES↗

Linking Deep and Shallow Convective Mass Fluxes via an Assumed Entrainment Distribution in CAM5‐CLUBB: Parameterization and Simulated Precipitation Variability

Abstract We modify the Zhang‐McFarlane (ZM) deep convection scheme in the Community Atmosphere Model version 5 to couple it with a unified parameterization for boundary‐layer turbulence and shallow convection, that is, Cloud Layers Unified by Binormals (CLUBB). By assuming a lognormal distribution of entrainment rate across the entire moist convective regimes, we link mass fluxes between shallow and deep convection, which are partitioned by the entrainment rate of the shallowest deep convective plume. Hence, a new deep convective closure is established which is coupled to the sub‐grid vertical motion variability in CLUBB. The convection feedback (or memory) effects are also considered to decrease the entrainment spectrum width and enhance the vertical velocity variability that further affect deep convection. Results show that the revised scheme improves the precipitation simulations in terms of the mean state and variability at various timescales, such as the alleviated double‐intertropical convergence zone and more realistic simulations of the seasonal variation of monsoon precipitation over East Asia, Madden‐Julian Oscillation, and precipitation diurnal phase propagations downstream of large terrains. The improvements are still seen in many aspects such as the mean‐state precipitation when turning off the convection feedback impacts in the revised scheme, emphasizing the benefits of using the modified mass‐flux closure. However, the convection feedbacks have considerable effects on the precipitation diurnal cycle simulations over regions with late‐afternoon precipitation peaks. Overall, the revised scheme provides a unified treatment for sub‐grid vertical motions across regimes of boundary‐layer turbulence, shallow convection, and deep convection, leading to better‐simulated precipitation at various timescales.

54 ENVIRONMENTAL SCIENCES↗

Global Temperature Responses to Large Tropical Volcanic Eruptions in Paleo Data Assimilation Products and Climate Model Simulations Over the Last Millennium

Abstract Large volcanic eruptions are one of the dominant perturbations to global and regional atmospheric temperatures on timescales of years to decades. Discrepancies remain, however, in the estimated magnitude and persistence of the surface temperature cooling caused by volcanic eruptions, as characterized by paleoclimatic proxies and climate models. We investigate these discrepancies in the context of large tropical eruptions over the Last Millennium using two state‐of‐the‐art data assimilation products, the Paleo Hydrodynamics Data Assimilation product (PHYDA) and the Last Millennium Reanalysis (LMR), and simulations from the National Center for Atmospheric Research Community Earth System Model‐Last Millennium Ensemble (NCAR CESM‐LME). We find that PHYDA and LMR estimate mean global and hemispheric cooling that is similar in magnitude and persistence once effects from eruptions occurring in short succession are removed. The estimates also compare well to Northern‐Hemisphere reconstructions based solely or partially on tree‐ring density, which have been proposed as the most accurate proxy estimates of surface cooling due to volcanism. All proxy‐based estimates also agree well with the magnitude of the mean cooling simulated by the CESM‐LME. Differences remain, however, in the spatial patterns of the temperature responses in the PHYDA, LMR, and the CESM‐LME. The duration of cooling anomalies also persists for several years longer in the PHYDA and LMR relative to the CESM‐LME. Our results demonstrate progress in resolving discrepancies between proxy‐ and model‐based estimates of temperature responses to volcanism, but also indicate these estimates must be further reconciled to better characterize the risks of future volcanic eruptions.

Tejedor, E.↗

Impact of vegetation physiology and phenology on watershed hydrology in a semiarid watershed in the Pacific Northwest in a changing climate

Concurrent with changes in climate conditions, increases in atmospheric composition and transport can affect land surface processes by regulating physiological and phenological responses of plants, and consequently perturb water, carbon, and nitrogen (N) dynamics at multiple scales. Howerver, the relationships bwtween plants physiological and phenological and land surface processes in a changing climate are less well known due to the limitation of representing palnts properties in traditional hydrological models. In this study, we use the Community Land Model version 5 (CLM5) to investigate how projected changes in large scale atmospheric conditions can affect hydrologic and biogeochemical dynamics in the Upper Columbia-Priest Rapids (UCPR) watershed, a typical semi-arid watershed located in the north-western United States dominated by croplands. Six numerical experiments considering climate change, CO2 concentration, and N deposition consistent with the Representative Concentration Pathways (RCPs) 4.5 and 8.5 are performed at 1-km resolution over the period of 2015-2099. Our results show that, compared with historical results (1985-2014), the elevated CO2 directly increases plants N uptake rates that provide additional leaf N to enhance photosynthesis by the end of 21st century (2070-2099). The higher GPP leads to changes in stomatal resistance and increased leaf area index (LAI) that affect evaporatranspiration (ET) and its components. As a result of stomatal responses to a higher CO2 concentration, ET over corn fields declines by 24.6 mm yr-1 under RCP4.5, and by 64.0 mm yr-1 under RCP8.5, respectively by the end of the century from its historical mean of 428 mm yr-1. Meanwhile, ET over soybean fields increases by 27.7 mm yr-1 under RCP4.5 and 9.98 mm yr-1 under RCP8.5 from its historical mean of 507 mm yr-1 as a result of CO2 fertilization as soybeans are able to uptake more N to increase LAI significantly. Furtherfore, the changes in physiological and phenological in corn and soybean can strongly modulate water budget terms over the watershed in response to expected changes in greenhouse gas concentrations and atmospheric composition. Additionally, increases in CO2 concentration and N deposition are beneficial for improving water use efficiency of corn and soybean, therefore lead to less water consumption per unit gram of carbon produced from these crops. Our study elucidates the need to jointly consider effects of surface climate and atmospheric composition and transport in understanding and predicting watershed hydrologic and biogeochemical responses, especially for an agricultural watershed.

Zhu, Bowen↗