README for data files accompanying Fiore et al. (2022)

This file is organized by Directory Names.

=================================================================

Emissions
	This directory includes 4 files containing global and regional total annual anthropgoenic and biomass burning (BB) emissions for NO or CO 
	CMIP6_CEDS_HIST_BB_{NO,CO}_regions.txt includes total annual biomass burning emissions from 1750 through 2014 (the full historical period) described in van Marle et al., 2017: https://doi.org/10.5194/gmd-10-3329-2017
	CMIP6_CEDS_HIST_{NO,CO}_regions.txt includes total annual anthropogenic emissions from 1750-2014 described in Hoesly et al. (2018): https://doi.org/10.5194/gmd-11-369-2018
	The global annual total emissions are shown in Figure 1, and trends over different time periods are given in Table S2.  In addition to the GLOBAL emissions (last column), these files also contain regional emissions (individual columns) for:
	NHML = 30-90N, 0-360E (biomass burning files only)
	MtnW = 30-46N, 115-105W
	WUS = 30-50N, 125-100W
	EUS = 25-50N, 95-60W
	NEUS = 30-50N, 65-90W
	USA = 24-50N, 125-65W
	NA = 15-55N, 125-60W
	EU = 25-65N, 10W -50E
	EA = 15-50N, 95-160E
	SA = 5-35N, 50-95E

==================================================================
Ensemble_LTtrends

trends_latlon_LT_wa6.ic3.csv		trends_latlon_LT_wa6_ic2.003.csv

These files contain the lower tropospheric ozone burden trends shown in Figure 4.  The file trends_latlon_LT_wa6_ic2.003.csv is for ensemble member #7;
trends_latlon_LT_wa6.ic3.csv is for ensemble member #15. 

Column 1, ‘Latitude’ = model latitude
Column 2, ‘Longitude’ = model longitude
Column 3, ‘Trend20’ = 20-year trend (1995-2014)
Column 4,  ‘Trend10’ = 10-year trend (2005-2014)

===================================================================
IAGOS

These files are used in Figures 3 and S4. 

o3_trends_IAGOS_CMIP_free_trop.csv  contains the free tropospheric trends
o3_trends_IAGOS_CMIP_trop_col.csv contains the trends calculated from the full tropospheric column. 

Both files contain the following columns:
Region: See Figure S1 and Figure 1 of Gaudel et al., 2020)
Model: See Table S1
Min Trend: Minimum trend simulated by any ensemble member
Mean Trend: Mean trend simulated by any ensemble member
Max Trend: Maximum trend simulated by any ensemble member

====================================================================
SurfaceSites

CMIP6_surface_sites_monthly.csv 

This file contains the monthly ozone concentrations sampled at the location of the 6 long-term surface monitoring sites, used in Figures 2, 5, and S3.  The file is structured as:

Column 1 = Model Name
Column 2 = Ensemble member
Column 3 = Site Name
Column 4 = Site Location
Remaining columns are the ozone values (monthly mean concentrations) from January 1950 through December 2014. 

These data were used to create the trends reported in the second file in this directory:

CMIP6_surface_sites_monthly_TrendOutput.csv

This file contains the fitted trend and related information for various time periods that is summarized in Table S3, and is structured as:
Model = 		Model Name. Note: bottom-most rows are trends calculated from observations, labeled ‘Observation’)
Scenario = 		Ensemble member 
Site = 			Site Name 
Country = 		Site Location 
Trend_(50_14) = 	Trend over full time period of CESM2-WACCM6 ensemble simulations 
2-sigma_(50-14) = 	2-sigma uncertainty on trend over full time period
P-value_(50-14)= 	p-value for 1950-2014 trend 

Trend_(71_14) = 	Trend over 1971-2014 (period for which Hohenpeissenberg, the site with the longest record, has observations)
2-sigma_(71-14) = 	2-sigma uncertainty on trend over 1971-2014
P-value_(71-14)= 	p-value for 1971-2014 trend 

