This dataset contains daily zonal stratospheric methane profile outputs between 1991-2021 simulated by the TOMCAT model.

The TOMCAT simulation is performed at T64L32 resolution that is similar to the one used in Dhomse et al., (2021, 2022) for 1991-2021 time period.

Collocated methane profiles are divided in five latitude bins: SH polar (90S-50S), SH mid-lat (70S-20S), tropics (40S-40N), NH mid-lat (20N-70N) and NH polar (50N-90N). Initially, differences are calculated for each zonal bins for 46 height levels (15km to 60km). Then separate XGBoost regression models are trained for the methane differences between TOMCAT and measurements at each level for a given latitude bin.

The same model is used for all day/night time (2 X11323 days) TOMCAT output sampled at 1.30 am and 1.30 pm local time at the equator. Thus bias corrections for a given model grid are added to the original TOMCAT day and night time profiles. Height resolved data are then interpolated on 28-pressure levels (300 - 0.1hPa). For overlapping latitude bins, we use averages and then calculate daily zonal mean values. For more details regarding this methodology see the associated presentation on Zenodo.

Dataset also includes two files containing daily mean zonal mean methane profiles on height (15-60 km) and pressure (300-0.1 hPa) levels:

- zmch4_TCOM_hlev_T2Dz_1991_2021.nc – height level data (15 to 60 km)
- zmch4_TCOM_plev_T2Dz_1991_2021.nc – pressure level data (300 to 0.1 hPa)

The exact cause of unusual methane variations during 1991-1994 is unknown, however some recent studies argue that it could be due to sudden changes in methane loss processes following Mount Pinatubo eruption as well as significant changes in methane emissions following collapse of the Soviet Union.