Skip to content

Code for Harris et al. "Satellite-Observed Vegetation Responses to Intraseasonal Precipitation Variability"

License

Notifications You must be signed in to change notification settings

ppharris/precipitation-VOD-ISV

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

55 Commits
 
 
 
 
 
 
 
 
 
 

Repository files navigation

precipitation-VOD-ISV

Code for Harris et al. "Satellite-Observed Vegetation Responses to Intraseasonal Precipitation Variability"

This repository contains code to reproduce the analysis and figures of Harris et al. (2022) Satellite-Observed Vegetation Responses to Intraseasonal Precipitation Variability. Geophysical Research Letters 49(15) e2022GL099635: https://doi.org/10.1029/2022GL099635.

Note that this analysis relies upon Fortran code containing an algorithm for the Lomb-Scargle periodogram, from:

Weedon et al. (2015). Evaluating the performance of hydrological models via cross-spectral analysis: case study of the Thames Basin, United Kingdom. Journal of Hydrometeorology, 16(1), 214--231.

This code is not included here, but it is freely available for to anyone to use from the Met Office Science Repository Service at https://code.metoffice.gov.uk/trac/lmed/browser#main/trunk/benchmarking (last access: 29 June 2022). Access requires registration for an account and this will be supported by a member of the JULES group. Requests for new accounts can be made by emailing [email protected] with details of the user’s name, email address, institution and purpose for requiring access.

Cross-spectral analysis

The code can be used to perform the cross-spectral analysis described in Section 2 of the paper as follows:

  • Obtain the Fortran code containing the Lomb-Scargle periodogram algorithm from Weedon et al. (2015) (csagan1.1.f from the Met Office repository described above).
  • Enable multiprocessing with csagan_multiprocess.py following steps in multiprocessing_setup.txt
  • Compile csagan (requires linking to netCDF libraries)
  • Run bash script csa_multiprocess_tiles. This will run the cross-spectral analysis for all pixels/seasons, saving the results in three latitude-band tiles for each season. The save directory is set in csagan_multiprocess.py. A single season takes approximately 8 hours to process and produces ~3--5 GB of saved data.
  • Run significant_coherent_intraseasonal_relationships.py for each tile ('tropics', 'northern', 'southern'), season and frequency-of-variability band. These variables have to be set at the top of the .py file before running. This computes and saves details of the coherent relationships that are 95% significant based on the 3-neighbour condition described in Section 2 of the paper.
  • Run lag_subplots_calc.py to prepare data for Fig 2.

Figures

Other files

About

Code for Harris et al. "Satellite-Observed Vegetation Responses to Intraseasonal Precipitation Variability"

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • Python 99.4%
  • Shell 0.6%