scieee AI-readable full text Open interactive document viewer

Peatland-VU-NUCOM (PVN 1.0): A peatland GHG emissions model using dynamic plant functional types

Tanya J. R. Lippmann

Abstract

This is the source code of the PVN model. For more information on how to use the model, please see the documentation included within the bitbucket repository (https://bitbucket.org/tlippmann/pvn_public/src/master/) or contact the authors. The model has been tested in a mac osx environment and may not work on other operating systems. Run 'make' in the terminal window to compile the code. To execute the model: './PVN'.

Full text

Peatland-VU model output files Some output files are always produced. If you are interested in other variables, you can change and increase the number of output files using the variable, ‘ProfileOutput’ in the params file using the syntax as follows: ProfileOutput = [ 4 5 7 18]; A positive flux is indicative of a gas moving from the ecosystem into the atmosphere. The SOM pools in the model are always in this order: 1: peat 2: dissolved organic matter manure 3: manure solids 4: root exudates 5: litter and roots 6: microbial biomass 7: humus; • The PFT output files are numbered (e.g. Methanefluxes_PFT1.dat) where the number of the PFT refers to the same order that the PFTs are written in the ‘PFTFiles’ variable in the ‘params’ file. The output of the PFT files follow the same structure as the files described below. i.e. Methanefluxes_PFT1.dat follows the same structure as methanefluxes.dat and Biomass_PFT1.dat follows the same structure at biomass.dat. Output files that are not per PFT represent the total ecosystem (or the sum of all PFTs.) Always output by the model: The following output files are produced by the model by default. The details are described below. • methanefluxes.dat • CO2reservoirs.dat • CO2layers.dat • biomass.dat • carbonbalance.dat • peat_masschange.dat • params The params file used to run the model simulation is copied to the output directory. • methanefluxes.dat and Methanefluxes_PFT*.dat The methane gas exchange from the soil surface to the atmosphere. The internal model variable, ‘TotalMethane’ is written to this output file. Units: mg CH4 m-2 hr-1 Columns: 1. Day number 2. Ebullition flux 3. Plant transported methane flux 4. Diffusion flux 5. N/A 6. Net CH4 flux • CO2reservoirs.dat and CO2reservoirs_PFT*.dat The internal model variable, ‘ReservoirTime’ is written to this output file. The CO2 fluxes from the following belowground sources: Units: mg CO2 m-2 hr-1 Columns: 1. Day number 2. Peat SOM 3. Dissolved manure SOM 4. Solid manure SOM 5. Root exudates SOM 6. Litter and roots SOM 7. Microbial biomass SOM 8. Humus SOM 9. CO2 produced by CH4 oxidation 10. CO2 from decomposition of surface litter 11. Total of the above columns • CO2layers.dat and CO2layers_PFT*.dat The internal model variable, ‘LayerTime’ is written to this output file. The file for each PFT is then the CO2 fluxes produced from the SOM produced by each PFT. The CO2 fluxes produced from each soil layer: Units: mg CO2 m-2 hr-1 Columns: 1. Day number 2. CO2 produced from all SOM pools in the topmost soil layer + CO2 produced by CH4 oxidation (from all soil layers) + CO2 from decomposition of surface litter 3. CO2 produced from all SOM pools in the 2nd topmost soil layer … X. CO2 produced from all SOM pools in the deepest soil layer • biomass.dat and Biomass_PFT*.dat The internal model variable, ‘BioMassRec’ is written to this output file. This file includes the output for a range of variables relating to aboveground plant activity: Units: see below Columns: 1. Day number 2. Mass of aboveground biomass [kg C m-2] 3. Net primary production (positive flux from atmosphere to land) [mg CO2 m-2 hr-1] 4. Plant respiration [mg CO2 m-2 hr-1] 5. Net Ecosystem Exchange (GPP – Reco) [mg CO2 m-2 hr-1]. Positive flux from land to atmosphere. 6. Reco (soil and plant respiration) [mg CO2 m-2 hr-1]. (Positive flux from land to atmosphere). 7. Litter mass [kg C m-2] • CBioMassPFT.dat The carbon biomass of each PFT is written to this output file. Units: kg C m-2 Columns: 1. Day number 2. carbon biomass of PFT 1 3. carbon biomass of PFT X… Final column: Total carbon biomass of all PFTs. • CoverFraction_PFT.dat The cover fraction of each PFT is written to this output file. Units: - Columns: 1. Day number 2. Cover fraction of PFT 1 3. Cover fraction of PFT X… Final column: Total cover fraction of all PFTs. • HeightGrowth_PFT.dat The height of each PFT is written to this output file. Units: m Columns: 1. Day number 2. Height of PFT 1 3. Height of PFT X… Final column: Summed height of all PFTs. • LAI.dat The LAI of each PFT is written to this output file. Units: m2 m-2 Columns: 1. Day number 2. LAI of PFT 1 3. LAI of PFT X… Final column: Summed LAI of all PFTs. • LightIntercepted.dat The fraction of light intercepted by each PFT is written to this output file. Units: - Columns: 1. Day number 2. fraction of light intercepted by PFT 1 3. fraction of light intercepted by PFT X… Final column: The sum of all light intercepted. • LitterLayer_PFT.dat The litter pool of each PFT is written to this output file. Units: kg C m-2 Columns: 1. Day number 2. Litter mass of PFT 1 3. Litter mass of PFT X… Final column: Summed litter of all PFTs. • MossGrowth_PFT.dat The height (or depth) of the living moss layer of each PFT is written to this output file. This value should be zero for vascular PFTs. Units: m Columns: 1. Day number 2. Depth of living moss layer of PFT 1 3. Depth of living moss layer of PFT X… Final column: Summed depth of all moss PFTs. • NPP_PFT.dat The net primary production (NPP) of each PFT is written to this output file. Notice the units differs from the other flux output. Units: kg CO2 m-2 day-1 Columns: 4. Day number 5. NPP of PFT 1 6. NPP of PFT X… Final column: Summed NPP of all PFTs. • PlantRespiration_PFT.dat The plant respiration of each PFT is written to this output file. Notice the units differs from the other flux output. Units: kg CO2 m-2 day-1 Columns: 7. Day number 8. plant respiration of PFT 1 9. plant respiration of PFT X… Final column: Summed plant respiration of all PFTs. Optional output files Here we provide the details of the optional output files. The number of output files can be defined in the params file using the parameter, ‘ProfileOutput’ using the following syntax: ProfileOutput = [ 4 5 7 18]; ProfileOutput = [1]; File: temperature.dat Daily soil temperature of each soil layer. First column is topmost soil layer. Final column is lowermost soil layer. Units: deg C ProfileOutput = [2]; File: moisture.dat Soil moisture (θ, water content per layer volume) Units: volume volume-1 ProfileOutput = [3]; File: methane.dat and Methane_PFT*.dat CH4 concentration in pore space per layer volume Units: μMol m-3 ProfileOutput = [4]; File: roots.dat and Roots_PFT*.dat Root mass per soil layer, starting at the top Units: kg C m-2 ProfileOutput = [5]; File: labileSOM.dat and labileSOM_PFT*.dat Sum of the SOM pools that constitute labile carbon, per soil layer, starting at the top. These include: liquid manure, solid manure, root exudates, litter and roots, microbial biomass Units: kg C m-2 ProfileOutput = [6]; File: ice.dat Mass of ice, per layer, starting at the top Units: kg m-2 ProfileOutput = [7]; File: watertable.dat Daily water table data. Metres below surface. Positive number if below surface. Only produced if water table is calculated internally by the model. Units: metres ProfileOutput = [8]; File: npp.dat Daily Net Primary Production (NPP) Units: kg C m-2 day-1 ProfileOutput = [9]; File: totalSOM.dat and TotalSOM_PFT*.dat The sum of all SOM pools, per soil layer, starting at the top. Units: kg C m-2 ProfileOutput = [11]; Deleted: filled pore space, Deleted: uMol Formatted: Superscript File: peat.dat and Peat_PFT*.dat Mass of the peat SOM pool, per soil layer, starting at the top Units: kg C m-2 ProfileOutput = [12]; File: liquid_manure.dat Mass of the liquid manure SOM pool, per soil layer, starting at the top Units: kg C m-2 ProfileOutput = [13]; File: solid_manure.dat Mass of the solid manure SOM pool, per soil layer, starting at the top Units: kg C m-2 ProfileOutput = [14]; File: rhizodep.dat and Rhizodep_PFT*.dat Mass of the root exudates SOM pool, per soil layer, starting at the top Units: kg C m-2 ProfileOutput = [15]; File: litter_roots.dat and Litter_roots_PFT*.dat Mass of the litter and root SOM pool, per soil layer, starting at the top Units: kg C m-2 ProfileOutput = [16]; File: microbes.dat and Microbes_PFT*.dat Mass of the microbial biomass SOM pool, per soil layer, starting at the top Units: kg C m-2 ProfileOutput = [17]; File: humus.dat and Humus_PFT*.dat Mass of the humus SOM pool, per soil layer, starting at the top Or mass of the PFTX humus SOM pool, per soil layer, starting at the top Units: kg C m-2