Investigation on Global Distribution of the Atmospheric Trapping Layer by Using Radio Occultation Dataset

: The trapping layer refers to the atmospheric layer with vertical gradient of atmospheric refractivity less than − 157 N-Units/km or vertical gradient of atmospheric modiﬁed refractivity 0 M-unit/km, which has a signiﬁcant impact on radar and radio communication systems. Based on COSMIC and other radio occultation data, we show the statistical characteristics of the global trapping layer during 2005–2020.The statistical results show that the occurrence rate of the trapping layers is mainly concentrated between 50 ◦ S and 50 ◦ N, and higher occurrences of the trapping layers with more than 50% mainly occur in the boundary area between ocean and land, such as the northwest coastal area of Mexico, the west coastal area of Africa, the Mediterranean Sea, the Red Sea and the Arabian Sea, and the northwest area of Australia, etc. The altitude of the trapping layer is lower near the land and increases with the distance away from the coastline. The intensity is mainly between 6 M-unit and 24 M-unit (an M-unit is the unit of atmospheric modiﬁed refractivity), and the average value in some regions is above 24 M-unit, such as in the Arabian Sea area. In addition, the thickness of the trapping layer is between 50 and 240 m, and is generally larger over the ocean than over the land. These results reveal that the generation


Introduction
Atmospheric super-refraction phenomena associated with sharp gradients in temperature and water vapor pressure frequently occur in the near-Earth atmospheric boundary layer, especially near coastal and ocean areas [1,2]. Ao et al. [3] indicated that there is an atmospheric trapping layer when the vertical gradient of atmospheric refractivity is dN/dh < −157 N-Units/km or vertical gradient of atmospheric modified refractivity is dM/dh < 0 M-unit/km, which significantly affects the performance of radio communication systems by changing the propagation path and attenuation characteristics of electromagnetic waves. The trapping layer is always embedded inside the duct, which produces information on the sharp layer and is related to the physical processes or conditions that form the duct [3,4]. The generation and distribution of the atmospheric trapping layer are closely related to the geography and meteorological water level conditions such as nocturnal radiation cooling, warm and dry land air passing over a cooler sea, sea-land breezes, monsoon and so on [5,6]. Studying the characteristics of the trapping layer and its influence on the propagation of electromagnetic waves is helpful to accurately assess the influence of the complex duct environment on the radio information system, and provide environmental analysis information and auxiliary decision-making methods for system design and development.
Statistical studies of the occurrence of the atmospheric trapping layer in different regions have been carried out in the last few decades [7][8][9][10]. By using 6-year European Centre for Medium-Range Weather Forecasts (ECMWF) data, Engeln [8] calculated the atmospheric refractivity and analyzed the seasonal variations in the occurrence rate, altitude, thickness, and gradient of the trapping layer at a height of 0~2.5 km. The observed results shown in Engeln [8] show that the distribution of the trapping layer occurrence has obvious seasonal variations, which are mainly affected by stratocumulus, fog, monsoon and boundary layer inversions in temperature and moisture [11][12][13][14]. However, ECMWF data have some limitations. The resolution of ECMWF data decreases with altitude, making trapping layer information from above 2.5 km is almost useless. The vertical resolution of ECMWF data is greater than 200 m where the height is over 1 km, and it can only detect some thicker trapping layers.
Atmospheric parameter profiles, such as atmospheric refractivity, temperature, pressure and humidity, can be measured by using radio occultation technology [15]. When the MicroLab-1 satellite was launched in 1995, Global Positioning System (GPS) radio occultation measurements were first applied to the study of global atmospheric distribution characteristics [16]. Subsequently, CHAMP [17] and SAC-C [18] equipped with GPS radio occultation receivers were successfully launched. Following the launch of COSMIC on 14 April 2006, approximately 2000 to 3000 atmospheric radio occultation profiles can now be obtained every day [19,20]. The large amount of COSMIC radio occultation data available is beneficial for studying the global distributions of the trapping layer occurrence rate. The vertical resolution of the COSMIC data is 100 m, which allows researchers to observe some of the finer structures of the trapping layer, which is helpful for investigating the generation mechanism of the trapping layer. It is noted that the number of RO profiles has declined greatly recently. Due to its significant advantages of high precision, high vertical resolution, and ability to perform all-weather observations [21][22][23], we firstly investigate the temporal and spatial distribution characteristics of the global trapping layer based on the GPS radio occultation dataset. COSMIC and other radio occultation data are used in this study. The data source and research methods are introduced in Section 2. In Section 3, we analyze the results of the temporal and spatial distribution of the trapping layer characteristics. Our conclusions are summarized in Section 4.

