Theodore J. Bohn

Multimodel Ensemble Reconstruction of Drought over the Continental United States

Abstract Retrospectively simulated soil moisture from an ensemble of six land surface/hydrological models was used to reconstruct drought events over the continental United States for the period 1920–2003. The simulations were performed at one-half-degree spatial resolution, using a common set of atmospheric forcing data and model-specific soil and vegetation parameters. Monthly simulated soil moisture was converted to percentiles using Weibull plotting position statistics, and the percentiles were then used to represent drought severities and durations. An ensemble method, based on an inverse mapping of the average of the individual model’s soil moisture percentiles, was also used to combine all models’ simulations. Major results are 1) all models and the ensemble reconstruct the known severe drought events during the last century. The spatial extents and severities of drought are plausible for the individual models although substantial among-model disparities exist. 2) The simulations are in more agreem...

Methane emissions from western Siberian wetlands: heterogeneity and sensitivity to climate change

The prediction of methane emissions from high-latitude wetlands is important given concerns about their sensitivity to a warming climate. As a basis for the prediction of wetland methane emissions at regional scales, we coupled the variable infiltration capacity macroscale hydrological model (VIC) with the biosphere–energy-transfer–hydrology terrestrial ecosystem model (BETHY) and a wetland methane emissions model to make large-scale estimates of methane emissions as a function of soil temperature, water table depth, and net primary productivity (NPP), with a parameterization of the sub-grid heterogeneity of the water table depth based on TOPMODEL. We simulated the methane emissions from a 100 km × 100 km region of western Siberia surrounding the Bakchar Bog, for a retrospective baseline period of 1980–1999 and have evaluated their sensitivity to increases in temperature of 0–5 °C and increases in precipitation of 0–15%. The interactions of temperature and precipitation, through their effects on the water table depth, played an important role in determining methane emissions from these wetlands. The balance between these effects varied spatially, and their net effect depended in part on sub-grid topographic heterogeneity. Higher temperatures alone increased methane production in saturated areas, but caused those saturated areas to shrink in extent, resulting in a net reduction in methane emissions. Higher precipitation alone raised water tables and expanded the saturated area, resulting in a net increase in methane emissions. Combining a temperature increase of 3 °C and an increase of 10% in precipitation to represent climate conditions that may pertain in western Siberia at the end of this century resulted in roughly a doubling in annual emissions.

A simplified, data-constrained approach to estimate the permafrost carbon–climate feedback

We present an approach to estimate the feedback from large-scale thawing of permafrost soils using a simplified, data-constrained model that combines three elements: soil carbon (C) maps and profiles to identify the distribution and type of C in permafrost soils; incubation experiments to quantify the rates of C lost after thaw; and models of soil thermal dynamics in response to climate warming. We call the approach the Permafrost Carbon Network Incubation–Panarctic Thermal scaling approach (PInc-PanTher). The approach assumes that C stocks do not decompose at all when frozen, but once thawed follow set decomposition trajectories as a function of soil temperature. The trajectories are determined according to a three-pool decomposition model fitted to incubation data using parameters specific to soil horizon types. We calculate litterfall C inputs required to maintain steady-state C balance for the current climate, and hold those inputs constant. Soil temperatures are taken from the soil thermal modules of ecosystem model simulations forced by a common set of future climate change anomalies under two warming scenarios over the period 2010 to 2100. Under a medium warming scenario (RCP4.5), the approach projects permafrost soil C losses of 12.2–33.4 Pg C; under a high warming scenario (RCP8.5), the approach projects C losses of 27.9–112.6 Pg C. Projected C losses are roughly linearly proportional to global temperature changes across the two scenarios. These results indicate a global sensitivity of frozen soil C to climate change (γ sensitivity) of −14 to −19 Pg C °C−1 on a 100 year time scale. For CH4 emissions, our approach assumes a fixed saturated area and that increases in CH4 emissions are related to increased heterotrophic respiration in anoxic soil, yielding CH4 emission increases of 7% and 35% for the RCP4.5 and RCP8.5 scenarios, respectively, which add an additional greenhouse gas forcing of approximately 10–18%. The simplified approach presented here neglects many important processes that may amplify or mitigate C release from permafrost soils, but serves as a data-constrained estimate on the forced, large-scale permafrost C response to warming.

On the causes of the shrinking of Lake Chad

