Skip to content
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.

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

Download and use ERA5 low and high veg cover fraction #1004

Open
kmdeck opened this issue Jan 31, 2025 · 0 comments
Open

Download and use ERA5 low and high veg cover fraction #1004

kmdeck opened this issue Jan 31, 2025 · 0 comments
Labels
enhancement New feature or request

Comments

@kmdeck
Copy link
Member

kmdeck commented Jan 31, 2025

Is your feature request related to a problem? Please describe.
ERA5 data includes a timeseries for high vegetation LAI and low vegetation LAI. We had been summing these, but they should be weighted by a coverage fraction. This can be downloaded from ERA5:https://cds.climate.copernicus.eu/datasets/reanalysis-era5-single-levels?tab=overview

Please link to any related issues/PRs.

Describe the solution you'd like

  1. Download these fields.
  2. These are independent of time, but when you download them, theyll have a time index. Remove this.
  3. Add to the era5 clima artifact 2008 and 40 years OR make a new artifact
  4. add to Climaland:
function prescribed_lai_era5(
    era5_lai_ncdata_path,
    era5_lai_cover_ncdata_path,
    surface_space,
    start_date;
    time_interpolation_method = LinearInterpolation(PeriodicCalendar()),
    regridder_type = :InterpolationsRegridder,
)
hv_cover = SpaceVaryingInput(era5_lai_cover_ncdata_path, surface_space, varname_hv)
lv_cover = SpaceVaryingInput(era5_lai_cover_ncdata_path, surface_space, varname_lv)
compose_function(lai_hv, lai_lv; hv_cover = hv_cover, lv_cover = lv_cover)
    return lai_hv* hv_cover + lai_lv * lv_cover # I dont think we need to normalize
end

    return TimeVaryingInput(
        era5_lai_ncdata_path,
        ["lai_hv", "lai_lv"],
        surface_space;
        start_date,
        regridder_type,
        method = time_interpolation_method,
        compose_function = compose_function,
    )
end
  1. Compare timeseries of global mean LAI in 2008 from MODIS and ERA5 with this fix. they should be similar.

Describe alternatives you've considered
We also downloaded MODIS for 2008, but this is only one year. On the other hand, we already have 40 years of ERA5 data.

Additional context
Add any other context or screenshots about the feature request here.

@kmdeck kmdeck added the enhancement New feature or request label Jan 31, 2025
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

1 participant