Skip to content

Calibrate soil co2 single site#1618

Open
AlexisRenchon wants to merge 2 commits intomainfrom
ar/em/calibrate_neon
Open

Calibrate soil co2 single site#1618
AlexisRenchon wants to merge 2 commits intomainfrom
ar/em/calibrate_neon

Conversation

@AlexisRenchon
Copy link
Copy Markdown
Member

@AlexisRenchon AlexisRenchon commented Jan 29, 2026

TO DO:

  • latest commit makes calibration model default looks good (like default run dk), but uses the netcdf LAI (not MODIS LAI) - which one do we want?
  • think about good values for RAI and SAI, check if NEE still looks reasonable
  • re-run calibration once mean params and prior spread looks good

wip

@AlexisRenchon AlexisRenchon force-pushed the ar/em/calibrate_neon branch 2 times, most recently from c7b5f68 to 76d30b9 Compare February 10, 2026 18:20
Comment thread experiments/calibrate_dk_sor/generate_observations.jl
Comment thread experiments/calibrate_dk_sor/generate_observations.jl Outdated
# Global noise variances (computed across all calibration years)
all_cal_mask =
(Date(first(cal_years), 1, 1) .<= flux_dates .<= Date(last(cal_years), 12, 31))
nee_var = mean(filter(!isnan, nee_uc_raw[all_cal_mask] .^ 2))
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

we may need to further scale nee_var by a constant to wait terms in the loss function equally. In our offline run with default parameters, the weighted RMSE for need is about 20x larger than the weighted RMSE for the heat fluxes.

julia> x = (lhf_model_at_obs .- qle_obs)./qle_unc;
julia> @show sqrt(mean(x.^2)), mean(x)
(6.310226345552107, -1.9822982854949414)

julia> x = (shf_model_at_obs .- qh_obs)./qh_unc;
julia> @show sqrt(mean(x.^2)), mean(x)
(5.226855236254443, 2.479323425387818)

julia> x = (nee_model_at_obs .- nee_obs)./nee_unc;
julia> @show sqrt(mean(x.^2)), mean(x)
(108.79570442024037, 19.221651412031065)

That is, we might want to scale the NEE uncertainty by 108/620.

Copy link
Copy Markdown
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I'll leave this open for consideration... I'll try running as is for now

Comment thread experiments/calibrate_dk_sor/model_interface.jl
Comment thread experiments/calibrate_dk_sor/model_interface.jl Outdated
Comment thread experiments/calibrate_dk_sor/model_interface.jl Outdated
Comment thread experiments/calibrate_dk_sor/run_calibration.jl Outdated
@AlexisRenchon AlexisRenchon force-pushed the ar/em/calibrate_neon branch 3 times, most recently from 86121bb to ea0290f Compare March 9, 2026 18:30
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants