Bibliography - Elena Shevliakova

1. Chaney, Nathaniel W., M H J Van Huijgevoort, Elena Shevliakova, Sergey Malyshev, P C D Milly, Paul P G Gauthier, and Benjamin N Sulman, June 2018: Harnessing Big Data to Rethink Land Heterogeneity in Earth System Models. Hydrology and Earth System Sciences, 22(6), DOI:10.5194/hess-22-3311-2018 .
   [ Abstract ]

   The continual growth in the availability, detail, and wealth of environmental data provides an invaluable asset to improve the characterization of land heterogeneity in Earth System models – a persistent challenge in macroscale models. However, due to the nature of these data (volume and complexity) and the computational constraints of macroscale models, until now these data have been underutilized for global applications. As a proof of concept, this study explores over a 1/4 degree (~ 25 km) grid cell in southeastern California how to effectively and efficiently harness these data in Earth System models. First, a novel hierarchical multivariate clustering approach (HMC) is used to summarize the high dimensional environmental data space into hydrologically interconnected representative clusters (i.e., tiles). These tiles and their associated properties are then used to parameterize the sub-grid heterogeneity of the Geophysical Fluid Dynamics Laboratory (GFDL) LM4-HB land model. To assess how this data-driven approach to assemble the model tiles impacts the simulated water, energy, and carbon cycles, model experiments are run using a series of different tile configurations assembled by HMC. The results over the 1/4 degree macroscale grid cell and the underlying 30-meter fine-scale grid in southeastern California show that: 1) the observed similarity over the landscape makes it possible to robustly account for the role of multi-scale heterogeneity in the macroscale states and fluxes with around 300 sub-grid land model tiles; 2) assembling the sub-grid tiles from observed data, at times, leads to noticeable differences in the macroscale water, energy, and carbon cycles; for example, explicit subsurface interactions between the tiles leads to a dampening of macroscale extremes; 3) connecting the fine-scale grid to the model tiles via HMC enables circumventing the classic scale discrepancies between the macroscale and field-scale estimates; this has potentially significant implications for the evaluation and application of Earth System models.

2. Kapnick, Sarah B., Xiaosong Yang, Gabriel A Vecchi, Thomas L Delworth, Richard G Gudgel, Sergey Malyshev, P C D Milly, Elena Shevliakova, Seth D Underwood, and S A Margulis, February 2018: Potential for western US seasonal snowpack prediction. Proceedings of the National Academy of Sciences, 115(6), DOI:10.1073/pnas.1716760115 .
   [ Abstract ]

   Mountain snowpack in the western United States provides a natural reservoir for cold season precipitation; variations in snowpack influence warm season water supply, wildfire risk, ecology, and industries like agriculture dependent on snow and downstream water availability. Efforts to understand snowpack variability have predominantly been focused on either weekly (weather) or decadal to centennial (climate variability and change) timescales. We focus on a timescale between these ranges by demonstrating that a global climate model suite can provide snowpack predictions 8 months in advance. The predictions from climate models outperform statistical methods from observations alone. Our results show that seasonal hydroclimate predictions are possible and highlight areas for future prediction system improvements.

3. Paulot, Fabien, Sergey Malyshev, T B Nguyen, J D Crounse, Elena Shevliakova, and Larry W Horowitz, in press: Representing sub-grid scale variations in nitrogen deposition associated with land use in a global Earth System Model: implications for present and future nitrogen deposition fluxes over North America. Atmospheric Chemistry and Physics Discussions. DOI:10.5194/acp-2018-572. July 2018.
   [ Abstract ]

   Reactive nitrogen (N) emissions have increased over the last 150 years as a result of greater fossil fuel combustion and food production. The resulting increase in N deposition can alter the function of ecosystems, but characterizing its ecological impacts remains challenging, in part because of uncertainties in model-based estimates of N dry deposition. Here, we leverage the tiled structure of the land component (LM3) of the Geophysical Fluid Dynamics Laboratory (GFDL) Earth System Model to represent the impact of physical, hydrological, and ecological heterogeneities on the surface removal of chemical tracers. We show that this framework can be used to estimate N deposition at more ecologically-relevant scales (e.g., natural vegetation, water bodies) than from the coarse-resolution global chemistry–climate model (GFDL-AM3). Focusing on North America, we show that the faster removal of N over forested ecosystems relative to cropland and pasture implies that coarse resolution estimates of N deposition from global models systematically underestimate N deposition to natural vegetation by 10 to 30% in the Central and Eastern US. Neglecting the subgrid scale heterogeneity of dry deposition velocities also results in an underestimate (overestimate) of the amount of reduced (oxidized) nitrogen deposited to water bodies. Overall, changes in land cover associated with human activities are found to slow down the removal of N from the atmosphere, causing a reduction in the dry oxidized, dry reduced, and total N deposition over the contiguous US of 8%, 26%, and 6%, respectively. We also find that the reduction in the overall rate of removal of N associated with land-use change tends to increase N deposition on the remaining natural vegetation and facilitate N export to Canada. We show that subgrid scale differences in the surface removal of oxidized and reduced nitrogen imply that near-term (2010–2050) changes in oxidized (−47%) and reduced (+40%) US N emissions will cause opposite changes in N deposition to water bodies (increase) and natural vegetation (decrease) in the Eastern US, with potential implications for acidification and ecosystems.

4. Rabin, S, Sergey Malyshev, B I Magi, Elena Shevliakova, and S W Pacala, March 2018: A fire model with distinct crop, pasture, and non-agricultural burning: Use of new data and a model-fitting algorithm for FINALv1. Geoscientific Model Development, 11(2), DOI:10.5194/gmd-11-815-2018 .
   [ Abstract ]

   This study describes and evaluates the Fire Including Natural & Agricultural Lands model (FINAL) which, for the first time, explicitly simulates cropland and pasture management fires separately from non-agricultural fires. The non-agricultural fire module uses empirical relationships to simulate burned area in a quasi-mechanistic framework, similar to past fire modeling efforts, but with a novel optimization method that improves the fidelity of simulated fire patterns to new observational estimates of non-agricultural burning. The agricultural fire components are forced with estimates of cropland and pasture fire seasonality and frequency derived from observational land-cover and satellite fire datasets. FINAL accurately simulates the amount, distribution, and seasonal timing of burned cropland and pasture over 2001–2009 (global totals: 0.434 × 106 and 2.02 × 106 km2 yr−1 modeled, 0.454 × 106 and 2.04 × 106 km2 yr−1 observed), but carbon emissions for cropland and pasture fire are overestimated (global totals: 0.297 PgC yr−1 and 0.712 PgC yr−1 modeled, 0.194 PgC yr−1 and 0.538 PgC yr−1 observed). The non-agricultural fire module underestimates global burned area (1.66 × 106 km2 yr−1 modeled, 2.44 × 106 km2 yr−1 observed) and carbon emissions (1.33 PgC yr−1 modeled, 1.84 PgC yr−1 observed). The spatial pattern of total burned area and carbon emissions is generally well reproduced across much of sub-Saharan Africa, Brazil, central Asia, and Australia, whereas the boreal zone suffers from underestimates. FINAL represents an important step in the development of global fire models, and offers a strategy for fire models to consider human-driven fire regimes on cultivated lands. At the regional scale, simulations would benefit from refinements in the parameterizations and improved optimization datasets.

5. Ward, Daniel S., Elena Shevliakova, Sergey Malyshev, and S Rabin, January 2018: Trends and variability of global fire emissions due to historical anthropogenic activities. Global Biogeochemical Cycles, 32(1), DOI:10.1002/2017GB005787 .
   [ Abstract ]

   Globally, fires are a major source of carbon from the terrestrial biosphere to the atmosphere, occurring on a seasonal cycle and with substantial interannual variability. To understand past trends and variability in sources and sinks of terrestrial carbon, we need quantitative estimates of global fire distributions. Here we introduce an updated version of the Fire Including Natural and Agricultural Lands model, version 2 (FINAL.2), modified to include multi-day burning and enhanced fire spread rate in forest crowns. We demonstrate that the improved model reproduces the interannual variability and spatial distribution of fire emissions reported in present day remotely sensed inventories. We use FINAL.2 to simulate historical (post-1700) fires and attribute past fire trends and variability to individual drivers: land use and land cover change, population growth, and lightning variability. Global fire emissions of carbon increase by about 10% between 1700 and 1900, reaching a maximum of 3.4 PgC yr-1 in the 1910s, followed by a decrease to about 5% below year 1700 levels by 2010. The decrease in emissions from the 1910s to the present day is driven mainly by land use change, with a smaller contribution from increased fire suppression due to increased human population, and is largest in Sub-Saharan Africa and South Asia. Interannual variability of global fire emissions is similar in the present day as in the early historical period, but present day wildfires would be more variable in the absence of land use change.

6. Zhao, Ming, J-C Golaz, Isaac M Held, Huan Guo, V Balaji, Rusty Benson, Jan-Huey Chen, Xi Chen, Leo J Donner, John P Dunne, Krista A Dunne, J W Durachta, Songmiao Fan, Stuart Freidenreich, Stephen T Garner, Paul Ginoux, Lucas Harris, Larry W Horowitz, John P Krasting, Amy R Langenhorst, Zhi Liang, Pu Lin, Shian-Jiann Lin, Sergey Malyshev, E Mason, P C D Milly, Yi Ming, Vaishali Naik, Fabien Paulot, David J Paynter, Peter Phillipps, Aparna Radhakrishnan, V Ramaswamy, Tom Robinson, M Daniel Schwarzkopf, Charles J Seman, Elena Shevliakova, Zhaoyi Shen, Hyeyum Hailey Shin, Levi G Silvers, R John Wilson, Michael Winton, Andrew T Wittenberg, Bruce Wyman, and Baoqiang Xiang, March 2018: The GFDL Global Atmosphere and Land Model AM4.0/LM4.0 - Part I: Simulation Characteristics with Prescribed SSTs. Journal of Advances in Modeling Earth Systems, 10(3), DOI:10.1002/2017MS001208 .
   [ Abstract ]

   In this two-part paper, a description is provided of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). This version, with roughly 100km horizontal resolution and 33 levels in the vertical, contains an aerosol model that generates aerosol fields from emissions and a “light” chemistry mechanism designed to support the aerosol model but with prescribed ozone. In Part I, the quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode – with prescribed sea surface temperatures (SSTs) and sea ice distribution – is described and compared with previous GFDL models and with the CMIP5 archive of AMIP simulations. The model's Cess sensitivity (response in the top-of-atmosphere radiative flux to uniform warming of SSTs) and effective radiative forcing are also presented. In Part II, the model formulation is described more fully and key sensitivities to aspects of the model formulation are discussed, along with the approach to model tuning.

7. Zhao, Ming, J-C Golaz, Isaac M Held, Huan Guo, V Balaji, Rusty Benson, Jan-Huey Chen, Xi Chen, Leo J Donner, John P Dunne, Krista A Dunne, J W Durachta, Songmiao Fan, Stuart Freidenreich, Stephen T Garner, Paul Ginoux, Lucas Harris, Larry W Horowitz, John P Krasting, Amy R Langenhorst, Zhi Liang, Pu Lin, Shian-Jiann Lin, Sergey Malyshev, E Mason, P C D Milly, Yi Ming, Vaishali Naik, Fabien Paulot, David J Paynter, Peter Phillipps, Aparna Radhakrishnan, V Ramaswamy, Tom Robinson, M Daniel Schwarzkopf, Charles J Seman, Elena Shevliakova, Zhaoyi Shen, Hyeyum Hailey Shin, Levi G Silvers, R John Wilson, Michael Winton, Andrew T Wittenberg, Bruce Wyman, and Baoqiang Xiang, March 2018: The GFDL Global Atmosphere and Land Model AM4.0/LM4.0 - Part II: Model Description, Sensitivity Studies, and Tuning Strategies. Journal of Advances in Modeling Earth Systems, 10(3), DOI:10.1002/2017MS001209 .
   [ Abstract ]

   In Part II of this two-part paper, documentation is provided of key aspects of a version of the AM4.0/LM4.0 atmosphere/land model that will serve as a base for a new set of climate and Earth system models (CM4 and ESM4) under development at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL). The quality of the simulation in AMIP (Atmospheric Model Intercomparison Project) mode has been provided in Part I. Part II provides documentation of key components and some sensitivities to choices of model formulation and values of parameters, highlighting the convection parameterization and orographic gravity wave drag. The approach taken to tune the model's clouds to observations is a particular focal point. Care is taken to describe the extent to which aerosol effective forcing and Cess sensitivity have been tuned through the model development process, both of which are relevant to the ability of the model to simulate the evolution of temperatures over the last century when coupled to an ocean model.

