Community Earth System Model (CESM) simulations disentangling CO2 effects - PI (4 of 4)
Data files
Jan 08, 2025 version files 194.05 GB
-
coupled_DEFmedslope_1xCO2_02.tar.gz
5.28 GB
-
coupled_DEFmedslope_1xCO2_03.tar.gz
45.84 GB
-
coupled_DEFmedslope_1xCO2_05.tar.gz
55.88 GB
-
coupled_DEFmedslope_1xCO2_SP_02.tar.gz
87.05 GB
-
README.md
18.18 KB
Abstract
It is widely recognized that water availability influences plant growth, but plants also exert a strong control on the water cycle. Plant transpiration constitutes about 60% of the total flux of water from the land to the atmosphere, and changes in vegetation coverage and ecosystem functioning can substantially alter land surface water and energy fluxes. A growing body of research has found that plant responses to increasing atmospheric CO2 can modify the water cycle. Rising CO2 concentrations at the land surface directly alter plant water use and land evapotranspiration, which previous studies have found can drive global-scale changes in precipitation, runoff, and streamflow. Concomitantly, the radiative effects of increasing CO2 modify atmospheric variables that affect plant functioning, such as temperature, vapor pressure deficit, and precipitation. Plant responses to these radiatively driven climate changes can further influence the hydrologic cycle. While multiple studies have found that plant responses to the physiological and radiative effects of CO2 can drive hydrologic impacts, there is limited consensus on the magnitude of plant responses’ overall influence. Significant uncertainty surrounds how plants respond to increasing CO2 and changing climate, as well as how a given plant change can impact the hydrologic cycle. This overall uncertainty is often discussed (e.g. IPCC 2021), but we still lack a systematic quantification of which processes are contributing the most to uncertainty in the hydrologic impact of plant responses to CO2. Ultimately, disentangling this uncertainty requires analysis of plants’ influence on the coupled biosphere-atmosphere system, because plant responses to increasing CO2 generate atmospheric feedbacks which can further modify land surface hydrology. In order to improve understanding of how plant responses to increasing atmospheric CO2 feed back on the global hydrologic cycle, we ran model experiments in the Community Earth System Model to systematically quantify the separate impacts of leaf-level and plant-level changes to different atmospheric drivers.
README: CESM simulations disentangling CO2 effects - PI
https://doi.org/10.5061/dryad.r4xgxd2nr
Related datasets
We ran a suite of model experiments to disentangle the physiological and radiative effects of CO2, and to identify how leaf area responses modulate those effects. The model output from this suite of simulations is divided among four repositories which isolate the physiological (PHYS), radiative (RAD), and full physiological+radiative (FULL) effects of CO2 compared to the preindustrial control (PI). This repository contains model output from the preindustrial control (PI) simulations.
Description of the data and file structure
This dataset is organized into four zipped files, each of which contains model output from a different simulation: coupled_DEFmedslope_1xCO2_02.tar.gz
coupled_DEFmedslope_1xCO2_03.tar.gz
, coupled_DEFmedslope_1xCO2_05.tar.gz
, and coupled_DEFmedslope_1xCO2_SP_02.tar.gz
. Each of these files contain model output from simulations where both the land and atmosphere experience preindustrial CO2 concentrations. coupled_DEFmedslope_1xCO2_02
is the 40 year preindustrial control spin up period, which was run with dynamic leaf area (i.e., BGC mode) - all other simulations (including the PHYS, RAD, and FULL simulations in related repositories) branched from the end of this simulation. coupled_DEFmedslope_1xCO2_03
and coupled_DEFmedslope_1xCO2_05
were run with dynamic leaf area (i.e., BGC mode), and are part of the same simulation. The coupled_DEFmedslope_1xCO2_03
case simulated years 41-94 (after branching from the 1xCO2 spin up at year 41), and the case coupled_DEFmedslope_1xCO2_05
continued this run in the same continuation for years 94-160. To analyze the full 41-160 year time series, output from coupled_DEFmedslope_1xCO2_03
and coupled_DEFmedslope_1xCO2_05
should be concatenated. The case coupled_DEFmedslope_1xCO2_SP_02
was run with fixed leaf area (i.e., SP mode), where the fixed leaf area was determined from the preindustrial leaf area climatology modeled in coupled_DEFmedslope_1xCO2_03
and coupled_DEFmedslope_1xCO2_05
. It is important to note that these preindustrial control simulations were run until year 160, while the simulations in the companion PHYS, RAD, and FULL datasets were run until year 200.
Within each of these zipped files, there are four folders: lnd
, atm
, rest
, and run_info
. lnd
and atm
contain model output in the form of postprocessed netcdf files from the land model and atmosphere model, respectively, where each file is a timeseries for a single variable. Each variable is described in the netcdf metadata. rest
contains the restart files from the end of the simulation, which can be used to continue the simulation for additional years, also in the form of netcdf files. Netcdf files can be opened and viewed using the xarray python package. run_info
documents how simulations were configured in more detail, and is likely only relevant to people already familiar with running the Community Earth System Model (CESM). The CESM user guide describes how to run CESM and how different inputs (e.g., namelist settings) are used when running the model. run_info
contains the shell script used to generate the simulation (e.g., coupled_BGC_DEFmedslope_1xCO2_spinup.sh
), the namelists
directory which contains files that document how each simulation was configured, and the SourceMods
directory which contains Fortran source files which document changes to the Community Earth System Model (CESM) code. All the files in run_info
can be opened in a text editor.
├── coupled_DEFmedslope_1xCO2_02.tar.gz
│ ├── lnd
│ │ └── tseries
│ │ └── month_1
│ │ ├── coupled_DEFmedslope_1xCO2_02.cam.h0.CLOUD.000101-004012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_02.cam.h0.FLDS.000101-004012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_02.cam.h0.FLNS.000101-004012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_02.cam.h0.FLNT.000101-004012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_02.cam.h0.FLUT.000101-004012.nc
│ │ └── ... (47 files total)
│ ├── rest
│ │ ├── 0004-01-01-00000
│ │ │ ├── coupled_DEFmedslope_1xCO2_02.cam.h0.0003-12.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_02.cam.i.0004-01-01-00000.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_02.cam.r.0004-01-01-00000.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_02.cam.rs.0004-01-01-00000.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_02.cice.r.0004-01-01-00000.nc
│ │ │ └── ... (16 files total)
│ │ └── 0041-01-01-00000
│ │ ├── coupled_DEFmedslope_1xCO2_02.cam.h0.0040-12.nc
│ │ ├── coupled_DEFmedslope_1xCO2_02.cam.i.0041-01-01-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_02.cam.r.0041-01-01-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_02.cam.rs.0041-01-01-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_02.cice.r.0041-01-01-00000.nc
│ │ └── ... (16 files total)
│ └── run_info
│ ├── coupled_BGC_DEFmedslope_1xCO2_spinup.sh
│ ├── namelists
│ │ ├── user_nl_cam
│ │ ├── user_nl_clm
│ │ └── user_nl_cpl
│ └── SourceMods
│ └── src.clm
│ ├── BareGroundFluxesMod.F90
│ ├── CanopyFluxesMod.F90
│ ├── PhotosynthesisMod.F90
│ ├── README
│ ├── UrbanFluxesMod.F90
│ └── WaterFluxType.F90
│
├── coupled_DEFmedslope_1xCO2_03.tar.gz
│ ├── atm
│ │ └── proc
│ │ └── tseries
│ │ ├── day_1
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h1.PRECT.00410101-00510101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h1.PRECT.00510102-00610101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h1.PRECT.00610102-00710101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h1.PRECT.00710102-00810101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h1.PRECT.00810102-00910101.nc
│ │ │ └── ... (18 files total)
│ │ ├── month_1
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h0.CLDHGH.004101-009012.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h0.CLDHGH.009101-009412.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h0.CLDLOW.004101-009012.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h0.CLDLOW.009101-009412.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h0.CLDMED.004101-009012.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h0.CLDMED.009101-009412.nc
│ │ │ └── ... (122 files total)
│ │ └── move_files.csh
│ ├── lnd
│ │ └── proc
│ │ └── tseries
│ │ ├── day_1
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.clm2.h2.BTRAN2.00410101-00910101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.clm2.h2.BTRAN2.00910102-00941231.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.clm2.h2.EFLX_LH_TOT.00410101-00910101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.clm2.h2.EFLX_LH_TOT.00910102-00941231.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_03.clm2.h2.FIRE.00410101-00910101.nc
│ │ │ └── ... (72 files total)
│ │ └── month_1
│ │ ├── coupled_DEFmedslope_1xCO2_03.clm2.h0.ACTUAL_IMMOB.004101-009012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_03.clm2.h0.ACTUAL_IMMOB.009101-009412.nc
│ │ ├── coupled_DEFmedslope_1xCO2_03.clm2.h0.AGNPP.004101-009012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_03.clm2.h0.AGNPP.009101-009412.nc
│ │ ├── coupled_DEFmedslope_1xCO2_03.clm2.h0.ALT.004101-009012.nc
│ │ └── ... (948 files total)
│ ├── rest
│ │ └── 0095-01-01-00000
│ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h0.0094-12.nc
│ │ ├── coupled_DEFmedslope_1xCO2_03.cam.h1.0094-01-02-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_03.cam.i.0095-01-01-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_03.cam.r.0095-01-01-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_03.cam.rs.0095-01-01-00000.nc
│ │ └── ... (17 files total)
│ └── run_info
│ ├── coupled_BGC_DEFmedslope_1xCO2.sh
│ ├── namelists
│ │ ├── user_nl_cam
│ │ ├── user_nl_clm
│ │ └── user_nl_cpl
│ └── SourceMods
│ └── src.clm
│ ├── BareGroundFluxesMod.F90
│ ├── CanopyFluxesMod.F90
│ ├── PhotosynthesisMod.F90
│ ├── README
│ ├── UrbanFluxesMod.F90
│ └── WaterFluxType.F90
│
├── coupled_DEFmedslope_1xCO2_05.tar.gz
│ ├── atm
│ │ └── proc
│ │ └── tseries
│ │ ├── day_1
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h1.PRECT.00950101-01050101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h1.PRECT.01050102-01150101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h1.PRECT.01150102-01250101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h1.PRECT.01250102-01350101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h1.PRECT.01350102-01450101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h1.PRECT.01450102-01550101.nc
│ │ │ └── ... (18 files total)
│ │ ├── month_1
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h0.CLDHGH.009501-014412.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h0.CLDHGH.014501-016012.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h0.CLDLOW.009501-014412.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h0.CLDLOW.014501-016012.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h0.CLDMED.009501-014412.nc
│ │ │ └── ... (122 files total)
│ │ └── move_files.csh
│ ├── lnd
│ │ └── proc
│ │ └── tseries
│ │ ├── day_1
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.clm2.h2.BTRAN2.00950101-01450101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.clm2.h2.BTRAN2.01450102-01601231.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.clm2.h2.EFLX_LH_TOT.00950101-01450101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.clm2.h2.EFLX_LH_TOT.01450102-01601231.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_05.clm2.h2.FIRE.00950101-01450101.nc
│ │ │ └── ... (72 files total)
│ │ └── month_1
│ │ ├── coupled_DEFmedslope_1xCO2_05.clm2.h0.ACTUAL_IMMOB.009501-014412.nc
│ │ ├── coupled_DEFmedslope_1xCO2_05.clm2.h0.ACTUAL_IMMOB.014501-016012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_05.clm2.h0.AGNPP.009501-014412.nc
│ │ ├── coupled_DEFmedslope_1xCO2_05.clm2.h0.AGNPP.014501-016012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_05.clm2.h0.ALT.009501-014412.nc
│ │ └── ... (948 files total)
│ ├── rest
│ │ └── 0201-01-01-00000
│ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h0.0200-12.nc
│ │ ├── coupled_DEFmedslope_1xCO2_05.cam.h1.0200-01-02-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_05.cam.i.0201-01-01-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_05.cam.r.0201-01-01-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_05.cam.rs.0201-01-01-00000.nc
│ │ └── ... (17 files total)
│ └── run_info
│ ├── coupled_BGC_DEFmedslope_1xCO2_continue.sh
│ ├── namelists
│ │ ├── user_nl_cam
│ │ ├── user_nl_clm
│ │ └── user_nl_cpl
│ └── SourceMods
│ └── src.clm
│ ├── BareGroundFluxesMod.F90
│ ├── CanopyFluxesMod.F90
│ ├── PhotosynthesisMod.F90
│ ├── README
│ ├── UrbanFluxesMod.F90
│ └── WaterFluxType.F90
│
├── coupled_DEFmedslope_1xCO2_SP_02.tar.gz
│ ├── atm
│ │ └── proc
│ │ └── tseries
│ │ ├── day_1
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h1.PRECT.00410101-00510101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h1.PRECT.00510102-00610101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h1.PRECT.00610102-00710101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h1.PRECT.00710102-00810101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h1.PRECT.00810102-00910101.nc
│ │ │ └── ... (18 files total)
│ │ └── month_1
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h0.CLDHGH.004101-009012.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h0.CLDHGH.009101-014012.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h0.CLDHGH.014101-016012.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h0.CLDLOW.004101-009012.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h0.CLDLOW.009101-014012.nc
│ │ │ └── ... (122 files total)
│ │ └── move_files.csh
│ ├── lnd
│ │ └── proc
│ │ └── tseries
│ │ ├── day_1
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.clm2.h2.ANSHA_LN.00410101-00910101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.clm2.h2.ANSHA_LN.00910102-01410101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.clm2.h2.ANSHA_LN.01410102-01601231.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.clm2.h2.ANSUN_LN.00410101-00910101.nc
│ │ │ ├── coupled_DEFmedslope_1xCO2_SP_02.clm2.h2.ANSUN_LN.00910102-01410101.nc
│ │ │ └── ... (72 files total)
│ │ └── month_1
│ │ ├── coupled_DEFmedslope_1xCO2_SP_02.clm2.h0.ANSHA_LN.004101-009012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_SP_02.clm2.h0.ANSHA_LN.009101-014012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_SP_02.clm2.h0.ANSHA_LN.014101-016012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_SP_02.clm2.h0.ANSUN_LN.004101-009012.nc
│ │ ├── coupled_DEFmedslope_1xCO2_SP_02.clm2.h0.ANSUN_LN.009101-014012.nc
│ │ └── ... (948 files total)
│ ├── rest
│ │ └── 0161-01-01-00000
│ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h0.0160-12.nc
│ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.h1.0160-01-02-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.i.0161-01-01-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.r.0161-01-01-00000.nc
│ │ ├── coupled_DEFmedslope_1xCO2_SP_02.cam.rs.0161-01-01-00000.nc
│ │ └── ... (17 files total)
│ └── run_info
│ ├── coupled_SP_DEFmedslope_1xCO2_v2.sh
│ ├── namelists
│ │ ├── user_nl_cam
│ │ ├── user_nl_clm
│ │ └── user_nl_cpl
│ └── SourceMods
│ └── src.clm
│ ├── BareGroundFluxesMod.F90
│ ├── CanopyFluxesMod.F90
│ ├── PhotosynthesisMod.F90
│ ├── README
│ ├── UrbanFluxesMod.F90
│ └── WaterFluxType.F90
Methods
We ran simulations in a partially coupled configuration of the Community Earth System Model version 2 (CESM2; git tag: cesm2_3_beta03) using an active atmosphere (CAM6; git tag cam6_3_017), an active land model (CLM5; git tag branch_tags/PPE.n08_ctsm5.1.dev030), and a slab ocean which is based on q-fluxes from preindustrial simulations of the full dynamic ocean model. The land was initialized with the spun-up land state from the default model parameterization which includes the carbon content of soil and vegetation pools.
We ran a suite of model experiments to disentangle the physiological and radiative effects of CO2, and to identify how leaf area responses modulate those effects. The model output from this suite of simulations is divided among four repositories which isolate the physiological (PHYS), radiative (RAD), and full physiological+radiative (FULL) effects of CO2 compared to the preindustrial control (PI). All simulations branched from a preindustrial control simulation at year 41. Simulations in the PHYS, RAD, and FULL repositories were run for 160 additional years (i.e., ending at year 200), including spin up. This repository contains model output from preindustrial control simulations, which were run for 120 additional years (i.e., ending at year 160), and the preindustrial controls were already largely spun up at year 41.. We ran two simulations: one where the leaf area was held fixed at preindustrial levels, and one where the leaf area was dynamically calculated by the model. Both simulations were equilibrium concentration-driven experiments where both the land and the atmosphere experienced preindustrial CO2 concentrations.