Next Article in Journal
Spatial and Temporal Characteristics of Drought Events in Southwest China over the Past 120 Years
Next Article in Special Issue
Comparative Study of the Atmospheric Gas Composition Detection Capabilities of FY-3D/HIRAS-I and FY-3E/HIRAS-II Based on Information Capacity
Previous Article in Journal
Adversarial Robustness Enhancement of UAV-Oriented Automatic Image Recognition Based on Deep Ensemble Models
Previous Article in Special Issue
Characterization of Bias in Fengyun-4B/AGRI Infrared Observations Using RTTOV
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Construction of the Global Reference Atmospheric Profile Database

1
Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
2
University of Chinese Academy of Sciences, Beijing 100101, China
3
College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 101408, China
4
College of Resources and Environment, Hefei University of Technology, Hefei 230009, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(12), 3006; https://doi.org/10.3390/rs15123006
Submission received: 30 April 2023 / Revised: 3 June 2023 / Accepted: 5 June 2023 / Published: 8 June 2023
(This article belongs to the Special Issue Advances in Infrared Observation of Earth’s Atmosphere II)

Abstract

:
Atmospheric profiles are important input parameters for atmospheric radiative transfer models and atmospheric parameter inversions. The construction of regionally representative reference atmospheric profiles can provide basic data for global atmospheric and environmental research. Most reference atmospheric profile databases commonly used lag behind in updating frequency. These databases usually have limited spatial and temporal resolution and differ greatly from the real atmospheric state. To present the real atmospheric state, this article constructs the Global Reference Atmospheric Profile Database (GRAP) based on ACE-FTS satellite products of 2021 and 2022, AIRS satellite products and ERA5 reanalysis data of 2022 u6sing a random forest regression model and a hierarchical mean algorithm. The radiance spectrum of FY-3E HIRAS-II using different profile databases was simulated and compared with the measured spectrum. The results show that GRAP spectral simulations fit better with the measured HIRAS-II spectrum. Comparing the CO2, CH4, O3 and N2O profiles of GRAP, AFGL, MIPAS, RTTOV and NDACC ground station profiles in equatorial, mid-latitude summer and polar winter, the results show that GRAP has high spatial and temporal resolution and better fits the current real atmospheric state. Comparing the temperature profiles of eight regions in China, the results illustrate that GRAP is a better representation of the state of the atmosphere in the Chinese region. GRAP can provide fundamental atmospheric data for radiative transfer studies and atmospheric parameter inversions.

1. Introduction

Atmospheric profiles describe the state of the atmosphere and fundamentally determine its optical properties. Atmospheric profile sample datasets and reference atmospheric profile databases are widely used in research on atmospheric radiative transfer models, atmospheric parameter inversions, simulations of the spectral properties of new satellite instruments and satellite data assimilation [1,2,3,4,5]. As the global atmospheric environment changes, the data on atmospheric profiles used for model development and instrument accuracy verification needs to be continuously updated. Therefore, the construction of a regionally representative global reference atmospheric profile database is of significant importance for the atmospheric environment and global change research.
Atmospheric profile sample datasets can be used to estimate the statistical properties of the background fields. Currently, there are several versions of atmospheric profile sample datasets commonly used internationally, such as TIGR (Thermodynamic Initial Guess Retrieval) [6], ECMWF (European Centre for Medium-Range Weather Forecast) 31L-SD, 50L-SD [6,7], 60L-SD [8,9,10], NESS-35 [11,12], NOAA88 [13,14], etc. Each sample dataset contains different atmospheric profile parameters, data sources and sample sizes due to their different intended applications. The TIGR is an atmospheric profile sample dataset created by the French Laboratory for Dynamical Meteorology, and there are currently five versions available. These profiles were selected from a large number of atmospheric samples from different periods around the world using topological methods. The ECMWF used the same methods as the TIGR and created the 31L-SD, 50L-SD, 60L-SD and 91-L short-range forecast atmospheric profile sample datasets. The ECMWF creates the NOAA88 atmospheric profile sample dataset of 7547 sounding profiles and the ECMWF-52 sample dataset with 52 atmospheric temperatures, humidity and ozone profiles in two atmospheric height level formats, 60 and 101 layers. The existing studies have analyzed these atmospheric sample datasets and found that only the TIGR-43 sample dataset contains one atmospheric profile located on Dachen Island, Zhejiang Province, China [15]. The other atmospheric sample datasets generally lack atmospheric samples that are representative of the Chinese region. To address this issue, Qi Chengli used the topological sampling method to establish the CRASD-1 and CRASD-2 sample datasets with characteristics specific to the Chinese region [15,16]. Due to the large latitude span, complex topography and diverse climate of China, using a single atmospheric profile to represent the whole Chinese region is unreasonable and may cause significant errors in research and analysis. Therefore, it is crucial to improve the Chinese regional atmospheric profiles.
The reference atmospheric profile databases are primarily used for the application performance evaluation and accuracy verification of satellite detectors, radiative transfer models and atmospheric inversion method models. The databases should include meteorological parameters such as air pressure, temperature, gas composition and profile distribution. For internationally popular atmospheric radiative transfer software, LOWTRAN [17], MODTRAN [18], LBLRTM [19], FASCODE [20] and RFM [21] are used in the six reference atmospheric profiles created by the US Air Force Geophysical Laboratory (AFGL), which are the tropical (15°N) atmosphere, the mid-latitude summer (45°N, July) atmosphere, the mid-latitude winter (45°N, January) atmosphere, the sub-polar summer (60°N, July) atmosphere, the sub-polar winter (60°N, January) atmosphere, and the 1976 US Standard Atmosphere [22,23]. The six reference atmospheric profiles take into account the changes of atmospheric parameters with latitude and season, but their spatial–temporal distribution only represents the summer and winter of latitude zones without considering the influence of longitude on atmospheric parameters or the seasonal changes in spring and autumn. Furthermore, the reference atmospheric profiles are updated less frequently. With the intensifying global climate change, atmospheric parameters such as global temperature, CO2, CH4 and O3 have undergone significant changes compared to previous ones, and the delay in updates could cause great errors in the application of studies applying these reference atmospheric profiles.
To address the issues of the existing reference atmospheric profile databases, such as long update periods, large spatial resolution and inadequate consideration of seasonal changes in spring and autumn, this article uses the ACE-FTS Level 2 Version 4.1 products in 2021 and 2022, the AIRS Support Level 2 Version 7 products, and ERA5 reanalysis data in 2022 to create the Global Reference Atmosphere Profile Database (GRAP) through the use of a random forest regression model and a stratified mean algorithm. The objective is to provide data support for research on global climate change and atmospheric component inversion.

