Anna M. Michalak

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.

Anthropogenic emissions of methane in the United States

Significance Successful regulation of greenhouse gas emissions requires knowledge of current methane emission sources. Existing state regulations in California and Massachusetts require ∼15% greenhouse gas emissions reductions from current levels by 2020. However, government estimates for total US methane emissions may be biased by 50%, and estimates of individual source sectors are even more uncertain. This study uses atmospheric methane observations to reduce this level of uncertainty. We find greenhouse gas emissions from agriculture and fossil fuel extraction and processing (i.e., oil and/or natural gas) are likely a factor of two or greater than cited in existing studies. Effective national and state greenhouse gas reduction strategies may be difficult to develop without appropriate estimates of methane emissions from these source sectors. 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.

Precision requirements for space-based XCO 2 data

Precision requirements are determined for space-based column-averaged CO_2 dry air mole fraction (X_(CO)_2) data. These requirements result from an assessment of spatial and temporal gradients in (X_(CO)_2) the relationship between (X_(CO)_2) precision and surface CO_2 flux uncertainties inferred from inversions of the (X_(CO)_2) data, and the effects of (X_(CO)_2) biases on the fidelity of CO_2 flux inversions. Observational system simulation experiments and synthesis inversion modeling demonstrate that the Orbiting Carbon Observatory mission design and sampling strategy provide the means to achieve these (X_(CO)_2) data precision requirements.

Inverse modeling estimates of the global nitrous oxide surface flux from 1998-2001

Northern Land, and Northern Oceans). We found that compared to our a priori estimate (from the International Geosphere-Biosphere Programme’s Global Emissions Inventory Activity), the a posteriori flux was much lower from 90� S–30� S and substantially higher from equator to 30� N. Consistent with these results, the a posteriori flux from the Southern Oceans region was lower than the a priori estimate, while Tropical Land and Tropical Ocean estimates were higher. The ratio of Northern Hemisphere to Southern Hemisphere fluxes was found to range from 1.9 to 5.2 (depending on the model setup), which is higher than the a priori ratio (1.5) and at the high end of previous estimates. Globally, ocean emissions contributed 26–36% of the total flux (again depending on the model setup), consistent with the a priori estimate (29%), though somewhat higher than some other previous estimates.

Maximum likelihood estimation of covariance parameters for Bayesian atmospheric trace gas surface flux inversions

[1] This paper introduces a Maximum Likelihood (ML) approach for estimating the statistical parameters required for the covariance matrices used in the solution of Bayesian inverse problems aimed at estimating surface fluxes of atmospheric trace gases. The method offers an objective methodology for populating the covariance matrices required in Bayesian inversions, thereby resulting in better estimates of the uncertainty associated with derived fluxes and minimizing the risk of inversions being biased by unrealistic covariance parameters. In addition, a method is presented for estimating the uncertainty associated with these covariance parameters. The ML method is demonstrated using a typical inversion setup with 22 flux regions and 75 observation stations from the National Oceanic and Atmospheric Administration-Climate Monitoring and Diagnostics Laboratory (NOAA-CMDL) Cooperative Air Sampling Network with available monthly averaged carbon dioxide data. Flux regions and observation locations are binned according to various characteristics, and the variances of the model-data mismatch and of the errors associated with the a priori flux distribution are estimated from the available data.

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.

A method for enforcing parameter nonnegativity in Bayesian inverse problems with an application to contaminant source identification

Received 28 May 2002; revised 18 September 2002; accepted 24 September 2002; published 14 February 2003. [1] When an inverse problem is solved to estimate an unknown function such as the hydraulic conductivity in an aquifer or the contamination history at a site, one constraint is that the unknown function is known to be everywhere nonnegative. In this work, we develop a statistically rigorous method for enforcing function nonnegativity in Bayesian inverse problems. The proposed method behaves similarly to a Gaussian process with a linear variogram (i.e., unrestricted Brownian motion) for parameter values significantly greater than zero. The method uses the method of images to define a prior probability density function based on reflected Brownian motion that implicitly enforces nonnegativity. This work focuses on problems where the unknown is a function of a single variable (e.g., time). A Markov chain Monte Carlo (MCMC) method, specifically, a highly efficient Gibbs sampler, is implemented to generate conditional realizations of the unknown function. The new method is applied to the estimation of the trichloroethylene (TCE) and perchloroethylene (PCE) contamination history in an aquifer at Dover Air Force Base, Delaware, based on concentration profiles obtained from an underlying aquitard. INDEX TERMS: 1831 Hydrology: Groundwater quality; 1869 Hydrology: Stochastic processes; 3260 Mathematical Geophysics: Inverse theory; KEYWORDS: stochastic inverse modeling, contaminant source identification, inference under constraints, Markov chain Monte Carlo (MCMC), Gibbs sampling, Bayesian inference

Observed suppression of ozone formation at extremely high temperatures due to chemical and biophysical feedbacks

Ground level ozone concentrations ([O3]) typically show a direct linear relationship with surface air temperature. Three decades of California measurements provide evidence of a statistically significant change in the ozone-temperature slope (ΔmO3-T) under extremely high temperatures (> 312 K). This ΔmO3-T leads to a plateau or decrease in [O3], reflecting the diminished role of nitrogen oxide sequestration by peroxyacetyl nitrates and reduced biogenic isoprene emissions at high temperatures. Despite inclusion of these processes in global and regional chemistry-climate models, a statistically significant change in ΔmO3-T has not been noted in prior studies. Future climate projections suggest a more frequent and spatially widespread occurrence of this ΔmO3-T response, confounding predictions of extreme ozone events based on the historically observed linear relationship.

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

Abstract 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.

Spatial and Temporal Trends in Lake Erie Hypoxia, 1987−2007

Hypoxic conditions, defined as dissolved oxygen (DO) concentrations below 2 mg/L, are a regular summertime occurrence in Lake Erie, but the spatial extent has been poorly understood due to sparse sampling. We use geostatistical kriging and conditional realizations to provide quantitative estimates of the extent of hypoxia in the central basin of Lake Erie for August and September of 1987 to 2007, along with their associated uncertainties. The applied geostatistical approach combines the limited in situ DO measurements with auxiliary data selected using the Bayesian Information Criterion. Bathymetry and longitude are found to be highly significant in explaining the spatial distribution of DO, while satellite observations of sea surface temperature and satellite chlorophyll are not. The hypoxic extent was generally lowest in the mid-1990s, with the late 1980s (1987, 1988) and the 2000s (2003, 2005) experiencing the largest hypoxic zones. A simple exponential relationship based on the squared average measured bottom DO explains 97% of the estimated variability in the hypoxic extent. The change in the observed maximum extent between August and September is found to be sensitive to the corresponding variability in the hypolimnion thickness.