/// <summary> /// Adds some biomass for a species to the WOOD litter pools at a site. /// </summary> public static void AddWoodLitter(double woodBiomass, ISpecies species, ActiveSite site) { if (woodBiomass > 0) { WoodLayer.PartitionResidue(woodBiomass, FunctionalType.Table[SpeciesData.FuncType[species]].WoodDecayRate, SpeciesData.WoodCN[species], SpeciesData.WoodLignin[species], LayerName.Wood, LayerType.Surface, site); } }
//--------------------------------------------------------------------- /// <summary> /// Kills coarse roots and add the biomass to the Wood Debris pool. /// </summary> public static void AddCoarseRootLitter(double abovegroundWoodBiomass, ICohort cohort, ISpecies species, ActiveSite site) { double coarseRootBiomass = CalculateCoarseRoot(cohort, abovegroundWoodBiomass); // Ratio above to below if (coarseRootBiomass > 0) { WoodLayer.PartitionResidue(coarseRootBiomass, FunctionalType.Table[SpeciesData.FuncType[species]].WoodDecayRate, SpeciesData.CoarseRootCN[species], SpeciesData.CoarseRootLignin[species], LayerName.CoarseRoot, LayerType.Soil, site); } }
/// <summary> /// Grows all cohorts at a site for a specified number of years. /// Litter is decomposed following growth. /// </summary> public static ISiteCohorts Run(ActiveSite site, int years, bool isSuccessionTimeStep) { ISiteCohorts siteCohorts = SiteVars.Cohorts[site]; IEcoregion ecoregion = PlugIn.ModelCore.Ecoregion[site]; for (int y = 0; y < years; ++y) { Year = y + 1; if (PlugIn.ModelCore.CurrentTime > 0 && Climate.Future_MonthlyData.ContainsKey(PlugIn.FutureClimateBaseYear + y + PlugIn.ModelCore.CurrentTime - years)) { EcoregionData.AnnualWeather[ecoregion] = Climate.Future_MonthlyData[PlugIn.FutureClimateBaseYear + y - years + PlugIn.ModelCore.CurrentTime][ecoregion.Index]; } SiteVars.ResetAnnualValues(site); if (y == 0 && SiteVars.FireSeverity != null && SiteVars.FireSeverity[site] > 0) { FireEffects.ReduceLayers(SiteVars.FireSeverity[site], site); } // Next, Grow and Decompose each month int[] months = new int[12] { 6, 7, 8, 9, 10, 11, 0, 1, 2, 3, 4, 5 }; if (OtherData.CalibrateMode) { months = new int[12] { 6, 7, 8, 9, 10, 11, 0, 1, 2, 3, 4, 5 } } ; for (MonthCnt = 0; MonthCnt < 12; MonthCnt++) { // Calculate mineral N fractions based on coarse root biomass if (MonthCnt == 0) { AvailableN.CalculateMineralNfraction(site); } //PlugIn.ModelCore.UI.WriteLine("SiteVars.MineralN = {0:0.00}, month = {1}.", SiteVars.MineralN[site], i); Month = months[MonthCnt]; SiteVars.MonthlyAGNPPcarbon[site][Month] = 0.0; SiteVars.MonthlyBGNPPcarbon[site][Month] = 0.0; SiteVars.MonthlyNEE[site][Month] = 0.0; SiteVars.MonthlyResp[site][Month] = 0.0; SiteVars.MonthlyStreamN[site][Month] = 0.0; SiteVars.SourceSink[site].Carbon = 0.0; SiteVars.TotalWoodBiomass[site] = Main.ComputeWoodBiomass((ActiveSite)site); //SiteVars.LAI[site] = Main.ComputeLAI((ActiveSite)site); double ppt = EcoregionData.AnnualWeather[ecoregion].MonthlyPrecip[Main.Month]; double monthlyNdeposition; if (EcoregionData.AtmosNintercept[ecoregion] != -1 && EcoregionData.AtmosNslope[ecoregion] != -1) { monthlyNdeposition = EcoregionData.AtmosNintercept[ecoregion] + (EcoregionData.AtmosNslope[ecoregion] * ppt); } else { monthlyNdeposition = EcoregionData.AnnualWeather[ecoregion].MonthlyNDeposition[Main.Month]; } if (monthlyNdeposition < 0) { throw new System.ApplicationException("Error: Nitrogen deposition input data are not present in climate library"); } EcoregionData.MonthlyNDeposition[ecoregion][Month] = monthlyNdeposition; EcoregionData.AnnualNDeposition[ecoregion] += monthlyNdeposition; SiteVars.MineralN[site] += monthlyNdeposition; //PlugIn.ModelCore.UI.WriteLine("Ndeposition={0},MineralN={1:0.00}.", monthlyNdeposition, SiteVars.MineralN[site]); double liveBiomass = (double)ComputeLivingBiomass(siteCohorts); double baseFlow, stormFlow; SoilWater.Run(y, Month, liveBiomass, site, out baseFlow, out stormFlow); // Calculate N allocation for each cohort AvailableN.SetMineralNallocation(site); if (MonthCnt == 11) { siteCohorts.Grow(site, (y == years && isSuccessionTimeStep), true); } else { siteCohorts.Grow(site, (y == years && isSuccessionTimeStep), false); } WoodLayer.Decompose(site); LitterLayer.Decompose(site); SoilLayer.Decompose(site); //...Volatilization loss as a function of the mineral n which // remains after uptake by plants. ML added a correction factor for wetlands since their denitrification rate is double that of wetlands //based on a review paper by Seitziner 2006. double volatilize = (SiteVars.MineralN[site] * EcoregionData.Denitrif[ecoregion]); // monthly value //PlugIn.ModelCore.UI.WriteLine("BeforeVol. MineralN={0:0.00}.", SiteVars.MineralN[site]); SiteVars.MineralN[site] -= volatilize; SiteVars.SourceSink[site].Nitrogen += volatilize; SiteVars.Nvol[site] += volatilize; SoilWater.Leach(site, baseFlow, stormFlow); SiteVars.MonthlyNEE[site][Month] -= SiteVars.MonthlyAGNPPcarbon[site][Month]; SiteVars.MonthlyNEE[site][Month] -= SiteVars.MonthlyBGNPPcarbon[site][Month]; SiteVars.MonthlyNEE[site][Month] += SiteVars.SourceSink[site].Carbon; } } ComputeTotalCohortCN(site, siteCohorts); return(siteCohorts); }