2. Materials and Methods

2.1. Data Sources

The data sources used in this article include ACE-FTS Level 2 Version 4.1 satellite products in 2021 and 2022, AIRS Support Level 2 Version 7 satellite products, and ERA5 reanalysis data in 2022.
(1)
ACE-FTS L2 products
The ACE-FTS instrument was launched on 12 August 2003 on board the SciSat-1 satellite. It has a spectral resolution of 0.02 cm−1, a vertical resolution of 1–2 km, a horizontal resolution of 500 km, a wavelength range of 750–4400 cm−1 (2.2–13.3 μm), and a high vertical resolution using occultation for atmospheric sounding [24]. ACE-FTS Level 2 Version 4.1 is a global dataset that includes pressure, temperature, and more than 40 atmospheric constituents such as CO2, CH4, H2O, O3 and N2O for the period 2004 to 2023.
(2)
AIRS L2 products
The Atmospheric Infrared Sounder (AIRS) on the EOS Aqua Spacecraft was launched on 4 May 2002. It orbits from one pole of the Earth to the other about fifteen times a day, covering the same region of the Earth twice a day. AIRS detects wavelengths in the range 650–2700 cm−1 (3.7–15.4 μm) and has a total of 2378 spectral channels [25]. AIRS can detect vertical profiles of atmospheric temperature and humidity, as well as the trace greenhouse gases CO2, CH4, O3 and CO.
(3)
ERA5 reanalysis data
ERA5 reanalysis data is the fifth generation of atmospheric reanalysis products produced by ECMWF, providing hourly data and monthly averages for many atmospheric, land surface and sea state parameters. ERA5 reanalysis data covers the time period from 1940 to the present, with daily ERA5 data updates currently 5 days behind real time. The data is stored in a globally gridded data format, GRIB and NetCDF, with a spatial resolution of 0.25° × 0.25°, vertical coverage from 1000 hPa to 1 hPa and a vertical resolution of 37 pressure layers [26].

2.2. Methods

In this article, considering the influence of time and space on the atmospheric state, GRAP is divided into January to December according to the month, and the globe is divided into 38 latitude zones and 14 longitude zones, each spanning 5° in latitude and 30° in longitude (2.5° for 0n, 0s, 90n, 90s and 15° for 0e, 0w, 180e, 180w). This division results in a total of 532 grids. GRAP includes two atmospheric state parameters as well as 59 atmospheric component parameters, which are detailed in Table 1. The flow chart of the methods of creating GRAP is shown in Figure 1.

2.2.1. Database Metadata

Information on the name, data format, spatial and temporal resolution, height stratification format and atmospheric parameters of the global reference atmospheric profile database is in Table 1.

2.2.2. Atmospheric Profile Samples Acquisition

The construction of GRAP necessitates meeting the demands for extensive spatial and temporal coverage as well as comprehensive atmospheric composition. Simultaneously, it is crucial to ensure that the reference atmospheric profiles accurately represent the local atmospheric conditions. The primary challenge lies in extracting realistic profile samples from a vast collection of atmospheric profiles.
Temperature, pressure and volume mixing ratio (vmr) profiles of different atmospheric constituents are derived from the ACE-FTS, AIRS and ERA5 datasets. These profiles are then aligned with a standardized global grid, considering the detection time and latitude/longitude information, resulting in a dataset of 61 atmospheric profiles covering the entire globe. The extracted atmospheric profiles encompass data from all seasons throughout the year for 365 days, providing comprehensive spatial coverage. Figure 2 illustrates the distribution of CH4 profile samples across the global grid, spanning from January to December.

2.2.3. Data Quality Control

To ensure the reliability of the satellite soundings, rigorous quality control measures are implemented to identify and eliminate erroneous profiles caused by detection errors, instrument malfunctions, cloud effects, or other exceptional circumstances. These quality control procedures aim to refine the original profile samples by excluding any incomplete or erroneous values.
The AIRS data comprises the volume mixing ratio (vmr) and the quality control mark (QC) of the atmospheric constituents. QC values are assigned as follows: 0 indicates the highest quality, 1 indicates good quality, and 2 indicates unusable. AIRS data include the vmr and the quality control mark (QC) of the atmospheric constituents. Profiles with a QC mark of 2 were excluded from the AIRS data.
The ACE-FTS data consists of volume mixing ratio (vmr) and vmr errors for different gaseous atmospheric constituents. Compute the error ratio for the volume mixing ratio and eliminate profiles with an error ratio exceeding 15%.
Since temperature, pressure and atmospheric composition vary with height with a certain regularity and the difference between adjacent heights is within a certain range [27], the vertical consistency is used to test the quality of the profiles and use Equation (1) to calculate the rate of vertical change of the atmosphere for each height layer.
d x h = x h + 1 x h H h + 1 H h
where d x h is the rate of vertical change of the profile in the height layer of h, x h is the profile value in the h layer and H h is the corresponding height layer.
The standard deviation of the vertical rate of change of the profiles, σ , as shown in Equation (2), was calculated. Moreover, the error profiles exceeding three σ were removed.
σ = i = 1 n x i μ 2 n
where x i is the rate of vertical change of the profile in layer i and μ is the average of x i across height layers of the profile.