Over the last 40 years, Lake Chad, once the sixth largest lake in the world, has decreased by more than 90% in area. In this study, we use a hydrological model coupled with a lake/wetland algorithm to simulate the effects of lake bathymetry, human water use, and decadal climate variability on the lakes level, surface area, and water storage. In addition to the effects of persistent droughts and increasing irrigation withdrawals on the shrinking, we find that the lakes unique bathymetry—which allows its division into two smaller lakes—has made it more vulnerable to water loss. Unfortunately the lakes split is favored by the 1952–2006 climatology. Failure of the lake to remerge with renewed rainfall in the 1990s following the drought years of the 1970s and 1980s is a consequence of irrigation withdrawals. Under current climate and water use, a full recovery of the lake is unlikely without an inter-basin water transfer. Breaching the barrier separating the north and south lakes would reduce the amount of supplemental water needed for recovery.

Uncertainties in North American Land Data Assimilation Systems over the Contiguous United States

AbstractThe Environmental Modeling Center (EMC) at the National Centers for Environmental Prediction (NCEP) and the University of Washington (UW) run parallel drought monitoring systems over the continental United States based on the North American Land Data Assimilation System (NLDAS). The NCEP system uses four land surface models (LSMs): Variable Infiltration Capacity (VIC), Noah, Mosaic, and Sacramento (SAC). The UW system uses VIC, SAC, Noah, and the Community Land Model (CLM). An assessment of differences in drought characteristics using both systems for the period 1979–2008 was performed. For soil moisture (SM) percentiles and runoff indices, differences are relatively small among different LSMs in the same system. However, the ensemble mean differences between the two systems are large over the western United States—in some cases exceeding 20% for SM and runoff percentile differences. These differences are most apparent after 2002 when the NCEP system transitioned to use the real-time North America...

Seasonal Hydrologic Forecasting: Do Multimodel Ensemble Averages Always Yield Improvements in Forecast Skill?

Abstract Multimodel techniques have proven useful in improving forecast skill in many applications, including hydrology. Seasonal hydrologic forecasting in large basins represents a special case of hydrologic modeling, in which postprocessing techniques such as temporal aggregation and time-varying bias correction are often employed to improve forecast skill. To investigate the effects that these techniques have on the performance of multimodel averaging, the performance of three hydrological models [Variable Infiltration Capacity, Sacramento/Snow-17, and the Noah land surface model] and two multimodel averages [simple model average (SMA) and multiple linear regression (MLR) with monthly varying model weights] are examined in three snowmelt-dominated basins in the western United States. These evaluations were performed for both simulating and forecasting [using the Ensemble Streamflow Prediction (ESP) method] monthly discharge, with and without monthly bias corrections. The single best bias-corrected mode...

Variability in the sensitivity among model simulations of permafrost and carbon dynamics in the permafrost region between 1960 and 2009

A significant portion of the large amount of carbon (C) currently stored in soils of the permafrost region in the Northern Hemisphere has the potential to be emitted as the greenhouse gases CO2 and CH4 under a warmer climate. In this study we evaluated the variability in the sensitivity of permafrost and C in recent decades among land surface model simulations over the permafrost region between 1960 and 2009. The 15 model simulations all predict a loss of near-surface permafrost (within 3 m) area over the region, but there are large differences in the magnitude of the simulated rates of loss among the models (0.2 to 58.8 × 103 km2 yr−1). Sensitivity simulations indicated that changes in air temperature largely explained changes in permafrost area, although interactions among changes in other environmental variables also played a role. All of the models indicate that both vegetation and soil C storage together have increased by 156 to 954 Tg C yr−1 between 1960 and 2009 over the permafrost region even though model analyses indicate that warming alone would decrease soil C storage. Increases in gross primary production (GPP) largely explain the simulated increases in vegetation and soil C. The sensitivity of GPP to increases in atmospheric CO2 was the dominant cause of increases in GPP across the models, but comparison of simulated GPP trends across the 1982–2009 period with that of a global GPP data set indicates that all of the models overestimate the trend in GPP. Disturbance also appears to be an important factor affecting C storage, as models that consider disturbance had lower increases in C storage than models that did not consider disturbance. To improve the modeling of C in the permafrost region, there is the need for the modeling community to standardize structural representation of permafrost and carbon dynamics among models that are used to evaluate the permafrost C feedback and for the modeling and observational communities to jointly develop data sets and methodologies to more effectively benchmark models. (Less)