8. Ballantyne, A P., William Smith, W Anderegg, P Kauppi, Jorge L Sarmiento, P P Tans, and Elena Shevliakova, et al., February 2017: Accelerating net terrestrial carbon uptake during the warming hiatus due to reduced respiration. Nature Climate Change, 7(2), DOI:10.1038/nclimate3204 .
   [ Abstract ]

   he recent ‘warming hiatus’ presents an excellent opportunity to investigate climate sensitivity of carbon cycle processes. Here we combine satellite and atmospheric observations to show that the rate of net biome productivity (NBP) has significantly accelerated from −0.007 ± 0.065 PgC yr−2 over the warming period (1982 to 1998) to 0.119 ± 0.071 PgC yr−2 over the warming hiatus (1998–2012). This acceleration in NBP is not due to increased primary productivity, but rather reduced respiration that is correlated (r = 0.58; P = 0.0007) and sensitive (γ = 4.05 to 9.40 PgC yr−1 per °C) to land temperatures. Global land models do not fully capture this apparent reduced respiration over the warming hiatus; however, an empirical model including soil temperature and moisture observations better captures the reduced respiration.

9. Clifton, O E., Arlene M Fiore, C E Morris, Sergey Malyshev, Larry W Horowitz, Elena Shevliakova, and Fabien Paulot, et al., January 2017: Interannual variability in ozone removal by a temperate deciduous forest. Geophysical Research Letters, 44(1), DOI:10.1002/2016GL070923 .
   [ Abstract ]

   The ozone (O3) dry depositional sink and its contribution to observed variability in tropospheric O3 are both poorly understood. Distinguishing O3 uptake through plant stomata versus other pathways is relevant for quantifying the O3 influence on carbon and water cycles. We use a decade of O3, carbon, and energy eddy covariance (EC) fluxes at Harvard Forest to investigate interannual variability (IAV) in O3 deposition velocities ( math formula). In each month, monthly mean math formula for the highest year is twice that for the lowest. Two independent stomatal conductance estimates, based on either water vapor EC or gross primary productivity, vary little from year to year relative to canopy conductance. We conclude that nonstomatal deposition controls the substantial observed IAV in summertime math formula during the 1990s over this deciduous forest. The absence of obvious relationships between meteorology and math formula implies a need for additional long-term, high-quality measurements and further investigation of nonstomatal mechanisms.

10. Findell, Kirsten L., Alexis Berg, P Gentine, John P Krasting, B R Lintner, Sergey Malyshev, J A Santanello, and Elena Shevliakova, October 2017: The impact of anthropogenic land use and landcover change on regional climate extremes. Nature Communications, 8, 989, DOI:10.1038/s41467-017-01038-w .
    [ Abstract ]

    Land surface processes modulate the severity of heat waves, droughts, and other extreme events. However, models show contrasting effects of land surface changes on extreme temperatures. Here, we use an earth system model from the Geophysical Fluid Dynamics Laboratory to investigate regional impacts of land use and land cover change on combined extremes of temperature and humidity, namely aridity and moist enthalpy, quantities central to human physiological experience of near-surface climate. The model’s near-surface temperature response to deforestation is consistent with recent observations, and conversion of mid-latitude natural forests to cropland and pastures is accompanied by an increase in the occurrence of hot-dry summers from once-in-a-decade to every 2–3 years. In the tropics,long time-scale oceanic variability precludes determination of how much of a small, but significant, increase in moist enthalpy throughout the year stems from the model’s novel representation of historical patterns of wood harvesting, shifting cultivation, and regrowth of secondary vegetation and how much is forced by internal variability within the tropical oceans.

11. Metcalf, C J., K S Walter, A Wesolowski, C O Buckee, and Elena Shevliakova, et al., August 2017: Identifying climate drivers of infectious disease dynamics: recent advances and challenges ahead. Proceedings of the Royal Society B, 284(1860), DOI:10.1098/rspb.2017.0901 .
    [ Abstract ]

    Climate change is likely to profoundly modulate the burden of infectious diseases. However, attributing health impacts to a changing climate requires being able to associate changes in infectious disease incidence with the potentially complex influences of climate. This aim is further complicated by nonlinear feedbacks inherent in the dynamics of many infections, driven by the processes of immunity and transmission. Here, we detail the mechanisms by which climate drivers can shape infectious disease incidence, from direct effects on vector life history to indirect effects on human susceptibility, and detail the scope of variation available with which to probe these mechanisms. We review approaches used to evaluate and quantify associations between climate and infectious disease incidence, discuss the array of data available to tackle this question, and detail remaining challenges in understanding the implications of climate change for infectious disease incidence. We point to areas where synthesis between approaches used in climate science and infectious disease biology provide potential for progress.

12. Sulman, Benjamin N., E R Brzostek, C Medici, and Elena Shevliakova, et al., August 2017: Feedbacks between plant N demand and rhizosphere priming depend on type of mycorrhizal association. Ecology Letters, 20(8), DOI:10.1111/ele.12802 .
    [ Abstract ]

    Ecosystem carbon (C) balance is hypothesised to be sensitive to the mycorrhizal strategies that plants use to acquire nutrients. To test this idea, we coupled an optimality-based plant nitrogen (N) acquisition model with a microbe-focused soil organic matter (SOM) model. The model accurately predicted rhizosphere processes and C–N dynamics across a gradient of stands varying in their relative abundance of arbuscular mycorrhizal (AM) and ectomycorrhizal (ECM) trees. When mycorrhizal dominance was switched – ECM trees dominating plots previously occupied by AM trees, and vice versa – legacy effects were apparent, with consequences for both C and N stocks in soil. Under elevated productivity, ECM trees enhanced decomposition more than AM trees via microbial priming of unprotected SOM. Collectively, our results show that ecosystem responses to global change may hinge on the balance between rhizosphere priming and SOM protection, and highlight the importance of dynamically linking plants and microbes in terrestrial biosphere models.

13. Evans, Stuart, Paul Ginoux, Sergey Malyshev, and Elena Shevliakova, November 2016: Climate-vegetation interaction and amplification of Australian dust variability. Geophysical Research Letters, 43(22), DOI:10.1002/2016GL071016 .
    [ Abstract ]

    Observations show that Australian dust activity varies by a factor of 4 on decadal timescales. General circulation models, however, typically fail to simulate this variability. Here we introduce a new dust parameterization into the NOAA/Geophysical Fluid Dynamics Laboratory climate model CM3 that represents land surface processes controlling dust sources including soil water and ice, snow cover, vegetation characteristics, and land type. In an additional novel step, we couple this new dust parameterization to the dynamic vegetation model LM3. In Australia, the new parameterization amplifies the magnitude and timescale of dust variability and better simulates the El Niño–Southern Oscillation-dust relationship by more than doubling its strength. We attribute these improvements primarily to the slow response time of vegetation to precipitation anomalies and show that vegetation changes account for approximately 50% of enhanced dust emission during El Niño events. The amplified dust leads to radiative forcing over Australia greater than −1 and −20 W/m2 at top of atmosphere and surface, respectively.

14. Kanter, D, X Zhang, D L Mauzerall, Sergey Malyshev, and Elena Shevliakova, September 2016: The importance of climate change and nitrogen use efficiency for future nitrous oxide emissions from agriculture. Environmental Research Letters, 11(9), DOI:10.1088/1748-9326/11/9/094003 .
    [ Abstract ]

    Nitrous oxide (N2O) is an important greenhouse gas and ozone depleting substance. Previous projections of agricultural N2O (the dominant anthropogenic source) show emissions changing in tandem, or at a faster rate than changes in nitrogen (N) consumption. However, recent studies suggest that the carbon dioxide (CO2) fertilization effect may increase plant N uptake, which could decrease soil N losses and dampen increases in N2O. To evaluate this hypothesis at a global scale, we use a process-based land model with a coupled carbon-nitrogen cycle to examine how changes in climatic factors, land-use, and N application rates could affect agricultural N2O emissions by 2050. Assuming little improvement in N use efficiency (NUE), the model projects a 24%–31% increase in global agricultural N2O emissions by 2040–2050 depending on the climate scenario—a relatively moderate increase compared to the projected increases in N inputs (42%–44%) and previously published emissions projections (38%–75%). This occurs largely because the CO2 fertilization effect enhances plant N uptake in several regions, which subsequently dampens N2O emissions. And yet, improvements in NUE could still deliver important environmental benefits by 2050: equivalent to 10 Pg CO2 equivalent and 0.6 Tg ozone depletion potential.

15. Lawrence, D, G C Hurtt, A Arneth, V Brovkin, K V Calvin, A Jones, C Jones, P Lawrence, N de Noblet-Ducoudre, J Pongratz, S I Seneviratne, and Elena Shevliakova, September 2016: The Land Use Model Intercomparison Project (LUMIP): Rationale and experimental design. Geoscientific Model Development, 9(9), DOI:10.5194/gmd-9-2973-2016 .
    [ Abstract ]

    Human land-use activities have resulted in large to the Earth surface, with resulting implications for climate. In the future, land-use activities are likely to expand and intensify further to meet growing demands for food, fiber, and energy. The Land Use Model Intercomparison Project (LUMIP) aims to further advance understanding of the impacts of land-use and land-cover change (LULCC) on climate, specifically addressing the questions: (1) What are the effects of LULCC on climate and biogeochemical cycling (past–future)? (2) What are the impacts of land management on surface fluxes of carbon, water, and energy and are there regional land-management strategies with promise to help mitigate and/or adapt to climate change? In addressing these questions, LUMIP will also address a range of more detailed science questions to get at process-level attribution, uncertainty, data requirements, and other related issues in more depth and sophistication than possible in a multi-model context to date. There will be particular focus on the separation and quantification of the effects on climate from land-use change relative to fossil fuel emissions, separation of biogeochemical from biogeophysical effects of land-use, the unique impacts of land-cover change versus land management change, modulation of land-use impact on climate by land-atmosphere coupling strength, and the extent that impacts of enhanced CO2 concentrations on plant photosynthesis are modulated by past and future land use. LUMIP involves three major sets of science activities: (1) development of an updated and expanded historical and future land-use dataset, (2) an experimental protocol for specific LUMIP experiments for CMIP6, and (3) definition of metrics and diagnostic protocols that quantify model performance, and related sensitivities, with respect to-. In this paper, we describe the LUMIP simulations that will formally be part of CMIP6. These experiments are explicitly designed to be complementary to experiments from the CMIP core, ScenarioMIP, and C4MIP. LUMIP includes a two-phase experimental design. Phase one features idealized coupled and land-only model experiments designed to advance process-level understanding of LULCC impacts on climate, as well as to quantify model sensitivity to potential land-cover and land-use change. Phase two experiments focus on quantification of the historic impact of land use and the potential for future land management decisions to aid in mitigation of climate change. This paper documents these simulations in detail, explains their rationale, outlines plans for analysis, and describes a new subgrid land-use tile data request (primary and secondary land, crops, pasture, urban). It is essential that modeling groups participating in LUMIP adhere to the experimental design as closely as possible.

16. Lee, M, Elena Shevliakova, Sergey Malyshev, P C D Milly, and P R Jaffe, July 2016: Climate variability and extremes, interacting with nitrogen storage, amplify eutrophication risk. Geophysical Research Letters, 43(14), DOI:10.1002/2016GL069254 .
    [ Abstract ]

    Despite 30 years of basin-wide nutrient-reduction efforts, severe hypoxia continues to be observed in the Chesapeake Bay. Here we demonstrate the critical influence of climate variability, interacting with accumulated nitrogen (N) over multidecades, on Susquehanna River dissolved nitrogen (DN) loads, known precursors of the hypoxia in the Bay. We used the process model LM3-TAN (Terrestrial and Aquatic Nitrogen), which is capable of capturing both seasonal and decadal-to-century changes in vegetation-soil-river N storage, and produced nine scenarios of DN-load distributions under different short-term scenarios of climate variability and extremes. We illustrate that after 1 to 3 yearlong dry spells, the likelihood of exceeding a threshold DN load (56 kt yr−1) increases by 40 to 65% due to flushing of N accumulated throughout the dry spells and altered microbial processes. Our analyses suggest that possible future increases in climate variability/extremes—specifically, high precipitation occurring after multiyear dry spells—could likely lead to high DN-load anomalies and hypoxia.

17. Li, Dan, Sergey Malyshev, and Elena Shevliakova, June 2016: Exploring historical and future urban climate in the Earth System Modeling framework: 1. Model development and evaluation. Journal of Advances in Modeling Earth Systems, 8(2), DOI:10.1002/2015MS000578 .
    [ Abstract ]

    A number of recent studies investigated impacts of Land-Use and Land-Cover Changes (LULCC) on climate with global Earth System Models (ESMs). Yet many ESMs are still missing a representation of the most extreme form of natural landscape modification – urban settlements. Moreover, long-term (i.e. decades to century) transitions between build-up and other land cover types due to urbanization and de-urbanization have not been examined in the literature. In this study we evaluate a new urban canopy model (UCM) that characterizes urban physical and biogeochemical processes within the sub-grid tiling framework of the Geophysical Fluid Dynamics Laboratory (GFDL) land model, LM3. The new model LM3-UCM is based on the urban canyon concept and simulates exchange of energy, water (liquid and solid), and carbon between urban land and the atmosphere. LM3-UCM has several unique features, including explicit treatment of vegetation inside the urban canyon and dynamic transition between urban, agricultural and unmanaged tiles. The model is evaluated using observational datasets collected at three urban sites: Marseille in France, Basel in Switzerland and Baltimore in the United States. It is found that LM3-UCM satisfactorily reproduces canyon air temperature, surface temperatures, radiative fluxes, and turbulent heat fluxes at the three urban sites. LM3-UCM can capture urban features in a computationally efficient manner and is incorporated into the land component of GFDL ESMs. This new capability will enable improved understanding of climate change effects on cities and the impacts of urbanization on climate.

18. Li, Dan, Sergey Malyshev, and Elena Shevliakova, June 2016: Exploring historical and future urban climate in the Earth System Modeling framework: 2. Impact of urban land use over the Continental United States. Journal of Advances in Modeling Earth Systems, 8(2), DOI:10.1002/2015MS000579 .
    [ Abstract ]

    Using a newly developed urban canopy model coupled to the Geophysical Fluid Dynamics Laboratory (GFDL) land model LM3 (LM3-UCM), this study examines the urban land use impacts over the Continental United States (CONUS) under the present-day climate and two future scenarios. Using natural (undisturbed) vegetation systems as references where no land use has occurred, the LM3-UCM simulations show that the spatial pattern of summer (June, July, and August) temperature differences between urban and natural vegetation systems is primarily controlled by the spatial pattern of differences in evapotranspiration, which further depends on the spatial distribution of precipitation. The magnitude of temperature differences generally increases as the summer precipitation amount increases and then levels off when the total summer precipitation amount exceeds 400 mm, which is broadly consistent with previous studies but with significant variability. In winter (December, January, February), the magnitude of temperature differences is more controlled by the building heating than the precipitation amount. At high latitudes where snow is an important factor in radiative balance, the magnitude is also affected by a larger net shortwave radiation input for urban areas due to the lower albedo of cities. Although both urban and natural vegetation temperatures increase as the climate warms, their increasing rates are different and hence their differences change with time. It is found that the multi-decadal trend of summer temperature difference is negligible. However, the winter temperature differences show a strong negative trend, which is caused by reduced building heating requirements under a warming climate.

19. Smith, N G., Sergey Malyshev, Elena Shevliakova, J Kattge, and J S Dukes, April 2016: Foliar temperature acclimation reduces simulated carbon sensitivity to climate. Nature Climate Change, 6(4), DOI:10.1038/nclimate2878 .
    [ Abstract ]

    Plant photosynthesis and respiration are the largest carbon fluxes between the terrestrial biosphere and the atmosphere1, and their parameterizations represent large sources of uncertainty in projections of land carbon uptake in Earth system models2, 3 (ESMs). The incorporation of temperature acclimation of photosynthesis and foliar respiration, commonly observed processes, into ESMs has been proposed as a way to reduce this uncertainty2. Here we show that, across 15 flux tower sites spanning multiple biomes at various locations worldwide (10° S–67° N), acclimation parameterizations4, 5 improve a model’s ability to reproduce observed net ecosystem exchange of CO2. This improvement is most notable in tropical biomes, where photosynthetic acclimation increased model performance by 36%. The consequences of acclimation for simulated terrestrial carbon uptake depend on the process, region and time period evaluated. Globally, including acclimation has a net effect of increasing carbon assimilation and storage, an effect that diminishes with time, but persists well into the future. Our results suggest that land models omitting foliar temperature acclimation are likely to overestimate the temperature sensitivity of terrestrial carbon exchange, thus biasing projections of future carbon storage and estimates of policy indicators such as the transient climate response to cumulative carbon emissions.

20. Ward, Daniel S., Elena Shevliakova, Sergey Malyshev, M Kuwata, and Andrew T Wittenberg, December 2016: Variability of fire emissions on interannual to multi-decadal timescales in two Earth System models. Environmental Research Letters, 11(12), DOI:10.1088/1748-9326/11/12/125008 .
    [ Abstract ]

    Connections between wildfires and modes of variability in climate are sought as a means for predicting fire activity on interannual to multi-decadal timescales. Several fire drivers, such as temperature and local drought index, have been shown to vary on these timescales, and analysis of tree-ring data suggests covariance between fires and climate oscillation indices in some regions. However, the shortness of the satellite record of global fire events limits investigations on larger spatial scales. Here we explore the interplay between climate variability and wildfire emissions with the preindustrial long control numerical experiments and historical ensembles of CESM1 and the NOAA/GFDL ESM2Mb. We find that interannual variability in fires is underpredicted in both Earth System models (ESMs) compared to present day fire emission inventories. Modeled fire emissions respond to the El Niño/southern oscillation (ENSO) and Pacific decadal oscillation (PDO) with increases in southeast Asia and boreal North America emissions, and decreases in southern North America and Sahel emissions, during the ENSO warm phase in both ESMs, and the PDO warm phase in CESM1. Additionally, CESM1 produces decreases in boreal northern hemisphere fire emissions for the warm phase of the Atlantic Meridional Oscillation. Through analysis of the long control simulations, we show that the 20th century trends in both ESMs are statistically significant, meaning that the signal of anthropogenic activity on fire emissions over this time period is detectable above the annual to decadal timescale noise. However, the trends simulated by the two ESMs are of opposite sign (CESM1 decreasing, ESM2Mb increasing), highlighting the need for improved understanding, proxy observations, and modeling to resolve this discrepancy.

21. Anderegg, W, C Schwalm, F Biondi, J J Camarero, G Koch, M Litvak, K Ogle, J D Shaw, and Elena Shevliakova, et al., July 2015: Pervasive drought legacies in forest ecosystems and their implications for carbon cycle models. Science, 349(6247), DOI:10.1126/science.aab1833 .
    [ Abstract ]

    The impacts of climate extremes on terrestrial ecosystems are poorly understood but important for predicting carbon cycle feedbacks to climate change. Coupled climate–carbon cycle models typically assume that vegetation recovery from extreme drought is immediate and complete, which conflicts with the understanding of basic plant physiology. We examined the recovery of stem growth in trees after severe drought at 1338 forest sites across the globe, comprising 49,339 site-years, and compared the results with simulated recovery in climate-vegetation models. We found pervasive and substantial “legacy effects” of reduced growth and incomplete recovery for 1 to 4 years after severe drought. Legacy effects were most prevalent in dry ecosystems, among Pinaceae, and among species with low hydraulic safety margins. In contrast, limited or no legacy effects after drought were simulated by current climate-vegetation models. Our results highlight hysteresis in ecosystem-level carbon cycling and delayed recovery from climate extremes.

22. Anderegg, W, A P Ballantyne, W Kolby Smith, J D Majkut, S Rabin, C Beaulieu, R A Birdsey, John P Dunne, R A Houghton, R B Myneni, Yude Pan, Jorge L Sarmiento, N Serota, and Elena Shevliakova, et al., December 2015: Tropical nighttime warming as a dominant driver of variability in the terrestrial carbon sink. Proceedings of the National Academy of Sciences, 112(51), DOI:10.1073/pnas.1521479112 .
    [ Abstract ]

    The terrestrial biosphere is currently a strong carbon (C) sink but may switch to a source in the 21st century as climate-driven losses exceed CO2-driven C gains, thereby accelerating global warming. Although it has long been recognized that tropical climate plays a critical role in regulating interannual climate variability, the causal link between changes in temperature and precipitation and terrestrial processes remains uncertain. Here, we combine atmospheric mass balance, remote sensing-modeled datasets of vegetation C uptake, and climate datasets to characterize the temporal variability of the terrestrial C sink and determine the dominant climate drivers of this variability. We show that the interannual variability of global land C sink has grown by 50–100% over the past 50 y. We further find that interannual land C sink variability is most strongly linked to tropical nighttime warming, likely through respiration. This apparent sensitivity of respiration to nighttime temperatures, which are projected to increase faster than global average temperatures, suggests that C stored in tropical forests may be vulnerable to future warming.

23. Malyshev, Sergey, Elena Shevliakova, Ronald J Stouffer, and S W Pacala, July 2015: Contrasting Local vs. Regional Effects of Land-Use-Change Induced Heterogeneity on Historical Climate: Analysis with the GFDL Earth System Model. Journal of Climate, 28(13), DOI:10.1175/JCLI-D-14-00586.1 .
    [ Abstract ]

    In this study we explore effects of land-use and land-cover change (LULCC) on surface climate using two ensembles of numerical experiments with the Geophysical Fluid Dynamics Laboratory (GFDL) comprehensive Earth System Model ESM2Mb. The experiments simulate historical climate with two different assumptions about LULCC: (1) no land use change with potential vegetation (PV) and (2) with the CMIP5 historical reconstruction of LULCC (LU). We used two different approached in the analysis: (1) we compare differences in LU and PV climates to evaluate the regional and global effects of LULCC, and (2) we characterize sub-grid climate differences among different land-use tiles within each grid cell in the LU experiment. Using the first method, we estimate the magnitude of LULCC effect to be similar to some previous studies. Using the second method we found a pronounced sub-grid signal of LULCC in near-surface temperature over majority of areas affected by LULCC. The signal is strongest on croplands, where it is detectable with 95% confidence over 68.5% of all non-glaciated land grid cells in June-July-August, compared to 8.3% in the first method. In agricultural areas, the sub-grid signal tends to be stronger than LU-PV signal by a factor of 1.3 in tropics in both summer and winter and by 1.5 in extra-tropics in winter. Our analysis for the first time demonstrates and quantifies the local, sub-grid scale LULCC effects with a comprehensive ESM and compares it to previous global and regional approaches.

24. Rabin, S, B I Magi, Elena Shevliakova, and S W Pacala, November 2015: Quantifying regional, time-varying effects of cropland and pasture on vegetation fire. Biogeosciences, 12(22), DOI:10.5194/bg-12-6591-2015 .
    [ Abstract ]

    The global extent of agriculture demands a thorough understanding of the ways it impacts the Earth system through both the modification of the physical and biological characteristics of the landscape as well as through emissions of greenhouse gases and aerosols. People use fire to manage cropland and pasture in many parts of the world, impacting both the timing and amount of fire. So far, much previous research into how these land uses affect fire regimes has either focused on individual small regions or global patterns at annual or decadal scales. Moreover, because pasture is not mapped globally at high resolution, the amount of fire associated with pasture has never been quantified as it has for cropland. The work presented here resolves the effects of agriculture – including pasture – on fire on a monthly basis for regions across the world, using globally gridded data on fire activity and land use at 0.25° resolution. The first global estimate of pasture-associated fire reveals that it accounts for over 40 % of annual burned area. Cropland, generally assumed to reduce fire occurrence, is shown to enhance or suppress fire at different times of year within individual regions. These results bridge important gaps in the understanding of how agriculture and associated management practices influence vegetation fire, enabling the next generation of vegetation and Earth system models more realistically incorporate these anthropogenic effects.

25. Weng, E S., Sergey Malyshev, J W Lichstein, C E Farrior, R Dybzinski, T Zhang, Elena Shevliakova, and S W Pacala, May 2015: Scaling from individuals to ecosystems in an Earth System Model using a mathematically tractable model of height-structured competition for light. Biogeosciences, 12(9), DOI:10.5194/bg-12-2655-2015 .
    [ Abstract ]

    The long-term and large scale dynamics of ecosystems are in large part determined by the performances of individual plants in competition with one another for light, water and nutrients. Woody biomass, a pool of carbon (C) larger than 50% of atmospheric CO2, exists because of height-structured competition for light. However, most of the current Earth System Models that predict climate change and C cycle feedbacks lack both a mechanistic formulation for height-structured competition for light and an explicit scaling from individual plants to the globe. In this study, we incorporate height-structured competition and explicit scaling from individuals to ecosystems into the land model (LM3) currently used in the Earth System Models developed by the Geophysical Fluid Dynamics Laboratory (GFDL). The height-structured formulation is based on the Perfect Plasticity Approximation (PPA), which has been shown to accurately scale from individual-level plant competition for light, water and nutrients to the dynamics of whole communities. Because of the tractability of the PPA, the coupled LM3–PPA model is able to include a large number of phenomena across a range of spatial and temporal scales, and still retain computational tractability, as well as close linkages to mathematically tractable forms of the model. We test a range of predictions against data from temperate broadleaved forests in the northern USA. The results show the model predictions agree with diurnal and annual C fluxes, growth rates of individual trees in the canopy and understory, tree size distributions, and species-level population dynamics during succession. We also show how the competitively optimal allocation strategy – the strategy that can competitively exclude all others – shifts as a function of the atmospheric CO2 concentration. This strategy is referred as an evolutionary stable strategy (ESS) in the ecological literature and is typically not the same as a productivity- or growth-maximizing strategy. Model simulations predict that C sinks caused by CO2 fertilization in forests limited by light and water will be down-regulated if allocation tracks changes in the competitive optimum. The implementation of the model in this paper is for temperate broadleaved forest trees, but the formulation of the model is general. It can be expanded to include other growth forms and physiologies simply by altering parameter values.

26. Hoffman, F, and Elena Shevliakova, et al., February 2014: Causes and Implications of Persistent Atmospheric Carbon Dioxide Biases in Earth System Models. Journal of Geophysical Research, 119(2), DOI:10.1002/2013JG002381 .
    [ Abstract ]

    The strength of feedbacks between a changing climate and future CO2 concentrations are uncertain and difficult to predict using Earth System Models (ESMs). We analyzed emission-driven simulations—in which atmospheric CO2 levels were computed prognostically—for historical (1850–2005) and future periods (RCP 8.5 for 2006–2100) produced by 15 ESMs for the Fifth Phase of the Coupled Model Intercomparison Project (CMIP5). Comparison of ESM prognostic atmospheric CO2 over the historical period with observations indicated that ESMs, on average, had a small positive bias in predictions of contemporary atmospheric CO2. Weak ocean carbon uptake in many ESMs contributed to this bias, based on comparisons with observations of ocean and atmospheric anthropogenic carbon inventories. We found a significant linear relationship between contemporary atmospheric CO2 biases and future CO2 levels for the multi-model ensemble. We used this relationship to create a contemporary CO2 tuned model (CCTM) estimate of the atmospheric CO2 trajectory for the 21st century. The CCTM yielded CO2 estimates of 600 ± 14 ppm at 2060 and 947 ± 35 ppm at 2100, which were 21 ppm and 32 ppm below the multi-model mean during these two time periods. Using this emergent constraint approach, the likely ranges of future atmospheric CO2, CO2-induced radiative forcing, and CO2-induced temperature increases for the RCP 8.5 scenario were considerably narrowed compared to estimates from the full ESM ensemble. Our analysis provided evidence that much of the model-to-model variation in projected CO2 during the 21st century was tied to biases that existed during the observational era, and that model differences in the representation of concentration-carbon feedbacks and other slowly changing carbon cycle processes appear to be the primary driver of this variability. By improving models to more closely match the long-term time series of CO2 from Mauna Loa, our analysis suggests uncertainties in future climate projections can be reduced.

27. Krasting, John P., John P Dunne, Elena Shevliakova, and Ronald J Stouffer, April 2014: Trajectory sensitivity of the transient climate response to cumulative carbon emissions. Geophysical Research Letters, 41(7), DOI:10.1002/2013GL059141 .
    [ Abstract ]

    The robustness of Transient Climate Response to cumulative Emissions (TCRE) is tested using an Earth System Model (Geophysical Fluid Dynamics Laboratory-ESM2G) forced with seven different constant rates of carbon emissions (2 GtC/yr to 25 GtC/yr), including low emission rates that have been largely unexplored in previous studies. We find the range of TCRE resulting from varying emission pathways to be 0.76 to 1.04°C/TtC. This range, however, is small compared to the uncertainty resulting from varying model physics across the Fifth Coupled Model Intercomparison Project ensemble. TCRE has a complex relationship with emission rates; TCRE is largest for both low (2 GtC/yr) and high (25 GtC/yr) emissions and smallest for present-day emissions (5–10 GtC/yr). Unforced climate variability hinders precise estimates of TCRE for periods shorter than 50 years for emission rates near or smaller than present day values. Even if carbon emissions would stop, the prior emissions pathways will affect the future climate responses.

28. Lee, M, Sergey Malyshev, Elena Shevliakova, P C D Milly, and D A Jaffe, October 2014: Capturing interactions between nitrogen and hydrological cycles under historical climate and land use: Susquehanna watershed analysis with the GFDL Land Model LM3-TAN. Biogeosciences, 11(20), DOI:10.5194/bg-11-5809-2014 .
    [ Abstract ]

    We developed a~process model LM3-TAN to assess the combined effects of direct human influences and climate change on Terrestrial and Aquatic Nitrogen (TAN) cycling. The model was developed by expanding NOAA's Geophysical Fluid Dynamics Laboratory land model LM3V-N of coupled terrestrial carbon and nitrogen (C-N) cycling and including new N cycling processes and inputs such as a~soil denitrification, point N sources to streams (i.e. sewage), and stream transport and microbial processes. Because the model integrates ecological, hydrological, and biogeochemical processes, it captures key controls of transport and fate of N in the vegetation-soil-river system in a comprehensive and consistent framework which is responsive to climatic variations and land use changes. We applied the model at 1/8° resolution for a study of the Susquehanna River basin. We simulated with LM3-TAN stream dissolved organic-N, ammonium-N, and nitrate-N loads throughout the river network, and we evaluated the modeled loads for 1986–2005 using data from 15 monitoring stations as well as a reported budget for the entire basin. By accounting for inter-annual hydrologic variability, the model was able to capture inter-annual variations of stream N loadings. While the model was calibrated with the stream N loads only at the last downstream station Marietta (40.02° N, 76.32° W), it captured the N loads well at multiple locations within the basin with different climate regimes, land use types, and associated N sources and transformations in the sub-basins. Furthermore, the calculated and previously reported N budgets agreed well at the level of the whole Susquehanna watershed. Here we illustrate how point and non-point N sources contribute to the various ecosystems are stored, lost, and exported via the river. Local analysis for 6 sub-basins showed combined effects of land use and climate on the soil denitrification rates, with the highest rates in the Lower Susquehanna sub-basin (extensive agriculture; Atlantic coastal climate) and the lowest rates in the West Branch Susquehanna sub-basin (mostly forest; Great Lakes and Midwest climate). In the re-growing secondary forests, most of the N from non-point sources was stored in the vegetation and soil, but in the agricultural lands most N inputs were removed by soil denitrification indicating that anthropogenic N applications could drive substantial increase of N2O emission, an intermediate of the denitrification process.

29. Lichstein, J W., Ni-Zhang Golaz, Sergey Malyshev, Elena Shevliakova, Tao Zhang, J Sheffield, R A Birdsey, Jorge L Sarmiento, and S W Pacala, June 2014: Confronting terrestrial biosphere models with forest inventory data. Ecological Applications, 24(4), DOI:10.1890/13-0600.1 .
    [ Abstract ]

    Efforts to test and improve terrestrial biosphere models (TBMs) using a variety of data sources have become increasingly common. However, geographically extensive forest inventories have been under-exploited in previous model-data fusion efforts. Inventory observations of forest growth, mortality, and biomass integrate processes across a range of time scales, including slow time-scale processes such as species turnover, that are likely to have important effects on ecosystem responses to environmental variation. However, the large number (thousands) of inventory plots precludes detailed measurements at each location, so that uncertainty in climate, soil properties, and other environmental drivers may be large. Errors in driver variables, if ignored, introduce bias into model-data fusion. We estimated errors in climate and soil drivers at U.S. Forest Inventory and Analysis (FIA) plots, and we explored the effects of these errors on model-data fusion with the Geophysical Fluid Dynamics Laboratory LM3V dynamic global vegetation model. When driver errors were ignored or assumed small at FIA plots, responses of biomass production in LM3V to precipitation and soil available water capacity appeared steeper than the corresponding responses estimated from FIA data. These differences became non-significant if driver errors at FIA plots were assumed large. Ignoring driver errors when optimizing LM3V parameter values yielded estimates for fine-root allocation that were larger than biometric estimates, which is consistent with the expected direction of bias. To explore if complications posed by driver errors could be circumvented by relying on intensive study sites where driver errors are small, we performed a power analysis. To accurately quantify the response of biomass production to spatial variation in mean annual precipitation within the eastern U.S. would require at least 40 intensive study sites, which is larger than the number of sites typically available for individual biomes in existing plot networks. Driver errors may be accommodated by several existing model-data fusion approaches, including hierarchical Bayesian methods and ensemble filtering methods; however, these methods are computationally expensive. We propose a new approach, in which the TBM functional response is fit directly to the driver-error-corrected functional response estimated from data, rather than to the raw observations.