Method Description
Radio waves are affected by atmospheric refraction during propagation, and refraction effects must be considered in the actual propagation process. The refractivity is determined by pressure, temperature and humidity [24]: where, p is atmospheric pressure (hPa), T is temperature (K), and e is water vapor pressure (hPa). The modified refractivity (M), corrected by the curvature of the Earth, is further given by Bean and Dutton [13].
where, R e is the average earth radius (6371 km), h is the altitude (m). The atmospheric refraction types are classified as subrefraction, standard refraction, super-refraction, and trapping layer (duct) according to the gradient of refractivity N or the gradient of modified refractivity M [3,13,25]. The relationship between N gradient and M gradient is shown in Table 1.
Trapping layer is defined as the layer with a gradient of refractivity dM/dh less than 0. In Figure 1, h m < h < h t indicates trapping layer, and h b < h < h t indicates elevated duct [3]. The atmospheric refraction types are classified as subrefraction, standard refraction, super-refraction, and trapping layer (duct) according to the gradient of refractivity N or the gradient of modified refractivity M [3,13,25]. The relationship between N gradient and M gradient is shown in Table 1. Trapping layer is defined as the layer with a gradient of refractivity dM dh ⁄ less than 0. In Figure 1, h < h < h indicates trapping layer, and h < h < h indicates elevated duct [3]. The steps for finding the trapping layer from each atmospheric refractivity profile are as follows: 1. Since the minimum effective height of each atmospheric profile is different, in order to ensure that the selected profile contains the effective information of the boundary layer and a sufficient number of samples, the lowest effective record height of the lower bound of the selected data is within 1 km from the ground. 2. Calculate the double-weighted average and standard deviation [26] of all refractivity profiles, and eliminate the profile that is greater than (less than) the average plus (minus) 5 times the standard deviation. 3. Using cubic spline interpolation to interpolate the profile to 5 km at the vertical height at intervals of 10 m. The formula of a cubic spline interpolation is where Ci = ai + bix + cix 2 + dix 3 (di ≠ 0), i = 1, …, n [8]. The error analysis of cubic spline interpolation has been presented in Figure 2. The comparative result indicates that the height distribution of refractivity is not affected by cubic spline interpolation. The interpolation results are shown in Figure 2. The red dotted line represents the original data value, and the solid blue line represents the interpolated data values. The results indicate that cubic spline interpolation has no effect on the data. The higher height resolution of The steps for finding the trapping layer from each atmospheric refractivity profile are as follows: 1.
Since the minimum effective height of each atmospheric profile is different, in order to ensure that the selected profile contains the effective information of the boundary layer and a sufficient number of samples, the lowest effective record height of the lower bound of the selected data is within 1 km from the ground.

2.
Calculate the double-weighted average and standard deviation [26] of all refractivity profiles, and eliminate the profile that is greater than (less than) the average plus (minus) 5 times the standard deviation.

3.
Using cubic spline interpolation to interpolate the profile to 5 km at the vertical height at intervals of 10 m. The formula of a cubic spline interpolation is where C i = a i + b i x + c i x 2 + d i x 3 (d i = 0), i = 1, . . . , n [8]. The error analysis of cubic spline interpolation has been presented in Figure 2. The comparative result indicates that the height distribution of refractivity is not affected by cubic spline interpolation. The interpolation results are shown in Figure 2. The red dotted line represents the original data value, and the solid blue line represents the interpolated data values.
The results indicate that cubic spline interpolation has no effect on the data. The higher height resolution of refractivity data can observe some of the finer structures of trapping layer, which is helpful for investigating the generation mechanism of trapping layer. 4.
Using the central difference method to calculate the vertical gradient of the refractivity at each height: where i the number of profile layers.

5.
Checking the value of M (i) from the bottom up. When M (i) < 0, we find the bottom of the trapping layer h m . When M (i) > 0, we find the top of the trapping layer h t . Here, only select the first trapping layer retrieved. In Figure 1, the height of the trapping layer is h t , the thickness of the trapping layer is h t − h m , and the intensity of the trapping layer is a m − a s . indicate that cubic spline interpolation has no effect on the data. The higher height resolution of refractivity data can observe some of the finer structures of trapping layer, which is helpful for investigating the generation mechanism of trapping layer. 4. Using the central difference method to calculate the vertical gradient of the refractivity at each height: where i the number of profile layers.

5.
Checking the value of M ′ (i) from the bottom up. When M ′ (i) < 0, we find the bottom of the trapping layer h m . When M ′ (i) > 0, we find the top of the trapping layer h t . Here, only select the first trapping layer retrieved. In Figure 1, the height of the trapping layer is h t , the thickness of the trapping layer is h t − h m , and the intensity of the trapping layer is a m − a s .
The processed refractivity data are binned within a 2 ∘ × 2 ∘ grid and are divided into four seasons according to the Northern Hemisphere, including spring (March, April, and May), summer (June, July, and August), autumn (September, October, and November), winter (December, January, and February), and four local times of LT00, LT06, LT12, LT18. We have statistically analyzed the temporal and spatial distribution of the global trapping layer occurrence probability, height, intensity, and thickness by using radio occultation data during 2005-2020.

Data Description
Radio occultation measurement has certain limitations. Some studies have shown that the atmospheric refractivity detected by radio occultation usually has a negative bias in the lower troposphere [27,28,29,30]. It is believed that tracking error is the main reason for refractivity deviation [31,32]. Another possible reason is the presence of the super- The processed refractivity data are binned within a 2 • × 2 • grid and are divided into four seasons according to the Northern Hemisphere, including spring (March, April, and May), summer (June, July, and August), autumn (September, October, and November), winter (December, January, and February), and four local times of LT00, LT06, LT12, LT18. We have statistically analyzed the temporal and spatial distribution of the global trapping layer occurrence probability, height, intensity, and thickness by using radio occultation data during 2005-2020.

Data Description
Radio occultation measurement has certain limitations. Some studies have shown that the atmospheric refractivity detected by radio occultation usually has a negative bias in the lower troposphere [27][28][29][30]. It is believed that tracking error is the main reason for refractivity deviation [31,32]. Another possible reason is the presence of the super-refractive layer that affects the non-uniformity of the Abelian inversion [28,33]. The successful application of Open Loop (OL) tracking technology can effectively reduce tracking errors [27,29,33], improve the detection ability of the bottom of the troposphere and the tracking ability of occultation events, and ameliorate the negative biases [34]. The Open Loop tracking technology was first verified on SAC-C and CHAMP satellites, and currently, it has been applied to radio occultation receivers such as COSMIC and METOP-A. The data used in this study is the level 2 wetPrf atmospheric refractivity profile data provided by UCAR's COSMIC Data Analysis and Archive Center (CDAAC), including temperature, air pressure, water vapor pressure, and refractivity (https://www.cosmic. ucar.edu/, accessed on 20 July 2020). Detailed data information is shown in Table 2.