2.2.4. Missing Data Filling

From Figure 1, it can be seen that there are missing atmospheric profiles in some of the grids, resulting in incomplete global data. In this article, the profiles in the missing grids are interpolated using interpolation in the space–time domain to fill the data. Interpolation is performed on the time domain based on the time series, using the data closest in time to the missing data to fill in. In the spatial domain, the inverse distance weight interpolation (IDW) is used as in Equation (3) [28]. using observations around the location of the interpolation point to fill in the missing positions.
Z 0 ^ = i = 0 n ( Z i Q i )
where Z 0 ^ is the estimated value at the point ( x 0 , y 0 ), Q i is the estimated weight coefficient of the interpolated point corresponding to the observed point, and n denotes the number of interpolated points. Q i is shown in Equation (4) as follows:
Q i = f ( d e j ) j = 1 n f ( d e j )
where n is the number of known observation points and f ( d e j ) denotes the weight function of the known distance d e j between the known observation points and the interpolated points. Equation (5) for f ( d e j ) is as follows:
f ( d e j ) = 1 d e j

2.2.5. Standardization of Atmospheric Profiles

The pressure levels of the different data sources depend on the effective sounding altitude of the instrument. The ACE-FTS data provide an altitude range of 0.5 to 149.5 km, corresponding to a pressure range of 1013 to 3.22 × 10−6 hpa, divided into 150 pressure levels. The AIRS data provide a pressure range of 1100 to 1.61 × 10−6 hpa, divided into 100 pressure levels. The ERA5 reanalysis data provides a pressure range of 1000 to 1 hPa, divided into 37 pressure levels. In this article, all profile samples are interpolated onto a uniform elevation grid. The profiles have an elevation range of 0–119 km and a vertical interval of 1 km over the entire height range (Table 2 gives the three data sources and the GRAP height range level) [29]. A non-linear relationship between height and sample contour values is constructed in each grid, and each profile is interpolated to a standard height grid using a spline function interpolation method.

2.2.6. Creation of the Global Reference Atmospheric Profile Database

The data sources for CH4, CO2, O3 and temperature profile samples are mainly from AIRS satellite data and ERA5 reanalysis data, which are large in number. The four atmospheric profile samples in each grid were fitted with a random forest regression model to obtain a standard profile representing that grid. Random Forest (RF) is an algorithm that uses multiple trees to train and predict a sample [30]. There are two advantages to using a random forest regression model:
(1)
The random forest model has a random nature in sample extraction and feature selection, and the algorithm is not prone to over-fitting;
(2)
When creating the random forest, an unbiased estimate of the Generalization Error is used, and the model has a strong generalization capability;
The random forest regression model construction process is as follows. Figure 3 shows the Random Forest Model Construction Flowchart.
(1)
Using the Bootstrap sampling method with put-back, n samples are randomly selected from the original dataset, and the samples that are not drawn (Out of Bag, OBB) form the test set;
(2)
Construct n decision trees, select m features from the training sample data, choose the best feature to split, and keep splitting each tree until all training samples at that node belong to the same class;
(3)
Repeat both steps (1) and (2), and finally form the generated multiple classification trees into a random forest regression model;
(4)
Integrate all the generated decision trees for prediction to obtain the final prediction results.
Other atmospheric component profile sample data sources are mainly from ACE-FTS products. The amount of data is small for all the profile samples in the same grid to adopt the stratified mean method, according to each height layer, to find the mean value of the profile samples to obtain a standard profile to represent the atmospheric parameters in this grid. Equation (6) for the stratified mean method is as follows:
x = i = 1 n x i n

3. GRAP-Based Simulation Validation

This study employs the RFM atmospheric radiative transfer model to simulate the location and absorption intensity of absorption spectral lines for CO2 and CH4 in both the full band and sensitive band ranges. The simulations select the Fengyun-3E HIRAS-II sample situated at mid- to low-latitudes (19.12°N, 98.65°E) in October. The parameters of the RFM were configured to match the surface temperature, surface reflectance, observation geometry, and spectral resolution (0.625 cm−1) of the HIRAS-II transit moment sampling image element. The spectrum data from HIRAS-II was captured on 10 October 2022, at 11:40 AM, with a spatial resolution of 14 km. The file name associated with the data is FY3E_HIRAS_GRAN_L1_20221010_1140_014KM_V0.HDF.
Figure 4 presents the comparison between the spectrum data obtained from HIRAS-II and the simulated spectrum of the RFM in the three atmospheric models across the full band range of 650–2550 cm−1. The results in Figure 4 indicate that the simulated spectrum of GRAP exhibits a much closer agreement with the measured spectrum from HIRAS-II.
The absorption characteristics of different gases vary across different spectral bands. Figure 5 illustrates the spectrum within the range of 650–760 cm−1, which corresponds to the strong absorption band of CO2. This band is influenced by interfering gases such as H2O, N2O, O3 and HNO3. Similarly, Figure 6 presents the spectrum within the range of 1200–1400 cm−1, which represents the strong absorption band of CH4. This band is affected by interfering gases such as H2O, N2O, CO2, CF4 and O3. By comparing the four spectrum curves in Figure 7 and Figure 8, it is evident that the measured HIRAS spectrum (represented by the red solid line) closely aligns with the simulated GRAP spectrum (represented by the green dashed line). The deviations between the three simulated spectrum curves and the HIRAS-II measured spectrum curves are calculated. Figure 7 presents the GRAP simulated CO2 absorption band spectrum within −10% to 12%, and Figure 8 presents the GRAP simulated CH4 absorption band spectrum within −30% to 25%. The simulated spectrum of GRAP exhibits smaller deviations compared to the simulated spectrum of AFGL and MIPAS. This finding indicates that the atmospheric profile values employed in the RFM model within GRAP exhibit better consistency with the true values of the current atmospheric state.