30. Milly, P C., Sergey Malyshev, Elena Shevliakova, Krista A Dunne, Kirsten L Findell, T Gleeson, Zhi Liang, Peter Phillipps, Ronald J Stouffer, and S C Swenson, October 2014: An enhanced model of land water and energy for global hydrologic and earth-system studies. Journal of Hydrometeorology, 15(5), DOI:10.1175/JHM-D-13-0162.1 .
    [ Abstract ]

    “LM3” is a new model of terrestrial water, energy, and carbon, intended for use in global hydrologic analyses and as a component of earth-system and physical-climate models. It is designed to improve upon the performance and extend the scope of the predecessor Land Dynamics (LaD) and LM3V models, by quantifying better the physical controls of climate and biogeochemistry and by relating more directly to components of the global water system that touch human concerns. LM3 includes multi-layer representations of temperature, liquid-water content, and ice content of both snow pack and macroporous soil/bedrock; topography-based description of saturated area and groundwater discharge; and transport of runoff to the ocean via a global river and lake network. Sensible heat transport by water mass is accounted throughout for a complete energy balance. Carbon and vegetation dynamics and biophysics are represented as in the model LM3V. In numerical experiments, LM3 avoids some of the limitations of the LaD model and provides qualitatively (though not always quantitatively) reasonable estimates, from a global perspective, of observed spatial and/or temporal variations of vegetation density, albedo, streamflow, water-table depth, permafrost, and lake levels. Amplitude and phase of annual cycle of total water storage are simulated well. Realism of modeled lake levels varies widely. The water table tends to be consistently too shallow in humid regions. Biophysical properties have an artificial step-wise spatial structure, and equilibrium vegetation is sensitive to initial conditions. Explicit resolution of thick (>100 m) unsaturated zones and permafrost is possible, but only at the cost of long (>>300 y) model spin-up times.

31. Subin, Zachary M., P C D Milly, Benjamin N Sulman, Sergey Malyshev, and Elena Shevliakova, in press: Resolving terrestrial ecosystem processes along a subgrid topographic gradient for an earth-system model. Hydrology and Earth System Sciences Discussions. DOI:10.5194/hessd-11-8443-2014. July 2014.
    [ Abstract ]

    Soil moisture is a crucial control on surface water and energy fluxes, vegetation, and soil carbon cycling. Earth-system models (ESMs) generally represent an areal-average soil-moisture state in gridcells at scales of 50–200 km and as a result are not able to capture the nonlinear effects of topographically-controlled subgrid heterogeneity in soil moisture, in particular where wetlands are present. We addressed this deficiency by building a subgrid representation of hillslope-scale topographic gradients, TiHy (Tiled-hillslope Hydrology), into the Geophysical Fluid Dynamics Laboratory (GFDL) land model (LM3). LM3-TiHy models one or more representative hillslope geometries for each gridcell by discretizing them into land model tiles hydrologically coupled along an upland-to-lowland gradient. Each tile has its own surface fluxes, vegetation, and vertically-resolved state variables for soil physics and biogeochemistry. LM3-TiHy simulates a gradient in soil moisture and water-table depth between uplands and lowlands in each gridcell. Three hillslope hydrological regimes appear in non-permafrost regions in the model: wet and poorly-drained, wet and well-drained, and dry; with large, small, and zero wetland area predicted, respectively. Compared to the untiled LM3 in stand-alone experiments, LM3-TiHy simulates similar surface energy and water fluxes in the gridcell-mean. However, in marginally wet regions around the globe, LM3-TiHy simulates shallow groundwater in lowlands, leading to higher evapotranspiration, lower surface temperature, and higher leaf area compared to uplands in the same gridcells. Moreover, more than four-fold larger soil carbon concentrations are simulated globally in lowlands as compared with uplands. We compared water-table depths to those simulated by a recent global model-observational synthesis, and we compared wetland and inundated areas diagnosed from the model to observational datasets. The comparisons demonstrate that LM3-TiHy has the capability to represent some of the controls of these hydrological variables, but also that improvement in parameterization and input datasets are needed for more realistic simulations. We found large sensitivity in model-diagnosed wetland and inundated area to the depth of conductive soil and the parameterization of macroporosity. With improved parameterization and inclusion of peatland biogeochemical processes, the model could provide a new approach to investigating the vulnerability of Boreal peatland carbon to climate change in ESMs.

32. Sulman, Benjamin N., R P Phillips, A Christopher Oishi, Elena Shevliakova, and S W Pacala, December 2014: Microbe-driven turnover offsets mineral-mediated storage of soil carbon under elevated CO₂. Nature Climate Change, 4(12), DOI:10.1038/nclimate2436 .
    [ Abstract ]

    The sensitivity of soil organic carbon (SOC) to changing environmental conditions represents a critical uncertainty in coupled carbon cycle–climate models1. Much of this uncertainty arises from our limited understanding of the extent to which root–microbe interactions induce SOC losses (through accelerated decomposition or ‘priming’2) or indirectly promote SOC gains (via ‘protection’ through interactions with mineral particles3, 4). We developed a new SOC model to examine priming and protection responses to rising atmospheric CO2. The model captured disparate SOC responses at two temperate free-air CO2 enrichment (FACE) experiments. We show that stabilization of ‘new’ carbon in protected SOC pools may equal or exceed microbial priming of ‘old’ SOC in ecosystems with readily decomposable litter and high clay content (for example, Oak Ridge5). In contrast, carbon losses induced through priming dominate the net SOC response in ecosystems with more resistant litters and lower clay content (for example, Duke6). The SOC model was fully integrated into a global terrestrial carbon cycle model to run global simulations of elevated CO2 effects. Although protected carbon provides an important constraint on priming effects, priming nonetheless reduced SOC storage in the majority of terrestrial areas, partially counterbalancing SOC gains from enhanced ecosystem productivity.

33. Todd-Brown, K E., and Elena Shevliakova, et al., April 2014: Changes in soil organic carbon storage predicted by Earth system models during the 21st century. Biogeosciences, 11(8), DOI:10.5194/bg-11-2341-2014 .
    [ Abstract ]

    Soil is currently thought to be a sink for carbon; however, the response of this sink to increasing levels of atmospheric carbon dioxide and climate change is uncertain. In this study, we analyzed soil organic carbon (SOC) changes from 11 Earth system models (ESMs) under the historical and high radiative forcing (RCP 8.5) scenarios between 1850 and 2100. We used a reduced complexity model based on temperature and moisture sensitivities to analyze the drivers of SOC losses. ESM estimates of SOC change over the 21st century (2090–2099 minus 1997–2006) ranged from a loss of 72 Pg C to a gain 253 Pg C with a multi-model mean gain of 63 Pg C. All ESMs showed cumulative increases in both NPP (15% to 59%) and decreases in SOC turnover times (15% to 28%) over the 21st century. Most of the model-to-model variation in SOC change was explained by initial SOC stocks combined with the relative changes in soil inputs and decomposition rates (R2 = 0.88, p<0.01). Between models, increases in decomposition rate were well explained by a combination of initial decomposition rate, ESM-specific Q10-factors, and changes in soil temperature (R2 = 0.80, p<0.01). All SOC changes depended on sustained increases in NPP with global change (primarily driven by increasing CO2) and conversion of additional plant inputs into SOC. Most ESMs omit potential constraints on SOC storage, such as priming effects, nutrient availability, mineral surface stabilization and aggregate formation. Future models that represent these constraints are likely to estimate smaller increases in SOC storage during the 21st century.

34. Dunne, John P., Jasmin G John, Elena Shevliakova, Ronald J Stouffer, John P Krasting, Sergey Malyshev, P C D Milly, Lori T Sentman, Alistair Adcroft, William F Cooke, Krista A Dunne, Stephen M Griffies, Robert Hallberg, Matthew J Harrison, Hiram Levy II, Andrew T Wittenberg, Peter Phillipps, and Niki Zadeh, April 2013: GFDL's ESM2 global coupled climate-carbon Earth System Models Part II: Carbon system formulation and baseline simulation characteristics. Journal of Climate, 26(7), DOI:10.1175/JCLI-D-12-00150.1 .
    [ Abstract ]

    We describe carbon system formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate good climate fidelity as described in Part I while incorporating explicit and consistent carbon dynamics. The two models differ almost exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. On land, both ESMs include a revised land model to simulate competitive vegetation distributions and functioning, including carbon cycling among vegetation, soil and atmosphere. In the ocean, both models include new biogeochemical algorithms including phytoplankton functional group dynamics with flexible stoichiometry. Preindustrial simulations are spun up to give stable, realistic carbon cycle means and variability. Significant differences in simulation characteristics of these two models are described. Due to differences in oceanic ventilation rates (Part I) ESM2M has a stronger biological carbon pump but weaker northward implied atmospheric CO2 transport than ESM2G. The major advantages of ESM2G over ESM2M are: improved representation of surface chlorophyll in the Atlantic and Indian Oceans and thermocline nutrients and oxygen in the North Pacific. Improved tree mortality parameters in ESM2G produced more realistic carbon accumulation in vegetation pools. The major advantages of ESM2M over ESM2G are reduced nutrient and oxygen biases in the Southern and Tropical Oceans.

35. Gerber, S, L O Hedin, S G Keel, S W Pacala, and Elena Shevliakova, October 2013: Land-use change and nitrogen feedbacks constrain the trajectory of the land carbon sink. Geophysical Research Letters, 40(19), DOI:10.1002/grl.50957 .
    [ Abstract ]

    Our understanding of Earth's carbon-climate system depends critically upon interactions between rising atmospheric CO₂, changing land use, and nitrogen limitation on vegetation growth. Using a global land model we show how these factors interact locally to generate the global land carbon sink over the past 200 years. Nitrogen constraints were alleviated by N₂ fixation in the tropics and by atmospheric nitrogen deposition in extra-tropical regions. Non-linear interactions between land-use change and land carbon and nitrogen cycling originated from three major mechanisms: (i) a sink foregone that would have occurred without land-use conversion; (ii) an accelerated response of secondary vegetation to CO₂ and nitrogen, and (iii) a compounded clearance loss from deforestation. Over time, these non-linear effects have become increasingly important and reduce the present-day net carbon sink by ~40% or 0.4 PgC yr^-1.

36. Jeong, S-J, David Medvigy, Elena Shevliakova, and Sergey Malyshev, January 2013: Predicting changes in temperate forest budburst using continental-scale observations and models. Geophysical Research Letters, 40(2), DOI:10.1029/2012GL054431 .
    [ Abstract ]

    A new framework for understanding the macro-scale variations in spring phenology is developed by using new data from the USA National Phenology Network. Changes in spring budburst for the U.S. are predicted by using Coupled Model Intercomparison Project phase 5 outputs. Macro-scale budburst simulations for the coming century indicate that projected warming leads to earlier budburst by up to 17 days. The latitudinal gradient of budburst becomes less pronounced due to spatially-varying sensitivity of budburst to climate change, even in the most conservative emissions scenarios. Currently existing inter-species differences in budburst date are predicted to become smaller, indicating the potential for secondary impacts at the ecosystem level. We expect that these climate-driven changes in phenology will have large effects on the carbon budget of U.S. forests and these controls should be included in dynamic global vegetation models.

37. Jones, C, and Elena Shevliakova, et al., July 2013: Twenty-First-Century Compatible CO2 Emissions and Airborne Fraction Simulated by CMIP5 Earth System Models under Four Representative Concentration Pathways. Journal of Climate, 26(13), DOI:10.1175/JCLI-D-12-00554.1 .
    [ Abstract ]

    The carbon cycle is a crucial earth system component affecting climate and atmospheric composition. The response of natural carbon uptake to CO2 and climate change will determine anthropogenic emissions compatible with a target CO2 pathway. For CMIP5 4 future Representative Concentration Pathways have been generated by Integrated Assessment Models and used as scenarios by state-of-the-art climate models, enabling quantification of compatible carbon emissions for the 4 scenarios by complex, process-based models. Here we present results from 15 such Earth System GCMs for future changes in land and ocean carbon storage and the implications for anthropogenic emissions. The results are consistent with the underlying scenarios, but show substantial model spread. Uncertainty in land carbon uptake due to differences among models is comparable with the spread across scenarios. Model estimates of historical fossil fuel emissions agree well with reconstructions and future projections for RCP2.6 and RCP4.5 are consistent with the IAMs. For high-end scenarios (6.0 and 8.5) GCMs simulate smaller compatible emissions than the IAMs, indicating a larger climate-carbon cycle feedback in the GCMs in these scenarios. For the RCP2.6 mitigation scenario an average reduction of 50% in emissions by 2050 from 1990 levels is required but with very large model spread (14-96%). The models also disagree on both the requirement for sustained negative emissions to achieve the RCP2.6 CO2 concentration and the success of this scenario to restrict global warming below 2°C. All models agree that the future airborne-fraction depends strongly on the emissions profile with higher airborne-fraction for higher emissions scenarios.