Probability of Trapping Occurrence
Here, we define the occurrence rate of the trapping layer as the number of the trapping layer occurrences divided by the total number of occultation events in each grid. The number of the radio occultation data and the trapping layer occurrence within the 2_deg by 2_deg grid can be seen from Figures 3 and 4. It is found that the occultation events can cover the world with a relative latitudinal dependence, as shown in Figure 3, and the trapping layer occurrence mainly focuses on the coastal regions, as shown in Figure 4. It should be noted that the latitudinal bands containing more data near 25 and 50 • in both the Northern Hemisphere and Southern Hemisphere in Figure 3 are due to the orbital configuration of the COSMIC satellites. Figure 5 shows the seasonal distribution of the occurrence rate of the trapping layer, where the white area indicates that there is no occultation event. It can be seen in Figure 5 that the regions with high occurrence rates are concentrated between 50 • S and 50 • N (except the Antarctic region), and are mainly found in the eastern areas of oceans and the western areas of landmasses. In addition, at the junction of continent and ocean, the occurrence rate of the trapping layer is higher, and decreases with increasing distance from the coast. ability of occultation events, and ameliorate the negative biases [34]. The Open Loop tracking technology was first verified on SAC-C and CHAMP satellites, and currently, it has been applied to radio occultation receivers such as COSMIC and METOP-A.
The data used in this study is the level 2 wetPrf atmospheric refractivity profile data provided by UCAR's COSMIC Data Analysis and Archive Center (CDAAC), including temperature, air pressure, water vapor pressure, and refractivity (https://www.cosmic.ucar.edu/ , accessed on 20 July 2020). Detailed data information is shown in Table 2.

