Atmospheric data used for calibrating the tropopause in global chemistry-climate or chemistry-transport models
Data files
Jan 20, 2025 version files 29.13 MB
-
O3sonde_ASC2.nc
5.73 MB
-
O3sonde_LAU2.nc
14.05 MB
-
O3sonde_WAL2.nc
9.07 MB
-
README.md
14.75 KB
-
tpp_FiguresTables.xlsx
264.63 KB
Abstract
We divide the atmosphere into distinct spheres based on their physical, chemical, and dynamical traits. In deriving chemical budgets and climate trends, which differ across spheres, we need clearly defined boundaries. Focusing on atmospheric mass and greenhouse gases, our primary spheres are the troposphere and stratosphere (~99.9 % by mass), and the boundary between them is the tropopause.
The standard method of locating the tropopause is based on the World Meteorological Organization’s lapse rate tropopause (LRT) defined from the vertical temperature gradient as observed by radiosonde balloons. Every global climate-weather model has one or more methods to calculate the LRT. Both involve subjective choices: expert judgment given meter-scale variability of sonde temperature profiles; methods for calculating gradients from km-thick model layers. Further, LRT and similar methods are consistent only in regions where gradients are primarily vertical (core tropics and midlatitudes) and fail in others (sub-tropical jets and polar regions).
Age-of-air tracers clock the effective time-distance from the tropopause, allowing unambiguous separation of stratosphere from troposphere in the chaotic jet regions. We apply a global model with synthetic tracer e90 (90-day e-folding), focusing on ozone and temperature structures about the tropopause using ozone sonde and satellite observations. We calibrate an observation-consistent tropopause for e90 using tropics-plus-midlatitudes and then apply it globally to calculate total tropospheric air-mass and tropopause ozone values. Such calibration can identify weak tropospheric mixing rates. The concept of calibrating an age-of-air tropopause can be readily applied to other models and even to observed age-of-air tracers like sulfur hexafluoride.
README: Tropopause data used for calibrating the tropopause in global chemistry-climate or chemistry-transport models
https://doi.org/10.5061/dryad.pk0p2ngzt
Description of the data and file structure
The data here focus on tropopause and near-tropopause profiles of P, T, O3 and synthetic age-of-air tracers from sondes, satellites and the UCI chemistry-transport model (CTM).
Files and variables
File: O3sonde_ASC2.nc
Description: Post-processed Ascension I. ozonesonde data. Cleaned to give a consistent ascent profile per text.
Format netcdf4
Dimensions: Levels = 1600
Sondes = 831
Variables:
SondeFile
Size: 831x1
Dimensions: Sondes
Datatype: string
Attributes:
id = 'Ascension O3 sonde files used 19980102:20220928'
rev = 'revised Nov 2024, version 2'
YYYY
Size: 831x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'Year of sonde YYYY'
MM
Size: 831x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'Month of sonde MM'
DD
Size: 831x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'Day of sonde DD'
ntop
Size: 831x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'topmost Level of sonde used here'
NB = 'ntop=0 means sonde not used, total of 44 out of 831'
NB2 = 'yoyo-ing removed, pressure now monotonic, Nov 2024'
pressure
Size: 1600x831
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'pressure (hPa) of measurement'
ozone
Size: 1600x831
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'ozone (mole fraction, ppb)'
temperature
Size: 1600x831
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'temperature (K) of measurement'
NB = 'Some sondes have O3 but no T, flt 614'
potential
Size: 1600x831
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'potential temperature (K) of measurement'
NB = 'Some sondes have O3 but no T, flt 614'
File: O3sonde_WAL2.nc
Description: Post-processed Wallops I. VA ozonesonde data. Cleaned to give a consistent ascent profile per text.
Format: netcdf4
Dimensions: Levels = 1600
Sondes = 1477
Variables:
SondeFile
Size: 1477x1
Dimensions: Sondes
Datatype: string
Attributes:
id = 'Wallops O3 sonde files used 19950725:20201110'
rev = 'revised Nov 2024, version 2'
YYYY
Size: 1477x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'Year of sonde YYYY'
MM
Size: 1477x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'Month of sonde MM'
DD
Size: 1477x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'Day of sonde DD'
ntop
Size: 1477x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'topmost Level of sonde used here'
NB = 'ntop=0 means sonde not used, total of 58 out of 1477'
NB2 = 'yoyo-ing removed, pressure now monotonic, Nov 2024'
pressure
Size: 1600x1477
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'pressure (hPa) of measurement'
ozone
Size: 1600x1477
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'ozone (mole fraction, ppb)'
temperature
Size: 1600x1477
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'temperature (K) of measurement'
NB = 'Some sondes have O3 but no T, flt 614'
potential
Size: 1600x1477
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'potential temperature (K) of measurement'
NB = 'Some sondes have O3 but no T, flt 614'
File: O3sonde_LAU2.nc
Description: Post-processed Lauder NZ ozonesonde data. Cleaned to give a consistent ascent profile per text.
Format: netcdf4
Dimensions: Levels = 1600
Sondes = 1973
Variables:
SondeFile
Size: 1973x1
Dimensions: Sondes
Datatype: string
Attributes:
id = 'Lauder O3 sonde files used 19860803:20210622'
rev = 'revised Nov 2024, version 2'
YYYY
Size: 1973x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'Year of sonde YYYY'
MM
Size: 1973x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'Month of sonde MM'
DD
Size: 1973x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'Day of sonde DD'
ntop
Size: 1973x1
Dimensions: Sondes
Datatype: int32
Attributes:
id = 'topmost Level of sonde used here'
NB = 'ntop=0 means sonde not used, total of 24 out of 1973'
NB2 = 'yoyo-ing removed, pressure now monotonic, Nov 2024'
pressure
Size: 1600x1973
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'pressure (hPa) of measurement'
ozone
Size: 1600x1973
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'ozone (mole fraction, ppb)'
temperature
Size: 1600x1973
Dimensions: Levels,Sondes
Datatype: single
Attributes:
id = 'temperature (K) of measurement'
potential
Size: 1600x1973
Dimensions: Levels,Sondes
Datatype: single
File: tpp_FiguresTables.xlsx
Description: Data sets used to plot Figures 1-9 in the manuscript. Each figure (or pair of) has its own page in the spreadsheet.
Fig01 tab
column variable units, notes number
A
B lower edge index # (1:1:58) 58
C eat-a P = eta-a + eta-b*Psrf (hPa) 58
D eta-b 58
E
F plot alt (km) 58
G
H latitude # (1:1:160) 160
I mid-pt (degrees N) 160
J e90 low tropopause model level 160
K e90 upp tropopause model level 160
L WMO upp tropopause model level 160
Fig02 tab
sonde variable units notes
A ASC sonde#600 z km
B ASC sonde#600 -dT/dz K/km 30 m difference
C ASC sonde#600 -dT/dz K/km 1 km difference
D ASC sonde#600 -dT/dz K/km 1 km average
E ASC sonde#600 dTHETA/dz K/km 30 m difference
F ASC sonde#600 dTHETA/dz K/km 1 km difference
G ASC sonde#600 dTHETA/dz K/km 1 km average
H
I WAL sonde#800 z km
J WAL sonde#800 -dT/dz K/km 30 m difference
K WAL sonde#800 -dT/dz K/km 1 km difference
L WAL sonde#800 -dT/dz K/km 1 km average
M WAL sonde#800 dTHETA/dz K/km 30 m difference
N WAL sonde#800 dTHETA/dz K/km 1 km difference
O WAL sonde#800 dTHETA/dz K/km 1 km average
P
Q LAU sonde#1100 z km
R LAU sonde#1100 -dT/dz K/km 30 m difference
S LAU sonde#1100 -dT/dz K/km 1 km difference
T LAU sonde#1100 -dT/dz K/km 1 km average
U LAU sonde#1100 dTHETA/dz K/km 30 m difference
V LAU sonde#1100 dTHETA/dz K/km 1 km difference
W LAU sonde#1100 dTHETA/dz K/km 1 km average
Fig03 tab see tables in texxt for tropopause criteria
extended rows for each season
column criteria/variable DJF # MAM # JJA # SON #
A Lat # 160 160 160 160
B PTG-tpp average upper troposphere model layer 160 160 160 160
C WMO-LRT average upper troposphere model layer 160 160 160 160
D e/90ppb average upper troposphere model layer 160 160 160 160
E e/80ppb average upper troposphere model layer 160 160 160 160
F e/70ppb average upper troposphere model layer 160 160 160 160
G WMO-alt average upper troposphere model layer 160 160 160 160
H
I PTG-tpp average lower troposphere model layer 160 160 160 160
J WMO-LRT average lower troposphere model layer 160 160 160 160
K e/90ppb average lower troposphere model layer 160 160 160 160
L e/80ppb average lower troposphere model layer 160 160 160 160
M e/70ppb average lower troposphere model layer 160 160 160 160
N WMO-alt average lower troposphere model layer 160 160 160 160
Fig04 tab
TAMF annual mean: rows 5:13
tropopause criteria source Global-upper Global-lower folds blank 60S60N-upper 60S60N-lower blank 25S25N-upper 25S25N-lower
A
B PTG-tpp CTM TAMF X X X X X X
C WMO-LRT CTM TAMF X X X X X X
D e/90ppb CTM TAMF X X X X X X X
E e/80ppb CTM TAMF X X X X X X X
F e/70ppb CTM TAMF X X X X X X X
G WMO-alt CTM TAMF X X X X X X
H
I as reported OMPS TAMF X X
Global
A day number (5:5:365) lines 17:89
B PTG-tpp CTM TAMF lines 17:89
C WMO-LRT CTM TAMF lines 17:89
D e/90ppb CTM TAMF lines 17:89
E e/80ppb CTM TAMF lines 17:89
F e/70ppb CTM TAMF lines 17:89
G WMO-alt CTM TAMF lines 17:89
60S60N
A CTM day number (5:5:365) lines 93:165
B PTG-tpp CTM TAMF lines 93:165
C WMO-LRT CTM TAMF lines 93:165
D e/90ppb CTM TAMF lines 93:165
E e/80ppb CTM TAMF lines 93:165
F e/70ppb CTM TAMF lines 93:165
G WMO-alt CTM TAMF lines 93:165
H OMPS month DJFMAMJJASON lines 93:104
I OMPS month CTM OMPS lines 93:104
25S25N
A CTM day number (5:5:365) lines 169:241
B PTG-tpp CTM TAMF lines 169:241
C WMO-LRT CTM TAMF lines 169:241
D e/90ppb CTM TAMF lines 169:241
E e/80ppb CTM TAMF lines 169:241
F e/70ppb CTM TAMF lines 169:241
G WMO-alt CTM TAMF lines 169:241
H OMPS month DJFMAMJJASON lines 169:180
I OMPS month CTM OMPS lines 169:180
Fig05 tab
A month JFMAMJJASOND tpp criteria
B WAL tpp P (hPa) CTM WMO-alt lines 6:17
C ASC tpp P (hPa) CTM WMO-alt lines 6:17
D LAU tpp P (hPa) CTM WMO-alt lines 6:17
E
F WAL tpp O3 (ppb) CTM WMO-alt lines 6:17
G ASC tpp O3 (ppb) CTM WMO-alt lines 6:17
H LAU tpp O3 (ppb) CTM WMO-alt lines 6:17
A month JFMAMJJASOND tpp criteria
B WAL tpp P (hPa) CTM e/90ppb lines 19:30
C ASC tpp P (hPa) CTM e/90ppb lines 19:30
D LAU tpp P (hPa) CTM e/90ppb lines 19:30
E
F WAL tpp O3 (ppb) CTM e/90ppb lines 19:30
G ASC tpp O3 (ppb) CTM e/90ppb lines 19:30
H LAU tpp O3 (ppb) CTM e/90ppb lines 19:30
A month JFMAMJJASOND tpp criteria
B WAL tpp P (hPa) OMPS lines 32:43
C ASC tpp P (hPa) OMPS lines 32:43
D LAU tpp P (hPa) OMPS lines 32:43
E
F WAL tpp O3 (ppb) OMPS lines 32:43
G ASC tpp O3 (ppb) OMPS lines 32:43
H LAU tpp O3 (ppb) OMPS lines 32:43
A month JFMAMJJASOND tpp criteria
B WAL tpp P (hPa) ACE lines 45:56
C ASC tpp P (hPa) ACE lines 45:56
D LAU tpp P (hPa) ACE lines 45:56
E
F WAL tpp O3 (ppb) ACE lines 45:56
G ASC tpp O3 (ppb) ACE lines 45:56
H LAU tpp O3 (ppb) ACE lines 45:56
A month JFMAMJJASOND tpp criteria
B WAL tpp P (hPa) sondes lines 59:70
C ASC tpp P (hPa) sondes lines 59:70
D LAU tpp P (hPa) sondes lines 59:70
E
F WAL tpp O3 (ppb) sondes lines 59:70
G ASC tpp O3 (ppb) sondes lines 59:70
H LAU tpp O3 (ppb) sondes lines 59:70
Fig06&07 tab
A month JFMAMJJASOND source
B WAL tpp P (hPa) sonde lines 6:17
C ASC tpp P (hPa) sonde lines 6:17
D LAU tpp P (hPa) sonde lines 6:17
E
F WAL tpp O3 (ppb) sonde lines 6:17
G ASC tpp O3 (ppb) sonde lines 6:17
H LAU tpp O3 (ppb) sonde lines 6:17
A month JFMAMJJASOND tpp criteria
B
C
D
E
F WAL tpp O3 (ppb) sonde +-2km lines 19:30
G ASC tpp O3 (ppb) sonde +-2km lines 19:30
H LAU tpp O3 (ppb) sonde +-2km lines 19:30
A month JFMAMJJASOND CTM criteria
B WAL tpp P (hPa) WMO-LRT lines 34:45
C ASC tpp P (hPa) WMO-LRT lines 34:45
D LAU tpp P (hPa) WMO-LRT lines 34:45
E
F WAL tpp O3 (ppb) WMO-LRT lines 34:45
G ASC tpp O3 (ppb) WMO-LRT lines 34:45
H LAU tpp O3 (ppb) WMO-LRT lines 34:45
A month JFMAMJJASOND CTM criteria
B WAL tpp P (hPa) WMO-alt lines 47:58
C ASC tpp P (hPa) WMO-alt lines 47:58
D LAU tpp P (hPa) WMO-alt lines 47:58
E
F WAL tpp O3 (ppb) WMO-alt lines 47:58
G ASC tpp O3 (ppb) WMO-alt lines 47:58
H LAU tpp O3 (ppb) WMO-alt lines 47:58
A month JFMAMJJASOND CTM criteria
B WAL tpp P (hPa) PTgrad lines 60:71
C ASC tpp P (hPa) PTgrad lines 60:71
D LAU tpp P (hPa) PTgrad lines 60:71
E
F WAL tpp O3 (ppb) PTgrad lines 60:71
G ASC tpp O3 (ppb) PTgrad lines 60:71
H LAU tpp O3 (ppb) PTgrad lines 60:71
A month JFMAMJJASOND CTM criteria
B WAL tpp P (hPa) e/70ppb lines 73:84
C ASC tpp P (hPa) e/70ppb lines 73:84
D LAU tpp P (hPa) e/70ppb lines 73:84
E
F WAL tpp O3 (ppb) e/70ppb lines 73:84
G ASC tpp O3 (ppb) e/70ppb lines 73:84
H LAU tpp O3 (ppb) e/70ppb lines 73:84
A month JFMAMJJASOND CTM criteria
B WAL tpp P (hPa) e/80ppb lines 86:97
C ASC tpp P (hPa) e/80ppb lines 86:97
D LAU tpp P (hPa) e/80ppb lines 86:97
E
F WAL tpp O3 (ppb) e/80ppb lines 86:97
G ASC tpp O3 (ppb) e/80ppb lines 86:97
H LAU tpp O3 (ppb) e/80ppb lines 86:97
A month JFMAMJJASOND CTM criteria
B WAL tpp P (hPa) e/90ppb lines 99:110
C ASC tpp P (hPa) e/90ppb lines 99:110
D LAU tpp P (hPa) e/90ppb lines 99:110
E
F WAL tpp O3 (ppb) e/90ppb lines 99:110
G ASC tpp O3 (ppb) e/90ppb lines 99:110
H LAU tpp O3 (ppb) e/90ppb lines 99:110
Fig08 tab
tpp O3 (ppb)
A index # lines 5:164
B latitude (deg N)
C DJF upper tpp PTG-tpp
D WMO-LRT
E e/90ppb
F e/80ppb
G e/70ppb
H WMO-alt
I
J lower tpp PTG-tpp
K WMO-LRT
L e/90ppb
M e/80ppb
N e/70ppb
O WMO-alt
P
Q MAM upper tpp PTG-tpp
R WMO-LRT
S e/90ppb
T e/80ppb
U e/70ppb
V WMO-alt
W
X lower tpp PTG-tpp
Y WMO-LRT
Z e/90ppb
AA e/80ppb
AB e/70ppb
AC WMO-alt
AD
AE JJA upper tpp PTG-tpp
AF WMO-LRT
AG e/90ppb
AH e/80ppb
AI e/70ppb
AJ WMO-alt
AK
AL lower tpp PTG-tpp
AM WMO-LRT
AN e/90ppb
AO e/80ppb
AP e/70ppb
AQ WMO-alt
AR
AS SON upper tpp PTG-tpp
AT WMO-LRT
AU e/90ppb
AV e/80ppb
AW e/70ppb
AX WMO-alt
AY
AZ lower tpp PTG-tpp
BA WMO-LRT
BB e/90ppb
BC e/80ppb
BD e/70ppb
BE WMO-alt
Fig09 tab
column variable range units number
A
B e90 60:2:98 ppb 20
C WAL Altitude (km) of e90 layer 20
D ASC Altitude (km) of e90 layer 20
E LAU Altitude (km) of e90 layer 20
F
G WAL Age (days per km) 20
H ASC Age (days per km) 20
I LAU Age (days per km) 20
J
K e90 61:2:99 ppb 20
L WAL Brunt-Väisälä: N^2 (1e-3 /sec^2) 20
M ASC Brunt-Väisälä: N^2 (1e-3 /sec^2) 20
N LAU Brunt-Väisälä: N^2 (1e-3 /sec^2) 20
O
P WAL -dln(theta)/dln(e90) (unitless) 20
Q ASC -dln(theta)/dln(e90) (unitless) 20
R LAU -dln(theta)/dln(e90) (unitless) 20
Code/software
Any netcdf reader will do for the .nc files. The Figure/Table data sets are given in a .xlsx file with separate page tabs for each Figure.
Access information
Data was collected from the following sources:
- OMPS
OMPS-NPP_LP-L2-O3-DAILY_v2.6
downloaded: 2023-04-01 via wget,
https://data.gesdisc.earthdata.nasa.gov/data/SNPP_OMPS_Level2/...
accepted profiles: 'O3Convergence' <10; ' O3Status' =2:7; ' O3Quality' =0;
'QMV' =0; 'ASI_PMCFlag' =0.
processing: 'TropopauseAltitude' reported using GMAO FP-IT T profiles
* ACE-FTS
ACEFTS_L2_v4p0_O3.nc (includes O3, p, T)
downloaded: 2023-12-17 from doi:10.20383/101.0291
https://databace.scisat.ca/level2/
processing: grid is 1 km altitude; drop all O3 <0 and below 6 km (k=1:6). For tropopause, find lowest troposphere level (LR > 2 K/km) with 2 stratosphere levels (LR ≤2) above it.
* ozonesondes
Wallops (WAL, 38°N), SHADOZV06; https://tropo.gsfc.nasa.gov/shadoz/Wallops.html
Ascension (ASC, 8°S) SHADOZV06: https://tropo.gsfc.nasa.gov/shadoz/Ascension.html
Lauder (LAU, 45°S) NIWA: https://woudc.org/data/explore.php?dataset=ozonesonde
downloaded: .dat's 2023-10-03 (WAL & ASC); csv's 2023-01-08 (LAU)
Methods
This dataset was collected from several sources as outlined in the table below
The CTM 4-D data is published in https://doi.org/10.5061/dryad.qbzkh18qq as netcdf files (CTRL data)
Table 1. Observational data sets |
||||
type |
source |
files |
profiles |
dates |
ozone sondes |
Wallops (WAL, 38°N), SHADOZV06 Ascension (ASC, 8°S) SHADOZV06 Lauder (LAU, 45°S) NIWA |
1,477 831 1,973 |
1,477/tot 831/tot 1,973/tot |
1995-2020 1998-2022 1986-2021 |
|
downloaded: .dat's 2023-10-03 (WAL & ASC); csv's 2023-01-08 (LAU) processing: see Table 3. |
|||
OMPS |
OMPS-NPP_LP-L2-O3-DAILY_v2.6 |
3,587 |
~1,160/day |
2014-2023 |
|
downloaded: 2023-04-01 via wget, https://data.gesdisc.earthdata.nasa.gov/data/SNPP_OMPS_Level2/... accepted profiles: 'O3Convergence' <10; ' O3Status' =2:7; ' O3Quality' =0; 'QMV' =0; 'ASI_PMCFlag' =0. processing: 'TropopauseAltitude' reported using GMAO FP-IT T profiles |
|||
ACE-FTS |
ACEFTS_L2_v4p0_O3.nc (includes O3, p, T) |
1 |
94,675/tot |
2004-2020 |
|
downloaded: 2023-12-17 from doi:10.20383/101.0291 https://databace.scisat.ca/level2/ processing: grid is 1 km altitude; drop all O3 <0 and below 6 km (k=1:6). For tropopause, find lowest troposphere level (LR > 2 K/km) with 2 stratosphere levels (LR ≤2) above it. |
The data was processed to identify the tropopause values of O3, T, potential temperature using several tropopause algorithms described in the table below
Table 2. CTM and the different tropopause definitions |
|
CTM grid cells |
The Chemistry-Transport Model here has a regular latitude-by-longitude 160x320 grid of resolution ~1.1º and 57 eta-coordinate levels defining the pressure edges of each layer based on the T159L60N80 grid of the ECMWF IFS model. The lowest 5 IFS L60 layers are combined into 2 CTM near-surface layers. The vertical resolution of model layers In the upper-troposphere/lower-stratosphere increases regularly from 0.8 km near 10 km altitude to 1.1 km at 17 km altitude. The CTM has layer numbers increasing upward. |
e90 tracer |
Synthetic chemical tracer, emitted uniformly everywhere from the surface, decaying with a 90-day e-fold time. Calculated in the CTM for every cell and every time step. Emissions are set to give an atmospheric mean abundance of 100 ppb. This techniques has been used to study O3 seasonality at the tpp (Prather et al., 2011) |
e/90ppbU |
Highest altitude model layer with e90 > 90 ppb is designated the uppermost tropospheric layer, and its upper boundary is the upper tropopause. The tropopause value for O3 and T is the mean value of the uppermost tropospheric layer. No interpolation is done because of the large change in dO3/dz. |
e/90ppbL |
Lowest altitude model layer with e90 ≤ 90 ppb is designated the lowest stratospheric layer, and its lower boundary is the lower tropopause. Often the same as e/90ppbU. |
e/80ppbU |
Highest altitude model layer with e90 > 80 ppb, as above |
e/80ppbL |
Lowest altitude model layer with e90 ≤ 80 ppb, as above |
e/70ppbU |
Highest altitude model layer with e90 > 70 ppb, as above |
e/70ppbL |
Lowest altitude model layer with e90 ≤ 70 ppb, as above |
Lapse Rate |
Lapse Rate (LR = -dT/dz, K/km) is calculated between 2 vertically aligned layers using the mean temperatures of each layer and the mid-layer altitude. The LR values and thresholds (LR <2 K/km = strat) apply to the upper boundary of the lower box used to calculate it. |
WMO-LRTU |
WMO LR Tropopause (Upper). Here we try to match the WMO definition as closely while following the logic of the PTG algorithm. The Lapse Rate (LR = -dT/dz, K/km) is calculated between 2 adjacent layers using the mean temperatures of each layer and the mid-layer altitude. Two LRs for layer k are calculated: LR1 = -[T(k+1)-T(k)]/[z(k+1)-z(k)]; and LR2 = -[T(k+2)-T(k)]/[z(k+2)-z(k)]. Thus, LR1 spans ~1 km above the mid-point of layer k and LR2 spans ~2 km above. All layers below 4.4 km are set as tropospheric (trop). All layers above 31 km are stratospheric (strat). Define strat layers as meeting the criteria LR1 <2 K/km. If the first (lowermost) strat layer has 2 trop layers immediately above, then relabel it trop and keep going. Having found the first strat layer that passes this test, we check that either LR2 from the current layer or LR1 from the layer above meets the <2 K/km criterion. If these conditions are met, we have the first strat layer and the tpp is the lower boundary of that layer. Now look for a 2nd tpp above: If there is an LR1 ≥ 2 K/km above (i.e., trop layer) then define the 2nd tpp (WMO-LRTU) as the lower boundary of the next strat layer. Otherwise, the 2nd tpp is the same as the 1st. |
WMO-LRTL |
WMO LR Tropopause (Lower). This is just the 1st tpp from WMO-LRTU. |
WMO-altU |
Alternate WMO Upper tpp. Similar to WMO-LRT but easier to implement with CTM grid. Identify trop layers from LR > 2 K/km. Set layers below 6 km altitude to trop. Set layers with θ > 500 K to strat. From the bottom, find the first trop layer with a gap of 2 strat layers (~2 km thick region) above it. This is the 2nd (Upper) tpp. Repeat looking for a trop layer with a gap of at least 1 strat layer above. This is the 1st (Lower) tpp. Both tpp can, and often are, the same, but strat-trop folds resolved by the CTM are easily identified. |
WMO-altL |
Alternate WMO Lower tpp. The 1st (Lower) tpp from WMO-altU search above. |
Potential Temperature Gradient |
The PTG (dθ/dz, K/km) is calculated like the LR between 2 vertically aligned layers using the mean θ of each layer and the mid-layer altitude. The PTG values and thresholds (PTG ≥ 10 K/km) apply to the upper boundary of the lower box used to calculate it. |
PTG-tppU |
PTG Tropopause (Upper). The algorithm here is similar to the original paper (Tinney et al., 2022), adapted to work efficiently with the CTM diagnostics. The PTG (dθ/dz, K/km) is calculated between 2 adjacent layers using the mean θ of each layer and the mid-layer altitude. Two PTGs for layer k are calculated: PTG1 = [θ(k+1)-θ(k)]/[z(k+1)-z(k)]; and PTG2 = [θ(k+2)-θ(k)]/[z(k+2)-z(k)]. Thus, PTG1 spans ~1 km above the mid-point of layer k and PTG2 spans ~2 km above. All layers below 4.4 km are set as trop. All layers above 31 km are strat. Define strat layers as meeting the criteria PTG ≥ 10 K/km. If the first (lowermost) strat layer has 2 trop layers immediately above, then relabel it trop and keep going. Having found the first strat layer that passes this test, we check that either PTG2 from the current layer or PTG1 from the layer above meets the ≥10 K/km criterion. If these conditions are met, we have the first strat layer and the 1st tpp is the lower boundary of that layer. Now look for a 2nd tpp: If we find a trop layer above the strat layer above the 1st tpp, then we search for the first strat layer above it using a revised threshold, PTG1 ≥ 15 K/km. This identifies the 2nd tpp (PTG-tppU) as the lower boundary of the next strat layer. Otherwise, the 2nd tpp is the same as the 1st. |
PTG-tppL |
PTG Tropopause (Lower). This is just the 1st tpp from the PTG-tppU search above. |
The ozonesonde data were preprocessed to obtain stable and consisten profiles of O3 and T vs. P as described in Table 3
Table 3. Ozonesonde processing and defining the tropopause |
|
find O3 NaNs |
Flagged (e.g., -999); negative; too large (>20 ppm); pressure (< 1 hPa) |
find T NaNs |
Flagged; T < 160 K |
collapse profile |
Remove: O3 NaN pts; yoyo (down-up) sections; descent part of profile; all p < 10 hPa. |
drop whole profiles |
With <25 O3 points; with top p > 65 hPa = ASC(44), LAU(24), WAL(58). |
reduce # pts |
Many profiles have 6000-8000 pts & resolution <10 m; create 4-pt (ASC & WAL) and 6-pt (LAU) averages to get 30-50 m spacing for p = 70-500 hPa. |
calculate 100-m gradients |
Use centered 2nd-order finite difference to get 100-m values for LR = -dT/dz and PTG = dθ/dz at each point. These are too noisy to identify tropopause, see Fig. 2 |
calculate stable 1-km gradients |
Average z, T, O3 for 20 points (~1 km) below/above each point; difference these to get a smoothed 1-km-average LR and PTG, see Figure 2. Due to the smoothed 1-km resolution profile, we do not apply the WMO criterion that gradients must be sustained above the threshold point. |
lower tropopause |
LRT: point below lowest LR ≤ 2 K/km, limited to 500 > p > 80 hPa PTG: point below lowest dθ/dz ≥ 10 K/km, ibid |
upper tropopause |
LRT: highest LR > 2 K/km, limited to 500 > p > 80 hPa; PTG: highest dθ/dz < 10 K/km, ibid; Thresholds are not changed for the upper LRT as in Tinney et al., 2022. |