CHANGELOG for Biome-BGCMuSo v7.0 2023-01-10 by Dóra HIDY -------------------------------- Changes from 7.0beta8 to 7.0beta9 -------------------------------- - DEVELOPMENT: output variables for litter CN ratio [2676-2679] - DEVELOPMENT: new output variables for summarizing N-mineralization and immobilization Gross mineralization flux [2213-2222] Actual immobilization flux [2223-2232] Potential immobilization flux [2233-2242] Summarized mineralization N flux from litter decomposition [2435-2444] Summarized mineralization N flux from soil decomposition [2445-2454] Summarized immobilization N flux from litter decomposition [2455-2464] Summarized immobilization N flux from soil decomposition [2465-2474] Summarized N input from environment [2475-2484] Total and cumulated variations of above mentioned output variables [2485-2497] -------------------------------- Changes from 7.0beta7 to 7.0beta8 -------------------------------- - BUG FXING: correction of pond water calculation in infiltANDpond.c - DEVELOPMENT: extension of log file with model estimated soil properties layer by layer (VWCsat, VWCfc, VWCwp, VWChw, PSIsat, PSIfc, PSIhw, soilB, bulkDensity, drainge coefficient, capillary fringe) - BUG FIXING: Previously dead stem and coarse woody biomass combustion proportion (DSCP and CWCP) were constant parameter in source code. From MuSo 7.0-b5 DSCP and CWCP can e set in EPC file, but model used the contstant values in source code. Now this bug has solved. - DEVELOPMENT: new output variables for summarizing soil and litter decomposition litrc_to_soilc_total [704] cwdc_to_litrc_total [710] itrn_to_soiln_total [2041] cwdn_to_litrn_total [2321] actIMMOBflux_total [2339] netMINERflux_total [2357] potIMMOBflux_total [2375] grossMINERflux_total [2391] litr1c_to_soil1c_total [968] -------------------------------- Changes from 7.0beta6 to 7.0beta7 -------------------------------- - DEVELOPMENT: extension of log file with aboveground litter and CWD-content - BUG FXING: correction of soil stress calculation function (error in ramp function calculation) - DEVELOPMENT: distinction of groundwater and capillary zone GWeff, CFeff variables shows groundwater zone and capillary zone - DEVELOPMENT: calculation of groundwater effect -------------------------------- Changes from 7.0beta5 to 7.0beta6 -------------------------------- - BUG FIXING: restart data contains also soil cover (aboveground litter and CWD) - DEVELOPMENT: rename of mulch.c into soilCover.c - DEVELOPMENT: new hydrological parameters relative soil water content (SAT-FCbase) [2614-2623] relative soil water content (FC-WP): [2624-2633] - DEVELOPMENT: implementation of control to remain balance in diffus_calc.c - limit of diffusion is the equalization of the difference in relative soil moisture - DEVELOPMENT: new process - flooding (similar to groundwater): User can set the first and the last days of flooding period (flooding_normal.txt, flooding_spinup.txt, flooding_trasient.txt), in which the water comes from above saturates the top soil layer and the pond water (if flooding height is greater than 0). new routine: flooding.c, flooding_init.c - DEVELOPMENT: new management type - CWD-extract: User can set the ratio of extract of aboveground coarse woody debris) new routines: CWDextract_init.c, CWDextract.c - DEVELOPMENT: new management type - mulching: User can set the amount and type of mulch which turns into the top litter or CWD layer) new routines: mulching_init.c, mulching.c - BUG FIXING: in case of groundwater, not only soil evaporation can be limited by evaporation limit, but also GWdischarge (pondANDrunoffD.c) -------------------------------- Changes from 7.0beta4 to 7.0beta5 -------------------------------- - DEVELOPMENT and BUG FIXING: separation of evapotranspiration limitation 1. potEVPandSUBLsurface limits soil water and pond water evaporation, plus snow water sublimation 2. potETcanopy limits canopy evaporation and transpiration 3. PET = potEVPandSUBLsurface + potETcanopy (ET must less than equal to PET) 4. ET_Elimit defines the maximum available energy - limits PET (ET_Elimit >= PET) change of subroutine's name: potEVPsurface_to_actEVPsurface.c (instead of potEVP_to_actEVP.c) cease of a subroutine: PriestleyTaylor.c new subroutine (merge of PriestleyTaylor.c and potEVP_calc): Elimit_and_PET.c new output variables: ET_Elimit[224] (instead of soilwTRPdemand_SUM - inactive variable) potEVPandSUBLsurface[212] - Potential evaporation (including sublimation) of surfac potETcanopy[243] - Potential evaporation and transpiration of canopy potEVPsurface[244] - Potential evaporation of surface - DEVELOPMENT: leap control (in leap years the last day is ignored istead of 29th February) in input data - DEVELOPMENT: new EPC parameters dead stem biomass combustion proportion - line 27 coarse woody biomass combustion proportion - line 28 - DEVELOPMENT: model error message in case of negative radiation data - BUG FIXING: aboveground litter and coarse woody debris into restart data structure (to avoid disuption between spinup and transient phase) - BUG FIXING: in canopy water calculation -------------------------------- Changes from 7.0beta3 to 7.0beta4 -------------------------------- - DEVELOPMENT: possibility of using varying fire mortality file (FM_spinup.txt, FM_transient.txt, FM_normal.txt) modification: the name of the varying whole plant mortality: WPM_spinup.txt, WPM_transient.txt, WPM_normal.txt - BUG FUXING: in drainage coefficient of tipping calculation -------------------------------- Changes from 7.0beta2 to 7.0beta3 -------------------------------- - DEVELOPMENT: new output variable for surface evaporation (sum of pond water and soil water evaporation) - surfaceEVP [239] - DEVELOPMENT: error message for 1-year simulation with truncated meteorological data in met_init.c minium simulation time: 1 year -------------------------------- Changes from 7.0beta1 to 7.0beta2 -------------------------------- - DEVELOPMENT: activating of Priestley-Taylor method for evapotranspiration calculation extra EPC-parameters (line 22): threshold temperature for PET-calculation using PT-method (-9999: no treshold data - only PM-method) - DEVELOPMENT: flag for the new optional intercepcion calculation in INI file line 60: (flag) interception calculation method (0: original - based on allLAI, 1: new - based on projLAI) - DEVELOPMENT: implementing of direct litter decomposition method extra SOI parameters (57-59): fraction of direct decomposition of labile/cellulose/lignin litter pool new output variables: litrc_to_release_total [1100], litrn_to_release_total[1770] - BUG FIXING: in drainage coefficient calculation in case of model estimation (tipping.c) - BUG FIXING: in CFD-calculation in groundwaterT_preproc.c - BUG FIXING: top soil layer must to be saturated if there is pond water (pondANDrunoffD.c, potEVP_to_actEVP.c) - BUG FIXING: canopy evaporation must be taken into account in the limitation calculation of potential evaporation (canopy_et.c, pondANDrunoffD.c) -------------------------------- Changes from 6.2 to 7.0beta1 -------------------------------- - BUG FIXING: in fertilization.c (water content of fertilizer goes into pond water) - DEVELOPMENT: implementation of groundwater calculation into Richards method - BUG FIXING: in flowering_heatstress.c (calculation of flowering heat stress) - DEVELOPMENT: after flowering heat stress, the carbon assimilation of the grain pool is limited (depending on the Closs due to flowering heat stress) changes in daily_allocation.c - MODIFICATION: resetting fruitC_HRV variable on the first day of the year (instead of the day of planting) - DEVELOPMENT: taking into account canopy water evaporation in calculation of potential evaporation (PET=evapPOT+transpPOT+cwe) - BUG FIXATION: potential transpiration can limit actual sum of transpiration (from version MuSo6.3) - DEVELOPMENT: in conditional irrigation, the number of the affected layer can be set by the User conditional irrigation is activated only if LAI>0.2 m2/m2 modified section of INI file: CONDITIONAL_MANAGEMENT_STRATEGIES 0 (flag) conditional mowing ? 0 - no, 1 - yes 0.0 (m2/m2) fixed value of the LAI before MOWING 0.0 (m2/m2) fixed value of the LAI after MOWING 0.0 (%) transported part of plant material after MOWING 0 (flag) conditional irrigation? 0 - no, 1 - based on SWCratio, 2 - based on SMSI) 4 (n) number of soil layers taking into account in cond. IRRIGATION (-9999: RZ) 0.6 (prop) starting point of cond. IRRIGATION (condIRGflag=1:SWCratio, condIRGflag=2:SMSI) 0.9 (prop) SWCratio of rootzone after cond. IRRIGATION 50 (kgH2O/m2) maximum amount of irrigated water - DEVELOPMENT: in senescence simulation three senescence mortality parameters are available in EPC file previous versions - aboveground, belowground, non-structured new versions - maximum of leaf senescence mortality coefficient, maximum of softstem/froot senescence mortality coefficient, multiplier for non-structured - BUG FIXATION: percolation of the bottom layer is zero if GW is available in soil column - DEVELOPMENT: testing option (if aerodyn. resistance== DATA_GAP -> no soil evaporation) - DEVELOPMENT: groundwater simulation capillary fringe can be estimated from PSIsat values modification of 2. groundwater calculation method percolation from capillary zone is limited - MODIFICATION: senescence is only occured (SMSI > 0): if leafC is greater than 0 - DEVELOPMENT: control to avoid oversaturated SWC pool in richards.c - DEVELOPMENT: adjustment of allocation parameters due to small error (<10-4) in the setting of allocation parameters (warning into logfile) - DEVELOPMENT: control of the groundwater data reading (to avoid memory allocation error) - DEVELOPMENT: new output variable (root index; 3127) - the ratio of live&dead root biomass to total live&dead biomass - DEVELOPMENT (only for expert users): in case of conditional managment a starting year can be set by the user - MODIFICATION: root index [3127] is calculated on the day of the harvest (similar to harvest index) new output variable: belowground ratio [3128] is the ratio of the root (woody: fine+coarse, non-woody: only fine root) and the total vegetation biomass on the given day - MODIFICATION: unit of leaching variables are kgN/m2 uniformly - BUG FIXATION: small corrections in Richards calculation - DEVELOPMENT: new output variables - layer-defined cumulated leached sminNH4 [2131-2140] and sminNO3 [2151-2160] - DEVELOPMENT: new output variables - cumulated C and N allocation [2587 and 2588] - MODIFICATION: cumulative SWC and N stress is resetted on the first day of year (instead of at the end of the vegetation period). Thanks to this it can be used as annual variable - BUG FIXING: in calculation of soil water stress using transpDemand method (soilstress_calculation.c) in layers where no root is available no root -> no transpiration -> no stress (instead of no root -> no transpiration -> full stress) - BUG FIXING: in calculation of available water (multilayer_transpiration - DEVELOPMENT: introduction of different groundwater calculation method for the soil hydrology calculations (Tipping and Richards) 3 new subroutines: groundwaterRICHARDS.c, groundwaterRICHARDSpostproc.c, groundwaterTIPPING.c in case of Tipping - "increasing FC values" groundwater method in case of Richards - "CF calculations" groundwater method for Richards: 12 soil layers in case of GW - normal layers, capillary zone and groundwater zone - DEVELOPMENT: thickness of capillary fringe is not an input data, but it is estimated by the model using PSIsat values of the different soil layers - DEVELOPMENT: calculation (and printing out) of net water flux (signed sum of percolation and diffusion) instead of percolation and diffusion separately (tipping.c, richards.c, output_map_init.c) - DEVELOPMENT: calculation (and printing out) of leaching (sminNH4, sminNO3, soilDOC and soilDON) based on net water flux instead of percolation and diffusion fluxes separately (multilayer_leaching.c, make_zero_flux_struct.c, output_map_init.c) - MODIFICATION: instead of soilw_from_GW - GWdischarge (groundwater_snk) - BUG FIXING: SWCstressLENGTH variable is not reinitialized on the first simulation year - MODIFICATION: cumulative water and N stress variable is resetted only on the first simulation year (in order to use as annual variable) - BUG FIXING: correction of N2O_Ceq calculation in summary.c (NGB calculation is affected also) - DEVELOPMENT: new output variable - cum_CH4 [3129] - BUG FIXING: correction of the logic of heat stress during anthesis - it causes mortality in flowering phenophase (not in the next one) - MODIFICATION: onday and offday output variables can be used as annual data (on last simulation day it contain the data referring to the actual year) - BUG FIXING: in initalization of soilwDiffus variable in case of SHCM_flag = 2 (tipping without diffusion - test version) - BUG FIXING: initialization of cumulative leaching variables on the first day of every simulation year - DEVELOPMENT: new routine (annVARinit.c) for initialization of annual/cumulative variables - DEVELOPMENT: new output variables for examinig cumulative N-immobilization and N-mineralization from different sources [2434-2443] - MODIFICATION: removal of the option of using the old version of the SOILf file (with extraSOIparameters.txt) - DEVELOPMENT: new SOI parameters: CN ratio of soil1, soi2, soil3 pool, estimation of total SOC-content to end spinup run - DEVELOPMENT: new input control in epc_init.c (SGS must be less than 365) - BUG FIXING: resetting of the onday/offday variables only at the end of the routine (phenphase.c) - BUG FIXING: in tipping.c - resetting soilw_diffus variables for all soil layers - DEVELOPMENT: new output variable - cum_nbp [3130] - BUG FIXING: correction of onday/offday output (annual variable) - BUG FIXING: correction remday_curgrowth and remday_transfer parameters - DEVELOPMENT: new SOIL parameters - decomposition limitation parameters (FP-WP ratio) to calculate rate constant scalar for soil water content - DEVELOPMENT: modified decomposition limitation parameters (FP-HW ratio) - BUG FIXING: corrected diffusion coefficient in Richards SWC calculation method - DEVELOPMENT: capillary fringe parameter can be set layer by layer (instead of line 76 -> line 98) - BUG FIXING: deleting an unnecessary semicolon in line 53 - DEVELOPMENT: logic of non-structured biomass simulation is modified. Non-stuctured biomass is part of the given plant part (e.q. leaf transfer pool is the part of the leaf, therefore in case of leaf mortality (senescence, management, etc), beside leaf actual pool, leaf transfer and storage pool decrease also. Due to this non-sturetured standing/cut-dow dead biomass pools and fluxes (STDBc_nsc, STDBn_nsc, CTDBc_nsc, CTDBn_nsc) ceased. - DEVELOPMENT: beside aboveground summary variables (live and dead aboveground biomass with/without non-stuctured biomass [3157-3168]), belowground summary variables (live and dead aboveground biomass with/without non-stuctured biomass [3145-3156]) are defined - DEVELOPMENT: Annual maximum of total living and living+dead aboveground/belowground biomass weight are defined [3141-3144]. - DEVELOPMENT: correction of the location of proprortion of NH4 flux of N-deposition parameter (instead of SITE block, NDEP_control block of INI file) - DEVELOPMENT: modification of tipping with more precise estimation of diffusuvity - BUG FIXING: in NGB calculations (Ceq) - BUG FIXING: in denitrification calculation- using SRlayer instead of SRtotal - DEVELOPMENT: mean annual air temperature and annual air temperature range is calculated based on input meteorological data instead of input parameter less parameter in SITE block of INI file new output variables: annTavgRA [77]: 10-year running average of mean annual air temperature annTrangeRA [78]: 10-year running average of mean annual air temperature range pa [79]: Atmospheric pressure - DEVELOPMENT: more structurated hydrological routines, new groundwater calculation method new output variables: GWmovchange [184-193] - Soil water change due to the movement of groundwater table in soil layer 1-10 GWmovchange [184-193] - Soil water change due to the movement of groundwater table in soil layer 1-10 2798] = &sprop->CFlayer; //%"Number of layer containing the upper layer REMOVED rountine: pondwANDrunoff.c NEW routines: diffusCalc.c - calculation of diffusion flux between two soil layers (used in tipping.c and groundwaterT_CFcharge.c) soilEVP_calc.c - calculate limitation of soil evaporation and update top soil water content infiltANDpond.c - calculation of pond water accumulation and potential infiltration pondANDrunoffD.c - state update of pond water, calculation of Dunnian runoff REMOVED/NEW routines: instead of groundwaterRICHARDS.c and groundwaterRICHARDSpostproc.c -> groundwaterR_preproc.c and groundwaterR_postproc.c instead of groundwaterTIPPING.c -> groundwaterT_preproc.c and groundwaterT_CFcharge.c instead of prcp_route.c -> prcpANDrunoffH.c instead of potential_evap.c -> potEVP_calc.c - DEVELOPMENT: irrigation routine - differentiation of above-canopy, below-canopy, and subsurface irrigation new column in irrigation file: the height of the irrigation (if positive - above the surface, if negative - below the surface) - DEVELOPMENT: fertilizing routine - water from fertilization goes into given soil layers fertilizing parameter), in case of oversaturation: pond water formation - DEVELOPMENT: new output variables for greenhouse gas balance analysis rootIndex [3127]: Ratio of fine root C content and vegetation C content on harvest day belowground_ratio [3128]: Ratio of root to vegetation C content cumCH4flux [3129]: Cumulative annual SUM of total CH4 flux cumNBP [3130]: Cumulative annual SUM of NBP N2OfluxCeq [3131]: Daily N2O flux in C eq cumN2OfluxCeq [3132]: Cumulative annual SUM N2O flux in C eq cumCflux_lateral [3133]: Cumulative lateral C flux livingBIOMabove, livingBIOMbelow, totalBIOMabove, totalBIOMbelow [3137-3140]: Living/total aboveground/belowground biomass annmax_livingBIOMabove, annmax_livingBIOMbelow, annmax_totalBIOMabove, annmax_totalBIOMbelow [3137-3140]: Annual maximum of living/total aboveground/belowground biomass LDbelowC_nw, LDbelowC_w, LDbelowCnsc_nw, LDbelowCnsc_w [3145-3148]: Living+dead belowground non-woody/woody biomass C content without/with non-structured carbohydrate LbelowC_nw, LbelowC_w, LbelowCnsc_nw, LbelowCnsc_w [3149-3152]: Living belowground non-woody/woody biomass C content without/with non-structured carbohydrate DbelowC_nw, DbelowC_w, DbelowCnsc_nw, DbelowCnsc_w [3145-3148]: Dead belowground non-woody/woody biomass C content without/with non-structured carbohydrate - DEVELOPMENT: new optional intercepcion calculation in prcpANDrunoff.c temporary solution: optional calculation method using inactive ETflag (in simulation control block of INI file) ETflag=0 - original calculation, ETflag=1 - new calculation