Probability of Trapping Occurrence
Here, we define the occurrence rate of the trapping layer as the number of the trapping layer occurrences divided by the total number of occultation events in each grid. The number of the radio occultation data and the trapping layer occurrence within the 2_deg by 2_deg grid can be seen from Figures 3 and 4. It is found that the occultation events can cover the world with a relative latitudinal dependence, as shown in Figure 3, and the trapping layer occurrence mainly focuses on the coastal regions, as shown in Figure 4. It should be noted that the latitudinal bands containing more data near 25 and 50° in both the Northern Hemisphere and Southern Hemisphere in Figure 3 are due to the orbital configuration of the COSMIC satellites.    Figure 5 shows the seasonal distribution of the occurrence rate of the trapping layer, where the white area indicates that there is no occultation event. It can be seen in Figure 5 that the regions with high occurrence rates are concentrated between 50° S and 50° N (except the Antarctic region), and are mainly found in the eastern areas of oceans and the western areas of landmasses. In addition, at the junction of continent and ocean, the occurrence rate of the trapping layer is higher, and decreases with increasing distance from the coast. Seasonal variations in the trapping layer are also obvious. In the spring, the trapping layer occurrence is higher in the Northern Hemisphere than in the Southern Hemisphere. The reason may be that the average sea temperature in the trade wind belts of the eastern ocean and the western continent is significantly lower than that of the land in the Northern Hemisphere in the spring. Trade winds cause convection of the warm and dry air on the land to the cold and wet sea surface, forming a strong inversion stratification [35].   Figure 5 shows the seasonal distribution of the occurrence rate of the trapping layer, where the white area indicates that there is no occultation event. It can be seen in Figure 5 that the regions with high occurrence rates are concentrated between 50° S and 50° N (except the Antarctic region), and are mainly found in the eastern areas of oceans and the western areas of landmasses. In addition, at the junction of continent and ocean, the occurrence rate of the trapping layer is higher, and decreases with increasing distance from the coast. Seasonal variations in the trapping layer are also obvious. In the spring, the trapping layer occurrence is higher in the Northern Hemisphere than in the Southern Hemisphere. The reason may be that the average sea temperature in the trade wind belts of the eastern ocean and the western continent is significantly lower than that of the land in the Northern Hemisphere in the spring. Trade winds cause convection of the warm and dry air on the land to the cold and wet sea surface, forming a strong inversion stratification [35]. Seasonal variations in the trapping layer are also obvious. In the spring, the trapping layer occurrence is higher in the Northern Hemisphere than in the Southern Hemisphere. The reason may be that the average sea temperature in the trade wind belts of the eastern ocean and the western continent is significantly lower than that of the land in the Northern Hemisphere in the spring. Trade winds cause convection of the warm and dry air on the land to the cold and wet sea surface, forming a strong inversion stratification [35].
The occurrence of trapping layers can exceed 50% in western Mexico, western North Africa, western South Africa, the Arabian Sea and the Gulf of Bengal in the spring. The areas of high trapping layer occurrence are consistent with those found in Engeln's research [8], but the occurrence of the trapping layer is nearly 100% in Engeln's work. Ao [3] pointed out that the trapping layer occurrence analyzed by the ECMWF was either very high or very low, which is inconsistent with the radiosonde analysis. The trapping layers are mainly affected by stratocumulus clouds, wind and cold ocean currents [8,36]. The trapping layers on the west coasts of America and Africa roughly correspond to the regions where stratocumulus clouds are prevalent [11]. Stratocumulus clouds can cause a sharp gradient of relative humidity and temperature, which is conducive to the generation of the trapping layer [11]. In addition, cold ocean currents are also associated with the trapping layers of these regions. The north-south coverage of the high trapping layer occurrence in the western Mexico sea area is roughly the same as the north-south coverage of the California cold current [37]. The western North African Ocean is affected by the Canary cold current [37], and the high occurrence of the trapping layer occurs in the range of 5 • N to 25 • N, 35 • W of western coastline of the African continent. Its north-south coverage is roughly the same as that covered by the Canary cold current. The western South African Ocean is also affected by Benguela cold current and the southeast trade wind [8,38]. When the warm-dry trade wind flows through the cold sea surface, an inversion layer is formed. The humidity in the layer drops sharply. The coverage areas of the trapping layer reach their maximum point in Arabian Sea and the Gulf of Bengal in the spring. In particular, the trapping layer occurrence rate of the Arabian Sea is can exceed 70%. It is closely related to the monsoon from the African continent.
In summer, affected by the seasonal shift of the equatorial convergence zone, the areas where there is high occurrence of the trapping layer shift northward to the Mediterranean, the Red Sea and other regions. The hot African wind and the cold and wet waters of the Mediterranean Sea and the Red Sea create favorable conditions for the emergence of the trapping layer [7]. At this time, due to the small air-sea temperature difference and the prevailing southwest monsoon in the Arabian Sea and the Bay of Bengal, the area of the trapping layer decreases sharply, as can been see in Engeln's paper. The occurrence of the trapping layer in the northwestern Australia is much higher in the autumn and winter than in the spring and summer. In the autumn and winter, affected by the southeast monsoon, a large area of the trapping layer appears in northwestern Australia. Figure 6 shows the distribution of the trapping layer's occurrence with local time variations. In general, the daily changes in the trapping layer are not obvious, and the locations with higher occurrence rates are basically fixed, mainly being concentrated on the west coast of Africa, the east side of the Atlantic Ocean, the Red Sea, the Mediterranean, the Arabian Sea, the west coast of the Americas, and the northwestern part of Australia. The trapping layer located in the Arabian Sea has obvious diurnal changes. At LT00, due to radiative cooling and the existence of a sea-land breeze [35], the trapping layer occurrence is significantly higher than that at other times.

Altitude of Trapping
The trapping layer altitude ℎ , has an important influence on the propagation of radio waves. Taking the average value of the trapping layer altitude in each grid, the seasonal distribution of the trapping layer altitude is shown in Figure 7.

Altitude of Trapping
The trapping layer altitude h t , has an important influence on the propagation of radio waves. Taking the average value of the trapping layer altitude in each grid, the seasonal distribution of the trapping layer altitude is shown in Figure 7. Figure 6. Distribution of the trapping layer occurrence rate with local time variations. The blank grids represent that there is no radio occultation event.