38. Shevliakova, Elena, Ronald J Stouffer, Sergey Malyshev, John P Krasting, G C Hurtt, and S W Pacala, October 2013: Historical warming reduced due to enhanced land carbon uptake. Proceedings of the National Academy of Sciences, 110(42), DOI:10.1073/pnas.1314047110 .
    [ Abstract ]

    Previous studies have demonstrated the importance of enhanced vegetation growth under future elevated atmospheric CO2 for 21st century climate warming. Surprisingly no study has completed an analogous assessment for the historical period, during which emissions of greenhouse gases increased rapidly and landuse changes (LUC) dramatically altered terrestrial carbon sources and sinks. Using the Geophysical Fluid Dynamics Laboratory comprehensive Earth System Model ESM2G and a reconstruction of the LUC, we estimate that enhanced vegetation growth has lowered the historical atmospheric CO2 concentration by 85 ppm, avoiding an additional 0.31 ± 0.06 °C warming. We demonstrate that without enhanced vegetation growth the total residual terrestrial carbon flux (i.e., the net land flux minus LUC flux) would be a source of 65–82 Gt of carbon (GtC) to atmosphere instead of the historical residual carbon sink of 186–192 GtC, a carbon saving of 251–274 GtC.

39. Yin, Lei, R Fu, Elena Shevliakova, and R Dickinson, December 2013: How well can CMIP5 simulate precipitation and its controlling processes over tropical South America? Climate Dynamics, 41(11-12), DOI:10.1007/s00382-012-1582-y .
    [ Abstract ]

    Underestimated rainfall over Amazonia was a common problem for the Coupled Model Intercomparison Project phase 3 (CMIP3) models. We investigate whether it still exists in the CMIP phase 5 (CMIP5) models and, if so, what causes these biases? Our evaluation of historical simulations shows that some models still underestimate rainfall over Amazonia. During the dry season, both convective and large-scale precipitation is underestimated in most models. GFDL-ESM2M and IPSL notably show more pentads with no rainfall. During the wet season, large-scale precipitation is still underestimated in most models. In the dry and transition seasons, models with more realistic moisture convergence and surface evapotranspiration generally have more realistic rainfall totals. In some models, overestimates of rainfall are associated with the adjacent tropical and eastern Pacific ITCZs. However, in other models, too much surface net radiation and a resultant high Bowen ratio appears to cause underestimates of rainfall. During the transition season, low pre-seasonal latent heat, high sensible flux, and a weaker influence of cold air incursions contribute to the dry bias. About half the models can capture, but overestimate, the influences of teleconnection. Based on a simple metric, HadGEM2-ES outperforms other models especially for surface conditions and atmospheric circulation. GFDL-ESM2M has the strongest dry bias presumably due to its overestimate of moisture divergence, induced by overestimated ITCZs in adjacent oceans, and reinforced by positive feedbacks between reduced cloudiness, high Bowen ratio and suppression of rainfall during the dry season, and too weak incursions of extratropical disturbances during the transition season.

40. Dunne, John P., Jasmin G John, Alistair Adcroft, Stephen M Griffies, Robert Hallberg, Elena Shevliakova, Ronald J Stouffer, William F Cooke, Krista A Dunne, Matthew J Harrison, John P Krasting, Sergey Malyshev, P C D Milly, Peter Phillipps, Lori T Sentman, Bonita L Samuels, Michael J Spelman, Michael Winton, Andrew T Wittenberg, and Niki Zadeh, October 2012: GFDL's ESM2 global coupled climate-carbon Earth System Models Part I: Physical formulation and baseline simulation characteristics. Journal of Climate, 25(19), DOI:10.1175/JCLI-D-11-00560.1 .
    [ Abstract ]

    We describe the physical climate formulation and simulation characteristics of two new global coupled carbon-climate Earth System Models, ESM2M and ESM2G. These models demonstrate similar climate fidelity as the Geophysical Fluid Dynamics Laboratory’s previous CM2.1 climate model while incorporating explicit and consistent carbon dynamics. The two models differ exclusively in the physical ocean component; ESM2M uses Modular Ocean Model version 4.1 with vertical pressure layers while ESM2G uses Generalized Ocean Layer Dynamics with a bulk mixed layer and interior isopycnal layers. Differences in the ocean mean state include the thermocline depth being relatively deep in ESM2M and relatively shallow in ESM2G compared to observations. The crucial role of ocean dynamics on climate variability is highlighted in the El Niño-Southern Oscillation being overly strong in ESM2M and overly weak ESM2G relative to observations. Thus, while ESM2G might better represent climate changes relating to: total heat content variability given its lack of long term drift, gyre circulation and ventilation in the North Pacific, tropical Atlantic and Indian Oceans, and depth structure in the overturning and abyssal flows, ESM2M might better represent climate changes relating to: surface circulation given its superior surface temperature, salinity and height patterns, tropical Pacific circulation and variability, and Southern Ocean dynamics. Our overall assessment is that neither model is fundamentally superior to the other, and that both models achieve sufficient fidelity to allow meaningful climate and earth system modeling applications. This affords us the ability to assess the role of ocean configuration on earth system interactions in the context of two state-of-the-art coupled carbon-climate models.

41. Jeong, S-J, David Medvigy, Elena Shevliakova, and Sergey Malyshev, March 2012: Uncertainties in terrestrial carbon budgets related to spring phenology. Journal of Geophysical Research, 117, G01030, DOI:10.1029/2011JG001868 .
    [ Abstract ]

    In temperate regions, the budburst date of deciduous trees is mainly regulated by temperature variation, but the exact nature of the temperature dependence has been a matter of debate. One hypothesis is that budburst date depends purely on the accumulation of warm temperature; a competing hypothesis states that exposure to cold temperatures is also important for budburst. In this study, variability in budburst is evaluated using 15 years of budburst data for 17 tree species at Harvard Forest. We compare two budburst hypotheses through reversible jump Markov Chain Monte Carlo. Then, we investigate how uncertainties in budburst date mapped into uncertainties in ecosystem carbon using GFDL's LM3 land model. For 15 of 17 species, we find that more complicated budburst models that account for a chilling period are favored over simpler models that do not include such dependence. LM3 simulations show that the choice of budburst model induces differences in the timing of carbon uptake commencement of ~11 days, in the magnitude of April-May carbon uptake of ~1.03 g C m^-2 day^-1, and in total ecosystem carbon stocks of ~2 kg C m^-2. While the choice of whether to include a chilling period in the budburst model strongly contributes to this variability, another important factor is how the species-dependent field data gets mapped into LM3's single deciduous plant functional type (PFT). We conclude budburst timing has a strong impact on simulated CO₂ fluxes, and uncertainty in the fluxes can be substantially reduced by improving the model's representation of PFT diversity.

42. Magi, B I., S Rabin, Elena Shevliakova, and S W Pacala, August 2012: Separating agricultural and non-agricultural fire seasonality at regional scales. Biogeosciences, 9(8), DOI:10.5194/bg-9-3003-2012 .
    [ Abstract ]

    The timing and length of burning seasons in different parts of the world depend on climate, land cover characteristics, and human activities. In this study, global fire data from satellite-based instruments are used in conjunction with global gridded distributions of agricultural land cover types (defined as the sum of cropland and pasture area) to separate the seasonality of agricultural burning practices from that of non-agricultural fire. The results presented in this study show that agricultural and non-agricultural land experience broadly different fire seasonality patterns that are not always linked to climate conditions. We highlight these differences on a regional basis, examining variations in both agricultural land cover and associated cultural practices to help explain our results. While we discuss two land cover categories, the methods can be generalized to derive seasonality for any number of land uses or cover types. This will be useful as global fire models evolve to be fully interactive with land use and land cover change in the next generation of Earth system models.

43. Donner, Leo J., Bruce Wyman, Richard S Hemler, Larry W Horowitz, Yi Ming, Ming Zhao, J-C Golaz, Paul Ginoux, Shian-Jiann Lin, M Daniel Schwarzkopf, John Austin, G Alaka, William F Cooke, Thomas L Delworth, Stuart Freidenreich, C Tony Gordon, Stephen M Griffies, Isaac M Held, William J Hurlin, Stephen A Klein, Thomas R Knutson, Amy R Langenhorst, Hyun-Chul Lee, Yanluan Lin, B I Magi, Sergey Malyshev, P C D Milly, Vaishali Naik, Mary Jo Nath, R Pincus, Jeff J Ploshay, V Ramaswamy, Charles J Seman, Elena Shevliakova, Joseph J Sirutis, William F Stern, Ronald J Stouffer, R John Wilson, Michael Winton, Andrew T Wittenberg, and Fanrong Zeng, July 2011: The dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component AM3 of the GFDL Global Coupled Model CM3. Journal of Climate, 24(13), DOI:10.1175/2011JCLI3955.1 .
    [ Abstract ]

    The Geophysical Fluid Dynamics Laboratory (GFDL) has developed a coupled general circulation model (CM3) for atmosphere, oceans, land, and sea ice. The goal of CM3 is to address emerging issues in climate change, including aerosol-cloud interactions, chemistry-climate interactions, and coupling between the troposphere and stratosphere. The model is also designed to serve as the physical-system component of earth-system models and models for decadal prediction in the near-term future, for example, through improved simulations in tropical land precipitation relative to earlier-generation GFDL models. This paper describes the dynamical core, physical parameterizations, and basic simulation characteristics of the atmospheric component (AM3) of this model. Relative to GFDL AM2, AM3 includes new treatments of deep and shallow cumulus convection, cloud-droplet activation by aerosols, sub-grid variability of stratiform vertical velocities for droplet activation, and atmospheric chemistry driven by emissions with advective, convective, and turbulent transport. AM3 employs a cubed-sphere implementation of a finite-volume dynamical core and is coupled to LM3, a new land model with eco-system dynamics and hydrology. Most basic circulation features in AM3 are simulated as realistically, or more so, than in AM2. In particular, dry biases have been reduced over South America. In coupled mode, the simulation of Arctic sea ice concentration has improved. AM3 aerosol optical depths, scattering properties, and surface clear-sky downward shortwave radiation are more realistic than in AM2. The simulation of marine stratocumulus decks and the intensity distributions of precipitation remain problematic, as in AM2. The last two decades of the 20th century warm in CM3 by .32°C relative to 1881-1920. The Climate Research Unit (CRU) and Goddard Institute for Space Studies analyses of observations show warming of .56°C and .52°C, respectively, over this period. CM3 includes anthropogenic cooling by aerosol cloud interactions, and its warming by late 20th century is somewhat less realistic than in CM2.1, which warmed .66°C but did not include aerosol cloud interactions. The improved simulation of the direct aerosol effect (apparent in surface clear-sky downward radiation) in CM3 evidently acts in concert with its simulation of cloud-aerosol interactions to limit greenhouse gas warming in a way that is consistent with observed global temperature changes.

44. Hurtt, G C., and Elena Shevliakova, et al., November 2011: Harmonization of land-use scenarios for the period 1500–2100: 600 years of global gridded annual land-use transitions, wood harvest, and resulting secondary lands. Climatic Change, 109(1-2), DOI:10.1007/s10584-011-0153-2 .
    [ Abstract ]

    In preparation for the fifth Assessment Report (AR5) of the Intergovernmental Panel on Climate Change (IPCC), the international community is developing new advanced Earth System Models (ESMs) to assess the combined effects of human activities (e.g. land use and fossil fuel emissions) on the carbon-climate system. In addition, four Representative Concentration Pathway (RCP) scenarios of the future (2005–2100) are being provided by four Integrated Assessment Model (IAM) teams to be used as input to the ESMs for future carbon-climate projections (Moss et al. 2010). The diversity of approaches and requirements among IAMs and ESMs for tracking land-use change, along with the dependence of model projections on land-use history, presents a challenge for effectively passing data between these communities and for smoothly transitioning from the historical estimates to future projections. Here, a harmonized set of land-use scenarios are presented that smoothly connects historical reconstructions of land use with future projections, in the format required by ESMs. The land-use harmonization strategy estimates fractional land-use patterns and underlying land-use transitions annually for the time period 1500–2100 at 0.5° × 0.5° resolution. Inputs include new gridded historical maps of crop and pasture data from HYDE 3.1 for 1500–2005, updated estimates of historical national wood harvest and of shifting cultivation, and future information on crop, pasture, and wood harvest from the IAM implementations of the RCPs for the period 2005–2100. The computational method integrates these multiple data sources, while minimizing differences at the transition between the historical reconstruction ending conditions and IAM initial conditions, and working to preserve the future changes depicted by the IAMs at the grid cell level. This study for the first time harmonizes land-use history data together with future scenario information from multiple IAMs into a single consistent, spatially gridded, set of land-use change scenarios for studies of human impacts on the past, present, and future Earth system.

45. Sentman, Lori T., Elena Shevliakova, Ronald J Stouffer, and Sergey Malyshev, October 2011: Time scales of terrestrial carbon response related to land-use application: Implications for initializing an earth system model. Earth Interactions, 15(30), DOI:10.1175/2011EI401.1 .
    [ Abstract ]

    The dynamic vegetation and carbon cycling component, LM3V, of the Geophysical Fluid Dynamics Laboratory (GFDL) prototype Earth System Model (ESM2.1), has been designed to simulate the effects of land use on terrestrial carbon pools, including secondary vegetation regrowth. Because of the long time scales associated with the carbon adjustment, special consideration is required when initializing the Earth System Model (ESM) when “historical” simulations are conducted. Starting from an equilibrated, preindustrial climate and potential vegetation state in an “offline” land only model (LM3V), estimates of historical land use are instantaneously applied in five experiments beginning in calendar years: 1500, 1600, 1700, 1750 and 1800. This application results in the land carbon pools experiencing an abrupt change – a “carbon shock”- and the secondary vegetation needs time to regrow into consistency with the harvesting history. We find that it takes approximately 100 years for the vegetation to recover from the carbon shock, while soils take at least 150 years to recover. The vegetation carbon response is driven primarily by land-use history, while the soil carbon response is affected by both land-use history and the geographic pattern of soil respiration rates. Based on these results, we recommend the application of historical land-use scenarios in 1700 to provide sufficient time for the land carbon in ESMs with secondary vegetation to equilibrate to adequately simulate carbon stores at the start of the historical integrations (i.e., 1860) in a computationally efficient manner.

46. Gerber, S, L O Hedin, M Oppenheimer, S W Pacala, and Elena Shevliakova, January 2010: Nitrogen cycling and feedbacks in a global dynamic land model. Global Biogeochemical Cycles, 24, GB1001, DOI:10.1029/2008GB003336 .
    [ Abstract ]

    Global anthropogenic changes in carbon (C) and nitrogen (N) cycles call for modeling tools that are able to address and quantify essential interactions between N, C, and climate in terrestrial ecosystems. Here we introduce a prognostic N cycle within the Princeton–Geophysical Fluid Dynamic Laboratory (GFDL) LM3V land model. The model captures mechanisms essential for N cycling and their feedbacks on C cycling: N limitation of plant productivity, the N dependence of C decomposition and stabilization in soils, removal of available N by competing sinks, ecosystem losses that include dissolved organic and volatile N, and ecosystem inputs through biological N fixation. Our model captures many essential characteristics of C-N interactions and is capable of broadly recreating spatial and temporal variations in N and C dynamics. The introduced N dynamics improve the model's short-term NPP response to step changes in CO2. Consistent with theories of successional dynamics, we find that physical disturbance induces strong C-N feedbacks, caused by intermittent N loss and subsequent N limitation. In contrast, C-N interactions are weak when the coupled model system approaches equilibrium. Thus, at steady state, many simulated features of the carbon cycle, such as primary productivity and carbon inventories, are similar to simulations that do not include C-N feedbacks.

47. Lichstein, J W., Ni-Zhang Golaz, Sergey Malyshev, Elena Shevliakova, Tao Zhang, J Sheffield, R A Birdsey, Jorge L Sarmiento, and S W Pacala, April 2010: Confronting terrestrial biosphere models with forest inventory data. Ecological Applications, 20(3), DOI:10.1890/13-0600.1 .
    [ Abstract ]

    Efforts to test and improve terrestrial biosphere models (TBMs) using a variety of data sources have become increasingly common. However, geographically extensive forest inventories have been under-exploited in previous model-data fusion efforts. Inventory observations of forest growth, mortality, and biomass integrate processes across a range of time scales, including slow time-scale processes such as species turnover, that are likely to have important effects on ecosystem responses to environmental variation. However, the large number (thousands) of inventory plots precludes detailed measurements at each location, so that uncertainty in climate, soil properties, and other environmental drivers may be large. Errors in driver variables, if ignored, introduce bias into model-data fusion. We estimated errors in climate and soil drivers at U.S. Forest Inventory and Analysis (FIA) plots, and we explored the effects of these errors on model-data fusion with the Geophysical Fluid Dynamics Laboratory LM3V dynamic global vegetation model. When driver errors were ignored or assumed small at FIA plots, responses of biomass production in LM3V to precipitation and soil available water capacity appeared steeper than the corresponding responses estimated from FIA data. These differences became non-significant if driver errors at FIA plots were assumed large. Ignoring driver errors when optimizing LM3V parameter values yielded estimates for fine-root allocation that were larger than biometric estimates, which is consistent with the expected direction of bias. To explore if complications posed by driver errors could be circumvented by relying on intensive study sites where driver errors are small, we performed a power analysis. To accurately quantify the response of biomass production to spatial variation in mean annual precipitation within the eastern U.S. would require at least 40 intensive study sites, which is larger than the number of sites typically available for individual biomes in existing plot networks. Driver errors may be accommodated by several existing model-data fusion approaches, including hierarchical Bayesian methods and ensemble filtering methods; however, these methods are computationally expensive. We propose a new approach, in which the TBM functional response is fit directly to the driver-error-corrected functional response estimated from data, rather than to the raw observations.

48. Hurtt, G C., and Elena Shevliakova, et al., June 2009: Harmonisation of global land-use scenarios for the period 1500-2100 for IPCC-AR5. iLEAPS Newsletter, 7, 6-8.
49. Shevliakova, Elena, S W Pacala, Sergey Malyshev, P C D Milly, and Lori T Sentman, et al., June 2009: Carbon cycling under 300 years of land use change: Importance of the secondary vegetation sink. Global Biogeochemical Cycles, 23, GB2022, DOI:10.1029/2007GB003176 .
    [ Abstract ]

    We have developed a dynamic land model (LM3V) able to simulate ecosystem dynamics and exchanges of water, energy, and CO2 between land and atmosphere. LM3V is specifically designed to address the consequences of land use and land management changes including cropland and pasture dynamics, shifting cultivation, logging, fire, and resulting patterns of secondary regrowth. Here we analyze the behavior of LM3V, forced with the output from the Geophysical Fluid Dynamics Laboratory (GFDL) atmospheric model AM2, observed precipitation data, and four historic scenarios of land use change for 1700–2000. Our analysis suggests a net terrestrial carbon source due to land use activities from 1.1 to 1.3 GtC/a during the 1990s, where the range is due to the difference in the historic cropland distribution. This magnitude is substantially smaller than previous estimates from other models, largely due to our estimates of a secondary vegetation sink of 0.35 to 0.6 GtC/a in the 1990s and decelerating agricultural land clearing since the 1960s. For the 1990s, our estimates for the pastures' carbon flux vary from a source of 0.37 to a sink of 0.15 GtC/a, and for the croplands our model shows a carbon source of 0.6 to 0.9 GtC/a. Our process-based model suggests a smaller net deforestation source than earlier bookkeeping models because it accounts for decelerated net conversion of primary forest to agriculture and for stronger secondary vegetation regrowth in tropical regions. The overall uncertainty is likely to be higher than the range reported here because of uncertainty in the biomass recovery under changing ambient conditions, including atmospheric CO2 concentration, nutrients availability, and climate.

50. Crevoisier, C, Elena Shevliakova, M Gloor, C Wirth, and S W Pacala, 2007: Drivers of fire in the boreal forests: Data constrained design of a prognostic model of burned area for use in dynamic global vegetation models. Journal of Geophysical Research, 112, D24112, DOI:10.1029/2006JD008372 .
    [ Abstract ]

    Boreal regions are an important component of the global carbon cycle because they host large stocks of aboveground and belowground carbon. Since boreal forest evolution is closely related to fire regimes, shifts in climate are likely to induce changes in ecosystems, potentially leading to a large release of carbon and other trace gases to the atmosphere. Prediction of the effect of this potential climate feedback on the Earth system is therefore important and requires the modeling of fire as a climate driven process in dynamic global vegetation models (DGVMs). Here, we develop a new data-based prognostic model, for use in DGVMs, to estimate monthly burned area from four climate (precipitation, temperature, soil water content and relative humidity) and one human-related (road density) predictors for boreal forest. The burned area model is a function of current climatic conditions and is thus responsive to climate change. Model parameters are estimated using a Markov Chain Monte Carlo method applied to on ground observations from the Canadian Large Fire Database. The model is validated against independent observations from three boreal regions: Canada, Alaska and Siberia. Provided realistic climate predictors, the model is able to reproduce the seasonality, intensity and interannual variability of burned area, as well as the location of fire events. In particular, the model simulates well the timing of burning events, with two thirds of the events predicted for the correct month and almost all the rest being predicted 1 month before or after the observed event. The predicted annual burned area is in the range of various current estimates. The estimated annual relative error (standard deviation) is twelve percent in a grid cell, which makes the model suitable to study quantitatively the evolution of burned area with climate.

51. Findell, Kirsten L., Elena Shevliakova, P C D Milly, and Ronald J Stouffer, July 2007: Modeled impact of anthropogenic land cover change on climate. Journal of Climate, 20(14), DOI:10.1175/JCLI4185.1 .
    [ Abstract ]

    Equilibrium experiments with the Geophysical Fluid Dynamics Laboratory’s climate model are used to investigate the impact of anthropogenic land cover change on climate. Regions of altered land cover include large portions of Europe, India, eastern China, and the eastern United States. Smaller areas of change are present in various tropical regions. This study focuses on the impacts of biophysical changes associated with the land cover change (albedo, root and stomatal properties, roughness length), which is almost exclusively a conversion from forest to grassland in the model; the effects of irrigation or other water management practices and the effects of atmospheric carbon dioxide changes associated with land cover conversion are not included in these experiments. The model suggests that observed land cover changes have little or no impact on globally averaged climatic variables (e.g., 2-m air temperature is 0.008 K warmer in a simulation with 1990 land cover compared to a simulation with potential natural vegetation cover). Differences in the annual mean climatic fields analyzed did not exhibit global field significance. Within some of the regions of land cover change, however, there are relatively large changes of many surface climatic variables. These changes are highly significant locally in the annual mean and in most months of the year in eastern Europe and northern India. They can be explained mainly as direct and indirect consequences of model-prescribed increases in surface albedo, decreases in rooting depth, and changes of stomatal control that accompany deforestation.

52. Delworth, Thomas L., Anthony J Broccoli, Anthony Rosati, Ronald J Stouffer, V Balaji, J A Beesley, William F Cooke, Keith W Dixon, John P Dunne, Krista A Dunne, J W Durachta, Kirsten L Findell, Paul Ginoux, Anand Gnanadesikan, C Tony Gordon, Stephen M Griffies, Richard G Gudgel, Matthew J Harrison, Isaac M Held, Richard S Hemler, Larry W Horowitz, Stephen A Klein, Thomas R Knutson, P J Kushner, Amy R Langenhorst, Hyun-Chul Lee, Shian-Jiann Lin, Jian Lu, Sergey Malyshev, P C D Milly, V Ramaswamy, J L Russell, M Daniel Schwarzkopf, Elena Shevliakova, Joseph J Sirutis, Michael J Spelman, William F Stern, Michael Winton, Andrew T Wittenberg, Bruce Wyman, Fanrong Zeng, and Rong Zhang, 2006: GFDL's CM2 Global Coupled Climate Models. Part I: Formulation and Simulation Characteristics. Journal of Climate, 19(5), DOI:10.1175/JCLI3629.1 .
    [ Abstract ]

    The formulation and simulation characteristics of two new global coupled climate models developed at NOAA's Geophysical Fluid Dynamics Laboratory (GFDL) are described. The models were designed to simulate atmospheric and oceanic climate and variability from the diurnal time scale through multicentury climate change, given our computational constraints. In particular, an important goal was to use the same model for both experimental seasonal to interannual forecasting and the study of multicentury global climate change, and this goal has been achieved. Two versions of the coupled model are described, called CM2.0 and CM2.1. The versions differ primarily in the dynamical core used in the atmospheric component, along with the cloud tuning and some details of the land and ocean components. For both coupled models, the resolution of the land and atmospheric components is 2° latitude × 2.5° longitude; the atmospheric model has 24 vertical levels. The ocean resolution is 1° in latitude and longitude, with meridional resolution equatorward of 30° becoming progressively finer, such that the meridional resolution is 1/3° at the equator. There are 50 vertical levels in the ocean, with 22 evenly spaced levels within the top 220 m. The ocean component has poles over North America and Eurasia to avoid polar filtering. Neither coupled model employs flux adjustments. The control simulations have stable, realistic climates when integrated over multiple centuries. Both models have simulations of ENSO that are substantially improved relative to previous GFDL coupled models. The CM2.0 model has been further evaluated as an ENSO forecast model and has good skill (CM2.1 has not been evaluated as an ENSO forecast model). Generally reduced temperature and salinity biases exist in CM2.1 relative to CM2.0. These reductions are associated with 1) improved simulations of surface wind stress in CM2.1 and associated changes in oceanic gyre circulations; 2) changes in cloud tuning and the land model, both of which act to increase the net surface shortwave radiation in CM2.1, thereby reducing an overall cold bias present in CM2.0; and 3) a reduction of ocean lateral viscosity in the extratropics in CM2.1, which reduces sea ice biases in the North Atlantic. Both models have been used to conduct a suite of climate change simulations for the 2007 Intergovernmental Panel on Climate Change (IPCC) assessment report and are able to simulate the main features of the observed warming of the twentieth century. The climate sensitivities of the CM2.0 and CM2.1 models are 2.9 and 3.4 K, respectively. These sensitivities are defined by coupling the atmospheric components of CM2.0 and CM2.1 to a slab ocean model and allowing the model to come into equilibrium with a doubling of atmospheric CO2. The output from a suite of integrations conducted with these models is freely available online (see http://nomads.gfdl.noaa.gov/). Manuscript received 8 December 2004, in final form 18 March 2005

53. Hurtt, G C., S Frolking, M G Fearon, B Moore, Elena Shevliakova, Sergey Malyshev, S W Pacala, and R A Houghton, 2006: The underpinnings of land-use history: three centuries of global gridded land-use transitions, wood-harvest activity, and resulting secondary lands. Global Change Biology, 12(7), DOI:10.1111/j.1365-2486.2006.01150.x .
    [ Abstract ]

    To accurately assess the impacts of human land use on the Earth system, information is needed on the current and historical patterns of land-use activities. Previous global studies have focused on developing reconstructions of the spatial patterns of agriculture. Here, we provide the first global gridded estimates of the underlying land conversions (land-use transitions), wood harvesting, and resulting secondary lands annually, for the period 1700–2000. Using data-based historical cases, our results suggest that 42–68% of the land surface was impacted by land-use activities (crop, pasture, wood harvest) during this period, some multiple times. Secondary land area increased 10–44 × 106 km2; about half of this was forested. Wood harvest and shifting cultivation generated 70–90% of the secondary land by 2000; permanent abandonment and relocation of agricultural land accounted for the rest. This study provides important new estimates of globally gridded land-use activities for studies attempting to assess the consequences of anthropogenic changes to the Earth's surface over time.

54. Anderson, Jeffrey L., V Balaji, Anthony J Broccoli, William F Cooke, Thomas L Delworth, Keith W Dixon, Leo J Donner, Krista A Dunne, Stuart Freidenreich, Stephen T Garner, Richard G Gudgel, C Tony Gordon, Isaac M Held, Richard S Hemler, Larry W Horowitz, Stephen A Klein, Thomas R Knutson, P J Kushner, Amy R Langenhorst, Ngar-Cheung Lau, Zhi Liang, Sergey Malyshev, P C D Milly, Mary Jo Nath, Jeff J Ploshay, V Ramaswamy, M Daniel Schwarzkopf, Elena Shevliakova, Joseph J Sirutis, Brian J Soden, William F Stern, Lori T Sentman, R John Wilson, Andrew T Wittenberg, and Bruce Wyman, December 2004: The New GFDL global atmosphere and land model AM2–LM2: Evaluation with prescribed SST simulations. Journal of Climate, 17(24), 4641-4673.
    [ Abstract ]

    for climate research developed at the Geophysical Fluid Dynamics Laboratory (GFDL) are presented. The atmosphere model, known as AM2, includes a new gridpoint dynamical core, a prognostic cloud scheme, and a multispecies aerosol climatology, as well as components from previous models used at GFDL. The land model, known as LM2, includes soil sensible and latent heat storage, groundwater storage, and stomatal resistance. The performance of the coupled model AM2–LM2 is evaluated with a series of prescribed sea surface temperature (SST) simulations. Particular focus is given to the model's climatology and the characteristics of interannual variability related to E1 Niño– Southern Oscillation (ENSO). One AM2–LM2 integration was performed according to the prescriptions of the second Atmospheric Model Intercomparison Project (AMIP II) and data were submitted to the Program for Climate Model Diagnosis and Intercomparison (PCMDI). Particular strengths of AM2–LM2, as judged by comparison to other models participating in AMIP II, include its circulation and distributions of precipitation. Prominent problems of AM2– LM2 include a cold bias to surface and tropospheric temperatures, weak tropical cyclone activity, and weak tropical intraseasonal activity associated with the Madden–Julian oscillation. An ensemble of 10 AM2–LM2 integrations with observed SSTs for the second half of the twentieth century permits a statistically reliable assessment of the model's response to ENSO. In general, AM2–LM2 produces a realistic simulation of the anomalies in tropical precipitation and extratropical circulation that are associated with ENSO.

55. Hurtt, G C., S W Pacala, and Elena Shevliakova, et al., February 2002: Projecting the future of the U.S. carbon sink. Proceedings of the National Academy of Sciences, 99, DOI:10.1073/pnas.012249999 .
    [ Abstract ]

    Atmospheric and ground-based methods agree on the presence of a carbon sink in the coterminous United States (the United States minus Alaska and Hawaii), and the primary causes for the sink recently have been identified. Projecting the future behavior of the sink is necessary for projecting future net emissions. Here we use two models, the Ecosystem Demography model and a second simpler empirically based model (Miami Land Use History), to estimate the spatio-temporal patterns of ecosystem carbon stocks and fluxes resulting from land-use changes and fire suppression from 1700 to 2100. Our results are compared with other historical reconstructions of ecosystem carbon fluxes and to a detailed carbon budget for the 1980s. Our projections indicate that the ecosystem recovery processes that are primarily responsible for the contemporary U.S. carbon sink will slow over the next century, resulting in a significant reduction of the sink. The projected rate of decrease depends strongly on scenarios of future land use and the long-term effectiveness of fire suppression.

56. Shevliakova, Elena, 2002: Modeling Potential Impacts of Climate Change on the Spatial Distribution of Vegetation in the United States with a Probabilistic Biogeography Approach In Wildlife Responses to Climate Change: North American Case Studies, Washington, DC, Island Press, 251-276.
57. Morgan, M G., L F Pitelka, and Elena Shevliakova, May 2001: Elicitation of Expert Judgments of Climate Change Impacts on Forest Ecosystems. Climatic Change, 49(3), DOI:10.1023/A:1010651300697 .
    [ Abstract ]

    Detailed interviews were conducted with 11 leading ecologists to obtain individual qualitative and quantitative estimates of the likely impact of a2 × [CO2] climate change on minimally disturbed forest ecosystems. Results display a much richer diversity of opinion than is apparent in qualitative consensus summaries, such as those of the IPCC.Experts attach different relative importance to key factors and processes such as soil nutrients, fire, CO2 fertilization, competition, and plant-pest-predator interactions. Assumptions and uncertainties about future fire regimes are particularly crucial. Despite these differences,most of the experts believe that standing biomass in minimally disturbed Northern forests would increase and soil carbon would decrease. There is less agreement about impacts on carbon storage in tropical forests. Estimates of migration rates in northern forests displayed a range of more than four orders of magnitude. Estimates of extinction rates and dynamic response show significant variation between experts. A series of questions about research needs found consensus on the importance of expanding observational and experimental work on ecosystem processes and of expanding regional and larger-scale observational, monitoring and modeling studies. Results of the type reported here can be helpful in performing sensitivity analysis in integrated assessment models, as the basis for focused discussions of the state of current understanding and research needs, and, if repeated over time, as a quantitative measure of progress in this and other fields of global change research.

58. Pacala, S W., G C Hurtt, Elena Shevliakova, Songmiao Fan, and Jorge L Sarmiento, et al., June 2001: Consistent land-and atmosphere-based US carbon sink estimates. Science, 292(5525), DOI:10.1126/science.1057320 .
    [ Abstract ]

    For the period 1980–89, we estimate a carbon sink in the coterminous United States between 0.30 and 0.58 petagrams of carbon per year (petagrams of carbon = 1015 grams of carbon). The net carbon flux from the atmosphere to the land was higher, 0.37 to 0.71 petagrams of carbon per year, because a net flux of 0.07 to 0.13 petagrams of carbon per year was exported by rivers and commerce and returned to the atmosphere elsewhere. These land-based estimates are larger than those from previous studies (0.08 to 0.35 petagrams of carbon per year) because of the inclusion of additional processes and revised estimates of some component fluxes. Although component estimates are uncertain, about one-half of the total is outside the forest sector. We also estimated the sink using atmospheric models and the atmospheric concentration of carbon dioxide (the tracer-transport inversion method). The range of results from the atmosphere-based inversions contains the land-based estimates. Atmosphere- and land-based estimates are thus consistent, within the large ranges of uncertainty for both methods. Atmosphere-based results for 1980–89 are similar to those for 1985–89 and 1990–94, indicating a relatively stable U.S. sink throughout the period.

59. Peterson, G, and Elena Shevliakova, et al., December 1997: Uncertainty, climate change, and adaptive management. Conservation Ecology, 1(2), Article 4.
    [ Abstract ]

    http://www.ecologyandsociety.org/vol1/iss2/art4/

Direct link to page: http://www.gfdl.noaa.gov/bibliography/results.php?author=1583