Trend_(73_14) = 	Trend over 1973-2014 (period when multiple sites (ABO, MLO, Hohenpeissenberg have observations)
2-sigma_(73-14) = 	2-sigma uncertainty on trend over 1973-2014
P-value_(73-14)= 	p-value for 1973-2014 trend 

Trend_(75_14) = 	Trend over 1975-2014 (period when 5 sites (not Cape Grim) have observations)
2-sigma_(75-14) = 	2-sigma uncertainty on trend over 1975-2014
P-value_(75-14)= 	p-value for 1975-2014 trend 

Trend_(82_14) = 	Trend over 1982-2014 (period when all sites have observations)
2-sigma_(82-14) = 	2-sigma uncertainty on trend over 1982-2014
P-value_(82-14)= 	p-value for 1982-2014 trend 

Note that ‘NA’ is used for the Observations in places where data is not available to calculate a trend for a time period (for example, no site has data for the full 1950-2014 simulation period). 

===================================================================
TropO3BurdenCols

ozone_lat_200hPa_all_20220526.csv

This is the main file for tropospheric ozone burdens used in Figures 1, S2, and to calculate the trends used in Table S2, as well as those shown in Figures 6 and 7 (but see also TropO3BurdenColsTrends below.  
For calculating the model trends for comparison with TOMS in the middle panel of Figure 6, we instead use the file 

ozone_lat_100hPa_all_20220526.csv

in which the entries for TROP_pr are calculated up to 100 hPa rather than to 200 hPa. Otherwise these files are identical.

Columns are:

Simulation: 	ensemble name
Year: 		ranges from 1950-2014
Season: 		includes ANNUAL, DJF, MAM, JJA, SON.

Latitude:		 includes a ‘global’ entry and then proceeds by 10 degree latitude, labeled by central latitude, from 85S (-85) to 85N (85).
 
Region: 		TROP_o3 uses the 150 ppb ozone tropopause
TROP_pr uses a 200 hPa or 100 hPa pressure surface for the tropopause. 
TROP_tr uses the WMO thermal lapse rate definition to identify the tropopause
TROP_lv	uses a constant vertical level with nominal 200 hPa pressure to define the tropopause
TROP_tr_DU	uses WMO thermal lapse rate but converts burden to Dobson Units (regional mean)
TROP_lv_DU	as for TROP_lv but burden converted to Dobson Units (regional mean)
UT	upper troposphere
MT	middle troposphere
LT 	lower troposphere

			
Ozone:    		This is the ozone burden (Tg) for all entries except for TROP_tr_DU and TROP_lv_DU which are provided in Dobson Units (DU). 


======================================================================
TropO3BurdensColsTrends

Contains two files
ozone_trend_lat_100hPa_all.csv	ozone_trend_lat_200hPa_all.csv

These trends reported in these files are calculated using the data provided in TropO3BurdensCols.  Please see above for the difference between the two files.

Columns are:

Simulation: 	ensemble name
Year: 		ranges from 1950-2014
Season: 		includes ANNUAL, DJF, MAM, JJA, SON.

Latitude:		 includes a ‘global’ entry and then proceeds by 10 degree latitude, labeled by central latitude, from 85S (-85) to 85N (85).
 
Region: 		TROP_o3 uses the 150 ppb ozone tropopause
TROP_pr uses a 200 hPa or 100 hPa pressure surface for the tropopause. 
TROP_tr uses the WMO thermal lapse rate definition to identify the tropopause
TROP_lv	uses a constant vertical level with nominal 200 hPa pressure to define the tropopause
TROP_tr_DU	uses WMO thermal lapse rate but converts burden to Dobson Units (regional mean)
TROP_lv_DU	as for TROP_lv but burden converted to Dobson Units (regional mean)
UT	upper troposphere
MT	middle troposphere
LT 	lower troposphere


TrendAll: 	Linear trend calculated over entire 1950-2014 simulation period.  This variable is used in Figure 6

Trend20: 	Linear trend over final 20 years of simulation (1995-2014). This variable is used in Figure 7.

Trend10: 	Linear trend over final 10 years of simulation (2005-2014).  This variable is used in Figure 6.

TrendTOMS: 	Linear trend over period of TOMS retrievals (1979-2005) shown in Figure 6.