Altitude of Trapping
The trapping layer altitude ℎ , has an important influence on the propagation of radio waves. Taking the average value of the trapping layer altitude in each grid, the seasonal distribution of the trapping layer altitude is shown in Figure 7. The average altitude of the trapping layer is lower near the coast of the mainland than on the sea. This phenomenon can also be seen in Engeln's results. With the distance away from the coastline, the altitude increases, such as is seen on the west coast of North America and the west coast of Africa. Similar situations also exist in continental seas, such as the Mediterranean Sea, the Persian Gulf, the Gulf of Aden, and the Red Sea [36]. The trapping layer on the land is mainly caused by radiative cooling, which mostly occurs on the surface or at low altitude [38]. Therefore, the trapping layer altitude is usually close to The average altitude of the trapping layer is lower near the coast of the mainland than on the sea. This phenomenon can also be seen in Engeln's results. With the distance away from the coastline, the altitude increases, such as is seen on the west coast of North America and the west coast of Africa. Similar situations also exist in continental seas, such as the Mediterranean Sea, the Persian Gulf, the Gulf of Aden, and the Red Sea [36]. The trapping layer on the land is mainly caused by radiative cooling, which mostly occurs on the surface or at low altitude [38]. Therefore, the trapping layer altitude is usually close to the land. As the distance into the sea increases, the temperature inversion gradually weakens, and the layer's altitude gradually increases [11,39]. The emergence of the trapping layer is due to the sharp gradient of humidity and temperature related to the PBL stratocumulus, and the increasing altitude of the trapping layer away from the west coastline is related to the trade wind cumulus [8]. On account of being near the equator in the Southern Hemisphere, the northwest sea area of Australia in winter, under the control of tropical depression and the strong water vapor in the updraft, it is easy for stratocumulus to form in the air, thus producing a higher trapping layer. In summer, controlled by high pressure, the trapping layer is mainly caused by subsidence inversion. Therefore, the trapping layer occurrence rate, range and altitudes are lower. One point of concern is that areas with the larger trapping layer altitudes begin to shift southward in the Southern Hemisphere in autumn and winter.
The local time distribution of the average trapping layer altitude is shown in Figure 8. From the analysis of the diurnal changes in altitude, it can be seen that the altitude of the trapping layer is much higher in the daytime on the ocean than at night on the land. At LT00, most of the land trapping layer altitude is less than 300 m. During the day, solar radiation produces dry convective motion, which leads to the growth of the boundary layer. Trapping layers at altitudes of more than 1.3 km mainly occur in areas with a high occurrence rate, such as western Mexico and the south of America.
The local time distribution of the average trapping layer altitude is shown in Figure  8. From the analysis of the diurnal changes in altitude, it can be seen that the altitude of the trapping layer is much higher in the daytime on the ocean than at night on the land. At LT00, most of the land trapping layer altitude is less than 300 m. During the day, solar radiation produces dry convective motion, which leads to the growth of the boundary layer. Trapping layers at altitudes of more than 1.3 km mainly occur in areas with a high occurrence rate, such as western Mexico and the south of America.

Intensity of Trapping
The trapping layer intensity refers to the difference in modified refractivity when the gradient dM dh < 0 ⁄ . As shown in Figure 1, the modified refractivity value corresponding to the top of the trapping layer h is a , and the modified refractivity value corresponding to the bottom of the trapping layer h is a . The value of a − a is the intensity of the trapping layer. The greater the trapping layer intensity is, the greater the influence on radio wave propagation is.
According to the seasonal division, the average value of the trapping layer intensity in each grid point of latitude and longitude is taken to obtain the seasonal distribution of the trapping layer intensity in Figure 9. The intensity of the trapping layer on land is mostly between 0M-unit and 15M-unit, which is generally lower than that of the sea.

