Warming may offset impact of precipitation changes on riverine nitrogen loading.

Climate change, especially in the form of precipitation and temperature changes, can alter the transformation and delivery of nitrogen on the land surface and to aquatic systems, impacting the trophic states of downstream water bodies. While the expected impacts of changes in precipitation have been explored, a quantitative understanding of the impact of temperature on nitrogen loading is lacking at landscape scales. Here, using several decades of nitrogen loading observations, we quantify how individual and combined future changes in precipitation and temperature will affect riverine nitrogen loading. We find that, contrary to recent decades, rising temperatures are likely to offset or even reverse previously reported impacts of future increases in total and extreme precipitation on nitrogen runoff across the majority of the contiguous United States. These findings highlight the multifaceted impacts of climate change on the global nitrogen cycle.

Widespread global increase in intense lake phytoplankton blooms since the 1980s.

Freshwater blooms of phytoplankton affect public health and ecosystem services globally1,2. Harmful effects of such blooms occur when the intensity of a bloom is too high, or when toxin-producing phytoplankton species are present. Freshwater blooms result in economic losses of more than US$4 billion annually in the United States alone, primarily from harm to aquatic food production, recreation and tourism, and drinking-water supplies3. Studies that document bloom conditions in lakes have either focused only on individual or regional subsets of lakes4-6, or have been limited by a lack of long-term observations7-9. Here we use three decades of high-resolution Landsat 5 satellite imagery to investigate long-term trends in intense summertime near-surface phytoplankton blooms for 71 large lakes globally. We find that peak summertime bloom intensity has increased in most (68 per cent) of the lakes studied, revealing a global exacerbation of bloom conditions. Lakes that have experienced a significant (P < 0.1) decrease in bloom intensity are rare (8 per cent). The reason behind the increase in phytoplankton bloom intensity remains unclear, however, as temporal trends do not track consistently with temperature, precipitation, fertilizer-use trends or other previously hypothesized drivers. We do find, however, that lakes with a decrease in bloom intensity warmed less compared to other lakes, suggesting that lake warming may already be counteracting management efforts to ameliorate eutrophication10,11. Our findings support calls for water quality management efforts to better account for the interactions between climate change and local hydrological conditions12,13.

Global patterns of drought recovery.

