Modeling of ecosystem dynamics and carbon cycling with coupled - - PowerPoint PPT Presentation

ENVIROMIS -2006, 1-8 July 2006, Tomsk, Russia

Modeling of ecosystem dynamics and carbon cycling with coupled climate model - LPJ-DGV model. .

Krupchatnikoff V.

( (ICMMG SB RAS/SibRHMRI, Novosibirsk, e-mail: vkrup@ommfao1.sscc.ru)

ACKNOWLEDGMENTS I would like to thank S. Sitch (PIK, Potsdam) who provide code of LPJ model. I am grateful to E. Volodin for providing the results of run with CM3.0 model. Finally, this effort is supported by RFBR grant number 05-05-64989

Introduction

n

Climate models are powerful tools for assessing future climate and interpreting observations. Model simulations can be diagnosed to assess contributing causes of particular phenomena. Climate is not controllable system. The models are only available tools for forecasting plausible climates in years and decades ahead. The biosphere plays pivotal role in the climate system. Land ecosystems are a dynamic component of the global carbon and water cycle. More than one third

• f CO2 in the atmosphere is exchange annually with the land
• biosphere. Land ecosystem processes also control exchange of energy

and momentum between atmosphere and surface. We must better understand the global evolution land ecosystem as main agent in Earth Climate System dynamics. The need has been driving the development

• f dynamic global vegetation models. Lund-Potsdam-Jena dynamic

global vegetation model coupled with INM/RAS climate model is described in this presentation.

Model Description

n

CGCM/INM RAS 5x4 horizontal resolution and 21-level vertical resolution

n

LSM/ICMMG SB RAS - biophysical and biochemical surface model

n

Lund-Potsdam-Jena dynamic global vegetation model (line-off)

n

Terrain-following vertical coordinate (21 σ-levels)

n

Semi-implicit formulation of integration in time

n

Energy conservation finite-difference schemes (5x 4) (Arakawa-Lamb,1981)

n

Convection (deep, middle, shallow)

n

Radiation (H2O, CO2, O3, CH4, N2O, O2; 18 spectral bands for SR and 10 spectral bands for LR)

n

PBL (5 σ-levels)

n

Gravity wave drag over irregular terrain

Coupled atmospheric - ocean model (INM/RAS):

Land surface model(ICM&MG/SB RAS):

n

Vegetation composition, structure

n

Radiative fluxes

n

Momentum and energy fluxes

n

Vegetation and ground temperature

n

Soil and lake temperature

n

Surface hydrology (snow, runoff, soil water, canopy water etc.)

n

CO2 emissions from terrestrial vegetation

n

CH4 emissions from natural wetlands

Leaf area index

Grid structure in land surface model

Global Net CO2 fluxes(mmol CO2/m^2 c), (coupled simulation)

Lund-Potsdam-Jena dynamic global vegetation model (S. Sitch et al., //Global ChangeBiology,(2003)9,161-185)

n

Plant functional types (10 PFT)

n

Annual vegetation and carbon dynamics

n

Penology

n

Production

• water availability
• photosynthesis
• respiration
• reproduction
• allocation
• mortality

n

Input data (monthly mean temperature, precipitation and cloud cover)

LPJ dynamic global vegetation model

• is able to reproduce carbon and water exchange with

atmosphere on seasonal time scale;

• is able to simulate spatial distributions of soil, litter

and vegetation pools and NPP, runoff within their accepted ranges and agree with observed patterns

• is able to reproduce global vegetation distribution in

general agreement with satellite derived maps of phenology and leaf type

• is able to evaluate seasonal cycle of CO2

Plant functional type

LPJ-simulated dominant PFT

Model were driven by CM3.0(INM) (scenario A2)

A2 Storyline and Scenario Family

The A2 scenario family represents a differentiated world. Compared to the A1

storyline it is characterized by lower trade flows, relatively slow capital stock turnover, and slower technological change. The A2 world "consolidates" into a series of economic regions.

Self-reliance in terms of resources and less emphasis on economic, social, and cultural interactions between regions are characteristic for this future. Economic growth is uneven and the income gap between now-industrialized and developing parts of the world does not narrow, unlike in the A1 and B1 scenario families. The A2 world has less international cooperation than the A1 or B1 worlds. People, ideas, and capital are less mobile so that technology diffuses more slowly than in the other scenario families. High-income but resource-poor regions shift toward advanced post-fossil technologies (renewables or nuclear), while low-income resource-rich regions generally rely on older fossil technologies.. As in other SRES storylines, the intention in this storyline is not to imply that the underlying dynamics of A2 are either good or bad.

INM-simulation cloudiness, 2010, dcld = cld(2040) – cld(2010) and dcld = cld(2080) – cld(2040)

INM-simulation of temperature 2010, dTemp = T(2040) - T(2010) and dTemp = T(2080) – T(2040)

INM-simulation of precipitation 2010, dPrec = Prec(2040) - Prec(2010) and dPrec = Prec(2080) – Prec(2040)

Increasing of CO2 concentration (ppm) according to scenario A2

Table of PFT

• TrBE - tropical broad-leaved evergreen
• TrBR - tropical broad-leaved raingreen
• TeNE – temperate needle – leaved evergreen
• TeBE - temperate broad – leaved evergreen
• TeBS - temperate broad – leaved summergreen
• BoNE – boreal needle – leaved evergreen
• BoNS - boreal needle – leaved summergreen
• BoBS - boreal broad – leaved summergreen
• TeH - temperate herbaceous
• TrH - tropical herbaceous

LPJ/INM(RAS)-simulation (scenario A2) Temperate(TeBS) (black), Boreal(BoNE) (green) forest dynamics for grid cell(60E,60N)

LPJ/INM(RAS)-simulation (scenario A2) Boreal(BoBS) forest (green) and Temperate(TeH) herbaceous (yellow) dynamics for grid cell(60E,60N)

LPJ/INM(RAS)-simulation (scenario A2) of carbon pools, (60E,60N) vegc – green, soilc – black, litterc - yellow

LPJ/INM(RAS)-simulation (scenario A2) flux C due to fire(gC/m^2),(60E,60N)

LPJ/INM(RAS)-simulation (scenario A2) of carbon pools(left) and flux C due to fire(gC/m^2),(20E,40N)

LPJ/INM(RAS)-simulation (scenario A2) Temperate(TeNE,TeBE, TeBS), forest and Temperate(TeH),Tropical(TrH) herbaceous dynamics for grid cell(20E,40N)

LPJ/INM(RAS)-simulation (scenario A2) of carbon pools (g C/m^2), (90E,60N)

LPJ/INM(RAS)-simulation (scenario A2) flux C due to fire(gC/m^2),(90E,60N)

LPJ/INM(RAS)-simulation (scenario A2) Annual grid cell NPP, NEP = Firec + RH - NPP, RH(gC/m^2/month) (90E,60N)

LPJ/INM(RAS)-simulation (scenario A2) Boreal (BoNE, BoBS) forest dynamics for grid cell(90E,60N)

LPJ/INM(RAS)-simulation (scenario A2) Temperate(TeH) herbaceous dynamics for grid cell(90E,60N)

Conclusion

LPJ model includes

• explicit representation of vegetation structure, dynamics and competition

between PFT over greed cell, soil biogeochemistry;

• LPJ run of modern climate are in agreement with other models;
• The study showed that simulated PFT after 2050 year (at ~ 500ppm)

begin to lose stability and reaches new balance.

“Those who have knowledge, do not predict, Those who predict, do not have knowledge.” (Lao Tzu)