Intensity of Trapping
The trapping layer intensity refers to the difference in modified refractivity when the gradient dM/dh < 0. As shown in Figure 1, the modified refractivity value corresponding to the top of the trapping layer h t is a s , and the modified refractivity value corresponding to the bottom of the trapping layer h m is a m . The value of a m − a s is the intensity of the trapping layer. The greater the trapping layer intensity is, the greater the influence on radio wave propagation is.
According to the seasonal division, the average value of the trapping layer intensity in each grid point of latitude and longitude is taken to obtain the seasonal distribution of the trapping layer intensity in Figure 9. The intensity of the trapping layer on land is mostly between 0 M-unit and 15 M-unit, which is generally lower than that of the sea. The trapping layer in the regions with a high occurrence rate usually has stronger intensity. In the eastern areas of oceans, the western coasts of continents and the Arabian Sea, and the northwest of Australia all illustrate this feature well. In the Arabian Sea, the occurrence rate and ranges of the trapping layer are the highest in spring, where a largescale trapping layer intensity >30M-unit appears. However, when this area is affected by The trapping layer in the regions with a high occurrence rate usually has stronger intensity. In the eastern areas of oceans, the western coasts of continents and the Arabian Sea, and the northwest of Australia all illustrate this feature well. In the Arabian Sea, the occurrence rate and ranges of the trapping layer are the highest in spring, where a largescale trapping layer intensity >30 M-unit appears. However, when this area is affected by the summer monsoon, the intensity reduces to 13~20 M-unit. In northwestern Australia, under the influence of the local pressure system, areas with intensity greater than 20 M-unit are more common in autumn and winter than in the other two seasons. Compared with the northwest and southwest regions of Africa, the intensity of the trapping layer in the northwest is over 24 M-unit in spring and summer, which is stronger than that in autumn and winter. In southwest Africa, on the other hand, it is stronger in autumn and winter than in spring and summer.
As it is affected by the intertropical convergence zone, the Northern Hemisphere has more areas with a stronger trapping layer in spring and summer than the Southern Hemisphere. In autumn and winter, the situation is reversed with the southward shift of the intertropical convergence zone. Figure 10 shows the local time distribution of the trapping layer intensity. The intensity in most areas does not exceed 12 M-unit, and in areas with high occurrence rates, the intensity of the trapping layer is above 20 M-unit. The trapping layer intensity in the Arabian Sea reached the maximum value of 45 M-unit in LT18. The emergence of this strong trapping layer is mainly related to the increase in the temperature gradient and the humidity gradient, both of which are affected by solar radiation and sea-land breezes.

Thickness of Trapping
The occurrence rate of the trapping layer is closely related to the inversion layers [40]. Several studies of the PBL have already found that the inversion thickness at the top of the PBL varies between about 100 and 400 m [41]. Figure 11 shows the average thickness of the trapping layer in each season. The results show that the average thickness varies from 80 to 240 m. The thickness on land is smaller than that over the ocean, usually below 120 m.

Thickness of Trapping
The occurrence rate of the trapping layer is closely related to the inversion layers [40]. Several studies of the PBL have already found that the inversion thickness at the top of the PBL varies between about 100 and 400 m [41]. Figure 11 shows the average thickness of the trapping layer in each season. The results show that the average thickness varies from 80 to 240 m. The thickness on land is smaller than that over the ocean, usually below 120 m.
The thickness is larger in the regions with a high occurrence rate of the trapping layer, too. The Arabian Sea has the most obvious seasonal variation, with the thickness exceeding 220 m in spring. Under the impact of the southwest monsoon in summer, the thickness decreases to about 100 m before increasing again in autumn and winter. In Australia, the thickness is obviously larger in autumn and winter than in spring and summer. These changes are consistent with the variation of the occurrence rate.

