All,
I calculated Soil C concentration (%SOC) from CLM5 SoilC_vr and soil bulk density outputs as:
%SOC = SoilC_vr (g/m3) / soil_bulk_density (g/m3) * 100
The calculated values make sense, however, I noticed that bulk density and vertically resolved soil C are not connected to each other in the CLM5 code. In fact, soil bulk density appears to be soil bulk density of wet soil:
soil_bulkdensity(c,j) = bd(c,j) + h2osoi_liq(c,j)/col%dz(c,j)
bd = (1._r8-soilstate_inst%watsat_col(c,lev))*2.7e3_r8
soilstate_inst%watsat_col(c,lev) = 0.489_r8 - 0.00126_r8*sand
bd = bulk density of dry soil material [kg/m^3]
h2osoi_liq = liquid soil water (kg/m2)
watsat = volumetric soil water at saturation
Therefore my question, if this is a reasonable approach to acquire %SOC? And if not, if there is another way to calculate it.
Thanks a lot.
I calculated Soil C concentration (%SOC) from CLM5 SoilC_vr and soil bulk density outputs as:
%SOC = SoilC_vr (g/m3) / soil_bulk_density (g/m3) * 100
The calculated values make sense, however, I noticed that bulk density and vertically resolved soil C are not connected to each other in the CLM5 code. In fact, soil bulk density appears to be soil bulk density of wet soil:
soil_bulkdensity(c,j) = bd(c,j) + h2osoi_liq(c,j)/col%dz(c,j)
bd = (1._r8-soilstate_inst%watsat_col(c,lev))*2.7e3_r8
soilstate_inst%watsat_col(c,lev) = 0.489_r8 - 0.00126_r8*sand
bd = bulk density of dry soil material [kg/m^3]
h2osoi_liq = liquid soil water (kg/m2)
watsat = volumetric soil water at saturation
Therefore my question, if this is a reasonable approach to acquire %SOC? And if not, if there is another way to calculate it.
Thanks a lot.