Frozen soil degradation and its effects on surface hydrology in the northern Tibetan Plateau

Frozen soil was simulated at six seasonally frozen and seven permafrost stations over the northern Tibetan Plateau using the Variable Infiltration Capacity (VIC) model for the period of 1962–2009. The VIC model resolved the seasonal cycle and temporal evolution of the observed soil temperatures and liquid soil moisture well. The simulated long-term changes during 1962–2009 indicated mostly positive trends for both soil temperature and soil moisture, and negative trends for soil ice content at annual and monthly time scales, although differences existed among the stations, soil layers, and seasons. Increases in soil temperature were due mainly to increases in daily air temperature maxima and internal soil heat conduction, while decreases in soil ice content were related to the warming of frozen soil. For liquid soil moisture, increases in the cold months can be attributed to increases in soil temperature and enhanced soil ice melt while changes in the warm months were the results of competition between positive precipitation and negative soil temperature effects. Precipitation and liquid soil moisture were strongly correlated with evapotranspiration and runoff but had various degrees of correlations with base flow during May–September. Seasonally frozen stations displayed longer and more active hydrological processes than permafrost stations. Slight enhancement of the surface hydrological processes at the study stations was indicated, due to the combined effects of precipitation changes, which were dominant, and frozen soil degradation.

Development and Evaluation of a Multi-Year Fractional Surface Water Data Set Derived from Active/Passive Microwave Remote Sensing Data

The sensitivity of Earth’s wetlands to observed shifts in global precipitation and temperature patterns and their ability to produce large quantities of methane gas are key global change questions. We present a microwave satellite-based approach for mapping fractional surface water (FW) globally at 25-km resolution. The approach employs a land cover-supported, atmospherically-corrected dynamic mixture model applied to 20+ years (1992–2013) of combined, daily, passive/active microwave remote sensing data. The resulting product, known as Surface WAter Microwave Product Series (SWAMPS), shows strong microwave sensitivity to sub-grid scale open water and inundated wetlands comprising open plant canopies. SWAMPS’ FW compares favorably (R2 = 91%–94%) with higher-resolution, global-scale maps of open water from MODIS and SRTM-MOD44W. Correspondence of SWAMPS with open water and wetland products from satellite SAR in Alaska and the Amazon deteriorates when exposed wetlands or inundated forests captured by the SAR products were added to the open water fraction reflecting SWAMPS’ inability to detect water underneath the soil surface or beneath closed forest canopies. Except for a brief period of drying during the first 4 years of observation, the inundation extent for the global domain excluding the coast was largely stable. Regionally, inundation in North America is advancing while inundation is on the retreat in Tropical Africa and North Eurasia. SWAMPS provides a consistent and long-term global record of daily FW dynamics, with documented accuracies suitable for hydrologic assessment and global change-related investigations.

Terrestrial ecosystem model performance in simulating productivity and its vulnerability to climate change in the northern permafrost region

Realistic projection of future climate-carbon (C) cycle feedbacks requires better understanding and an improved representation of the C cycle in permafrost regions in the current generation of Earth system models. Here we evaluated 10 terrestrial ecosystem models for their estimates of net primary productivity (NPP) and responses to historical climate change in permafrost regions in the Northern Hemisphere. In comparison with the satellite estimate from the Moderate Resolution Imaging Spectroradiometer (MODIS; 246±6gCm-2yr-1), most models produced higher NPP (309±12gCm-2yr-1) over the permafrost region during 2000-2009. By comparing the simulated gross primary productivity (GPP) with a flux tower-based database, we found that although mean GPP among the models was only overestimated by 10% over 1982-2009, there was a twofold discrepancy among models (380 to 800gCm-2yr-1), which mainly resulted from differences in simulated maximum monthly GPP (GPPmax). Most models overestimated C use efficiency (CUE) as compared to observations at both regional and site levels. Further analysis shows that model variability of GPP and CUE are nonlinearly correlated to variability in specific leaf area and the maximum rate of carboxylation by the enzyme Rubisco at 25°C (Vcmax_25), respectively. The models also varied in their sensitivities of NPP, GPP, and CUE to historical changes in climate and atmospheric CO2 concentration. These results indicate that model predictive ability of the C cycle in permafrost regions can be improved by better representation of the processes controlling CUE and GPPmax as well as their sensitivity to climate change. (Less)