4. Comparison of Reference Profiles and Discussion

This article compares and analyzes the profiles of CO2, CH4, O3 and N2O in different latitude zones in the GRAP database, the AFGL database [22], the MIPAS database [31], the RTTOV database [32], and the profiles of CH4, O3 and N2O of the NDACC ground station in different latitude zones. It further selects and analyzes the temperature reference profiles of eight different grids in the Chinese region.

4.1. Comparison of Equatorial Reference Profiles

This study compares various atmospheric profiles in the equatorial climate zone, including the reference profiles of GRAP for the 0n0e grid in July, the equatorial reference profiles of the AFGL, the MIPAS equatorial reference profiles, the RTTOV reference profiles, and the profiles measured in July 2021 at the Izaña ground station in Tenerife, Spain, which is located at the equator and affiliated with the Network for the Detection of Atmospheric Composition Change (NDACC).
The CO2 profiles in Figure 9a exhibit overall consistency among the four databases. However, below 80 km altitude, the CO2 profile values of AFGL are 330 ppmv, whereas the CO2 profile values of MIPAS, RTTOV and GRAP are approximately 370 ppmv, 400 ppmv and 418 ppmv, respectively. There is a significant degree of difference between the four reference profiles. Between 80 and 120 km, the three reference CO2 profiles of GRAP, MIPAS and AFGL experience a rapid decrease, whereas the CO2 profile of GRAP remains higher than the other two profiles. Figure 9b presents the comparison of the CH4 profiles, showing that the shape of the GRAP CH4 profile resembles the other three profiles. However, its values are consistently higher from 0 to 120 km, peaking at approximately 1.95 ppmv in the troposphere. The primary focus of CH4 is in the troposphere, and the comparison between the four profiles and the observed profiles from the NDACC ground station reveals the smallest difference between the GRAP profile and the NDACC observations in the troposphere. Figure 9c indicates minimal differences in O3 profile values between GRAP and MIPAS, AFGL, NDACC, with the peak O3 concentrations occurring at 28–32 km. In contrast, the differences between the RTTOV reference atmospheric profiles and the other three reference profiles are more pronounced. Figure 9d illustrates that the N2O profile values of GRAP, RTTOV and NDACC are slightly larger than those of AFGL and MIPAS within the 0–15 km range, while MIPAS values are larger within the 15–45 km range, with N2O concentration reaching approximately zero above 50 km.
The four atmospheric profiles were utilized to calculate their corresponding total column density, and the resulting values are presented in Table 3. According to the WMO Greenhouse Gas Bulletin (No. 18, 2022) [33], published by the World Meteorological Organization (WMO), annual average global atmospheric concentrations of major greenhouse gases reached new highs in 2021. These are 415.7 ± 0.2 ppmv for CO2, 1.908 ± 0.002 ppmv for CH4 and 0.3345 ± 0.0001 ppmv for N2O, which are 149%, 262% and 124% of pre-industrial (pre-1750) levels, respectively. The total column density measurements of atmospheric constituents at the equatorial Ascension Island ground station from the World Data Centre for Greenhouse Gases (WDCGG) were selected for validation, yielding values of 415.53 ppmv for CO2, 1.868 ppmv for CH4, and 0.334 ppmv for N2O. In summary, the atmospheric parameter profile values of GRAP demonstrate closer agreement with the current atmospheric state compared to AFGL, RTTOV and MIPAS.

4.2. Comparison of Reference Northern Hemisphere Mid-Latitude Summer Profiles

This study compares atmospheric profiles for summer in the mid-latitude climatic zone. The selected profiles include reference profiles of GRAP for July located in the US region (40n90w) and China region (40n120e), mid-latitude summer reference profiles of the AFGL, mid-latitude daytime reference profiles of the MIPAS, the RTTOV reference profiles, and NDACC measured profiles for July 2021 at the Boulder ground station site in Boulder, CO, United States.
Figure 10a illustrates the comparison of CO2 profiles, revealing the oscillation of CO2 profile values of GRAP in the 18–75 km altitude range in both the China and US regions, with a maximum difference of 15 ppmv. However, above 75 km, the CO2 values are lower in the US region compared to the China region. Figure 10b presents the comparison of CH4 profiles. Between 0 and 20 km, there is a slight difference between the CH4 profiles of GRAP in the US and China regions. The profile values of GRAP demonstrate closer agreement with the measured values at the NDACC ground station than the other three reference profiles. The CH4 profile values of RTTOV significantly exceed those of the other three atmospheric models at 20–55 km. Above 55 km, the CH4 profile of GRAP remains relatively constant and is notably higher than the other three atmospheric profiles. Figure 10c displays the comparison of O3 profiles, indicating two peaks in the 25–55 km and 85–100 km ranges. The O3 profiles of GRAP in the US and China regions differ by 0.2 ppmv, and the maximum difference among the six O3 profile values is 2 ppmv. Figure 10d presents the comparison of N2O profiles, showing that the six profiles exhibit a similar trend with differences primarily located below 60 km. The N2O profile of GRAP exhibits the least difference from the measured profile at the NDACC ground station.
The total column density values for the six atmospheric constituents are presented in Table 4. According to the China Greenhouse Gas Bulletin No. 11 report [34], observations from the China Meteorological Administration’s Waliguan National Atmospheric Background Station in 2021 indicate annual average atmospheric concentrations of CO2, CH4 and N2O as 417.0 ± 0.2 ppmv, 1.965 ± 0.0006 ppmv, and 0.3351 ± 0.0001 ppmv, respectively. These values are comparable to the same period in the northern hemisphere mid-latitudes, although slightly higher than the global mean. It is important to note that, despite being in the same latitudinal zones, the atmospheric conditions can vary across different longitude zones. In terms of data source, the GRAP is more recent, ensuring that the atmospheric parameter profiles align more closely with current atmospheric conditions when compared to the other three reference profiles.

4.3. Comparison of Reference Polar Winter Profiles

This study compares reference profiles for the polar winter climatic zone. The selected profiles for comparison include the reference atmospheric profiles of the 90n30e grid of GRAP in January, the reference polar winter profiles of the AFGL, the reference polar winter profiles of the MIPAS, the RTTOV reference profiles, and the profiles measured in March 2022 at the NDACC ground station in Ny Ålesund, Norway, which is located in the polar regions.
Figure 11a presents the comparison of the CO2 reference profiles, which exhibit a similar trend to the equatorial and mid-latitude regions. Figure 11b presents the comparison of the CH4 reference profiles. The CH4 profile of GRAP exhibits higher values compared to the other three reference profiles, with a concentration peak at 15 km reaching nearly 2 ppmv. Furthermore, the GRAP values in the troposphere closely align with the measured profile values in the NDACC. Figure 11c presents the comparison of the O3 profiles. The GRAP profile represents the lowest O3 concentration, reaching as low as 4.7 ppmv, while the RTTOV profile exhibits the highest O3 concentration, reaching up to 7.2 ppmv. The O3 profile concentration from AFGL is in closer agreement with the measured profile from NDACC. Figure 11d illustrates the comparison of the N2O profiles. The trends in the five N2O profiles align closely with the equatorial and mid-latitude regions, with the N2O profile values of the GRAP and NDACC measured profiles being highly similar.

4.4. Comparison of Reference Atmospheric Temperature Profiles in China

Table 5 presents the characteristics of the selected Chinese regional reference atmospheric temperature profiles. Temporally, representative months for each season, namely January for winter, April for spring, July for summer, and October for autumn, were chosen in the Chinese region. Spatially, the profiles encompass various regions in China, including the northern, central, southern, northeastern, northwestern and southwestern parts of the country.
Figure 12 presents the comparison of reference temperature profiles for eight regions in China. Analysis of the figure reveals variations in surface temperatures (0 km) among all regions, with a maximum difference of 30 K occurring in winter, a minimum difference of only 15 K in summer, and approximately 20 K differences in both spring and autumn. In the troposphere, temperatures decrease with altitude, exhibiting the smallest variations in summer and a nearly parallel decreasing temperature profile. However, during spring, autumn and winter, a significant inflection point is observed around 10 km in Harbin, Urumqi, Beijing and Shanghai, where the rate of temperature decrease diminishes and even shows a rising trend. In the stratosphere, temperatures increase with altitude, displaying substantial differences between summer and winter. Harbin, Urumqi and Beijing exhibit significantly higher temperatures compared to the other five regions, while spring and autumn temperatures exhibit less disparity. In the mesosphere, temperatures decrease with altitude, reaching values of 160–190 K at the mesosphere’s upper boundary. Finally, in the thermosphere, temperatures rise rapidly with altitude, peaking at 440 K.
The analysis of temperature profile characteristics in various regions of China reveals that the temperature in the troposphere decreases as the latitude of the region increases. During spring, autumn and winter, the troposphere’s upper boundary decreases with latitude in eight regions, reaching approximately 15 km in Harbin, Urumqi and Beijing. In summer, the troposphere’s upper boundary is around 18 km in all seven regions except Urumqi, where it stands at 16 km. The studies by Qi Chengli reveal significant seasonal and regional variations in the vertical distribution of temperature profiles in China [15,16]. Comparing the reference temperature profiles of GRAP for China with Qi Chengli’s studies and the actual situation in China yields similar results, suggesting that latitudinal variations can lead to substantial disparities in tropospheric temperatures. Therefore, further subdivision of the atmospheric reference profile based on latitude holds considerable significance.

5. Conclusions

Most reference atmospheric profile databases commonly used have limited spatial and temporal resolution. These classification criteria fail to meet the research requirements for global regions with intricate topography and diverse climates. To address this, this study constructs the GRAP. The random forest regression model and the stratified mean method are adopted to process the ACE-FTS L2 products, the AIRS L2 products and the ERA5 reanalysis data. The data is divided into monthly intervals from January to December and spatially organized into 532 grids with dimensions of 5° × 30°, creating a comprehensive global atmospheric profile reference database. Leading to the following conclusions:
(1)
GRAP provides extensive coverage on a global scale, presenting a comprehensive composition of the atmosphere that accurately presents its current state.
(2)
Four atmospheric reference profile databases were used as input parameters for the RFM radiative transfer model to simulate the FY-3E HIRAS-II absorption spectrum and compare it with the measured spectrum. The results illustrate that the spectral simulation of GRAP as an input parameter is a better fit for the measured HIRAS-II spectrum.
(3)
The atmospheric profiles of the four reference atmospheric profile databases were compared to the measured atmospheric profiles from NDACC and the column total concentrations measured by WDCGG. The findings indicate substantial updates in the gas components of GRAP compared to the other three databases. Notably, the four greenhouse gases (CO2, CH4, O3 and N2O) of GRAP demonstrate better alignment with the current atmospheric conditions.
(4)
Comparing the reference temperature profiles of GRAP for eight distinct regions of China reveals that these profiles effectively capture the climatic conditions. The fine spatial and temporal grids enable GRAP to achieve superior regional representativeness compared to previous reference atmospheric profile databases. Consequently, GRAP exhibits enhanced regional representation capabilities.

Author Contributions

Conceptualization, X.L.; Methodology, Y.G.; Software, Y.G.; Formal analysis, Y.G.; Investigation, X.Z., W.L. and W.F.; Data curation, Y.G.; Writing—original draft, Y.G.; Writing—review & editing, Y.G., X.L., T.C. and S.L.; Supervision, X.L.; Project administration, X.L.; Funding acquisition, X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China under grant number 2022YFF0606400, the Basic Strengthening Program Technical Field Foundation, grant number 2023KJC-Y-0193, and Unit Commissioned Projects, grant number 2023KJC-Y-0032.