Thickness of Trapping
The occurrence rate of the trapping layer is closely related to the inversion layers [40]. Several studies of the PBL have already found that the inversion thickness at the top of the PBL varies between about 100 and 400 m [41]. Figure 11 shows the average thickness of the trapping layer in each season. The results show that the average thickness varies from 80 to 240 m. The thickness on land is smaller than that over the ocean, usually below 120 m. The thickness is larger in the regions with a high occurrence rate of the trapping layer, too. The Arabian Sea has the most obvious seasonal variation, with the thickness exceeding 220 m in spring. Under the impact of the southwest monsoon in summer, the thickness The trapping layer thickness on the northwest coast of Africa tends to shift northward in summer, and the thickness saw little overall variation; the average thickness is around 200 m. A large-scale lower trapping layer thickness appears over the Asian continent in summer, mainly due to the surface temperature inversion caused by radiative cooling.
In general, the trapping layer thickness seasonal variation is not significant; only a few areas have obvious changes. In spring and summer, the thickness of the Northern Hemisphere is greater than that of the Southern Hemisphere, but the opposite is true in autumn and winter.
The local time variations distribution of the trapping layer thickness is shown in Figure 12. Large areas of the trapping layer with an average thickness of more than 210 m appeared in the western sea areas of South America at LT12, and the average thickness in this area was between 150 m and 180 m at other times. The thickness of the areas in the northwest of Australia are higher in the daytime than at night, while the thickness of the central and northern African continent is greater at LT18 than at LT06. decreases to about 100 m before increasing again in autumn and winter. In Australia, the thickness is obviously larger in autumn and winter than in spring and summer. These changes are consistent with the variation of the occurrence rate. The trapping layer thickness on the northwest coast of Africa tends to shift northward in summer, and the thickness saw little overall variation; the average thickness is around 200 m. A large-scale lower trapping layer thickness appears over the Asian continent in summer, mainly due to the surface temperature inversion caused by radiative cooling.
In general, the trapping layer thickness seasonal variation is not significant; only a few areas have obvious changes. In spring and summer, the thickness of the Northern Hemisphere is greater than that of the Southern Hemisphere, but the opposite is true in autumn and winter.
The local time variations distribution of the trapping layer thickness is shown in Figure 12. Large areas of the trapping layer with an average thickness of more than 210m appeared in the western sea areas of South America at LT12, and the average thickness in this area was between 150m and 180m at other times. The thickness of the areas in the northwest of Australia are higher in the daytime than at night, while the thickness of the central and northern African continent is greater at LT18 than at LT06. The local time variations of the trapping layer's thickness in other regions are not obvious

Conclusions
We have used many years of historical occultation observation data and post-processing data to analyze and study the temporal and spatial distribution characteristics of

Conclusions
We have used many years of historical occultation observation data and post-processing data to analyze and study the temporal and spatial distribution characteristics of the occurrence probability, altitude, intensity, and thickness of the trapping layer. The following conclusions can be drawn: In general, the regions with high occurrence rates of the trapping layer are mainly concentrated in the eastern areas of the ocean and the western areas of the continents between 50 • S and 50 • N. Among them, the western oceans of Mexico, the western coast of Africa, the Arabian Sea, Red Sea, Mediterranean Sea and the sea areas of northwest Australia have a more than 50% occurrence rate of the trapping layer. The high trapping layer occurrence rate usually corresponds to relatively large altitude, intensity, and thickness.
There are geographical and seasonal differences in the distribution of the trapping layer characteristics in the Northern and Southern Hemispheres. The altitude, intensity and thickness of the trapping layer over the ocean are greater than that over the land. Generally, the altitude of the trapping layer is lower in the boundary areas between the oceans and the continents. The increased altitude of the trapping layer away from the western coastline is related to the trade wind cumulus. In spring and summer, the occurrence rate and range of the trapping layer in the Northern Hemisphere are higher than those in the Southern Hemisphere, but in autumn and winter, there are more occurrences in Southern Hemisphere than in the Northern Hemisphere.
The geographical and seasonal differences in the distribution of the trapping layer are affected by many climatic factors. Climatology shows that the areas with high occurrence rates of the trapping layers in western Mexico, the coastal areas of western Africa, the Red Sea, the Arabian Sea, the Mediterranean, and the Bay of Bengal are mainly affected by trade wind cumulus, ocean currents, equatorial convergence zone deviation, radiative cooling and other factors.
The diurnal variation of the occurrence rate, altitude, intensity and thickness is not obvious except for in a few areas. The main reasons for daily changes can be roughly divided into the following categories: firstly, boundary layer growth and boundary inversion caused by strong solar radiation heating; and secondly, boundary layer changes influenced by sea-land breezes.