Drought, a recurring phenomenon with major impacts on both human and natural systems, is the most widespread climatic extreme that negatively affects the land carbon sink. Although twentieth-century trends in drought regimes are ambiguous, across many regions more frequent and severe droughts are expected in the twenty-first century. Recovery time-how long an ecosystem requires to revert to its pre-drought functional state-is a critical metric of drought impact. Yet the factors influencing drought recovery and its spatiotemporal patterns at the global scale are largely unknown. Here we analyse three independent datasets of gross primary productivity and show that, across diverse ecosystems, drought recovery times are strongly associated with climate and carbon cycle dynamics, with biodiversity and CO2 fertilization as secondary factors. Our analysis also provides two key insights into the spatiotemporal patterns of drought recovery time: first, that recovery is longest in the tropics and high northern latitudes (both vulnerable areas of Earth's climate system) and second, that drought impacts (assessed using the area of ecosystems actively recovering and time to recovery) have increased over the twentieth century. If droughts become more frequent, as expected, the time between droughts may become shorter than drought recovery time, leading to permanently damaged ecosystems and widespread degradation of the land carbon sink.

India's Riverine Nitrogen Runoff Strongly Impacted by Monsoon Variability.

Agricultural intensification in India has increased nitrogen pollution, leading to water quality impairments. The fate of reactive nitrogen applied to the land is largely unknown, however. Long-term records of riverine nitrogen fluxes are nonexistent and drivers of variability remain unexamined, limiting the development of nitrogen management strategies. Here, we leverage dissolved inorganic nitrogen (DIN) and discharge data to characterize the seasonal, annual, and regional variability of DIN fluxes and their drivers for seven major river basins from 1981 to 2014. We find large seasonal and interannual variability in nitrogen runoff, with 68% to 94% of DIN fluxes occurring in June through October and with the coefficient of variation across years ranging from 44% to 93% for individual basins. This variability is primarily explained by variability in precipitation, with year- and basin-specific annual precipitation explaining 52% of the combined regional and interannual variability. We find little correlation with rising fertilizer application rates in five of the seven basins, implying that agricultural intensification has thus far primarily impacted groundwater and atmospheric emissions rather than riverine runoff. These findings suggest that riverine nitrogen runoff in India is highly sensitive to projected future increases in precipitation and intensification of the seasonal monsoon, while the impact of projected continued land use intensification is highly uncertain.

The terrestrial biosphere as a net source of greenhouse gases to the atmosphere.

The terrestrial biosphere can release or absorb the greenhouse gases, carbon dioxide (CO2), methane (CH4) and nitrous oxide (N2O), and therefore has an important role in regulating atmospheric composition and climate. Anthropogenic activities such as land-use change, agriculture and waste management have altered terrestrial biogenic greenhouse gas fluxes, and the resulting increases in methane and nitrous oxide emissions in particular can contribute to climate change. The terrestrial biogenic fluxes of individual greenhouse gases have been studied extensively, but the net biogenic greenhouse gas balance resulting from anthropogenic activities and its effect on the climate system remains uncertain. Here we use bottom-up (inventory, statistical extrapolation of local flux measurements, and process-based modelling) and top-down (atmospheric inversions) approaches to quantify the global net biogenic greenhouse gas balance between 1981 and 2010 resulting from anthropogenic activities and its effect on the climate system. We find that the cumulative warming capacity of concurrent biogenic methane and nitrous oxide emissions is a factor of about two larger than the cooling effect resulting from the global land carbon dioxide uptake from 2001 to 2010. This results in a net positive cumulative impact of the three greenhouse gases on the planetary energy budget, with a best estimate (in petagrams of CO2 equivalent per year) of 3.9 ± 3.8 (top down) and 5.4 ± 4.8 (bottom up) based on the GWP100 metric (global warming potential on a 100-year time horizon). Our findings suggest that a reduction in agricultural methane and nitrous oxide emissions, particularly in Southern Asia, may help mitigate climate change.

Global vegetation biomass production efficiency constrained by models and observations.

Plants use only a fraction of their photosynthetically derived carbon for biomass production (BP). The biomass production efficiency (BPE), defined as the ratio of BP to photosynthesis, and its variation across and within vegetation types is poorly understood, which hinders our capacity to accurately estimate carbon turnover times and carbon sinks. Here, we present a new global estimation of BPE obtained by combining field measurements from 113 sites with 14 carbon cycle models. Our best estimate of global BPE is 0.41 ± 0.05, excluding cropland. The largest BPE is found in boreal forests (0.48 ± 0.06) and the lowest in tropical forests (0.40 ± 0.04). Carbon cycle models overestimate BPE, although models with carbon-nitrogen interactions tend to be more realistic. Using observation-based estimates of global photosynthesis, we quantify the global BP of non-cropland ecosystems of 41 ± 6 Pg C/year. This flux is less than net primary production as it does not contain carbon allocated to symbionts, used for exudates or volatile carbon compound emissions to the atmosphere. Our study reveals a positive bias of 24 ± 11% in the model-estimated BP (10 of 14 models). When correcting models for this bias while leaving modeled carbon turnover times unchanged, we found that the global ecosystem carbon storage change during the last century is decreased by 67% (or 58 Pg C).

Long-Term Changes in Precipitation and Temperature Have Already Impacted Nitrogen Loading.

Increases in nitrogen loading over the past several decades have led to widespread water quality impairments across the U.S. Elevated awareness of the influence of climate variability on nitrogen loading has led to several studies investigating future climate change impacts on water quality. However,  it remains unclear whether long-term climate impacts can already be observed in the historical record. Here, we quantify long-term trends in total nitrogen loading over the period 1987-2012 across the contiguous U.S. and attribute these trends to long-term changes in nitrogen inputs and climatic variables. We find that annual precipitation, extreme springtime precipitation, and springtime temperature are key drivers of trends in historical loading in most regions. These decadal climate trends have either  amplified or offset loading trends expected from nitrogen inputs alone. We also find that rising temperatures have been insufficient to offset precipitation-induced loading increases, suggesting that future increases in temperature under climate change may have limited potential to counteract loading increases expected as a result of anticipated changes in precipitation. This work demonstrates the important role of decadal climate variability in long-term nitrogen loading, emphasizing the need to consider climate change risks when designing and monitoring nutrient reduction programs.

Long-Term Phosphorus Loading and Springtime Temperatures Explain Interannual Variability of Hypoxia in a Large Temperate Lake.

Anthropogenic eutrophication has led to the increased occurrence of hypoxia in inland and coastal waters around the globe. While low dissolved oxygen conditions are known to be driven primarily by nutrient loading and water column stratification, the relative importance of these factors and their associated time scales are not well understood. Here, we explore these questions for Lake Erie, a large temperate lake that experiences widespread annual summertime hypoxia. We leverage a three-decade data set of summertime hypoxic extent (1985-2015) and examine the role of seasonal and long-term nutrient loading, as well as hydrometeorological conditions. We find that a linear combination of decadal total phosphorus loading from tributaries and springtime air temperatures explains a high proportion of the interannual variability in average summertime hypoxic extent (R2 = 0.71). This result suggests that the lake responds primarily to long-term variations in phosphorus inputs, rather than springtime or annual loading as previously assumed, which is consistent with internal phosphorus loading from lake sediments likely being an important contributing mechanism. This result also demonstrates that springtime temperatures have a substantial impact on summertime hypoxia, likely by impacting the timing of onset of thermal stratification. These findings imply that management strategies based on reducing tributary phosphorus loading would take several years to reap full benefits, and that projected future increases in temperatures are likely to exacerbate hypoxia in Lake Erie and other temperate lakes.

Elliptic Cylinder Airborne Sampling and Geostatistical Mass Balance Approach for Quantifying Local Greenhouse Gas Emissions.

In this study, we explore observational, experimental, methodological, and practical aspects of the flux quantification of greenhouse gases from local point sources by using in situ airborne observations, and suggest a series of conceptual changes to improve flux estimates. We address the major sources of uncertainty reported in previous studies by modifying (1) the shape of the typical flight path, (2) the modeling of covariance and anisotropy, and (3) the type of interpolation tools used. We show that a cylindrical flight profile offers considerable advantages compared to traditional profiles collected as curtains, although this new approach brings with it the need for a more comprehensive subsequent analysis. The proposed flight pattern design does not require prior knowledge of wind direction and allows for the derivation of an ad hoc empirical correction factor to partially alleviate errors resulting from interpolation and measurement inaccuracies. The modified approach is applied to a use-case for quantifying CH4 emission from an oil field south of San Ardo, CA, and compared to a bottom-up CH4 emission estimate.

Regional atmospheric CO2 inversion reveals seasonal and geographic differences in Amazon net biome exchange.

Understanding tropical rainforest carbon exchange and its response to heat and drought is critical for quantifying the effects of climate change on tropical ecosystems, including global climate-carbon feedbacks. Of particular importance for the global carbon budget is net biome exchange of CO2 with the atmosphere (NBE), which represents nonfire carbon fluxes into and out of biomass and soils. Subannual and sub-Basin Amazon NBE estimates have relied heavily on process-based biosphere models, despite lack of model agreement with plot-scale observations. We present a new analysis of airborne measurements that reveals monthly, regional-scale (~1-8 × 10(6)  km(2) ) NBE variations. We develop a regional atmospheric CO2 inversion that provides the first analysis of geographic and temporal variability in Amazon biosphere-atmosphere carbon exchange and that is minimally influenced by biosphere model-based first guesses of seasonal and annual mean fluxes. We find little evidence for a clear seasonal cycle in Amazon NBE but do find NBE sensitivity to aberrations from long-term mean climate. In particular, we observe increased NBE (more carbon emitted to the atmosphere) associated with heat and drought in 2010, and correlations between wet season NBE and precipitation (negative correlation) and temperature (positive correlation). In the eastern Amazon, pulses of increased NBE persisted through 2011, suggesting legacy effects of 2010 heat and drought. We also identify regional differences in postdrought NBE that appear related to long-term water availability. We examine satellite proxies and find evidence for higher gross primary productivity (GPP) during a pulse of increased carbon uptake in 2011, and lower GPP during a period of increased NBE in the 2010 dry season drought, but links between GPP and NBE changes are not conclusive. These results provide novel evidence of NBE sensitivity to short-term temperature and moisture extremes in the Amazon, where monthly and sub-Basin estimates have not been previously available.

A multi-year estimate of methane fluxes in Alaska from CARVE atmospheric observations.

Methane (CH4) fluxes from Alaska and other arctic regions may be sensitive to thawing permafrost and future climate change, but estimates of both current and future fluxes from the region are uncertain. This study estimates CH4 fluxes across Alaska for 2012-2014 using aircraft observations from the Carbon in Arctic Reservoirs Vulnerability Experiment (CARVE) and a geostatistical inverse model (GIM). We find that a simple flux model based on a daily soil temperature map and a static map of wetland extent reproduces the atmospheric CH4 observations at the state-wide, multi-year scale more effectively than global-scale, state-of-the-art process-based models. This result points to a simple and effective way of representing CH4 flux patterns across Alaska. It further suggests that contemporary process-based models can improve their representation of key processes that control fluxes at regional scales, and that more complex processes included in these models cannot be evaluated given the information content of available atmospheric CH4 observations. In addition, we find that CH4 emissions from the North Slope of Alaska account for 24% of the total statewide flux of 1.74 ± 0.44 Tg CH4 (for May-Oct.). Contemporary global-scale process models only attribute an average of 3% of the total flux to this region. This mismatch occurs for two reasons: process models likely underestimate wetland area in regions without visible surface water, and these models prematurely shut down CH4 fluxes at soil temperatures near 0°C. As a consequence, wetlands covered by vegetation and wetlands with persistently cold soils could be larger contributors to natural CH4 fluxes than in process estimates. Lastly, we find that the seasonality of CH4 fluxes varied during 2012-2014, but that total emissions did not differ significantly among years, despite substantial differences in soil temperature and precipitation; year-to-year variability in these environmental conditions did not affect obvious changes in total CH4 fluxes from the state.

Precipitation Dominates Interannual Variability of Riverine Nitrogen Loading across the Continental United States.

Excessive nitrogen loading to waterways leads to increased eutrophication and associated water quality impacts. An understanding of the regional and interannual variability in nitrogen loading and associated drivers is necessary for the design of effective management strategies. Here we develop a parsimonious empirical model based on net anthropogenic nitrogen input, precipitation, and land use that explains 68% of the observed variability in annual total nitrogen flux (QTN) (76% of ln(QTN)) across 242 catchment years. We use this model to present the first spatially and temporally resolved estimates of QTN for all eight-digit hydrologic unit (HUC8) watersheds within the continental United States (CONUS), focusing on the period 1987-2007. Results reveal high spatial and temporal variability in loading, with spatial variability primarily driven by nitrogen inputs, but with interannual variability and the occurrence of extremes dominated by precipitation across over three-quarters of the CONUS. High interannual variability and its correlation with precipitation persist at large aggregated scales. These findings point to a fundamental challenge in managing regions with high nutrient loading, because these regions also exhibit the strongest interannual variability and because the impact of changes in management practices will be modulated by meteorological variability and climatic trends.

Anthropogenic emissions of methane in the United States.

This study quantitatively estimates the spatial distribution of anthropogenic methane sources in the United States by combining comprehensive atmospheric methane observations, extensive spatial datasets, and a high-resolution atmospheric transport model. Results show that current inventories from the US Environmental Protection Agency (EPA) and the Emissions Database for Global Atmospheric Research underestimate methane emissions nationally by a factor of ∼1.5 and ∼1.7, respectively. Our study indicates that emissions due to ruminants and manure are up to twice the magnitude of existing inventories. In addition, the discrepancy in methane source estimates is particularly pronounced in the south-central United States, where we find total emissions are ∼2.7 times greater than in most inventories and account for 24 ± 3% of national emissions. The spatial patterns of our emission fluxes and observed methane-propane correlations indicate that fossil fuel extraction and refining are major contributors (45 ± 13%) in the south-central United States. This result suggests that regional methane emissions due to fossil fuel extraction and processing could be 4.9 ± 2.6 times larger than in EDGAR, the most comprehensive global methane inventory. These results cast doubt on the US EPA's recent decision to downscale its estimate of national natural gas emissions by 25-30%. Overall, we conclude that methane emissions associated with both the animal husbandry and fossil fuel industries have larger greenhouse gas impacts than indicated by existing inventories.

Record-setting algal bloom in Lake Erie caused by agricultural and meteorological trends consistent with expected future conditions.

In 2011, Lake Erie experienced the largest harmful algal bloom in its recorded history, with a peak intensity over three times greater than any previously observed bloom. Here we show that long-term trends in agricultural practices are consistent with increasing phosphorus loading to the western basin of the lake, and that these trends, coupled with meteorological conditions in spring 2011, produced record-breaking nutrient loads. An extended period of weak lake circulation then led to abnormally long residence times that incubated the bloom, and warm and quiescent conditions after bloom onset allowed algae to remain near the top of the water column and prevented flushing of nutrients from the system. We further find that all of these factors are consistent with expected future conditions. If a scientifically guided management plan to mitigate these impacts is not implemented, we can therefore expect this bloom to be a harbinger of future blooms in Lake Erie.

Global patterns and controls of soil organic carbon dynamics as simulated by multiple terrestrial biosphere models: Current status and future directions.

Soil is the largest organic carbon (C) pool of terrestrial ecosystems, and C loss from soil accounts for a large proportion of land-atmosphere C exchange. Therefore, a small change in soil organic C (SOC) can affect atmospheric carbon dioxide (CO2) concentration and climate change. In the past decades, a wide variety of studies have been conducted to quantify global SOC stocks and soil C exchange with the atmosphere through site measurements, inventories, and empirical/process-based modeling. However, these estimates are highly uncertain, and identifying major driving forces controlling soil C dynamics remains a key research challenge. This study has compiled century-long (1901-2010) estimates of SOC storage and heterotrophic respiration (Rh) from 10 terrestrial biosphere models (TBMs) in the Multi-scale Synthesis and Terrestrial Model Intercomparison Project and two observation-based data sets. The 10 TBM ensemble shows that global SOC estimate ranges from 425 to 2111 Pg C (1 Pg = 1015 g) with a median value of 1158 Pg C in 2010. The models estimate a broad range of Rh from 35 to 69 Pg C yr-1 with a median value of 51 Pg C yr-1 during 2001-2010. The largest uncertainty in SOC stocks exists in the 40-65°N latitude whereas the largest cross-model divergence in Rh are in the tropics. The modeled SOC change during 1901-2010 ranges from -70 Pg C to 86 Pg C, but in some models the SOC change has a different sign from the change of total C stock, implying very different contribution of vegetation and soil pools in determining the terrestrial C budget among models. The model ensemble-estimated mean residence time of SOC shows a reduction of 3.4 years over the past century, which accelerate C cycling through the land biosphere. All the models agreed that climate and land use changes decreased SOC stocks, while elevated atmospheric CO2 and nitrogen deposition over intact ecosystems increased SOC stocks-even though the responses varied significantly among models. Model representations of temperature and moisture sensitivity, nutrient limitation, and land use partially explain the divergent estimates of global SOC stocks and soil C fluxes in this study. In addition, a major source of systematic error in model estimations relates to nonmodeled SOC storage in wetlands and peatlands, as well as to old C storage in deep soil layers.

Assessing biophysical controls on Gulf of Mexico hypoxia through probabilistic modeling.

A mechanistic model was developed to predict midsummer bottom-water dissolved oxygen (BWDO) concentration and hypoxic area on the Louisiana shelf of the northern Gulf of Mexico, USA (1985-2011). Because of its parsimonious formulation, the model possesses many of the benefits of simpler, more empirical models, in that it is computationally efficient and can rigorously account for uncertainty through Bayesian inference. At the same time, the model incorporates important biophysical processes such that its parameterization can be informed by field-measured biological and physical rates. The model is used to explore how freshwater flow, nutrient load, benthic oxygen demand, and wind velocity affect hypoxia on the western and eastern sections of the shelf, delineated by the Atchafalaya River outfall. The model explains over 70% of the variability in BWDO on both shelf sections, and outperforms linear regression models developed from the same input variables. Model results suggest that physical factors (i.e., wind and flow) control a larger portion of the year-to-year variability in hypoxia than previously thought, especially on the western shelf, though seasonal nutrient loads remain an important driver of hypoxia, as well. Unlike several previous Gulf hypoxia modeling studies, results do not indicate a temporal shift in the system's propensity for hypoxia formation (i.e., no regime change). Results do indicate that benthic oxygen demand is a substantial BWDO sink, and a better understanding of the long-term dynamics of this sink is required to better predict how the size of the hypoxic zone will respond to proposed reductions in nutrient loading.

Record-breaking Lake Erie hypoxia during 2012 drought.

Hypoxia has been observed in the central basin of Lake Erie for decades. To understand the impact of various controlling factors, we analyze a record of hypoxic extents for Lake Erie for 1985­2012 and develop a parsimonious model of their interannual variability. We find that the 2012 North American drought and accompanying low tributary discharge was associated with a record-breaking hypoxic event in Lake Erie, whereas a record-setting harmful algal bloom in 2011 was likely associated with only mild hypoxia. River discharge and the timing of nutrient input therefore impact western basin bloom growth and central basin oxygen demand in distinct ways that merit further investigation. Overall, April to June tributary discharge, May to July soluble reactive phosphorus loading, July wind stress, and June northwesterly wind duration explain 82% of the interannual variability of hypoxia, and discharge alone explains 39%, indicating that meteorological factors need to be considered in the development of nutrient management strategies, especially as both extreme precipitation events and droughts become more frequent under a changing climate.

Detecting fossil fuel emissions patterns from subcontinental regions using North American in situ CO2 measurements.

The ability to monitor fossil fuel carbon dioxide (FFCO2) emissions from subcontinental regions using atmospheric CO2 observations remains an important but unrealized goal. Here we explore a necessary but not sufficient component of this goal, namely, the basic question of the detectability of FFCO2 emissions from subcontinental regions. Detectability is evaluated by examining the degree to which FFCO2 emissions patterns from specific regions are needed to explain the variability observed in high-frequency atmospheric CO2 observations. Analyses using a CO2 monitoring network of 35 continuous measurement towers over North America show that FFCO2 emissions are difficult to detect during nonwinter months. We find that the compounding effects of the seasonality of atmospheric transport patterns and the biospheric CO2 flux signal dramatically hamper the detectability of FFCO2 emissions. Results from several synthetic data case studies highlight the need for advancements in data coverage and transport model accuracy if the goal of atmospheric measurement-based FFCO2 emissions detection and estimation is to be achieved beyond urban scales. KEY POINTS: Poor detectability of fossil fuel CO2 emissions from subcontinental regionsDetectability assessed via attribution of emissions patterns in atmospheric dataLoss in detectability due to transport modeling errors and biospheric signal.