Data Availability Statement

Not applicable.

Acknowledgments

The AIRS Support Level 2 Version 7 products used in this research were downloaded from the NASA Goddard Space Flight Center Earth Science Data and Information Service, the ACE-FTS Level 2 Version 4.1 products were downloaded from the Atmospheric Chemistry Experiment Science Operations Centre, Department of Chemistry, University of Waterloo, the ERA5 reanalysis data were downloaded from the Copernicus Climate Change Service Climate Data Repository, and HIRAS-II L1 data from the National Centre for Meteorological Satellites (NCMMS) Feng Yun Satellite Remote Sensing Data Service Network, the NDACC ground station profiles data used in this study were downloaded from the NDACC Rapid Delivery (RD) Data Access, for which we would like to express our sincere thanks!

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Shi, G.Y.; Dai, T.; Xu, N. Recent progress in the study of atmospheric CO2 concentration by satellite remote sensing. Adv. Earth Sci. 2010, 25, 7–13. [Google Scholar]
  2. Li, X.Y.; Xu, J.; Cheng, T.H.; Shi, H.L.; Zhang, X.Y.; Ge, S.L.; Wang, H.M.; Zhu, S.Y.; Miao, J.; Luo, Q. Monitoring Trace Gases over the Antarctic Using Atmospheric Infrared Ultraspectral Sounder Onboard GaoFen-5: Algorithm Description and First Retrieval Results of O3, H2O, and HCl. Remote Sens. 2019, 11, 1991. [Google Scholar] [CrossRef] [Green Version]
  3. Wang, H.M.; Li, X.Y.; Chen, L.F.; Wang, Y.P.; Zhang, Y.; Zou, M.M.; Yu, C.; Zhu, S.Y. Sensitivity analysis of terahertz proximity detection of O3 and HCl. J. Remote Sens. 2017, 21, 653–664. [Google Scholar]
  4. Bauer, P.; Auligne, T.; Bell, W.; Geer, A.; Guidard, V.; Heilliette, S.; Kazumori, M.; Kim, M.; Liu, E.; McNally, A.; et al. Satellite Cloud and Precipitation Assimilation at Operational NWP Centres. Q. J. R. Meteorol. Soc. 2011, 137, 1934–1951. [Google Scholar] [CrossRef] [Green Version]
  5. Wei, H.L.; Dai, C.M. Research of atmospheric transfer correction in radiance measurement: Atmospheric radiative transfer model and the analysis of key atmospheric parameters. Infrared Laser Eng. 2014, 43, 884–890. [Google Scholar] [CrossRef]
  6. Chevallier, F.; Chedin, A.; Cheruy, F.; Morcrette, J. TIGR-like Atmospheric-Profile Databases for Accurate Radiative-Flux Computation. Q. J. R. Meteorol. Soc. 2000, 126, 777–785. [Google Scholar] [CrossRef]
  7. Simmons, A.J.; Gibson, J.K. ERA-40 Project Report Series No. 1; European Centre for Medium Range Weather Forecasts: Reading, UK, 2000. [Google Scholar]
  8. Matricardi, M.; Saunder, R.W. A Fast Radiative Transfer Model for Simulation of IASI Radiances; Technical Memoranda; ECMWF: Reading, UK, 2000; 23p. [Google Scholar] [CrossRef]
  9. Tjemkes, S.; Schmetz, J. Synthetic Satellite Radiances Using the Radiance Sampling Method. J. Geophys. Res. Atmos. 1997, 102, 1807–1818. [Google Scholar] [CrossRef]
  10. Chevallier, F. Sampled Database of 60-Level Atmospheric Profiles from the ECMWF Analyses; ECMWF: Reading, UK, 2002. [Google Scholar]
  11. Saunders, R.W.; Matricardi, M.; Brunel, P. A Fast Radiative Transfer Model for Assimilation of Satellite Radiance Observations—RTTOV-5; ECMWF: Reading, UK, 1999. [Google Scholar]
  12. Matricardi, M.; Chevallier, F.; Kelly, G.; Thepaut, J. An Improved General Fast Radiative Transfer Model for the Assimilation of Radiance Observations. Q. J. R. Meteorol. Soc. 2004, 130, 153–173. [Google Scholar] [CrossRef]
  13. Chevallier, F. TIGR-Like Sampled Databases of Atmospheric Profiles from the ECMWF 50-Level Forecast Model; European Centre for Medium Range Weather Forecasts: Reading, UK, 1999. [Google Scholar]
  14. Eyre, J. A Fast Radiative Transfer Model for Satellite Sounding Systems; Tech. Memo 176; ECMWF: Reading, UK, 1991. [Google Scholar]
  15. Qi, C.L.; Dong, C.H.; Zhang, W.J. Characterization of a typical regional atmospheric sample base in China based on radiative transfer model. Meteorol. Sci. Technol. 2004, 32, 357–362. [Google Scholar] [CrossRef]
  16. Qi, C.L.; Liu, H.; Ma, G.; Zhang, P.; Wu, X.B. A selection method for a sample library of typical atmospheric profiles in the Chinese region. J. Appl. Meteorol. 2010, 21, 70–75. [Google Scholar] [CrossRef]
  17. Kneizys, F.X. Users Guide to LOWTRAN 7; Air Force Geophysics Laboratory: Hanscom AFB, MA, USA, 1988; Volume 88. [Google Scholar]
  18. Berk, A.; Anderson, G.P.; Bernstein, L.S.; Acharya, P.K.; Dothe, H.; Matthew, M.W.; Adler-Golden, S.M.; Chetwynd, J.H., Jr.; Richtsmeier, S.C.; Pukall, B. MODTRAN4 Radiative Transfer Modeling for Atmospheric Correction; SPIE: Bellingham, WA, USA, 1999; Volume 3756, pp. 348–353. [Google Scholar]
  19. Turner, D.; Tobin, D.; Clough, S.; Brown, P.; Ellingson, R.; Mlawer, E.; Knuteson, R.; Revercomb, H.; Shippert, T.; Smith, W.; et al. The QME AERI LBLRTM: A Closure Experiment for Downwelling High Spectral Resolution Infrared Radiance. J. Atmos. Sci. 2004, 61, 2657–2675. [Google Scholar] [CrossRef]
  20. Isaacs, R.G.; Wang, W.C.; Worsham, R.D.; Goldenberg, S. Multiple scattering Lowtran and Fascode models. Appl. Opt. 1987, 26, 1272–1281. [Google Scholar] [CrossRef] [PubMed]
  21. Dudhia, A. The Reference Forward Model (RFM). J. Quant. Spectrosc. Radiat. Transf. 2017, 186, 243–253. [Google Scholar] [CrossRef] [Green Version]
  22. Anderson, G.P.; Chetwynd, J.H.; She, E.P. AFGL Atmospheric Constituent Profiles (0–120 km); Air Force Geophysics Laboratory: Hanscom AFB, MA, USA, 1986. [Google Scholar]
  23. NOAA-S/T 76-1526; National Aeronautics and Space Administration; United States Air Force; United States Weather Bureau. U.S. Standard At-mosphere (1976). U.S. Government Publishing Office: Washington, DC, USA, 1976.
  24. Boone, C.D.; Bernath, P.F.; Cok, D.; Jones, S.C.; Steffen, J. Version 4 Retrievals for the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) and Imagers. J. Quant. Spectrosc. Radiat. Transf. 2020, 247, 106939. [Google Scholar] [CrossRef]
  25. Pagano, T.; Chahine, M.; Fetzer, E. The Atmospheric Infrared Sounder (AIRS) on the NASA Aqua Spacecraft: A General Remote Sensing Tool for Understanding Atmospheric Structure, Dynamics and Composition. In Remote Sensing of Clouds and the Atmosphere XV; Picard, R., Schafer, K., Comeron, A., Eds.; SPIE: Bellingham, WA, USA, 2010; Volume 7827. [Google Scholar]
  26. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. The ERA5 Global Reanalysis. Q. J. R. Meteorol. Soc. 2020, 146, 1999–2049. [Google Scholar] [CrossRef]
  27. Fu, X.S.; Tan, J.G. Quality control methods for ground-based microwave radiometer sounding data. J. Appl. Meteorol. 2017, 28, 209–217. [Google Scholar] [CrossRef]
  28. Li, H.T.; Shao, Z.D. A review of spatial interpolation analysis algorithms. Comput. Syst. Appl. 2019, 28, 1–8. [Google Scholar]
  29. Sofieva, V.; Rahpoe, N.; Tamminen, J.; Kyrola, E.; Kalakoski, N.; Weber, M.; Rozanov, A.; von Savigny, C.; Laeng, A.; von Clarmann, T.; et al. Harmonized Dataset of Ozone Profiles from Satellite Limb and Occultation Measurements. Earth Syst. Sci. Data 2013, 5, 349–363. [Google Scholar] [CrossRef] [Green Version]
  30. Segal, M.R. Machine Learning Benchmarks and Random Forest Regression; UCSF Center for Bioinformatics and Molecular Biostatistics: San Francisco, CA, USA, 2004; Available online: https://escholarship.org/uc/item/35x3v9t4 (accessed on 5 April 2023).
  31. Remedios, J.J.; Leigh, R.J.; Waterfall, A.M.; Moore, D.P.; Sembhi, H.; Parkes, I.; Greenhough, J.; Chipperfield, M.P.; Hauglustaine, D. MIPAS reference atmospheres and comparisons to V4.61/V4.62 MIPAS level 2 geophysical data sets. Atmos. Chem. Phys. Discuss. 2007, 7, 9973–10017. [Google Scholar] [CrossRef] [Green Version]
  32. Saunders, R.; Hocking, J.; Rundle, D.; Rayer, P.; Havemann, S.; Matricardi, M.; Geer, A.; Lupu, C.; Brunel, P.; Vidot, J. RTTOV-12 Science and Validation Report; EUMETSAT: Darmstadt, Germany, 2017. [Google Scholar]
  33. World Meteorological Organization (WMO). The WMO Greenhouse Gas Bulletin (2022) No. 18. Available online: https://library.wmo.int/index.php?lvl=notice_display&id=22149 (accessed on 11 April 2023).
  34. China Meteorological Administration. The China Greenhouse Gas Bulletin No. 11 Report. Available online: https://www.cma.gov.cn/zfxxgk/gknr/qxbg/202301/t20230119_5274988.html (accessed on 11 April 2023).
