Skip to content

Hector V3.5.0 Calibration #789

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
kdorheim merged 13 commits into from
Oct 16, 2025
Merged

Hector V3.5.0 Calibration #789

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
13 commits
Select commit Hold shift + click to select a range
1ab145f
update v number
kdorheim Sep 5, 2025
c332249
saving wip
kdorheim Sep 18, 2025
0a659cc
save new diff value
kdorheim Sep 18, 2025
18846f2
update beta
kdorheim Sep 18, 2025
6dc42b2
update q10 rh
kdorheim Sep 18, 2025
5b77bd1
update pi control
kdorheim Sep 23, 2025
eb0b43e
update parameter values
kdorheim Oct 9, 2025
2db162b
update natural emissions
kdorheim Oct 9, 2025
0b257e2
update documenation
kdorheim Oct 9, 2025
f2d6536
add pi control test
kdorheim Oct 9, 2025
de05d05
update old new test data due to behavior changes
kdorheim Oct 9, 2025
445ad1e
update pi control params
kdorheim Oct 10, 2025
d63445e
adding some more documenation
kdorheim Oct 16, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
2 changes: 1 addition & 1 deletion DESCRIPTION
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
Package: hector
Title: The Hector Simple Climate Model
Version: 3.4.9
Version: 3.4.95
Comment thread
kdorheim marked this conversation as resolved.
Authors@R: c(person("Kalyn", "Dorheim",
email = "kalyn.dorheim@pnnl.gov",
role = c("aut", "cre"),
Expand Down
2 changes: 1 addition & 1 deletion inst/include/h_util.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -26,7 +26,7 @@
* \brief The model version number to be included in logs and outputs.
* \note Manually update the git tag to match this.
*/
#define MODEL_VERSION "3.4.9"
#define MODEL_VERSION "3.4.95"

#define OUTPUT_DIRECTORY "output/"

Expand Down
12 changes: 6 additions & 6 deletions inst/input/hector_picontrol.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -81,8 +81,8 @@ CO2_constrain[1745]=277.15 ;

; CO2 and temperature effects on the carbon cycle
; these are global values, can optionally specify biome-specific ones as above
beta=0.55 ; CO2 fertilization factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details
q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see Dorheim et al. in prep for details
beta=0.65 ; CO2 fertilization factor (unitless), calibrated to historical observations see PR 789 for details
q10_rh=1.2 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see PR 789 for details

; Optional biome-specific warming factors
; by default, assume 1.0 (i.e., warms as fast as the globe)
Expand All @@ -107,10 +107,10 @@ SV[1745]=0.185745674 ; volcanic radiative forcing time series
[CH4]
M0=731.41 ; 731.41 ppb preindustrial methane IPCC AR6 Table 7.SM.1, the CH4 forcing equations is calibrated to a M0 of 731.41 ppb
Tsoil=120 ; lifetime of soil sink (years) Myhre et al., 2013, doi: 10.1017/cbo9781107415324.018
Tstrat=150 ; lifetime of tropospheric sink (years) Myhre et al., 2013, doi: 10.1017/cbo9781107415324.018
Tstrat=150 ; lifetime of tropospheric sink (years) Myhre et al., 2013, doi: 10.1017/cbo9781107415324.018
UC_CH4=2.78 ; Tg(CH4)/ppb unit conversion between emissions and concentrations
CH4N[1750]=223 ; Natural CH4 emissions (Tgrams) constant over simluation
CH4_emissions[1745]=337.7313 ; emissions time series
CH4N=csv:tables/picontrol.csv ; Natural CH4 emissions (Tgrams) constant over simluation
CH4_emissions[1745]=0 ;

;------------------------------------------------------------------------
[OH]
Expand Down Expand Up @@ -172,7 +172,7 @@ RF_misc[1750]=0
enabled=1

S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6)
diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations see Dorheim et al. in prep for details
diff=1.042 ; ocean heat diffusivity, cm2/s calibrated to historical observations see PR 789 for details
qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6)

; Optional global temperature constraint
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp119.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,10 +79,10 @@ luc_emissions=csv:tables/ssp119_emiss-constraints_rf.csv
luc_uptake=csv:tables/ssp119_emiss-constraints_rf.csv

; CO2 and temperature effects on the carbon cycle
beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
beta=0.65 ; CO2 fertilization factor (unitless), calibrated to historical observations see PR 789 for details
;pf.beta=0.53
;nonpf.beta=0.53
q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
q10_rh=1.2 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see PR 789 for details
;pf.q10_rh=1.76
;nonpf.q10_rh=1.76

Expand Down Expand Up @@ -152,7 +152,7 @@ TN2O0=132 ; initial lifetime of N2O, years

; An example of setting a time series by individual values
; The model will interpolate as necessary between them
N2O_natural_emissions[1750]=9.72 ; default natural emissions based on a best fit to IPCC AR6 N2O RF
N2O_natural_emissions=csv:tables/core_inputs.csv ; default natural emissions based on a best fit to IPCC AR6 N2O RF
N2O_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
;N2O_constrain=csv:tables/ssp119_emiss-constraints_rf.csv ; N2O concentration constraint

Expand Down Expand Up @@ -180,7 +180,7 @@ RF_misc[1750]=0
;------------------------------------------------------------------------
[temperature]
S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6)
diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details
diff=1.042 ; ocean heat diffusivity, cm2/s calibrated to historical observations see PR 789 for details
qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6)

; Optional global temperature constraint
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp126.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,10 +79,10 @@ luc_emissions=csv:tables/ssp126_emiss-constraints_rf.csv
luc_uptake=csv:tables/ssp126_emiss-constraints_rf.csv

; CO2 and temperature effects on the carbon cycle
beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
beta=0.65 ; CO2 fertilization factor (unitless), calibrated to historical observations see PR 789 for details
;pf.beta=0.53
;nonpf.beta=0.53
q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
q10_rh=1.2 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see PR 789 for details
;pf.q10_rh=1.76
;nonpf.q10_rh=1.76

Expand Down Expand Up @@ -152,7 +152,7 @@ TN2O0=132 ; initial lifetime of N2O, years

; An example of setting a time series by individual values
; The model will interpolate as necessary between them
N2O_natural_emissions[1750]=9.72 ; default natural emissions based on a best fit to IPCC AR6 N2O RF
N2O_natural_emissions=csv:tables/core_inputs.csv ;
N2O_emissions=csv:tables/ssp126_emiss-constraints_rf.csv ; emissions time series
;N2O_constrain=csv:tables/ssp126_emiss-constraints_rf.csv ; N2O concentration constraint

Expand Down Expand Up @@ -180,7 +180,7 @@ RF_misc[1750]=0
;------------------------------------------------------------------------
[temperature]
S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6)
diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details
diff=1.042 ; ocean heat diffusivity, cm2/s calibrated to historical observations see PR 789 for details
qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6)

; Optional global temperature constraint
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp245.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,10 +79,10 @@ luc_emissions=csv:tables/ssp245_emiss-constraints_rf.csv
luc_uptake=csv:tables/ssp245_emiss-constraints_rf.csv

; CO2 and temperature effects on the carbon cycle
beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
beta=0.65 ; CO2 fertilization factor (unitless), calibrated to historical observations see PR 789 for details
;pf.beta=0.53
;nonpf.beta=0.53
q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
q10_rh=1.2 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see PR 789 for details
;pf.q10_rh=1.76
;nonpf.q10_rh=1.76

Expand Down Expand Up @@ -152,7 +152,7 @@ TN2O0=132 ; initial lifetime of N2O, years

; An example of setting a time series by individual values
; The model will interpolate as necessary between them
N2O_natural_emissions[1750]=9.72 ; default natural emissions based on a best fit to IPCC AR6 N2O RF
N2O_natural_emissions=csv:tables/core_inputs.csv ;
N2O_emissions=csv:tables/ssp245_emiss-constraints_rf.csv ; emissions time series
;N2O_constrain=csv:tables/ssp245_emiss-constraints_rf.csv ; N2O concentration constraint

Expand Down Expand Up @@ -180,7 +180,7 @@ RF_misc[1750]=0
;------------------------------------------------------------------------
[temperature]
S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6)
diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details
diff=1.042 ; ocean heat diffusivity, cm2/s calibrated to historical observations see PR 789 for details
qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6)

; Optional global temperature constraint
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp370.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,10 +79,10 @@ luc_emissions=csv:tables/ssp370_emiss-constraints_rf.csv
luc_uptake=csv:tables/ssp370_emiss-constraints_rf.csv

; CO2 and temperature effects on the carbon cycle
beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
beta=0.65 ; CO2 fertilization factor (unitless), calibrated to historical observations see PR 789 for details
;pf.beta=0.53
;nonpf.beta=0.53
q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
q10_rh=1.2 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see PR 789 for details
;pf.q10_rh=1.76
;nonpf.q10_rh=1.76

Expand Down Expand Up @@ -152,7 +152,7 @@ TN2O0=132 ; initial lifetime of N2O, years

; An example of setting a time series by individual values
; The model will interpolate as necessary between them
N2O_natural_emissions[1750]=9.72 ; default natural emissions based on a best fit to IPCC AR6 N2O RF
N2O_natural_emissions=csv:tables/core_inputs.csv ;
N2O_emissions=csv:tables/ssp370_emiss-constraints_rf.csv ; emissions time series
;N2O_constrain=csv:tables/ssp370_emiss-constraints_rf.csv ; N2O concentration constraint

Expand Down Expand Up @@ -180,7 +180,7 @@ RF_misc[1750]=0
;------------------------------------------------------------------------
[temperature]
S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6)
diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details
diff=1.042 ; ocean heat diffusivity, cm2/s calibrated to historical observations see PR 789 for details
qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6)

; Optional global temperature constraint
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp434.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,10 +79,10 @@ luc_emissions=csv:tables/ssp434_emiss-constraints_rf.csv
luc_uptake=csv:tables/ssp434_emiss-constraints_rf.csv

; CO2 and temperature effects on the carbon cycle
beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
beta=0.65 ; CO2 fertilization factor (unitless), calibrated to historical observations see PR 789 for details
;pf.beta=0.53
;nonpf.beta=0.53
q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
q10_rh=1.2 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see PR 789 for details
;pf.q10_rh=1.76
;nonpf.q10_rh=1.76

Expand Down Expand Up @@ -152,7 +152,7 @@ TN2O0=132 ; initial lifetime of N2O, years

; An example of setting a time series by individual values
; The model will interpolate as necessary between them
N2O_natural_emissions[1750]=9.72 ; default natural emissions based on a best fit to IPCC AR6 N2O RF
N2O_natural_emissions=csv:tables/core_inputs.csv ;
N2O_emissions=csv:tables/ssp434_emiss-constraints_rf.csv ; emissions time series
;N2O_constrain=csv:tables/ssp434_emiss-constraints_rf.csv ; N2O concentration constraint

Expand Down Expand Up @@ -180,7 +180,7 @@ RF_misc[1750]=0
;------------------------------------------------------------------------
[temperature]
S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6)
diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details
diff=1.042 ; ocean heat diffusivity, cm2/s calibrated to historical observations see PR 789 for details
qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6)

; Optional global temperature constraint
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp460.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,10 +79,10 @@ luc_emissions=csv:tables/ssp460_emiss-constraints_rf.csv
luc_uptake=csv:tables/ssp460_emiss-constraints_rf.csv

; CO2 and temperature effects on the carbon cycle
beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
beta=0.65 ; CO2 fertilization factor (unitless), calibrated to historical observations see PR 789 for details
;pf.beta=0.53
;nonpf.beta=0.53
q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
q10_rh=1.2 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see PR 789 for details
;pf.q10_rh=1.76
;nonpf.q10_rh=1.76

Expand Down Expand Up @@ -152,7 +152,7 @@ TN2O0=132 ; initial lifetime of N2O, years

; An example of setting a time series by individual values
; The model will interpolate as necessary between them
N2O_natural_emissions[1750]=9.72 ; default natural emissions based on a best fit to IPCC AR6 N2O RF
N2O_natural_emissions=csv:tables/core_inputs.csv ;
N2O_emissions=csv:tables/ssp460_emiss-constraints_rf.csv ; emissions time series
;N2O_constrain=csv:tables/ssp460_emiss-constraints_rf.csv ; N2O concentration constraint

Expand Down Expand Up @@ -180,7 +180,7 @@ RF_misc[1750]=0
;------------------------------------------------------------------------
[temperature]
S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6)
diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details
diff=1.042 ; ocean heat diffusivity, cm2/s calibrated to historical observations see PR 789 for details
qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6)

; Optional global temperature constraint
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp534-over.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,10 +79,10 @@ luc_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv
luc_uptake=csv:tables/ssp534-over_emiss-constraints_rf.csv

; CO2 and temperature effects on the carbon cycle
beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
beta=0.65 ; CO2 fertilization factor (unitless), calibrated to historical observations see PR 789 for details
;pf.beta=0.53
;nonpf.beta=0.53
q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
q10_rh=1.2 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see PR 789 for details
;pf.q10_rh=1.76
;nonpf.q10_rh=1.76

Expand Down Expand Up @@ -152,7 +152,7 @@ TN2O0=132 ; initial lifetime of N2O, years

; An example of setting a time series by individual values
; The model will interpolate as necessary between them
N2O_natural_emissions[1750]=9.72 ; default natural emissions based on a best fit to IPCC AR6 N2O RF
N2O_natural_emissions=csv:tables/core_inputs.csv ;
N2O_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv ; emissions time series
;N2O_constrain=csv:tables/ssp534-over_emiss-constraints_rf.csv ; N2O concentration constraint

Expand Down Expand Up @@ -180,7 +180,7 @@ RF_misc[1750]=0
;------------------------------------------------------------------------
[temperature]
S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6)
diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details
diff=1.042 ; ocean heat diffusivity, cm2/s calibrated to historical observations see PR 789 for details
qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6)

; Optional global temperature constraint
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp585.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -79,10 +79,10 @@ luc_emissions=csv:tables/ssp585_emiss-constraints_rf.csv
luc_uptake=csv:tables/ssp585_emiss-constraints_rf.csv

; CO2 and temperature effects on the carbon cycle
beta=0.53 ; CO2 fertilization factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
beta=0.65 ; CO2 fertilization factor (unitless), calibrated to historical observations see PR 789 for details
;pf.beta=0.53
;nonpf.beta=0.53
q10_rh=1.76 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations Dorheim et al. 2024 for details
q10_rh=1.2 ; Heterotrophic respiration temperature sensitivity factor (unitless), calibrated to historical observations see PR 789 for details
;pf.q10_rh=1.76
;nonpf.q10_rh=1.76

Expand Down Expand Up @@ -152,7 +152,7 @@ TN2O0=132 ; initial lifetime of N2O, years

; An example of setting a time series by individual values
; The model will interpolate as necessary between them
N2O_natural_emissions[1750]=9.72 ; default natural emissions based on a best fit to IPCC AR6 N2O RF
N2O_natural_emissions=csv:tables/core_inputs.csv ;
N2O_emissions=csv:tables/ssp585_emiss-constraints_rf.csv ; emissions time series
;N2O_constrain=csv:tables/ssp585_emiss-constraints_rf.csv ; N2O concentration constraint

Expand Down Expand Up @@ -180,7 +180,7 @@ RF_misc[1750]=0
;------------------------------------------------------------------------
[temperature]
S=3.0 ; equilibrium climate sensitivity for 2xCO2 degC (A.4.4 of IPCC AR6)
diff=2.38 ; ocean heat diffusivity, cm2/s calibrated to historical observations Dorheim et al. 2024 for details
diff=1.042 ; ocean heat diffusivity, cm2/s calibrated to historical observations see PR 789 for details
qco2=3.75 ; 2×CO2 RF (7.3.2 of IPCC AR6)

; Optional global temperature constraint
Expand Down
Loading