CTSM FATES integration

Goals Notes and milestones for CTSM-FATES integration

v0 integration: Goal Sept 1, 2020

• Reduced complexity FATES [Rosie Fisher]

  • Create a reduced complexity versions of FATES for testing, development and application including FATES-SP and fixed biogeography (with and without competition).
  • Marked as completed, but maybe more work needs to be done here?

• LULCC implementation in FATES [Charlie Koven]

  • Beta version of LULCC applied in historical, transient simulations

    • initially uses area weighted harvest (CLM5 currently uses C mass harvested)
    • complicated to investigate how biases may influence primary and secondary forest in FATES
    • FATES want to know just the fraction of vegetated land unit for each column that's being harvested (which is directly provided from LUH2).
      • Do we back out the preprocessing LUH data that's being provided on land use time series?
      • Area harvest needs to be on land_use_timeseries file.
    • FATES disturbance / harvest happens daily, but harvest data from CTSM provided annually.
      • Currently FATES harvest just happens 1x/year.
    • FATES marks trees that are harvested, surviving trees that are moved to degraded (secondary) land. Additional collateral damage happens from treefall and roadbuilding.
      • LUH2 assumes 30% slash, collateral damage from harvest

  • More work needed to investigate potential issues with forest age classes being tracked in FATES

  • Also questions about using area harvest vs. biomass harvested.

  • Current implementation uses area weights and seems to produce sensible changes in secondary forest area

  • Charlie has questions for Bill and Erik if the implementation is correct.

• Rebase FATES to CTSM MAIN [Gregory Lemieux]

  • FATES needs to migrate to MAIN, not the released version of CTSM.
  • Need conversations about maintaining code & integration between two large and active development communities

• Coupled model smoke test [Keith Oleson]

  • Need appropriate compsets [Erik]
  • Initially just run CAM-CTSM case, as chemistry likely won't work.
  • Worth running chemistry, so we can see what else we'll need (e.g. dry deposition).

• Performance, FATES costs

  • Dave thinks we should have cost estimates for FATES. Rosie notes that we don’t currently do a formal analysis.
  • Greg notes that we note the memcomp and tutcomp warnings in the regression tests only right now. Ryan has done spot checking of cost (notably with PARTEH refactor).
  • Erik notes that the cost estimate is going to be somewhat difficult due to dynamic nature of FATES. Patches and cohorts spawn and merge initially (mostly an issue during spin up from bare ground).
  • Smart load balancing as a way to address performance issues, but Bill notes dynamic load balancing would be challenging.

v1 integration: Goal Jan 1, 2021

• Code refactoring to allow FATES to run on it's own land unit [Erik]

  • Erik will scope out what this will look like, but the idea would be to have FATES run on seperate land units from big leaf CTSM & Crops to make FATES-Crop simulations possible

• Running FATES with the crop model

  • requires code refactoring (above)
  • If crops are running with nitrogen, and natural veg without, need to think about how to initialize nitrogen content in newly-initializing crop columns with dynamic land use. Charlie: FATES does have a nitrogen model that could be turned on.
  • Need to consider fire related to LULCC, clearning, agricultural fires, etc (see next task)?

• Fire & human-fire interactions [Jackie]

  • Lots of additional complexity WRT fire, LULCC & crops.
  • There is a framework to add in population & GDP; Jackie is planning to incorporate them in the next year.
  • May need to add some capabilities for global fire: allow spatially-varying parameters, or have some of the existing parameters depend on other variables
  • Significant discussion on ignition sources related to landuse change & deforestation forest
    • Peter noted roughly 1 of 3 Pg C y-1 fire budget that’s currently handled ‘implicitly’.
  • Jackie is cautiously optimistic based on the range of regions over which this is being tested.

• Parameterization & Individual PFT performance investigations.

  • Calibrate dynamic vegetation model: Ideally with a relatively simple PFT set. If individual PFTs are performing sensibly, this should be easier.
  • Drought deciduous phenology: This remains challenging, as per CLM5. Some clear structural errors need addressing as well as parametric issues.
  • C4 grass: This is historically over productive in CLM. Anthony Walker made a start at looking at this in MAAT.
  • Shrub introduction: Allow shrubs to tolerate harsher growing conditions v.s allometric limits to height etc. Maybe collaboration with EMERALD/NGEE-arctic/?
  • N.b. This is a place where domain experts can potentially take charge of, e.g. of the testing of a new PFT type. Need to think of a way of making those studies publishable. Sensitivity tests? Regional simulations?
  • Also FATES makes this easier than big-leaf mode in that you can run with one-PFT at a time.
  • Jackie highlighted the importance of coexistence for this target. Noted Polly is working on this.

• Seed carbon conservation?

• Shared code base for parameterizations, functions and routines between FATES and CTSM. This should likely live in the FATES repo, but would need to be accessed by CTSM for big-leaf simulations (e.g., phenology & photosynthesis code).

  • Avoid redundancies in code and allow for smoother development in the future
  • Afford opportunities for multi-hypothesis testing in both FATES and big-leaf
  • Requires SE attention on both CTSM and FATES side

• Test and decide whether HYDRO would be on in default

• Isotopes

  • Carbon: the functionality is basically in place for this
  • Water: if plant hydraulics are off, FATES doesn’t handle any water pools / fluxes, but just hands conductance back to the host land model. If plant hydraulics are on, then this will need some work
    • Whether to turn on plant hydraulics is still an open question

• Need for dry deposition to work with FATES. And maybe fire emissions passed to the atmosphere.

• Benchmarking, additional FATES specific benchmarks: NOTE, there are additional discussion on benchmarks in the workshop meeting notes. This short list is intended highlights priorities that are of greatest utility for global, transient simulations.

  • PFT Distributions
    • Sandy Harrison: kappa statistic, simulated vs. observed statistic (Monserud & Leemans 1992)
    • Charlie: Regarding PFT distributions, part of what they’ve been pushing back on is the idea of canonical PFTs. So separating into trees & grasses is fine, but want to be a little cautious about trying to match canonical PFTs.
  • Abby: idea of comparing distribution of traits rather than PFT distributions per se (e.g., isohydricity / anisohydricity, etc.)
  • Tree height data, tree density data
    • Peter: in CLM5-SP, we specify canopy height from ICESat, could be used for benchmarking instead
    • Lidar-based global canopy structure and height distribution (GEDI) https://gedi.umd.edu/ GEDI science team also has a large database of forest metrics for modeling biomass - they are usually willing to share
    • Emulators of direct satellite observations (NASA MAP project?). This would also allow comparison of forest change after disturbance from LandTrender or other satellite obs.
    • S. S. Saatchi, N. L. Harris, S. Brown, M. Lefsky, E. T. A. Mitchard, W. Salas, B. R. Zutta, W. Buermann, S. L. Lewis, S. Hagen, S. Petrova, L. White, M. Silman, and A. Morel. Benchmark map of forest carbon stocks in tropical regions across three continents. Proceedings of the National Academy of Sciences of the United States of America, 108(24):9899–9904, June 2011.

• Hillslope interaction? Probably not an initial target. The main challenge here is the interaction between hillslopes and LULCC.