Green roofs are promoted to provide ecosystem services and to mitigate climate change in urban areas. This is largely due to their supposed benefits for biodiversity, rainwater management, evaporative cooling, and carbon sequestration. One scientific challenge is quantifying the various contributions of green roofs using reliable methods. In this context, the green roof module already running in the Town Energy Balance urban canopy model for water and energy exchanges was improved by implementing CO2 fluxes and carbon sequestration potential. This parameterization combines the Interactions between Soil–Biosphere–Atmosphere (ISBA) photosynthesis, biomass, and soil respiration module with the green roof module in order to quantify the net CO2 amount emitted or fixed by the green roof over a time period. Measurement data from an extensive non-irrigated Sedum green roof located at Berlin Brandenburg (BER) Airport in Germany from 2016 to 2020 are used to calibrate and evaluate the parameterization. Based on the 5 years of measurements, a sensitivity analysis is conducted to quantify the significance of selected parameterization parameters for the photosynthesis process, followed by a calibration over the most important parameters and an evaluation. Results show good agreement of the simulated leaf area index and CO2 fluxes with in situ observations, with good diurnal, seasonal, and inter-annual variability, although the model tends to be overly responsive on the day-to-day variability. The model effectively reproduces net ecosystem exchange, which provides a reliable estimation of the annual carbon sequestration. These results are encouraging in terms of quantifying the potential of carbon sequestration by green roofs and open up the possibility of applying the new parameterization on a city-wide scale to evaluate green roof scenarios.
Agence Nationale de la RechercheANR-22-CE92-0001-01Région Occitanie Pyrénées-Méditerranée00087377/21012846Introduction
Green roofs refer to roofs with a vegetated surface on top of a growing layer. They are mainly divided into two categories: extensive green roofs, which are made with shallow substrate and low-profile plant species and do not necessarily require irrigation, and intensive green roofs, which can support shrubs and trees and require irrigation and deeper substrate. Recent studies have investigated the various advantages of both green roof types. Different types of impacts and interactions have been studied: the effects on buildings' energy saving at building and city scales and under different climates ; the effects on water runoff quantity and quality ; interactions with plants, animals, and abiotic environments ; and the benefits for air quality .
A further advantage is the reduction in carbon emissions in two ways, as pointed out by . First, an indirect reduction occurs due to less carbon being emitted as a result of building energy savings, depending on the energy generation source. In addition, a direct reduction comes from the carbon sequestration by the soil and vegetation on the green roof. From this perceptive, recent studies, like the work of , estimated direct carbon sequestration in containers with dry-weight measurement for Sedum acre, Frankenia thymifolia, and Vinca major species at 38, 565, and 166 gCm-2yr-1, respectively. In addition, indirect carbon sequestration was quantified by modelling at 7680, 7222, and 6393 gCm-2yr-1, respectively, for the same species. , also with dry-weight measurements, estimated direct carbon sequestration by an extensive green roof for irrigated Cynodon dactylon, Festuca arundinacea, Zoysia matrella, irrigated Sedum aizoon, and non-irrigated Sedum aizoon at 690, 751, 671, 459, and 336 gCm-2yr-1, respectively. They also evaluated the payback of time for a non-irrigated extensive Sedum aizoon green roof ranging between 8.5 and 14.0 years, suggesting that CO2 sequestration could be a real positive effect of green roofs. However, as highlighted by , most of the studies quantifying potential carbon sequestration rely on short-term observations.
In addition, in order to really determine the feasibility and interest of such installations, impact studies must be conducted considering all interactions on a city-wide scale and over long enough time periods to cover a variety of climatic conditions. This illustrates the relevance of representing green roof fluxes, including CO2 exchanges, in appropriate urban land surface models, such as the Town Energy Balance TEB; model, which can be run in various configurations: with or without complex atmospheric coupling; from the street level to the neighbourhood, city, or regional scale; and for specific meteorological events or seasonal (and even multi-annual) time periods.
The TEB model already includes an extensive green roof module for addressing heat, energy, and water exchanges . The model has already been applied in the Paris region to quantify the benefits of green roofs in terms of urban cooling, improved thermal comfort, and energy savings . Here we improve the existing model by adding the calculation of CO2 fluxes for green roofs. The aim of this article is to provide a full description and evaluation of the added CO2 fluxes in TEB's extensive green roof module in order to have a model that is able to quantify the CO2 sequestration of green roofs with its annual and inter-annual variations. This is done by reusing the model called Interactions between Soil–Biosphere–Atmosphere ISBA; , which represents CO2 fluxes of the soil and vegetation. But, in order to fit with the specificity of green roofs (shallow soil and Sedum vegetation), a new parameterization of ISBA is achieved with calibration and evaluation on experimental data collected for several years over an extensive green roof with Sedum species on top of an airport car park in Berlin, Germany.
The experimental data are presented in Sect. . Subsequently, the TEB model, the TEB-GREENROOF module, and the implementation of CO2 fluxes are developed in Sect. , followed by a description of the model configuration and numerical setup for the case study in Sect. . Then, the sensitivity analysis, calibration, and evaluation of the new TEB-GREENROOF version including CO2 fluxes with the observational data are presented in Sect. . Finally, Sect. is a discussion about the behaviour of the Sedum simulated with the new parameterization, the diurnal cycle of CO2 fluxes, and the quantification of the amount of carbon fixed by the green roof.
Instrumented green roof for model development and evaluation
The modelling is informed and evaluated by comparison with continuous observations collected on a green roof site for several years by the Technische Universität Braunschweig . The site is a non-irrigated extensive green roof of 8600 m2 (see Appendix A) constructed in May 2012. It is located on the flat roof of an 18 m high car park at Berlin Brandenburg Airport in Germany (referred to as BER; 52.37° N, 13.51° E; altitude of 61 m above sea level). The green roof is composed of four layers: (1) a vegetation layer made up mainly of Sedum (Sedum floriferum `Weihenstephaner Gold', Sedum album) with herbaceous plants (Allium schoenoprasum, Trifolium sp.); (2) a 0.09 m deep substrate layer composed of a mix of expanded shale, pumice, and compost; (3) a 0.003 m thick protection mat; and (4) a 0.05 m thick insulation layer. It is supported on a 1.60 cm layer of ferroconcrete. Gardeners provide basic maintenance of the roof vegetation approximately once a year.
The site is equipped with a 3D ultrasonic anemometer and an open-path infrared gas analyser at 1.15 m above roof level to determine the net CO2 fluxes and the turbulent latent (LE) and sensible (H) heat fluxes using the eddy-covariance technique. The site is also equipped with radiometers on top of the green roof horizontally to measure both downwelling (received by the green roof) and upwelling (reflected and emitted by the green roof) components of long-wave (LW↓, LW↑) and short-wave (SW↓, SW↑) radiation. Air temperature and humidity are measured at 2 m above the green roof surface by an HMP155 probe, and precipitation data are gathered from the nearby German Weather Service station at Berlin Schönefeld Airport (BSCH; ID: 00427). Probes are also placed in the substrate layer to measure the soil temperature and water content at different depths (0.025, 0.05, and 0.075 m). The fraction of vegetation cover (Fveg) and the leaf area index (LAI) are estimated occasionally (≃10 times a year) through photograph analysis taken at 10 different locations randomly selected on the roof. For details, the reader is referred to .
Evolution of meteorological conditions over the period 2016–2020: daily maximal temperature (TX, in °C), daily minimal temperature (TN, in °C), and daily cumulative precipitation (black bars, in mm). The data come from the German weather station at Berlin Schönefeld Airport (BSCH; ID: 00427).
The site has been in operation since mid-2014. This study focuses on the period 2016–2020, for which a very comprehensive dataset is available. According to the Köppen climate classification , Berlin is located in a region of temperate oceanic climate (Cfb). Winters are cold, and summers are warm and humid. The rainfall pattern indicates moderate rainfall throughout the year annual average of 591 mm according to, with climatological maxima in June and July. Snowfall is typical from December to March. Figure presents the monthly values for the daily minimum and maximum temperatures (TN and TX, respectively) and precipitation for the period of interest. The 5 years selected show contrasting meteorological conditions with a noticeable inter-annual variability in temperature and precipitation. The year 2018 was a particularly dry year, with total rainfall well below normal (only 66 mm recorded in summer compared to 157 mm on average for summer for the entire period), and also slightly warmer (average summer TX of 26.0 °C, compared to 25.0 °C on average for all the period). The year 2019 was also slightly drier and, above all, warmer than average (120 mm of precipitation and average TX of 26.4 °C in summer). In contrast, 2017 was wetter and colder in summer (268 mm of precipitation and an average TX of 23.5 °C).
Model description and implementation Overview of the TEB urban canopy model
The Town Energy Balance (TEB) model is integrated in the open-access SURFEX land surface modelling platform together with other dedicated surface models such as ISBA for vegetation and natural soils and FLake for inland waters .
The TEB model is a surface scheme developed by to represent heat, water, and momentum exchanges between urban surfaces and the atmosphere and to compute street-level microclimatic conditions. TEB can be run at street or city scale with meteorological forcings or coupled to an atmospheric model due to its simplified geometry.
In TEB, the urban geometry is represented by the concept of street canyon (Oke, 1987). This hypothesis considers that urban areas can be roughly represented as a single road between two facing buildings of the same height and infinite length and with flat roofs. With this geometry, the model computes the radiation, energy, and water balances on each surface of the street canyon (wall, road, and roof) and aggregates the fluxes to simulate the exchanges between the overall urban canopy layer and the atmosphere above. To better describe the heterogeneity of the urban environment in TEB, the interactions between urban vegetation (i.e. ground vegetation and street trees) and built-up surfaces are now represented within the street canyon . To model the functioning of urban vegetation, TEB relies on the Interactions between Soil–Biosphere–Atmosphere (ISBA) model.
Initial version of the TEB-GREENROOF module
The TEB-GREENROOF module was developed in TEB to allow for the simulation of extensive green roofs. It is also based on the physics of ISBA in order to describe the soil and vegetation layers of the green roof and model its hydrological and thermal performances in interacting with the building on which it is installed, as well as exchanges with the atmosphere above. Figure shows the different fluxes simulated by the original version of TEB-GREENROOF for the hydrological and thermal processes (panels a and b). At surface level, the thermal balance is estimated according to the following equation:
cs∂TS∂t=Q*-H-LE-G0,
where Q* is the net radiation (Wm-2), H is the sensible heat flux (Wm-2), LE is the latent heat flux (Wm-2), G0 is the surface ground heat flux (Wm-2), cs is the surface soil heat capacity (JK-1m-3), and Ts is the surface temperature (K).
(a) TEB-GREENROOF hydrological processes: water gain from precipitation (RAIN), surface runoff (RUNOFF), drainage (K), total drainage (Ktotal), vertical water fluxes (F), ground evaporation (Eg), evaporation of water intercepted by vegetation (EV), and vegetation transpiration (Etr), all in kgm-2s-1. (b) TEB-GREENROOF thermal processes: net radiation (Q*), latent heat flux (LE), sensible heat flux (H), and ground heat flux (G), all in Wm-2. (c) TEB-GREENROOF CO2 processes: gross primary production (GPP) by photosynthesis, leaf respiration (Rleaf), and soil respiration (Rsoil).
The thermal balance connects to the water balance through the latent heat flux (LE). The evaporation on the vegetated surface E is the sum of the evaporation of the soil (Eg) and the evaporation of the vegetation (Eveg), E=Eg+Eveg. The evaporation due to vegetation is also split into direct evaporation EV from the fraction δ of the foliage covered by intercepted water and transpiration Etr of the remaining part of the leaves (Eveg=Ev+Etr). Eg, EV, and Etr are calculated with the following ISBA parameterization:
2Eg=(1-Fveg)⋅ρaCHVa(αqsat(Ts)-qa),3Ev=Fveg⋅ρaδRa(qsat(Ts)-qa),4Etr=Fveg⋅ρa1-δRa+Rs(qsat(Ts)-qa),
where Fveg is the fraction of vegetation, ρa is the air density (kgm-3), CH is the turbulent exchange coefficient (kgm-1s-1), Va is the wind speed (ms-1), qsat(Ts) (kgkg-1) is the saturated specific humidity at the temperature Ts, qa is the atmospheric specific humidity at the lowest atmospheric level (kgkg-1), and α is a coefficient that depends on soil moisture. Ra=(CHVa)-1 is the aerodynamic resistance and Rs the canopy surface resistance that depend on atmospheric factors, soil water content, and LAI, which are modelled later with the implementation of CO2 fluxes.
The soil of the green roof is described using the multi-layer diffusion version ISBA-DF , which discretizes the soil into different vertical layers. In each soil layer, the evolution equation of soil temperature is given by the following expression:
cg∂T∂t=∂G∂z+ϕ,
where T is the soil temperature (K), G is the vertical ground heat flux (Jm-2s-1), cg is the soil heat capacity (JK-1m-3), ϕ is a latent heat source/sink resulting from phase transformation of soil water (Jm-3s-1), and z is the soil depth (m).
For the evolution of soil liquid water and soil ice volumetric content (m3m-3), the equations are
6∂wl∂t=-∂F∂z-ϕLfρw-Slρw(wmin≤wl≤wsat),7∂wi∂t=ϕLfρw(0≤wi≤wsat-wmin),
where F is the vertical water flux (ms-1), Lf is the latent heat of water fusion (Jkg-1), ρw is the liquid water density (kgm-3), and Sl represents external sources/sinks for liquid water (kgm-3s-1).
In the soil, the hydrological and thermal balances are coupled through the effective thermal properties of each soil layer j, which change over time as a function of the soil water content. This is the case for both the layer-averaged soil heat capacity (cgj, in JK-1m-3) and thermal conductivity (λj, in kgms-3K-1) for layer j, according to
8cgj=(1-wsatj)cdryj+wljcl+wijci,9λj=Keλsatj+(1-Ke)λdryj,
where wlj and wij are the soil liquid water content and soil ice content in layer j (m3m-3). cdryj is the heat capacity of the dry-soil matrix for layer j, and cl and ci are the heat capacities of liquid water and ice, respectively. Ke is the Kersten number (that is, a dimensionless parameter representing the normalized thermal conductivity as a function of the degree of saturation only), and λsatj and λdryj are the saturated and dry-soil thermal conductivities for layer j, respectively.
Implementation of <inline-formula><mml:math id="M143" display="inline"><mml:mrow class="chem"><mml:msub><mml:mi mathvariant="normal">CO</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> fluxes in TEB-GREENROOF
Similar to the formulation of radiative, thermal, and hydrological exchanges, CO2 fluxes in TEB-GREENROOF are adapted from an existing module in ISBA, and the CO2 fluxes are modelled in TEB-GREENROOF by activating, with some adaptations, an existing module of ISBA designed to represent C fluxes in vegetation and soils.
The main fluxes are plant photosynthesis (or gross primary production, referred to as GPP, in gCm-2yr-1) and ecosystem respiration (RECO, in gCm-2yr-1), which combine leaf respiration (Rleaf, in gCm-2yr-1 from Eqs. B10–B11) and soil heterotrophic respiration (Rsoil, in gCm-2yr-1) (see Fig. c). The difference between these two large fluxes is defined as the net ecosystem exchange, NEE=RECO-GPP. Photosynthesis and leaf respiration are modelled with the A-gs parameterization , based on the assimilation scheme proposed by Jacobs (see Appendix B). This semi-empirical model has the advantage of being simple and suited for both C3 and C4, the most prominent physiological plant groups based on different photosynthetic pathways. The approach chosen considers that light and CO2 atmospheric concentration are the two limiting factors impacting the net photosynthetic rate.
The goal of this study is to adapt some aspects of the present parameterization in order to model Berlin Brandenburg Airport's extensive green roof by representing the behaviour of the dominant Sedum species (which is widely used for extensive green roofs). Sedum are facultative C3–CAM (crassulacean acid metabolism) , meaning that they behave like a C3 plant when they are well watered and switch to a CAM photosynthesis pathway when they lack water. With CAM behaviour, the stomata of the plant open only at night to fix CO2. The photosynthesis takes place during the day, requiring light energy, but the leaf stomata remain closed to prevent water loss through evapotranspiration. Nonetheless, the flux measurements carried out on the green roof at Berlin Brandenburg Airport did not reveal any periods with CAM behaviour, which would be characterized by the absorption of CO2 during the night. Consequently, the question of modelling the photosynthesis mechanisms specific to this particular species in the ISBA model (which does not include CAM behaviour like most large-scale vegetation models) was not investigated further in this study. However, it still requires adaptation of the model with a specific parameterization to match the behaviour of Sedum on a shallow substrate.
First, the maximum net CO2 assimilation (Am,max, in gm-2s-1) and the mesophyll conductance (gmes*, in ms-1) use inhibition functions (see Appendix Eqs. A2–A3) with reference temperature adjusted to better represent optimum temperature. This was done to differentiate optimum assimilation between C3 and C4 plants. However, an analysis of the observed CO2 fluxes measured by eddy covariance did not show any high- or low-temperature inhibition of photosynthesis in the range of temperatures observed. Because we could not find data triggering high- and low-temperature inhibition for these species, we decided to set the inhibition functions to 1, which modifies the equations for Am,max and gmes:
10Am,max(Ts)=Am,max(25)⋅Q10Ts-2510,11gmes*(Ts)=gmes*(25)⋅Q10Ts-2510,
where Am,max(25) is the maximum net CO2 assimilation at 25 °C in mgm-2s-1; Ts is the leaf skin temperature (K), here considered the first layer of soil temperature; and gmes*(25) is the mesophyll conductance at 25 °C. Q10 is the response function defined as a proportional increase in a parameter for a 10° increase.
When soil moisture content drops, plants tend to close their stomata to limit water loss. This is described empirically by a soil water stress function (F2) that simply models the stomata closing and opening when the plant lacks water according to the regulation of the mesophyll conductance:
gmes=gmes*(Ts)⋅F2,
with the normalized soil water stress factor estimated as
F2=∑j=1NΥj⋅wj-wwiltjwfcj-wwiltj,
where N is the number of soil layers, and Υk is the root fraction in layer j. wwiltj is the wilting point, wfcj the field capacity, and wj the soil water content in layer j. As the substrate of the extensive green roof is shallow, thresholds F2min and F2max are prescribed for F2 to prevent photosynthesis cut or variations that are too excessive.
The respiration of the soil (in kgm-2s-1) is estimated with the simple respiration scheme:
Rsoil=4.4×10-8×(13.5+5.4⋅LAI)14⋅w10⋅e0.069⋅(Ts10-25),
where LAI is the leaf area index (m2m-2), w10 is the weighted soil volumetric water content between 0–10 cm depth (%), and Ts10 is the weighted soil temperature between 0–10 cm depth (°C). Since the soil on the green roof is very shallow, water content thresholds w10min and w10max are prescribed during the calibration to prevent extreme values.
Furthermore, the A-gs CO2 assimilation scheme can be either forced by a prescribed LAI or coupled to a biomass scheme, making the LAI evolve over time. In the model, LAI is the variable that reflects the vegetation evolution (i.e. the CO2 accumulation). Making the LAI evolve dynamically during the simulation, rather than prescribing LAI values imposed on the model, enables us to represent the impact of meteorological conditions on vegetation over the long term and thus to simulate more realistic CO2 accumulation scenarios over longer periods.
We use the ISBA NIT biomass scheme that only represents above-ground biomass, separated into three different reservoirs: a leaf biomass reservoir, an above-ground stem biomass reservoir, and a residual reservoir.
Each biomass reservoir (Bi, expressed in kg of dry matter (DM) m−2) follows the same time evolution equation, the principle of which can be expressed as follows :
ΔBiΔt=Ai-Di-Ri,
with Ai being the biomass input and Di and Ri being biomass loss due to mortality, allocation to other reservoirs, and respiration. The input flux for leaf biomass is the carbon assimilation due to photosynthesis (expressed in kgCm-2d-1):
Aleaf=An,day⋅Δt.
After a daily biomass iteration, the LAI is calculated from the leaf biomass and the specific leaf area (SLA; in m2 per kgDM) as follows:
LAI=SLA⋅Bleaf.
Except for crops, ISBA assumes a constant vegetation cover with time. However, extensive green roof vegetation can be patchy, and the vegetation cover may vary in time. To take this heterogeneity and variability into account, the fraction of vegetation cover is estimated by fitting the equation on the estimation of LAI and vegetation cover:
Fveg=1-ea⋅LAI,
where a is a coefficient set on the basis of observations.
All the vegetation parameters required for the CO2 flux modelling are listed in Table . Standard values are available in ISBA for each of these parameters, for both C3 and C4 plants . Since it is challenging to find appropriate ecophysiological data to derive these parameter values for Sedum, we chose to calibrate some of them in Sect. 5.
Configuration of TEB simulationAtmospheric conditions
In this study, TEB is run in an offline simulation configuration over 5 years, from 2016 to 2021. It is applied to a single grid point, i.e. for an average urban canyon whose characteristics are representative of the study site, in particular the properties of the green roof. In offline mode, the time evolution of atmospheric conditions over the urban canyon must be provided to the TEB model at both the given altitude and the time step. The data required are the above urban canopy air temperature, humidity, pressure, and CO2 concentration; the wind speed and direction; precipitation; and the short-wave (direct and scattered) and long-wave incoming radiation.
Here, the meteorological forcing is built with local observations described in Sect. and provided to the model with a time resolution of 30 min. The specific humidity is calculated with the relative humidity, air pressure, and air temperature measurements. The partitioning of the incoming global short-wave radiation into scattered and direct components is made based on the method of . Solid and liquid precipitation rates are determined by disaggregating the total precipitation, first according to the rain/snow precipitation classification directly provided from the nearby German Weather Service station at Berlin Schönefeld Airport (BSCH; ID: 00427). When classification data are missing or indicate rain and snow during the same day, disaggregation is done using the air temperature recorded on the roof by applying a threshold of 273.15 K.
TEB model configuration
The TEB model is run on a 1D grid, with the urban canyon dimensions and properties defined as the average parameters of the measuring site. The height of the building is set to 18 m, which is consistent with the airport parking lot characteristics. The forcing height is determined according to the sensor location on the green roof, i.e. 19.15 m above ground level (a.g.l.) for wind speed and direction and 20 ma.g.l. for air temperature and humidity. TEB runs with a simplified surface-boundary-layer scheme , allowing for the vertical discretization of the atmosphere in the canyon (from the ground to the forcing height) into six layers of microclimate variables. The fraction of building is set to 0.6 and fully covered by green roofs using the TEB-GREENROOF module, the fraction of road is 0.2, and the fraction of low vegetation is 0.2 with the GARDEN module . The road is discretized into two layers, corresponding to a surface layer of artificial coating and a basement layer of natural soil. The parking lot is a concrete building with large openings on the facades. The walls are defined in the model as having two layers of concrete, and the roof consists of a surface layer of concrete and a second layer of insulation.
TEB-GREENROOF input parameters for green roof substrate and vegetation. Geometric parameters are defined according to the soil depth of the site; surface and thermal properties are based on the previous study of , with an adjustment to soil thermal conductivity. For hydrological properties, Bcoef is set based on ; the wilting point, field capacity, and porosity profile are derived from soil water content measurements; the saturated soil matrix is defined on the basis of ex situ analysis; and saturated hydraulic conductivity is based on soil manufacturer documentation.
TypeParameterUnitValuesGeometryNumbers of numerical soil layer(–)6Layer thickness (from top to bottom)cm(0.3, 1.9, 1.9, 1.9, 1.5, 1.5)SurfaceAlbedo of bare soil(–)0.154Emissivity of bare soil(–)0.83Albedo of vegetation(–)0.2Emissivity of vegetation(–)0.83Roughness length for momentumm0.03ThermalDry-soil thermal conductivityWm-1K-10.15propertiesDry-soil heat capacityJm-3K-12 000 000HydrologicalWilting pointm3m-30.001propertiesField capacitym3m-30.20Porosity profilem3m-30.57Matrix potential at saturationm-0.1Bcoef coefficient for water retention curve(–)4.0Hydraulic conductivity at saturationms-12.183×10-3
Description of the input parameters for the calculation of CO2 fluxes in TEB-GREENROOF. All the parameter values tested for the sensitivity analysis, pre-calibration, and final calibration stages are listed. The data in square brackets define the ranges of values tested, and the data in brackets define the pairs of values tested.
ParameterDescriptionUnitSensitivity analysis/pre-calibrationCalibrationBest configA-gs photosynthesis parameters f0*Value of f if there is no saturation deficit (with no soil water stress)(–)[0.45;0.935](0.45,0.55)0.45ϵ0Maximum initial quantum use efficiencykgCO2J-1PAR[12.6×10-9;18.7×10-9](14×10-9,17×10-9)14×10-9Γ(25)CO2 compensation concentrationppmv[2.52;49.5](2.8,55)55Am,max(25)Maximum net CO2 assimilation at 25 °Cmgm-2s-1[1.53×10-6;2.42×10-6](1.53×10-6,2.42×10-6)2.2×10-6gmes*(25)Mesophyll conductance at 25 °Cms-1[10-3;4.0×10-2](0.001,0.002)0.001gcCuticular conductancems-1[0;0.0001]/0.00010.00010.0001DmaxMaximum saturation deficit of atmosphere tolerated by vegetationkgH2O per kg air[0.3;0.6]/0.30.30.3Response to drought parameters F2minMinimum normalized soil water stress factor(–)0.150.150.15F2maxMaximum normalized soil water stress factor(–)[0.2;1](0.35,0.75)0.75Biomass parameters τMMaximal lifespan of leavesd/(75,150)75SLASpecific leaf aream2 per kgDM/12.912.9Respiration parameters w10minThe weighted soil volumetric water content between 0–10 cm depth for respiration(–)0.050.050.05w10maxThe weighted soil volumetric water content between 0–10 cm depth for respiration(–)0.150.150.15aCoefficient for estimating vegetation coverage(–)/-0.35-0.35
Besides the photosynthesis parameters (described in Table ), the input parameters of the TEB-GREENROOF module, including the description of the green roof and the thermal and hydrological properties, are described in Table . The green roof soil is discretized into six layers with the same composition for a total depth of 0.09 m. With regard to hydrological properties, the potential of the soil matrix at saturation is defined based on ex situ analysis. The empirical coefficient for the water retention curve (Bcoef) varies according to soil and connects the water potential to the water content in a soil matrix derived by . It is set based on the results of a previous case study referenced in . The hydraulic conductivity at saturation is based on the soil manufacturer's documentation. The soil porosity profile, field capacity, and wilting point are set directly based on measurements of soil water content on the green roof. The definition of the thermal properties is based on the initial calibration proposed by for a different instrumented extensive green roof: the same value of dry-soil heat capacity is applied, but the dry-soil thermal conductivity is slightly increased to better model the heat conduction in the soil according to measurements (not presented here). The vegetation albedo and emissivity are taken from the previous study of .
Definition of green roof <inline-formula><mml:math id="M288" display="inline"><mml:mrow class="chem"><mml:msub><mml:mi mathvariant="normal">CO</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> flux parameters
Without information on the specific characteristics of the current green roof to be simulated, the input parameters for the photosynthesis model must be defined. The initialization of these parameters is done following three successive steps:
Sensitivity analysis on the main parameters of the A-gs photosynthesis model. This step is based on a standalone and very low computing time version of A-gs (detailed in Appendix C). Forced by predefined microclimate conditions, a very large number of simulations are carried out to test a wide range of parameter values and clarify which parameters are the most significant and have the greatest effect on the calculation of CO2 fluxes.
Pre-calibration. This step consists of running additional simulations, again using the standalone version of A-gs but with a focus on the parameters selected in the previous step in order to narrow the range of plausible values (Appendix E).
Calibration of the selected parameters according to the values identified in the previous step. This time, several combinations of plausible parameter values are tested by running the full TEB configuration considering dynamic LAI. Based on these simulations, the best configuration can be identified.
Sensitivity analysis of ISBA–A-gs parametersSobol index method
In order to assess which ISBA–A-gs parameters are the most predominant in modelling carbon uptake, a global sensitivity analysis is performed on the parameters listed in Table . The approach used is the Sobol method that can be applied for either linear or non-linear models. Two indices are considered here. The first one, Si, is the first-order global sensitivity of output Y to a single selected input parameter Xi, i.e. the variance of Y when Xi is the only parameter that varies. Si is normalized by the variance of Y to obtain a score between 0 (no sensitivity) and 1 (full sensitivity). The second index considered, STi, also known as the total-order index, calculates the variance of output Y to Xi when Xi varies with the other parameters in every possible combination (Xi varies solely, then Xi varies with each parameter, then three parameters including Xi vary together, etc.). As for Si, STi is normalized by the variance of Y with a score between 0 (no sensitivity) and 1.
To compute these two indices for each A-gs parameter, we use a Monte Carlo approach developed by . This involves working with samplings that adequately span the space of possible values for each parameter. In this study, the sampling is performed with Latin hypercube sampling , which is a stratified sampling method aiming to spread the sample points evenly across all possible values. This is done using the R package “lhs” with the geneticLHS method , which samples using a genetic algorithm to maximize the mean distance from each point to all other points. The complete sensitivity analysis process is summarized in Fig. .
Flowchart of sensitivity analysis. This chart is for the first order only and is the same for the total effect index, except for the calculation of the Sobol index box.
Application to the photosynthesis A-gs model
To implement this sensitivity study on the input parameters of the A-gs photosynthesis model, a simplified modelling configuration is developed. The sensitivity analysis is carried out with the standalone A-gs model (Appendix C) on scores of simulated CO2 fluxes. The mean absolute error (MAE) and root-mean-square error (RMSE) are computed over the 5 years of observations selected at a temporal scale of 30 min. Among all input parameters of the A-gs model, eight parameters (Table ) are not set on observations and are selected for the sensitivity analysis. This includes gmes*(25) and Am,max(25) presented in Sect. 3.3, with the cuticular conductance (gc), the maximum saturation deficit of the atmosphere tolerated by vegetation (Dmax), the maximum initial quantum use efficiency (ϵ0), the CO2 compensation concentration (Γ(25)), the value of the f factor when there is no saturation deficit (f0*), and the maximum normalized soil water stress factor (F2max). The range of values for part of the parameters is based on a previous sensitivity analysis from with reference values from the literature: gmes*(25) from , gc from , and Dmax from . For the other parameters, the ranges are set according to the default values defined for types C3 and C4 in ISBA–A-gs, with a ±10 % margin: ϵ0, Γ(25), and Am,max(25) from and f0* from . The parameters and associated value ranges are listed in Table . The distribution within the fixed range of each parameter is defined according to a uniform distribution.
The size of each sampling matrix (the number of different values for one parameter) is set to 5000, which implies a total of 50 000 simulations based on the Monte Carlo method. Here, the output variables of interest for which the Sobol indices are calculated are the mean absolute error (MAE) and the root-mean-square error (RMSE) of the net CO2 fluxes simulated by A-gs compared with the net CO2 flux measurements collected on the instrumented green roof.
Sensitivity analysis results
The results of the sensitivity analysis on the A-gs input parameters for its implementation in the TEB-GREENROOF model for CO2 calculation are presented in Fig. for both MAE and RMSE. The values obtained for the Sobol index make it possible to first identify the parameters with no influence on the net carbon flux scores, i.e. Dmax and gc. As a result, they are set to default values, as listed in Table . All other parameters are retained for the calibration with the same range for the pre-calibration. Among them, gmes*(25) is the most sensitive one, followed by Γ(25) and ϵ0, with similar results for MAE and RMSE. As a result, these three input parameters have the greatest influence on the performance of the A-gs model in calculating CO2 fluxes and thus require more meticulous calibration (presented in the next section). The parameters Am,max(25), f0*, and F2max have a smaller but non-negligible impact so that they are also retained for the pre-calibration.
Comparison of the Sobol first-order (red) and total-order (blue) indexes calculated for the eight variables in the sensitivity analysis, for both RMSE (left bar) and MAE (right bar).
Calibration of A-gs parametersMethod
The model is calibrated based on the observation dataset of net CO2 fluxes divided into two distinct time periods: a 4-year calibration period (from 2016 to 2019) and a 1-year evaluation period (2020). Pre-calibration (Appendix E) is done to reduce the range of values for the six photosynthesis parameters to be calibrated (gmes*(25), Γ(25)ϵ0, Am,max(25), f0*, and F2max).
An ensemble of TEB-GREENROOF simulations is then run according to the full configuration presented in Sect. (and considering a dynamic LAI). Each simulation represents a different combination of parameter values. For each of the six photosynthesis parameters, two values are tested based on the pre-calibration results (Table ). Since TEB-GREENROOF simulations are performed with a dynamic LAI, additional parameters related to the biomass evolution need to be calibrated, i.e. the maximal lifespan of leaves (τM) and the specific leaf area (SLA). For Sedum album, the value of SLA is set to 12.9 m2 per kgDM according to the TRY database , which brings together the different plant trait databases worldwide (the last version of TRY contains more than 1 million trait records for 6.24 million individual plants). The parameter τM is tested for the two values of 75 and 150 d, with 150 being the default value in ISBA for C3 and C4 plants and 75 being half of this value. In total, 128 combinations of parameters are simulated with TEB-GREENROOF. To identify the best configuration, the statistical scores of the 128 experiments calculated over 2016–2019 are compared by crossing the mean absolute error (MAE) of net CO2 fluxes and the root-mean-square error (RMSE) of LAI.
Calibration results
Figure illustrates the outcomes of all simulations in terms of MAE for the net CO2 fluxes and RMSE for LAI, with details of the values of each parameter in all simulations. In accordance with the findings of the sensitivity analysis, the mesophyll conductance at 25 °C is identified as the most influential parameter. In Fig. c nearly all simulations with gmes*(25)=0.001ms-1 demonstrate better performances than simulations with gmes*(25)=0.002ms-1 for both net CO2 fluxes and LAI. The maximal lifespan of leaves (τM) is also found to have a significant impact. However, as evidenced by the two distinct curves (Fig. a), this parameter primarily affects the representation of LAI, with systematically better RMSE for τM=75 d than for 150 d. For Γ(25) and ϵ0 (Fig. d and f), the two values tested for each parameter give quite comparable performances. It is nonetheless possible to identify a single value, considered the best, set to 55 ppmv for Γ(25) and 1.4 kgCO2PAR for ϵ0. For the other parameters, namely Am,max(25), f0*, and F2max (Fig. b, e, and g), it is more difficult to conclude on the values leading to the best configuration. Indeed, the differences in performances between the simulations performed with the two values of Am,max(25) are very low, which corresponds well with the sensitivity analysis results. For f0* and F2max, the simulations with the best RMSE for LAI do not correspond to the simulations with the best MAE for CO2 fluxes. The best configuration is highlighted in Fig. a, with the corresponding parameter values listed in Table . It is selected because it gives the best simulation for LAI and still close to the best simulation for FCO2 MAE.
Comparison of the performances of the 128 simulations combining the RMSE of daily LAI (x axis) and the MAE of 30 min average instantaneous net CO2 fluxes (y axis). Panel (a) highlights the simulation selected as the best (black dot) compared to the other simulations (grey dots). For each of the other panels corresponding to one of the variables to be calibrated (b–h), the red/blue colour code distinguishes the performance of the two values tested for that variable.
Evaluation of the TEB-GREENROOF
The final simulation (i.e. the evaluation simulation) is based on the best configuration obtained from the calibration stage, with all parameters listed in Table . It is analysed and evaluated here in relation to observations of LAI and CO2 fluxes.
Dynamic LAI modelling
The modelling of the evolving LAI using the biomass model is presented in Fig. and compared with the LAI retrieved from observations. The observed LAI was estimated by capturing the variation in the green chromatic information (green fraction) in the RGB space from photographs taken at 10 different randomly selected locations on the roof, approximately once a month . Since observation-based LAI is estimated from the analysis of discontinuous photographs, it is difficult to assess the accuracy of the model in detail. Nevertheless, some interesting results stand out from the comparison. The temporal evolution of the LAI reveals a good representation of the overall seasonal dynamics for the evaluation year 2020, with a clear identification of the growing period starting around April and ending in September. During this period, the annual maximum LAI is found in July, reaching 3.08 m2m-2 for modelling and 3.71 m2m-2 for estimation. For the calibration period, the modelled annual maximum LAI reaches 2.56, 1.84, and 2.61 m2m-2 for 2016, 2018, and 2019, respectively, which corresponds with the variation in the annual maximum observed LAI at 2.55, 2.49, and 3.11 m2m-2 for the same years. Note that comparison cannot be done for the year 2017 since there is no estimation of LAI after mid-July. With regard to the senescence period, a minimum threshold of LAI is prescribed in the model at 0.5 m2m-2 based on the estimated LAI between 2016–2020 (in particular, according to winter values in 2018–2019; see Fig. ), which still matches reasonably well with what is estimated in 2020. Finally, it is noteworthy that the model is able to simulate inter-annual variability in LAI. For the particularly dry year of 2018, the average LAI simulated during the growing season (April–September) is much weaker (0.9 m2m-2) than that of other years (1.6–1.7 m2m-2).
Daily evolution of the modelled LAI (black line) during the calibration (2016–2019) and evaluation (2020) periods compared to the estimated on-site LAI values (red dots), in m3m-3.
Monthly evolution of daily variation in modelled (a) and observed (b)CO2 fluxes (FCO2) for the evaluation period 2020. The colour range varies from green for negative CO2 fluxes (photosynthesis) to red for positive CO2 fluxes (respiration).
Net <inline-formula><mml:math id="M393" display="inline"><mml:mrow class="chem"><mml:msub><mml:mi mathvariant="normal">CO</mml:mi><mml:mn mathvariant="normal">2</mml:mn></mml:msub></mml:mrow></mml:math></inline-formula> flux modelling
The monthly evolution of the net CO2 flux diurnal cycle simulated and observed over the year 2020 is represented in Fig. . The model provides a good representation of the amplitude of the net CO2 fluxes for the evaluation period. Both the annual cycle and the diurnal cycle are close to the observations. The model is able to capture the net CO2 flux seasonal variations in response to variation in climatic forcing. In August and late April, both the model and the observations catch less photosynthesis due to a drier period of stress for the vegetation. Conversely, the model also reproduces periods in July when photosynthesis is enhanced. However, the model remains excessively responsive to meteorological fluctuations, particularly following precipitation events, when soil water content is high, leading to an overestimation of photosynthesis. This is illustrated after the rain of 13 June 2020, when the daily average (between 06:00 and 18:00 UTC) of the difference between the observation and the model is 3.35, 4.17, 2.93, and 2.42 µmolm-2s-1 for the 14th, 15th, 16th, and 17th, respectively, leading to an overestimation of the flux for the 4 d after the rain event. When this trend reverses, the difference between the observations and the model is then 0.89, -1.57, and -2.35µmolm-2 for the 18th, 19th, and 20th, respectively. Furthermore, outside the growing season, the photosynthesis simulated by the model is close to 0, whereas observations indicate that photosynthesis continues even during winter. This can be seen by taking the average over the meteorological winter period of the daily minimum FCO2. The observed value is 2.53 µmolm-2s-1, while the modelled value is only 0.40 µmolm-2s-1.
The diurnal cycle of the net CO2 fluxes averaged monthly over the 5-year period is illustrated in Fig. . Similarly to what is shown in Fig. , outside the growth period (between November and March), the model does not simulate the diurnal CO2 cycle, which is still noticeable in the observations and reflects weak but still active photosynthesis. However, the amplitude of the net CO2 fluxes is on average quite well represented during the growing period, especially between June and October. During the day, the CO2 fluxes are well represented, although the model tends to be overly responsive compared to the observations, resulting in a greater standard deviation for the modelling. The modelled respiration at night is in close agreement with the measurements, as can be seen between 20:00 and 06:00 UTC, and follows the annual variation well, being greater during the growing period and close to 0 outside the growing period.
Comparison of the modelled (black line) and observed (red line) diurnal cycles of net CO2 fluxes (gs-1m-2) averaged monthly over the 5-year period. The transparent ranges indicate 1 standard deviation.
Comparison of the modelled (black line) and observed (dashed black line) monthly averaged daily NEE (in gCm-2d-1) for the period 2016–2020. For simulation only, the monthly averages of both daily ecosystem respiration (red line) and daily GPP (green) are also presented.
Comparison of modelled and observed annual NEE in gCm-2yr-1 for each year from 2016 to 2020. By definition, a negative NEE corresponds to a quantity of carbon fixed by the green roof.
The net ecosystem exchange (NEE; see Sect. 3.3) quantifies the net sequestration of CO2 if the NEE is negative or the net CO2 emissions if the NEE is positive. Figure represents the observed and modelled daily NEE averaged monthly for the 5 successive years from 2015 to 2020. The modelled NEE is in close agreement with the measurements for all years. It follows the observed annual variation, with a slightly positive NEE from October to March and negative NEE for the rest of the year. This means that the green roof acts as a net carbon emitter in autumn and winter but as a net carbon sink in spring and summer. The partitioning on the modelling between GPP and RECO in Fig. demonstrates that both GPP and RECO follow similar inter-annual variations, being greater in summer than in winter. However, since no partitioning is available between GPP and RECO on observations during this study, it is not possible to confirm the accuracy of the model processes. At the annual scale (Table ), the model shows that the green roof is a net carbon sink, in accordance with measurements, and quantifies the amount of carbon fixed by the green roof within the range of error of the measurement estimated at 16 gCm-2yr-1 fairly well. Indeed, the observed and modelled annual NEE for the evaluation year 2020 is -168 and -166gCm-2yr-1, respectively. The model is also able to capture the inter-annual variations in NEE, which are partly governed by changes in weather conditions, especially in the year 2018, which stands out as a dry year compared with normal years and shows a lower annual NEE than the other years in both observations and the simulation (-85 and -51gCm-2yr-1, respectively). Inversely, the wetter year 2017 presents an observed NEE of -178gCm-2yr-1 that is substantially greater than for all other years. This is partly seen in the modelling results, where the NEE is high and reaches -163gCm-2yr-1 in 2017, but it is in the same range than for the year 2020 (-166gCm-2yr-1).
Discussion
The work presented in the previous section has enabled us to characterize, for the first time, the parameters of the photosynthesis and vegetation growth model for Sedum in a soil–vegetation–atmosphere transfer model. We now discuss in more detail the behaviour of the net CO2 fluxes from green roofs in the TEB model and the perspectives that evolving vegetation opens up for the modelling of green roofs.
Sedum model response to micro-meteorological conditions
The calibration proposed in this study allows for the representation of Sedum behaviour when they do not use a CAM photosynthesis pathway. Here we investigate the response of the new Sedum parameterization to the environmental variables driving photosynthesis (soil temperature, water content, and radiation) compared to the standard C3 and C4 parameterizations. Figure shows the response curves of GPP to the three environmental variables, Ts, PAR, and VWC, in each panel. For each response curve corresponding to one variable, the other environmental variables are fixed. For comparison, the observations are selected within the range of the most extreme values of the three curves (the dashed and dotted lines). The observed GPP displayed in the three figures is estimated from the observation of net CO2 fluxes on the BER green roof site subtracted by the modelled RECO of the best simulation. All observations were conducted within the specified range of each figure, with variables influencing photosynthesis varying between the most and the least favourable values for photosynthesis. For the three curves, the Sedum parameterization performs better than the ISBA C3 and ISBA C4 parameterizations, with a significantly lower photosynthetic rate. The response to volumetric water content (Fig. a) demonstrates that the threshold values for the normalized water stress factor (F2min and F2max) and the weighted soil volumetric water content between 0–10 cm depth for soil respiration (w10min and w10max) were appropriately selected. This is visible on the plateau reached at a VWC of 0.056 m3m-3 under average micro-meteorological conditions (straight line), which correspond well with the observations. With regard to the photosynthetically active radiation (Fig. b), the Sedum parameterization models correctly represent the increase in GPP but also the threshold above which the GPP no longer increases with photosynthetically active radiation. This response is not noted with the ISBA C3 and ISBA C4 parameterizations. For the response to temperature (Fig. c), the range of temperatures observed does not allow determination of the maximum GPP achievable at higher temperatures, beyond which the GPP decreases with temperature.
Evolution of the modelled (lines) GPP with (a) volumetric water content (VWC; m3m-3), (b) photosynthetically active radiation (PAR; Wm-2), and (c) temperature (Ts; °C), for the parameterizations of Sedum, ISBA C3, and ISBA C4 and the comparison against selected observations (dots).
Dynamic vegetation
The use of a vegetation growth model makes it possible to simulate the explicit evolution of LAI. LAI is a key parameter for vegetation, including green roofs, which drives and sizes physiological processes and water, energy, and gas exchanges. For the same micro-meteorological conditions, fluxes can vary considerably depending on LAI and its evolution over time. Therefore, LAI has a direct impact on the energy and water balance through the calculation of evapotranspiration and on the carbon balance through the calculation of GPP. However, quantifying LAI and its evolution over time accurately is a challenge, given that both can vary considerably over the course of a year and from one year to the next, depending on the variability of climatic conditions. At the present study site in particular, LAI varies between 0.05 and 3.71 m2m-2, with significant inter-annual variations as well.
Most green roof models use LAI as an input prescribed variable. The first approach for short-term simulations of several days is to define a fixed LAI value for the duration of the simulation e.g. and compare simulations with different prescribed LAIs to find the best configuration and quantify the sensitivity of green roof performance to LAI variation, as done by and . Another approach is to use on-site LAI measurement, as done in the works of , with direct measurement of LAI or with indirect estimation methods, as in , where LAI is calculated with measurements of soil evapotranspiration. show that taking into account the seasonal evolutive LAI during the simulation, rather than applying a fixed LAI, leads to a significantly more relevant energy balance, resulting in a better quantification of the reduction in building temperature and energy consumption, in particular for long-term simulations. With this aim, they use a model to represent vegetation seasonal changes based on temperature and fixed maximum and minimum LAI using a modified Dickinson semi-empirical equation. The work presented here includes more climatic and environmental factors that impact vegetation growth, namely the soil water content, atmospheric pressure, temperature, humidity, photoactive radiation, and CO2 concentration. In this way, the model can represent and estimate the biomass accumulation or loss due to meteorological changes. In contrast to other works, it enables us to more accurately represent long-term changes in the performance of green roofs, in particular due to the development of vegetation.
Conclusions
A new parameterization for the net CO2 flux calculation has been implemented in the TEB-GREENROOF model using the ISBA–A-gs photosynthesis model with a biomass module and an empirical parameterization of ecosystem respiration. The modelling was informed by observations on an extensive non-irrigated Sedum green roof located at BER Airport. The sensitivity analysis results showed that the main parameters driving CO2 fluxes on the green roof are the mesophyll conductance at 25 °C (gmes*(25)), the CO2 compensation concentration at 25 °C (Γ(25)), and the maximum initial quantum use efficiency (ϵ0), which needed to be calibrated more carefully. The results after calibration showed that the model performed well in reproducing the diurnal cycle during the growing period and its dependence on the variations in meteorological conditions and soil water conditions. Outside the growing period, the model struggled in simulating the weak photosynthesis process that seems to persist based on observations. Nonetheless, as CO2 fluxes remain very low during this period, the impact on the overall quantification of carbon sequestration is limited, leading to a good estimation of the annual NEE and its inter-annual variations. This work also ultimately allowed us to characterize, for the first time, the photosynthesis and growth parameters appropriate for modelling Sedum with the ISBA soil vegetation atmosphere transfer model.
Future development will need to include comparison between the ecophysiological parameters calibrated here and on-site measurements of green roof plant photosynthesis parameters. To study the full carbon cycle of the green roof, the management and maintenance of green roofs need to be addressed, especially with respect to the biomass that can be removed by gardeners. Similarly, the carbon sequestered in the soil needs to be quantified with on-site sampling. On the modelling side, the growth module needs to be improved, and a soil organic carbon module needs to be added in order to be able to perform longer-term studies.
In addition to thermal and hydrological effects, the short-term carbon sequestration can now be added to impact studies in order to have a full picture of the impact of green roofs on the fluxes at city scale and under different climate events and conditions. However, the green roof effect should be looked at not only in terms of fluxes, but also with regard to other effects that cannot be modelled in models like SURFEX, such as the effect on individuals, and biodiversity should also be taken into account, thus requiring researchers to cross-reference the results of different approaches.
Photograph of the green roof experimental plot located on top of the car park of Berlin Brandenburg Airport (Germany).
A-gs equation
The following equations are used in ISBA–A-gs in SURFEX v9. The model uses an empirical light response function of net assimilation (An) to combine the effects of light and CO2 as limiting factors. When light is not limiting, Am is limited by a maximum photosynthetic rate Am,max:
Am=Am,max1-e-gmes*(Ci-Γ)Am,max,
where Ci is the internal leaf CO2 concentration in kgCO2 per kgair. Am,max is the maximum net CO2 assimilation in mgm-2s-1, gmes* is the mesophyll conductance in ms-1, and Γ is the CO2 concentration compensation point in ppmv, defined according to the following equations:
B2Am,max(Ts)=Am,max(25)⋅Q10Ts-2510(1+e0.3(T1-Ts))⋅(1+e0.3(Ts-T2)),B3gmes*(Ts)=gmes*(25)⋅Q10Ts-2510(1+e0.3(T1-Ts))⋅(1+e0.3(Ts-T2)),B4Γ(Ts)=Γ(25)⋅Q10Ts-2510,
where Am,max(25) is the maximum net CO2 assimilation at 25 °C in mgm-2s-1, gmes*(25) is the mesophyll conductance at 25° in ms-1, Ts is the leaf skin temperature in °C, and T1 and T2 are reference temperatures in °C.
The internal CO2 concentration depends directly on the atmospheric CO2 concentration (Eq. ) and is controlled by the air humidity (Eq. ):
Ci=f⋅Cs+(1-f)Γ,
where Cs is the atmospheric CO2 concentration in kgCO2 per kgair, and f is a coupling factor estimated via
f=f0*⋅1-DsDmax*+fmin⋅DsDmax*,
where Dmax* is the maximum specific humidity deficit of the air tolerated by vegetation in kgH2O per kgair, Ds is the leaf-to-air saturation deficit in kgH2O per kgair, and f0* is the value of the f factor for Ds=0gkg-1 and fmin given by
fmin=gcgc+gmes*,
where gc is the cuticular conductance in ms-1.
The CO2 assimilation is then limited by the photosynthetically active radiation via
An=(Am+Rd)⋅(1-e-ϵ*IaAm+Rd)-Rd,
where ϵ is the initial quantum use efficiency in kgCO2J-1PAR, Ia is the photosynthetically active radiation, and Rd is the dark respiration (Rd=Am9). ϵ is estimated with the following (Eq. ):
ϵ=ϵ0⋅Ci-ΓCi+2Γ,
where ϵ0 is the maximum quantum use efficiency in kgCO2J-1PAR.
Plant respiration
Depending on whether the model is run with forced LAI or with interactive vegetation, plant respiration is estimated in two different ways. First, when the model is run with forced LAI, plant respiration Rleaf is calculated as
Rleaf=Am9⋅LAI.
When the model is run with interactive vegetation, plant respiration is estimated, and an A term for the respiration of the above-ground stem biomass reservoir is added:
Rleaf=Am9⋅LAI+ηRBsQ10Ts-25Δt,
where Bs (expressed in kilograms of dry matter (DM) m−2) is the above-ground stem biomass reservoir, Ts is the leaf temperature in °C, ηR is a respiration rate fixed at 1 % d−1 (s-1), and Q10= 2.0.
TEB-GREENROOF input parameters
TEB-GREENROOF input parameters for roads, walls, and roofs.
In order to reduce the cost of calculation, the A-gs component of the ISBA–A-gs model has been rewritten and implemented in R, allowing for the rapid execution of simulations without the necessity of computing each step individually. Indeed, the separated A-gs model is forced with computed ISBA–A-gs input from a reference TEB-GREENROOF simulation with forced LAI. This approach enabled the simulation to run every time step simultaneously; however, this approach removed the retroactive effect of the CO2 fluxes. The LAI monthly evolution for the reference simulation was constructed using the monthly average of the spline interpolation of the punctually estimated values of LAI at the green roof case study site.
Description of the numerical setup implemented for running the A-gs photosynthesis model in a standalone configuration.
Pre-calibration
Before the calibration with the full TEB configuration, pre-calibration is done in order to reduce the plausible range of values and combinations of parameters for the calibration simulations. During the pre-calibration, the identified sensitive parameters (gmes*(25), Γ(25)ϵ0, Am,max(25), f0*, and F2max) are modified by conducting multiple simulations on the forced A-gs model, with the same parameter range and sampling method as those used in the sensitivity analysis but only on the sensitive parameters. A total of 50 000 simulations were conducted on the specified range of parameters. The intersection of the 0.1 % best simulations for the three scores, root-mean-square error (RMSE), mean absolute error (MAE), and coefficient of determination (r2), was retained and is displayed in Fig. . The results showed that for the 0.1 % intersection for the three scores, the average value of gmes*(25) for the simulations was 0.0019 ms-1. Thus, it was decided to test for the two values 0.001 and 0.002 ms-1. For f0*, the average was 0.56, but since the lower range was at 0.45, the values selected were 0.45 and 0.55. For Γ(25), the average was 43.1 ppmv, and the standard deviation was 8.86 ppmv, but, like for f0*, the values reached the upper range set at 49.5 ppmv, so the two values selected were set to 45 and 55 ppmv. For ϵ0, Am,max(25), and F2max, there was no clearly highlighted values. Consequently, the pair of values defined for ϵ0 and for Am,max(25) was chosen based on the ISBA–A-gs C3 and C4 values. For F2max, the values selected were 0.35 and 0.75.
Pre-calibration results. On the left is the distribution for RMSE, MAE, and r2.
Code and data availability
TEB is part of the software SURFEX from the CNRM open-source website https://opensource.umr-cnrm.fr under the CeCILL Free Software License Agreement v1.0 license. The version with net CO2 flux modelling for green roofs is available on 10.5281/zenodo.14289462. The experimental data that are used for the calibration and evaluation were provided by Stephan Weber from the Technische Universität Braunschweig; it is necessary to contact Stephan Weber directly.
Author contributions
AM did the model development, calibration, and validation and wrote the paper. CdM, AL, BB, VM, and AL supervised the project, gave their expertise on modelling, and reviewed the paper. SW provided experimental data, gave expertise on these data, and reviewed the paper.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
Acknowledgements
This research was funded by the Agence Nationale de la Recherche (ANR) under the ANR-22-CE92-0001-01 GREENVELOPES project and by the German Research Foundation (DFG) under project 505703010. The thesis work of Aurélien Mirebeau under the supervision of Cécile de Munck and Aude Lemonsu received a doctoral research grant from the Occitanie region (no. 00087377/21012846) and from Météo-France (no. 2130C0013).
Financial support
This research has been supported by the Agence Nationale de la Recherche (grant no. ANR-22-CE92-0001-01) and the Région Occitanie Pyrénées-Méditerranée (grant no. 00087377/21012846).
Review statement
This paper was edited by Cynthia Whaley and reviewed by two anonymous referees.
ReferencesAouade, G., Jarlan, L., Ezzahar, J., Er-Raki, S., Napoly, A., Benkaddour, A., Khabba, S., Boulet, G., Garrigues, S., Chehbouni, A., and Boone, A.: Evapotranspiration partition using the multiple energy balance version of the ISBA-A-gs land surface model over two irrigated crops in a semi-arid Mediterranean region (Marrakech, Morocco), Hydrol. Earth Syst. Sci., 24, 3789–3814, 10.5194/hess-24-3789-2020, 2020.Ascione, F., Bianco, N., De' Rossi, F., Turni, G., and Vanoli, G. P.: Green roofs in European climates. Are effective solutions for the energy savings in air-conditioning?, Appl. Energ., 104, 845–859, 10.1016/j.apenergy.2012.11.068, 2013. Bonan, G. B.: Ecological Climatology: Concepts and Applications, in: 3rd edn., Cambridge University Press, New York, ISBN 978-1-107-04377-0, 978-1-107-61905-0, 2016.Boone, A., Masson, V., Meyers, T., and Noilhan, J.: The Influence of the Inclusion of Soil Freezing on Simulations by a Soil–Vegetation–Atmosphere Transfer Scheme, J. Appl. Meteorol., 39, 1544–1569, 10.1175/1520-0450(2000)039<1544:TIOTIO>2.0.CO;2, 2000.Calvet, J.-C.: Investigating soil and atmospheric plant water stress using physiological and micrometeorological data, Agr. Forest Meteorol., 103, 229–247, 10.1016/S0168-1923(00)00130-1, 2000.Calvet, J.-C. and Soussana, J.-F.: Modelling CO2-enrichment effects using an interactive vegetation SVAT scheme, Agr. Forest Meteorol., 108, 129–152, 10.1016/S0168-1923(01)00235-0, 2001.Calvet, J.-C., Noilhan, J., Roujean, J.-L., Bessemoulin, P., Cabelguenne, M., Olioso, A., and Wigneron, J.-P.: An interactive vegetation SVAT model tested against data from six contrasting sites, Agr. Forest Meteorol., 92, 73–95, 10.1016/S0168-1923(98)00091-4, 1998.Calvet, J.-C., Rivalland, V., Picon-Cochard, C., and Guehl, J.-M.: Modelling forest transpiration and CO2 fluxes–response to soil moisture stress, Agr. Forest Meteorol., 124, 143–156, 10.1016/j.agrformet.2004.01.007, 2004.Carnell, R.: lhs: Latin Hypercube Samples, https://cran.r-project.org/web/packages/lhs/index.html (last access: 12 June 2024), 2022.Clapp, R. B. and Hornberger, G. M.: Empirical equations for some soil hydraulic properties, Water Resour. Res., 14, 601–604, 10.1029/WR014i004p00601, 1978.CNRM – Centre National de Recherches Météorologiques: CNRM Open Source, https://opensource.umr-cnrm.fr (last access: 6 June 2025), 2025.Cook-Patton, S. C. and Bauerle, T. L.: Potential benefits of plant diversity on vegetated roofs: A literature review, J. Environ. Manage., 106, 85–92, 10.1016/j.jenvman.2012.04.003, 2012.Currie, B. A. and Bass, B.: Estimates of air pollution mitigation with green plants and green roofs using the UFORE model, Urban Ecosyst., 11, 409–422, 10.1007/s11252-008-0054-y, 2008.Decharme, B., Boone, A., Delire, C., and Noilhan, J.: Local evaluation of the Interaction between Soil Biosphere Atmosphere soil multilayer diffusion scheme using four pedotransfer functions, J. Geophys. Res.-Atmos., 116, D20126, 10.1029/2011JD016002, 2011. Del Barrio, E. P.: Analysis of the green roofs cooling potential in buildings, Energ. Buildings, 27, 179–193, 1998.de Munck, C., Lemonsu, A., Masson, V., Le Bras, J., and Bonhomme, M.: Evaluating the impacts of greening scenarios on thermal comfort and energy and water consumptions for adapting Paris city to climate change, Urban Climate, 23, 260–286, 10.1016/j.uclim.2017.01.003, 2018.de Munck, C. S., Lemonsu, A., Bouzouidja, R., Masson, V., and Claverie, R.: The GREENROOF module (v7.3) for modelling green roof hydrological and energetic performances within TEB, Geosci. Model Dev., 6, 1941–1960, 10.5194/gmd-6-1941-2013, 2013.Erbs, D. G., Klein, S. A., and Duffie, J. A.: Estimation of the diffuse radiation fraction for hourly, daily and monthly-average global radiation, Sol. Energy, 28, 293–302, 10.1016/0038-092X(82)90302-4, 1982.Garisoain, R.: Évolution du cycle du carbone des tourbières pyrénéennes dans un contexte de changement climatique global: observation et modélisation, PhD thesis, Institut National Polytechnique de Toulouse – INPT, https://theses.hal.science/tel-05031426 (last access: 23 May 2025), 2023.Gibelin, A., Calvet, J., Roujean, J., Jarlan, L., and Los, S. O.: Ability of the land surface model ISBA-A-gs to simulate leaf area index at the global scale: Comparison with satellites products, J. Geophys. Res.-Atmos., 111, 2005JD006691, 10.1029/2005JD006691, 2006.Goudriaan, J., van Laar, H. H., van Keulen, H., and Louwerse, W.: Photosynthesis, CO2 and Plant Production, in: Wheat Growth and Modelling, edited by: Day, W. and Atkin, R. K., Springer US, Boston, MA, 107–122, 10.1007/978-1-4899-3665-3_10, 1985.Hamdi, R. and Masson, V.: Inclusion of a Drag Approach in the Town Energy Balance (TEB) Scheme: Offline 1D Evaluation in a Street Canyon, J. Appl. Meteorol. Clim., 47, 2627–2644, 10.1175/2008JAMC1865.1, 2008.Heusinger, J. and Weber, S.: Extensive green roof CO2 exchange and its seasonal variation quantified by eddy covariance measurements, Sci. Total Environ., 607–608, 623–632, 10.1016/j.scitotenv.2017.07.052, 2017a.Heusinger, J. and Weber, S.: Surface energy balance of an extensive green roof as quantified by full year eddy-covariance measurements, Sci. Total Environ., 577, 220–230, 10.1016/j.scitotenv.2016.10.168, 2017b.Jacobs, C. M. J.: Direct impact of atmospheric CO2 entrichment on regional transpiration, https://www.researchgate.net/publication/40182734_Direct_Impact_of_Atmospheric_CO2_Enrichment_on_Regional_Transpiration (last access: 17 June 2024), 1994.Kattge, J., Bönisch, G., Díaz, S., Lavorel, S., Prentice, I. C., Leadley, P., Tautenhahn, S., and Werner, G. D. A.: TRY plant trait database – enhanced coverage and open access, Glob. Change Biol., 26, 119–188, 10.1111/gcb.14904, 2020.Konopka, J., Heusinger, J., and Weber, S.: Extensive Urban Green Roof Shows Consistent Annual Net Uptake of Carbon as Documented by 5 Years of Eddy-Covariance Flux Measurements, J. Geophys. Res.-Biogeo., 126, e2020JG005879, 10.1029/2020JG005879, 2021. Köppen, W.: Versuch einer Klassifikation der Klimate, vorzugsweise nach ihren Beziehungen zur Pflanzenwelt, Geogr. Z., 6, 657–679, 1900.Kuronuma, T., Watanabe, H., Ishihara, T., Kou, D., Toushima, K., Ando, M., and Shindo, S.: CO2 Payoff of Extensive Green Roofs with Different Vegetation Species, Sustainability-Basel, 10, 2256, 10.3390/su10072256, 2018.Lazzarin, R. M., Castellotti, F., and Busato, F.: Experimental measurements and numerical modelling of a green roof, Energ. Buildings, 37, 1260–1267, 10.1016/j.enbuild.2005.02.001, 2005.Lemonsu, A., Masson, V., Shashua-Bar, L., Erell, E., and Pearlmutter, D.: Inclusion of vegetation in the Town Energy Balance model for modelling urban green areas, Geosci. Model Dev., 5, 1377–1393, 10.5194/gmd-5-1377-2012, 2012.Li, Y. and Babcock, R. W.: Green roofs against pollution and climate change. A review, Agron. Sustain. Dev., 34, 695–705, 10.1007/s13593-014-0230-9, 2014.Lorenz, J. M., Kronenberg, R., Bernhofer, C., and Niyogi, D.: Urban Rainfall Modification: Observational Climatology Over Berlin, Germany, J. Geophys. Res.-Atmos., 124, 731–746, 10.1029/2018JD028858, 2019.Masson, V.: A Physically-Based Scheme For The Urban Energy Budget In Atmospheric Models, Bound.-Lay. Meteorol., 94, 357–397, 10.1023/A:1002463829265, 2000.Masson, V. and Seity, Y.: Including Atmospheric Layers in Vegetation and Urban Offline Surface Schemes, J. Appl. Meteorol. Clim., 48, 1377–1397, 10.1175/2009JAMC1866.1, 2009.Masson, V., Le Moigne, P., Martin, E., Faroux, S., Alias, A., Alkama, R., Belamari, S., Barbu, A., Boone, A., Bouyssel, F., Brousseau, P., Brun, E., Calvet, J.-C., Carrer, D., Decharme, B., Delire, C., Donier, S., Essaouini, K., Gibelin, A.-L., Giordani, H., Habets, F., Jidane, M., Kerdraon, G., Kourzeneva, E., Lafaysse, M., Lafont, S., Lebeaupin Brossier, C., Lemonsu, A., Mahfouf, J.-F., Marguinaud, P., Mokhtari, M., Morin, S., Pigeon, G., Salgado, R., Seity, Y., Taillefer, F., Tanguy, G., Tulet, P., Vincendon, B., Vionnet, V., and Voldoire, A.: The SURFEXv7.2 land and ocean surface platform for coupled or offline simulation of earth surface variables and fluxes, Geosci. Model Dev., 6, 929–960, 10.5194/gmd-6-929-2013, 2013.Mirebeau, A.: Model changes : extensive green roof CO2 exchanges in the TEB urban canopy model, Zenodo [code], 10.5281/zenodo.14289462, 2024. Mironov, D., Heise, E., Kourzeneva, E., Ritter, B., Schneider, N., and Terzhevik, A.: Implementation of the lake parameterisation scheme flake into the numerical weather prediction model cosmo, Boreal Environ. Res., 15, 218–230, 2010.Noilhan, J. and Mahfouf, J.-F.: The ISBA land surface parameterisation scheme, Global Planet. Change, 13, 145–159, 10.1016/0921-8181(95)00043-7, 1996.Noilhan, J. and Planton, S.: A Simple Parameterization of Land Surface Processes for Meteorological Models, Mon. Weather Rev., 117, 536–549, 10.1175/1520-0493(1989)117<0536:ASPOLS>2.0.CO;2, 1989.Norman, J. M., Garcia, R., and Verma, S. B.: Soil surface CO2 fluxes and the carbon budget of a grassland, J. Geophys. Res.-Atmos., 97, 18845–18853, 10.1029/92JD01348, 1992.Ouldboukhitine, S.-E., Belarbi, R., Jaffal, I., and Trabelsi, A.: Assessment of green roof thermal behavior: A coupled heat and mass transfer model, Build. Environ., 46, 2624–2631, 10.1016/j.buildenv.2011.06.021, 2011.Redon, E., Lemonsu, A., and Masson, V.: An urban trees parameterization for modeling microclimatic variables and thermal comfort conditions at street level with the Town Energy Balance model (TEB-SURFEX v8.0), Geosci. Model Dev., 13, 385–399, 10.5194/gmd-13-385-2020, 2020.Redon, E. C., Lemonsu, A., Masson, V., Morille, B., and Musy, M.: Implementation of street trees within the solar radiative exchange parameterization of TEB in SURFEX v8.0, Geosci. Model Dev., 10, 385–411, 10.5194/gmd-10-385-2017, 2017.Sailor, D. J.: A green roof model for building energy simulation programs, Energ. Buildings, 40, 1466–1478, 10.1016/j.enbuild.2008.02.001, 2008.Saltelli, A., Annoni, P., Azzini, I., Campolongo, F., Ratto, M., and Tarantola, S.: Variance based sensitivity analysis of model output. Design and estimator for the total sensitivity index, Comput. Phys. Commun., 181, 259–270, 10.1016/j.cpc.2009.09.018, 2010.Seyedabadi, M. R., Eicker, U., and Karimi, S.: Plant selection for green roofs and their impact on carbon sequestration and the building carbon footprint, Environmental Challenges, 4, 100119, 10.1016/j.envc.2021.100119, 2021.Shafique, M., Xue, X., and Luo, X.: An overview of carbon sequestration of green roofs in urba, areas, Urban For. Urban Gree., 47, 126515, 10.1016/j.ufug.2019.126515, 2020. Sobol, I. M.: Sensitivity Estimates for Nonlinear Mathematical Model, Math. Modelling Comput. Exp., 1, 407–414, 1993.Sobol, I. M.: Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates, Math. Comput. Simulat., 55, 271–280, 10.1016/S0378-4754(00)00270-6, 2001.Stein, M.: Large Sample Properties of Simulations Using Latin Hypercube Sampling, Technometrics, 29, 143–151, 10.2307/1269769, 1987. Stocki, R.: A method to improve design reliability using optimal Latin hypercube sampling, Computer Assisted Mechanics and Engineering Sciences, 12, 393–411, 2005.Tan, T., Kong, F., Yin, H., Cook, L. M., Middel, A., and Yang, S.: Carbon dioxide reduction from green roofs: A comprehensive review of processes, factors, and quantitative methods, Renew. Sust. Energ. Rev., 182, 113412, 10.1016/j.rser.2023.113412, 2023.Virk, G., Jansz, A., Mavrogianni, A., Mylona, A., Stocker, J., and Davies, M.: Microclimatic effects of green and cool roofs in London and their impacts on energy use for a typical office building, Energ. Buildings, 88, 214–228, 10.1016/j.enbuild.2014.11.039, 2015.Wang, M., Yu, H., Liu, Y., Lin, J., Zhong, X., Tang, Y., Guo, H., and Jing, R.: Unlock city-scale energy saving and peak load shaving potential of green roofs by GIS-informed urban building energy modelling, Appl. Energ., 366, 123315, 10.1016/j.apenergy.2024.123315, 2024.Winter, K.: Ecophysiology of constitutive and facultative CAM photosynthesis, J. Exp. Bot., 70, 6495–6508, 10.1093/jxb/erz002, 2019.Zheng, X., Zou, Y., Lounsbury, A. W., Wang, C., and Wang, R.: Green roofs for stormwater runoff retention: A global quantitative synthesis of the performance, Resour. Conserv. Recy., 170, 105577, 10.1016/j.resconrec.2021.105577, 2021.Zhou, L. W., Wang, Q., Li, Y., Liu, M., and Wang, R. Z.: Green roof simulation with a seasonally variable leaf area index, Energ. Buildings, 174, 156–167, 10.1016/j.enbuild.2018.06.020, 2018.