Figure 1. Flow chart of the methods of creating GRAP.
Figure 1. Flow chart of the methods of creating GRAP.
Remotesensing 15 03006 g001
Figure 2. Sample numbers of CH4 profiles in the global from January to December.
Figure 2. Sample numbers of CH4 profiles in the global from January to December.
Remotesensing 15 03006 g002aRemotesensing 15 03006 g002b
Figure 3. Random forest regression model construction process.
Figure 3. Random forest regression model construction process.
Remotesensing 15 03006 g003
Figure 4. Comparison of the simulated spectrum and the measured spectrum of HIRAS-II in the 650–2550 cm−1 band (full absorption band).
Figure 4. Comparison of the simulated spectrum and the measured spectrum of HIRAS-II in the 650–2550 cm−1 band (full absorption band).
Remotesensing 15 03006 g004
Figure 5. Comparison of the simulated spectrum and the measured spectrum of HIRAS-II in the 650–760 cm−1 band (CO2 absorption band).
Figure 5. Comparison of the simulated spectrum and the measured spectrum of HIRAS-II in the 650–760 cm−1 band (CO2 absorption band).
Remotesensing 15 03006 g005
Figure 6. Comparison of the simulated spectrum and the measured spectrum of HIRAS-II in the 1200–1400 cm−1 band (CH4 absorption band).
Figure 6. Comparison of the simulated spectrum and the measured spectrum of HIRAS-II in the 1200–1400 cm−1 band (CH4 absorption band).
Remotesensing 15 03006 g006
Figure 7. Deviation of the simulated and measured spectra of HIRAS-II in the 650–760 cm−1 band.
Figure 7. Deviation of the simulated and measured spectra of HIRAS-II in the 650–760 cm−1 band.
Remotesensing 15 03006 g007
Figure 8. Deviation of the simulated and measured spectra of HIRAS-II in the 1200–1400 cm−1 band.
Figure 8. Deviation of the simulated and measured spectra of HIRAS-II in the 1200–1400 cm−1 band.
Remotesensing 15 03006 g008
Figure 9. Comparison of standard equatorial profiles.
Figure 9. Comparison of standard equatorial profiles.
Remotesensing 15 03006 g009
Figure 10. Comparison of standard summer profiles in the northern hemisphere at mid-latitudes.
Figure 10. Comparison of standard summer profiles in the northern hemisphere at mid-latitudes.
Remotesensing 15 03006 g010
Figure 11. Comparison of standard polar winter profiles.
Figure 11. Comparison of standard polar winter profiles.
Remotesensing 15 03006 g011
Figure 12. Comparison of regional atmospheric temperature standard profiles in China.
Figure 12. Comparison of regional atmospheric temperature standard profiles in China.
Remotesensing 15 03006 g012aRemotesensing 15 03006 g012b
Table 1. Metadata information for the GRAP.
Table 1. Metadata information for the GRAP.
KeywordsDetailed Information
NameGlobal Reference Atmospheric Profile Database, GRAP
Data storage format.txt
Time resolutionJanuary–December, one-month interval
(MON1, MON2, MON3, …, MON11, MON12)
Spatial resolutionLatitude zones: 5° interval. n for north, s for south.
(0n, 0s, 90n, 90s latitude zones are 2.5°)
Longitude zones: 30° interval. e for eastern, w for western.
(0e, 0w, 180e, 180w longitude zones are 15°)
Height level0–119 km at 1km intervals, divided into 120 level
Atmospheric parametersPressure, Temperature, CO2, CH4, O2, NO, N2O, O3, SO2, NH3, SF6, CO, N2, HF, HBr, CF4, NO2, HI, OCS, H2CO, H2O2, C2H2, C2H6, PH3, COF2, H2S, CFCl3, CF2Cl2, CClF3, CHCl2F, OH, CHClF2, ClO, C2Cl3F3, C2Cl2F4, C2ClF5, CCl4, ClONO2, N2O5, HNO4, BrO, CH3Cl, CH3CN, CH3OH, H2O, HCl, HCN, HNO3, HCOOH, HOCl, CCl2F2, HO2, CCl3F, COCl2, COClF, pan(CH3C(O)OONO2), CHF3, HO2NO2, HCFC141b, HCFC142b, HFC134a
Table 2. Data sources and GRAP height range level.
Table 2. Data sources and GRAP height range level.
Data SourcesHeight RangePressure RangeNumber of Layers
ACE-FTS0.5–149.5 km1013–3.22 × 10−6 hpa150
AIRS-1100–1.61 × 10−6 hpa100
ERA5-1000–1 hPa37
GRAP0–119 kmDifferences between grids120
Table 3. Equatorial profiles correspond to total column density (units are ppmv).
Table 3. Equatorial profiles correspond to total column density (units are ppmv).
Atmospheric CompositionGRAPRTTOVMIPASAFGLNDACCWDCGGWMO
CO2415.25401.54368.03330-415.53415.7
CH41.8751.8071.7441.6481.9131.8681.908
O30.3690.3010.3040.3460.374--
N2O0.3230.3190.3110.3070.3260.3340.3345
Table 4. Northern Hemisphere mid-latitude summer profiles correspond to total column density (units are ppmv).
Table 4. Northern Hemisphere mid-latitude summer profiles correspond to total column density (units are ppmv).
Atmospheric CompositionGRAP_40n120eGRAP_40n90wRTTOVMIPASAFGLNDACCCGGB
CO2415.15415.46401.54368.03330-417
CH41.8781.8591.8071.7261.6481.8811.965
O30.4510.3780.3010.3750.3920.469-
N2O0.3280.3270.3190.3040.2970.3220.3351
Table 5. Location of the Regional Reference Atmospheric Temperature Profile in China.
Table 5. Location of the Regional Reference Atmospheric Temperature Profile in China.
Profile NumberLocation NameLatitude and LongitudeGrid Location
1Harbin45.75N, 126.63E45n120e
2Urumqi43.77N, 87.68E45n90e
3Beijing39.92N, 116.42E40n120e
4Shanghai34.50N, 121.4335n120e
5Lhasa29.60N, 91.00E30n90e
6Kunming25.05N, 102.73E25n90e
7Guangzhou23.13N, 113.27E25n120e
8Sanya18.25N, 109.5E20n120e
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, Y.; Li, X.; Cheng, T.; Li, S.; Zhang, X.; Lu, W.; Fang, W. Construction of the Global Reference Atmospheric Profile Database. Remote Sens. 2023, 15, 3006. https://doi.org/10.3390/rs15123006

AMA Style

Guo Y, Li X, Cheng T, Li S, Zhang X, Lu W, Fang W. Construction of the Global Reference Atmospheric Profile Database. Remote Sensing. 2023; 15(12):3006. https://doi.org/10.3390/rs15123006

Chicago/Turabian Style

Guo, Yuhang, Xiaoying Li, Tianhai Cheng, Shenshen Li, Xinyuan Zhang, Wenjing Lu, and Weifang Fang. 2023. "Construction of the Global Reference Atmospheric Profile Database" Remote Sensing 15, no. 12: 3006. https://doi.org/10.3390/rs15123006

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop