1. Introduction
The planetary boundary layer (PBL) is the low atmosphere close to the Earth's surface, is directly affected by the forced action of the Earth's surface, with a height up to several kilometers, and is the main area of human activities[
1]. The Earth's surface exchanges heat, momentum, water vapor, and chemicals with the free atmosphere through the PBL. Since the PBL involves many processes, such as convective activity, turbulent mixing, low-layer cloud formation, pollutant diffusion and surface energy budget, the PBL plays a key role in many aspects, such as weather, climate and air pollution[
2]. The planetary boundary layer height (PBLH) is an important parameter for characterizing the structural characteristics of the PBL. It determines the vertical degree of turbulent mixing, convective transport and pollutant diffusion in the PBL and strongly influences the development and evolution of convective activity [
3,
4,
5]. PBLH is one of the most important parameters in weather forecasting, climate and air quality forecasting models.
Currently, there are many datasets and methods to obtain PBLH. Among them, radiosonde datasets are one of the most widely used datasets. The radiosonde can obtain the vertical profiles of physical quantities such as temperature, humidity, wind direction and wind speed and diagnose PBLH through the parcel method, potential temperature gradient method, relative humidity gradient method, bulk Richardson number method, and so on[
6]. In addition to radiosonde data, remote sensing data can also be used to estimate the PBLH, including lidar[
7,
8], wind profiler radar[
9,
10], ceilometer[
11,
12], etc. Radar can achieve continuous observations, thus obtaining high-temporal-resolution data. With the development of satellite remote sensing technology, GPS radio occultation data are also used in PBL studies[
13,
14], as are cloud-aerosol lidar data with orthogonal polarization[
15,
16]. Although the resolution is relatively low, they can expand the spatial coverage of observations. In addition to the above sounding data, reanalysis data and model data are also commonly used in PBLH studies[
17,
18,
19]. They have the advantages of high temporal and spatial resolutions and can realize long-term and large-scale PBLH climatology studies.
Climatology research on the PBLH began with the work of Holzworth[
20], who first used radiosonde data to obtain the characteristics of the height of the mixed layer in the United States. Later, Liu and Liang[
21] studied the diurnal cycle climatology of the PBLH using radiosonde data from 14 stations around the world. Seidel et al.[
17] used radiosonde data, reanalysis data and two model datasets to analyze the diurnal and seasonal variations in the PBLH in the United States and Europe and obtained more climate characteristics of the PBLH. At present, many scholars have conducted research on the PBLH in different regions of China through different datasets and methods [
16,
22,
23,
24]. Zhao et al.[
25] studied the characteristics of the PBLH in Northwest China using radiosonde data from 2015 to 2016. Guo et al.[
23] studied the PBLH climatology in China based on radiosonde data from 2011 to 2015. However, due to the difficulty of PBL observation, it is difficult to obtain long-term, high-temporal-resolution observational data, resulting in most studies being confined to specific locations or datasets with short durations and low resolutions. Therefore, it is still necessary to use reanalysis data for long-term PBLH climatology studies. Seidel et al.[
17] showed that reanalysis data such as ERA-Interim and NCAR-CAM5 could reproduce well the spatial and seasonal distributions of the PBLH. Zhao et al.[
26] studied the variation in the PBLH in arid and semiarid regions of East Asia over 110 years through ERA-20C reanalysis data. von Engeln et al.[
18] studied global PBLH climatology from 1990 to 2009 using ERA-Interim reanalysis data. ERA5 is the latest reanalysis data of the European Centre for Medium-Range Weather Forecasts (ECMWF). Altabakash et al.[
27] studied the PBLH of the Korean Peninsula using ERA5 data and noted that the PBLH obtained from ERA5 data is very consistent with the PBLH obtained from radiosonde and GPS radio occultation data. We use PBLH data from ERA5 to carry out our research. First, we use radiosonde data to verify ERA5 data and then use ERA5 data to study the climatology of the PBLH, which will provide some reference for the application of ERA5 data and the study of the PBLH in the future.
In recent years, the variation trend of the PBLH has been widely studied by scholars[
24,
28,
29]. A study by Yang et al.[
30] used ground-based lidar data and noted that the PBLH in Hong Kong showed a downward trend from 2004 to 2009. Zhang et al.[
28] used radiosonde data from 25 sites in Europe from 1973 to 2010 to reveal the increasing trend of PBLH and attributed it to the decrease in near-surface relative humidity and the increase in near-surface temperature. A recent study showed that the PBLH in China experienced dramatic mutation in 2004, but the trend varied from place to place[
24]. Obviously, previous studies have shown that the variation trend of the PBLH may vary due to different data, methods or study locations. Therefore, it is necessary to study the long-term variation trend of the PBLH.
Located in the eastern coastal area of China, Jiangsu Province is China's second largest economic province with a dense population and rapid industrial development and urbanization. The PBL is the main area of human activities. The meteorological conditions and atmospheric environment in the PBL have an important influence on human life and health. PBLH determines the atmospheric environmental capacity and is one of the important factors affecting the transportation and diffusion of pollutants. At present, there have been few studies on the PBLH in Jiangsu. Based on the ECMWF ERA5 reanalysis dataset, this paper studies the climatic characteristics of the PBLH in Jiangsu and its relationship with meteorological factors, analyzes the physical mechanisms affecting the development of the PBL, and reveals the long-term variation trend and causes of the PBLH.
The rest of the paper is arranged as follows:
Section 2 introduces the data and methods used;
Section 3 verifies the PBLH based on ERA5;
Section 4 analyzes the climatic characteristics of the PBLH in Jiangsu and its relationship with meteorological factors and analyzes the long-term trend and causes of the PBLH; and
Section 5 presents a summary and discussion.
2. Data and Methods
2.1. Data
2.1.1. ERA5 Reanalysis Data
ERA5 reanalysis data are the fifth-generation atmospheric reanalysis data of the global climate from the European Centre for Medium-Range Weather Forecasts (ECMWF) and the latest generation reanalysis data of the ECMWF. They assimilate several satellite and radiosonde datasets, including various observational datasets from the World Meteorological Organization's Global Telecommunication System (GTS). The ERA5 data cover the entire globe, with a horizontal resolution of 0.25° x 0.25° grid, vertical resolution of 37 standard pressure layers, and hourly temporal resolution. In this paper, ERA5 boundary layer height data are used to study the climatic characteristics of the planetary boundary layer height in Jiangsu. The selected time range is from 1979 to 2019, a total of 41 years. In addition, to study the relationship between the PBLH and meteorological factors, data such as 2 m temperature, 10 m wind speed and near-surface relative humidity (replaced by relative humidity at 1000 hPa) from ERA5 are used, as well as surface pressure and 700 hPa temperature data, which are used to calculate lower tropospheric stability (defined as the potential temperature difference between 700 hPa and the surface[
31]).
2.1.2. Radiosonde Data
We used radiosonde data in Jiangsu to verify the PBLH based on ERA5. The selected radiosonde data are obtained from the radiosonde network of the China Meteorological Administration (CMA). Since 2002, the CMA has been upgrading its radiosonde system by using an L-band radar sounding system that combines a secondary wind-measuring radar and a GTS1 digital electronic radiosonde. The system mainly detects meteorological elements such as wind direction, wind speed, temperature and pressure and provides high vertical resolution (second-level) profile data. By 2010, the L-band sounding system had been applied to 120 radiosonde stations across China. The radiosonde data used in this paper are the L-band conventional sounding data (twice daily: 08:00, 20:00, LST) at Xuzhou Station, Sheyang Station and Nanjing Station in Jiangsu from January 1, 2010 to December 31, 2019, and include temperature, pressure, wind direction, wind speed and other conventional meteorological elements. The locations of the three stations are shown in
Figure 1. In total, we made use of 21588 profiles at the three stations, including 3600 profiles at 08:00 (LST) and 3600 profiles at 20:00 (LST) at Xuzhou Station, 3595 profiles at 0800 (LST) and 3600 profiles at 2000 (LST) at Sheyang Station, and 3591 profiles at 0800 (LST) and 3602 profiles at 2000 (LST) at Nanjing Station.
2.2. PBLH diagnosis method
There are many methods to diagnose the PBLH, and the commonly used methods include the potential temperature gradient method, parcel method, relative humidity gradient method, specific humidity gradient method, refractive index gradient method, and Richardson number method. Seidel et al.[
6,
17] compared and evaluated different diagnostic methods and confirmed that the Richardson number method is more suitable for the diagnosis of PBLH in large volume datasets, and this method is suitable for both the stable boundary layer (SBL) and convective boundary layer (CBL). Therefore, in this paper, we choose the Richardson number method to diagnose the PBLH. The Richardson number method determines the PBLH at which the bulk Richardson number first reaches the critical value. Referring to the study of Sicard et al.[
32], the formula for calculating the bulk Richardson number (Ri) is:
where g is the gravitational acceleration, z0 is the surface altitude, θ is the potential temperature, and u and v are the zonal wind and meridional wind components, respectively. Previous theories and studies have shown that laminar flow is unstable when the bulk Richardson number Ri is less than the critical value (~ 0.25) [
1]. Guo et al.[
23] determined the critical value as 0.25 in their study and noted that the uncertainty caused by the selection of the critical value (0.2, 0.25, 0.3) is very small. Therefore, in this paper, the critical value is determined as 0.25, and the height when the bulk Richardson number Ri first reaches 0.25 is determined as the PBLH.
2.3. Trend analysis
Following the method used by Guo et al.[
24], Sen’s slope estimator[
33] to calculate the trend is selected here. The trend is determined as the median of the slopes between all paired values, which is calculated as follows:
where
and
represent the ith and jth values in the time series of X, respectively. The median is the 50% value of the total datasets. Compared with linear regression, Sen's slope is less sensitive to outliers. In addition, the nonparametric Mann–Kendall statistical test[
34,
35] is used to determine the 95% confidence level.
3. Data validation
The accuracy of some general variables in ERA5 data has been confirmed in some studies[
27,
36], but the quality of the PBLH in these data has not been verified. To use this dataset to study the characteristics of the PBLH in Jiangsu, the first problem we need to solve is to evaluate the accuracy of the PBLH. In this paper, the PBLH diagnosed by radiosonde data from Xuzhou, Sheyang and Nanjing stations is used to evaluate the PBLH based on ERA5. To compare the two kinds of PBLH, first, the ERA5 data are sampled at nine grid points centered on the radiosonde position, and then the arithmetic mean is calculated. Second, to maintain consistency with the time of radiosonde data, ERA5 data are sampled at 08:00 (LST) and 20:00 (LST) from 2010 to 2019.
Figure 2 shows the scatter plots of the ERA5 PBLH and the PBLH diagnosed by the radiosonde data from the Xuzhou, Sheyang and Nanjing stations. At Xuzhou Station, the correlation coefficient between the ERA5 and radiosonde datasets is 0.81 at 08:00 (
Figure 2a), and the PBLHs of the two datasets have very good agreement. At 20:00 (
Figure 2d), the correlation coefficient between these two is 0.71, and the root mean square error is 198 m. Compared with the radiosonde data, ERA5 generally underestimates the PBLH by 84 m on average. The scatter plots of Sheyang Station show that the correlation coefficient between the ERA5 and radiosonde datasets at 08:00 (
Figure 2b) is 0.77, and the root mean square error is 208 m. Compared with radiosonde data, ERA5 generally overestimates the PBLH by 127 m on average. At 20:00 (
Figure 2e), the PBLHs based on the two kinds of datasets are very consistent, and the correlation coefficient reaches 0.82. At Nanjing Station, the correlation coefficient between the ERA5 and radiosonde datasets at 08:00 (
Figure 2C) is 0.85, and the root mean square error is 127 m. Compared with the radiosonde data, ERA5 slightly overestimates the PBLH by 62 m on average. At 20:00 (
Figure 2f), the correlation coefficient between these two is 0.79, and the root mean square error is 156 m. ERA5 slightly underestimates the PBLH by 31 m on average. All the above correlation coefficients pass the 99% significance test. The scatter plots of the PBLHs between the radiosonde and ERA5 datasets demonstrate the close correspondence between these two kinds of PBLHs and provide confidence that the PBLH based on ERA5 can be used to study the climatology of the PBLH over Jiangsu, China.
4. Results
4.1. Climatological pattern of the PBLH
As a gross climatological overview,
Figure 3 shows the spatial distribution of the mean PBLH at different times in Jiangsu. As shown in the figure, the PBLH was highest in the daytime (
Figure 3(b)), ranging from 700 to 1200 m. At night (
Figure 3(d)), the PBLH was lowest, ranging from 100 to 400 m. The morning (
Figure 3(a)) and evening (
Figure 3(c)) were in the transitional period between the daytime CBL and nighttime SBL, with PBLHs between 200 and 500 m. In the daytime, the PBLH was higher in the north and lower in the south. The highest PBLH was in the northwest, reaching 1200 m, and the lowest PBLH was in Hongze Lake in the western part and Tai Lake in the southern part of Jiangsu, at approximately 700 m. The PBLH in the eastern coastal areas was generally lower than that in the inland areas. At night, morning and evening, the PBLH was lower in the north and higher in the south. The lowest PBLH was in the northwest, and the highest PBLH was in the two lakes (Hongze Lake in the west and Tai Lake in the south) and the eastern coastal area, with the highest PBLH reaching approximately 500 m.
4.2. Seasonal variation
Figure 4 shows the seasonal distribution of the PBLH at different times in Jiangsu. As shown in the figure, in the daytime (14:00 LST), the PBLH ranged from 200 m to 2000 m, with a median of approximately 1000 m. In spring, the PBLH could reach more than 2000 m. The PBLH was usually highest in spring, followed by fall and summer, and lowest in winter. At night (02:00 LST), the PBLH was below 600 m, and the median was approximately 200 m. At night, morning (08:00, LST) and evening (20:00, LST), the PBLH generally followed a seasonal variation, that is, summer > spring > fall > winter, and the seasonal variation was more obvious in the morning.
Figure 5 shows the spatial distribution of the seasonal mean PBLH at different times, revealing the spatial variation characteristics not shown in
Figure 4.
Figure 5 (a1-a4) shows that at 08:00 (LST, the same below), the seasonal variation in the PBLH in inland areas of Jiangsu was higher in spring and summer than in fall and winter. The average PBLH was approximately 300-500 m in spring and summer and 200-400 m in fall and winter. However, the seasonal variation in the PBLH in coastal areas was obviously opposite to that in inland areas, showing that the PBLH in autumn and winter was higher than that in spring and summer, and the average PBLH in autumn and winter was approximately 400-700 m. This may be because the PBLH in coastal areas is affected by the ocean and has the characteristics of the height of the marine atmospheric boundary layer; that is, lower in spring and summer and higher in autumn and winter[
37,
38]. PBLH showed obvious southeast‒northwest gradient changes in the four seasons. The PBLH was the highest in the southeast, reaching over 400 m in summer, the lowest in the northwest, and only approximately 200 m in winter. This may be because the southeast receives stronger solar radiation and the surface temperature is higher[
23].
At 14:00 (
Figure 5(b1-b4)), the seasonal variation in the PBLH in inland areas of Jiangsu was ranked as spring > fall > summer > winter. The average PBLH in spring was approximately 1000-1400 m, and that in winter was approximately 900-1100 m. Affected by the ocean, the PBLH in coastal areas was still higher in fall and winter than in spring and summer. The spatial distribution of the PBLH in spring and summer was opposite to that at 08:00; that is, the PBLH in the northwest was higher than that in the southeast. In fall and winter, the highest PBLH appeared in central Jiangsu. Guo et al.[
23] pointed out that the PBLH is higher in dry areas during the daytime. Compared with the southeastern part, the northwestern part of Jiangsu is drier, which means that in addition to solar radiation, humidity and other factors play an important role in controlling the spatial distribution of the PBLH during the daytime. At the same time, it can be noted that the PBLH in the two lakes was relatively low in all four seasons, especially in winter, when the average PBLH was below 700 m. This is because the high specific heat capacity of water causes the lake to warm up slowly and heat the atmosphere weakly, which then leads to the low PBLH.
At 20:00 (
Figure 5(c1-c4)) and 02:00 (
Figure 5(d1-d4)), the seasonal and spatial variables of the PBLH were generally similar to those at 08:00. However, since the SBL began to appear at night, the PBLH decreased, the average PBLH was approximately 100-500 m, and the seasonal variables were less obvious.
4.3. Diurnal variation
The high-resolution hourly ERA5 data allow us to obtain more accurate diurnal variations in the PBLH than twice-daily radiosonde data.
Figure 6 shows the diurnal cycle of the PBLH in the four seasons in Jiangsu. As shown in the figure, in all seasons, the PBLH had obvious diurnal variation characteristics, and the difference between different seasons was small. The diurnal cycle of PBLH was as follows: at night, the ground released longwave radiation, the surface cooled down, the near-surface layer was relatively stable, the PBLH decreased slowly reaching the lowest point at 06:00 (LST), and the PBLH was approximately 200 m00m; after sunrise, the solar radiation increased and heated the ground, the CBL began to develop, and the PBLH gradually increased; in the afternoon, the PBLH reached its maximum, approximately 900-1200 m200m; then, with the decrease of solar radiation, the turbulence intensity in CBL gradually attenuated, and the PBLH decreased rapidly; after sunset, the SBL began to form, and the PBLH remained low at night, approximately 200-300 m00m. At the same time, the decrease in the PBLH in the afternoon was more rapid than the increase in the PBLH in the morning, which is consistent with the findings of Liu and Liang. At night, the seasonal variation in the PBLH was very small, and the PBLH in summer was slightly higher than that in other seasons, which is consistent with the results of Gu et al.[
38] In the daytime, especially in the afternoon (12:00-17:00 LST), the PBLH showed obvious seasonal variables, but the PBLH during summer was not the highest. In the afternoon, the PBLH was the highest in spring, followed by the PBLH in fall and summer, and the lowest was in winter, indicating that the PBLH in the daytime is related to other influencing factors besides solar radiation, which is analyzed in
Section 4.4.
Figure 7 (left) shows the diurnal variation amplitude of the PBLH in Jiangsu, calculated by the daily maximum minus the minimum. As shown in the figure, the diurnal variation amplitude of the PBLH has obvious spatiotemporal characteristics. The amplitude in spring was the largest and decreased gradually from north to south. In spring, the strongest diurnal variation was in the northern area of Jiangsu, and the average diurnal variation amplitude was more than 1300 m. The diurnal variation was the weakest in the southern and coastal areas of Jiangsu, with an average amplitude below 700 m. The diurnal variation amplitude of fall and summer was lower than that of spring. The largest amplitude in fall was in central Jiangsu. The spatial distribution of the diurnal variation amplitude in summer was similar to that in spring, but the intensity was much weaker than that in spring. The diurnal variation amplitude of winter was the smallest. It is worth noting that the diurnal variation amplitude of the PBLH in coastal areas was smaller than that in inland areas, especially in spring and summer. At the same time, the diurnal variation amplitude of the PBLH in Hongze Lake and Tai Lake was also lower. This is because the high specific heat capacity of water leads to smaller diurnal variations in water temperature, which in turn leads to smaller diurnal variations in the PBLH.
Figure 7 (right) shows the time when the daily PBLH reached its peak. As shown in the figure, in most areas, the PBLH reached its peak at 15:00 (LST) in spring and winter and at 14:00 (LST) in summer and fall. This difference corresponds to the time of sunrise. After sunrise, the ground heats up, which is the determining factor in the development of the PBL. Due to late sunrise in winter and spring, it is later for the ground to accumulate enough heat to reach the peak. Furthermore, the time when the PBLH reached its peak in the eastern area of Jiangsu was usually earlier than that in the western area, which was also because the time of sunrise in the eastern area was earlier.
4.4. Affecting factors for PBLH diurnal variation
In addition to solar radiation, the PBLH also has a significant relationship with meteorological factors [
17]. According to a study by Gu et al.[
38], the diurnal cycle of the PBLH is significantly positively correlated with the near-surface temperature and negatively correlated with the relative humidity. Zhang et al.[
22] noted that higher surface temperature and lower relative humidity lead to higher sensible heat flux and lower latent heat flux, resulting in deeper convection and a higher PBLH. To further understand the influence of meteorological factors on the diurnal cycle of the PBLH in Jiangsu and analyze the reason for the seasonal variation in the diurnal cycle, we study the diurnal cycles of several meteorological factors in the four seasons and calculate the correlations between them and the PBLH. The selected meteorological factors include 2 m temperature (T2m), 10 m wind speed (WS), near-surface relative humidity (RH) and lower tropospheric stability (LTS). LTS is defined as the potential temperature difference between 700 hPa and the surface, which can represent the thermodynamic state of the low troposphere. A larger LTS usually reflects a more stable low troposphere[
39,
40].
Figure 8 shows the diurnal cycles of meteorological factors in the four seasons. Comparing the diurnal cycle of the PBLH (
Figure 6), we find that the diurnal cycles of T2m and WS have a positive relationship with the diurnal cycle of the PBLH. A higher T2m and stronger WS correspond to a higher PBLH, indicating that stronger surface heating and surface forcing contribute to the development of the PBL. On the other hand, the diurnal cycles of RH and LTS have obvious inverse relationships with the diurnal cycle of the PBLH. A lower RH and lower LTS correspond to a higher PBLH because the lower the RH is, the lower the latent heat flux, the more sensible is the heat flux, and the stronger is the development of the PBL. Moreover, the lower the LTS is, the more conducive it is to convective activity, and the easier it is to develop the PBL.
Figure 6 also shows that at night, the PBLH in summer was slightly higher than that in other seasons. This may have been due to the higher T2m (
Figure 8(a)) and lower LTS (
Figure 8(d)) in summer. However, during the daytime, the highest PBLH occurred in spring, not in summer.
Figure 8 (b, c) shows that in the daytime, the WS of spring was highest and the RH was lowest, while in summer, the WS was lower and the RH was higher. Therefore, it can be concluded that the seasonal variation in daytime PBLH is more related to WS and RH. PBLH was the highest in spring, which may have been due to the high WS in spring[
41,
42]. High WS enhances low-layer wind shear and upward forcing, thus enhancing turbulence and increasing the PBLH. The lower PBLH in summer may be because Jiangsu is in the plum rain season with high RH and more precipitation[
43]. Water vapor absorbs more heat, which reduces the heating efficiency of the surface to the air, thus inhibiting the development of the PBL.
We calculate the correlation coefficients between the PBLH and meteorological factors based on the diurnal cycles in the four seasons (
Table 1), with bold data indicating significant correlations at the 95% or higher confidence level. The results show that the diurnal cycle of the PBLH is significantly correlated with the diurnal cycle of meteorological factors in the four seasons. Specifically, the PBLH is significantly positively correlated with the T2m and WS and significantly negatively correlated with the RH and LTS. This is consistent with the research results of Gu et al.[
38] These correlations indicate that on the daily scale, higher T2m and WS and lower RH and LTS will correspond to higher PBLH.
4.5. Long-term trend of PBLH over Jiangsu and affecting factors
In addition to the climatic characteristics and influencing factors of the PBLH, the long-term variation trend of the PBLH is another issue that needs attention. We calculate the variation trend of the annual mean PBLH from 1979 to 2019 in Jiangsu. The results show that the daytime PBLH had a significant increasing trend (
Figure 9), while the PBLH at night had no significant variation trend (not shown), which is consistent with the results of Zhang et al.[
28] and Gu et al.[
38]
Figure 9a shows the time series of the daytime (14:00 LST) annual mean PBLH in Jiangsu. As shown in the figure, the daytime PBLH had an increasing trend, especially since the 1990s, when the PBLH had obviously increased.
Figure 9b shows the trend of the PBLH in Jiangsu. The results show that the PBLH in most areas of Jiangsu had a significant growth trend, and the fastest growing areas were in the central and southern areas of Jiangsu, with the highest trend reaching 3 m•yr
-1.
To explore the reason for the increase in the daytime PBLH, we analyze the variation trend of the impact factors from 1979 to 2019. The results show that daytime T2m, RH and LTS had significant variation trends (Figs. 10, 11), but WS had no significant variation trend (not shown).
Figure 10 shows the time series of the daytime (14:00 LST) annual mean T2m, RH and LTS and the correlation coefficients with the annual mean PBLH. The annual mean PBLH has a significantly positive correlation with T2m, with a correlation coefficient of 0.56 (
Figure 10a), and a significantly negative correlation with RH and LTS, with correlation coefficients of -0.9 and -0.82, respectively (
Figure 10b, c). At the same time, we also find that in the past 41 years, T2m had obviously increased, and RH and LTS had obviously decreased, especially since the 1990s. This is consistent with the findings of Xia et al.[
44] and Chen et al.[
45] Therefore, we can calculate the variation trend of meteorological factors in Jiangsu (
Figure 11). The results show that T2m had a significantly increasing trend, and RH and LTS had significantly decreasing trends. At the same time, we find that the spatial distribution of the variation trend of meteorological factors was similar to that of the PBLH; that is, the T2m in the southern area increased faster than that in the northern area, and the RH and LTS also decreased faster, which may be the reason for the faster growth of the PBLH in the south.
Finally, we calculate the correlation coefficients between the annual mean PBLH and T2m and RH and LTS in Jiangsu (
Figure 12). The results show that the annual mean PBLH has a significantly positive correlation with T2m and a significantly negative correlation with RH and LTS. Moreover, the areas with the strongest correlation correspond to the areas with the fastest growing trend of PBLH. Therefore, we can conclude that the increase in daytime PBLH may be due to the increase in T2m and the decrease in RH and LTS. This is consistent with the above results: a higher T2m, lower RH and lower LTS contribute to a higher PBLH.
5. Conclusion and Discussion
Based on the hourly reanalysis data from ERA5 from 1979 to 2019, the climatology of the PBLH in Jiangsu, China, is studied in this paper.
The PBLH based on ERA5 is verified by radiosonde data firstly, and the results show that the PBLH based on ERA5 fits very well with the PBLH diagnosed by radiosonde data, with correlation coefficients exceeding 0.7 at all three stations. This shows that it is valid to use the PBLH based on ERA5 data in the Jiangsu region.
In Jiangsu, the daytime PBLH was highest, and the average PBLH was approximately 700-1200 m, which was higher in the north and lower in the south. At night, the PBLH was lowest, and the average PBLH was between 100 and 400 m, lower in the north and higher in the south. PBLH exhibited complex spatiotemporal variation. In the daytime, the PBLH in inland areas was highest in spring, followed by fall and summer, and lowest in winter. A higher wind speed in spring produces stronger ground forcing, which increases the PBLH. Higher relative humidity in summer increases the latent heat flux and decreases the sensible heat flux, leading to a lower PBLH. At night, the seasonal variation in the PBLH was less obvious. Affected by the ocean, the seasonal variation in the PBLH in coastal areas was opposite to that in inland areas, higher in fall and winter and lower in spring and summer.
The PBLH had an obvious diurnal cycle, usually reaching its peak at 14:00 (LST) or 15:00 (LST). The diurnal variation in the PBLH was largest in spring and smallest in winter. The diurnal cycle of the PBLH was significantly positively correlated with the near-surface temperature and wind speed and significantly negatively correlated with the near-surface relative humidity and lower tropospheric stability. The daytime PBLH had a significantly increasing trend in most areas. The increase in PBLH can be attributed to the increase in near-surface temperature and the decrease in near-surface relative humidity and lower tropospheric stability.
The results of this paper will provide the basis for the application of ERA5 reanalysis data in eastern China. Furthermore, the study of the PBLH in Jiangsu will also help us to understand the characteristics of the PBL and the transportation and diffusion of pollutants and provide a reference for weather forecasting, climate and air quality prediction.
Author Contributions
Conceptualization, Xiang Li; Methodology, Xiang Li and Yiwen Dong; Formal Analysis, Yiwen Dong; Investigation, Yiwen Dong; Data Curation, Jiawei Yao; Writing – Original Draft Preparation, Xiang Li and Yiwen Dong; Writing – Review & Editing, Yunfei Zhang and Zhen Shi.
Acknowledgments
This work is supported by the National Key R&D Program of China (2019YFC1408004), the National Natural Science Foundation of China (42205064) and the Scientific research project of China Three Gorges Corporation (Grant. Nos. 202103460). We also thank anonymous reviewers for their comments, which improved the paper.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Stull, R.B. An Introduction to Boundary Layer Meteorology; Springer Netherlands: 1988.
- Bakas, N.A.; Fotiadi, A.; Kariofillidi, S. Climatology of the Boundary Layer Height and of the Wind Field over Greece. Atmosphere 2020, 11. [Google Scholar] [CrossRef]
- Davies, F.; Middleton, D.R.; Bozier, K.E. Urban air pollution modelling and measurements of boundary layer height. Atmospheric Environment 2007, 41, 4040–4049. [Google Scholar] [CrossRef]
- Quan, J.; Gao, Y.; Zhang, Q.; Tie, X.; Cao, J.; Han, S.; Meng, J.; Chen, P.; Zhao, D. Evolution of planetary boundary layer under different weather conditions, and its impact on aerosol concentrations. Particuology 2013, 11, 34–40. [Google Scholar] [CrossRef]
- Yang, T.; Wang, Z.; Zhang, W.; Gbaguidi, A.; Sugimoto, N.; Wang, X.; Matsui, I.; Sun, Y. Technical note: Boundary layer height determination from lidar for improving air pollution episode modeling: Development of new algorithm and evaluation. Atmospheric Chemistry and Physics 2017, 17, 6215–6225. [Google Scholar] [CrossRef]
- Seidel, D.J.; Ao, C.O.; Li, K. Estimating climatological planetary boundary layer heights from radiosonde observations: Comparison of methods and uncertainty analysis. Journal of Geophysical Research 2010, 115. [Google Scholar] [CrossRef]
- Hennemuth, B.; Lammert, A. Determination of the Atmospheric Boundary Layer Height from Radiosonde and Lidar Backscatter. Boundary-Layer Meteorology 2006, 120, 181–200. [Google Scholar] [CrossRef]
- Tucker, S.C.; Senff, C.J.; Weickmann, A.M.; Brewer, W.A.; Banta, R.M.; Sandberg, S.P.; Law, D.C.; Hardesty, R.M. Doppler Lidar Estimation of Mixing Height Using Turbulence, Shear, and Aerosol Profiles. Journal of Atmospheric and Oceanic Technology 2009, 26, 673–688. [Google Scholar] [CrossRef]
- Angevine, W.M.; White, A.B.; Avery, S.K. Boundary-layer depth and entrainment zone characterization with a boundary-layer profiler. Boundary-Layer Meteorology 1994, 68, 375–385. [Google Scholar] [CrossRef]
- Bianco, L.; Wilczak, J.M. Convective Boundary Layer Depth: Improved Measurement by Doppler Radar Wind Profiler Using Fuzzy Logic Methods. Journal of Atmospheric and Oceanic Technology 2002, 19, 1745–1758. [Google Scholar] [CrossRef]
- Eresmaa, N.; Karppinen, A.; Joffre, S.M.; Räsänen, J.; Talvitie, H. Mixing height determination by ceilometer. Atmospheric Chemistry & Physics Discussions, 2005. [Google Scholar]
- Kamp, D.; Mckendry, I. Diurnal and Seasonal Trends in Convective Mixed-Layer Heights Estimated from Two Years of Continuous Ceilometer Observations in Vancouver, BC. Boundary-Layer Meteorology 2010, 137, 459–475. [Google Scholar] [CrossRef]
- Ratnam, M.V.; Basha, S.G. A robust method to determine global distribution of atmospheric boundary layer top from COSMIC GPS RO measurements. Atmospheric Science Letters 2010, 11, 216–222. [Google Scholar] [CrossRef]
- Ao, C.O.; Waliser, D.E.; Chan, S.K.; Li, J.L.; Tian, B.; Xie, F.; Mannucci, A.J. Planetary boundary layer heights from GPS radio occultation refractivity and humidity profiles. Chi O. Ao;Duane E. Waliser;Steven K. Chan;Jui-Lin Li;Baijun Tian;Feiqin Xie;Anthony J. Mannucci 2012, 117. [Google Scholar] [CrossRef]
- Su, T.; Li, J.; Li, C.; Xiang, P.; Lau, A.K.-H.; Guo, J.; Yang, D.; Miao, Y. An intercomparison of long-term planetary boundary layer heights retrieved from CALIPSO, ground-based lidar, and radiosonde measurements over Hong Kong. Journal of Geophysical Research: Atmospheres 2017, 122, 3929–3943. [Google Scholar] [CrossRef]
- Zhang, W.; Guo, J.; Miao, Y.; Liu, H.; Zhang, Y.; Li, Z.; Zhai, P. Planetary boundary layer height from CALIOP compared to radiosonde over China. Atmospheric Chemistry and Physics 2016, 16, 9951–9963. [Google Scholar] [CrossRef]
- Seidel, D.J.; Zhang, Y.; Beljaars, A.; Golaz, J.C.; Jacobson, A.R.; Medeiros, B. Climatology of the planetary boundary layer over the continental United States and Europe. Journal of Geophysical Research Atmospheres 2012, 117. [Google Scholar] [CrossRef]
- von Engeln, A.; Teixeira, J. A Planetary Boundary Layer Height Climatology Derived from ECMWF Reanalysis Data. Journal of Climate 2013, 26, 6575–6590. [Google Scholar] [CrossRef]
- Si, Y.; Li, S.; Chen, L.; Yu, C.; Wang, Z.; Wang, Y.; Wang, H. Validation and Spatiotemporal Distribution of GEOS-5–Based Planetary Boundary Layer Height and Relative Humidity in China. Advances in Atmospheric Sciences 2018, 35, 479–492. [Google Scholar] [CrossRef]
- Holzworth, G.C. Estimates of Mean Maximum Mixing Depths in the Contiguous United States. Monthly Weather Review 1964, 92, 235–242. [Google Scholar] [CrossRef]
- Liang, X.-Z.; Liu, S. Observed Diurnal Cycle Climatology of Planetary Boundary Layer Height. Journal of Climate 2010, 23, 5790–5809. [Google Scholar] [CrossRef]
- Zhang, Y.; Zhang, S.; Huang, C.; Huang, K.; Gong, Y.; Gan, Q. Diurnal variations of the planetary boundary layer height estimated from intensive radiosonde observations over Yichang, China. Science China Technological Sciences 2014, 57, 2172–2176. [Google Scholar] [CrossRef]
- Guo, J.; Miao, Y.; Zhang, Y.; Liu, H.; Li, Z.; Zhang, W.; He, J.; Lou, M.; Yan, Y.; Bian, L.; et al. The climatology of planetary boundary layer height in China derived from radiosonde and reanalysis data. Atmospheric Chemistry and Physics 2016, 16, 13309–13319. [Google Scholar] [CrossRef]
- Guo, J.; Li, Y.; Cohen, J.B.; Li, J.; Chen, D.; Xu, H.; Liu, L.; Yin, J.; Hu, K.; Zhai, P. Shift in the Temporal Trend of Boundary Layer Height in China Using Long-Term (1979–2016) Radiosonde Data. Geophysical Research Letters 2019, 46, 6080–6089. [Google Scholar] [CrossRef]
- Zhao, C.; Li, Y.; Liu, Y.; Zhou, G. The variation characteristics of planetary boundary layer height in Northwest China: Based on radiosonde and ERA-Interim reanalysis data. Plateau Meteorology 2019, 38, 1181–1193. [Google Scholar]
- ZHAO, Y.; MAO, W.; ZHANG, K.; MA, Y.; LI, H.; ZHANG, W. Climatic Variations in the Boundary Layer Height of Arid and Semiarid Areas in East Asia and North Africa. Journal of the Meteorological Society of Japan. Ser. II 2017, 95. [Google Scholar] [CrossRef]
- Allabakash, S.; Lim, S. Climatology of Planetary Boundary Layer Height-Controlling Meteorological Parameters Over the Korean Peninsula. Remote Sensing 2020, 12. [Google Scholar] [CrossRef]
- Zhang, Y.; Seidel, D.J.; Zhang, S. Trends in Planetary Boundary Layer Height over Europe. Journal of Climate 2013, 26, 10071–10076. [Google Scholar] [CrossRef]
- Li, J.; Chu, Y.; Li, X.; Dong, Y. Long-term trends of global maximum atmospheric mixed layer heights derived from radiosonde measurements. Environmental Research Letters 2020, 15. [Google Scholar] [CrossRef]
- Yang, D.; Li, C.; Lau, K.H.; Li, Y. Long-term measurement of daytime atmospheric mixing layer height over Hong Kong. Journal of Geophysical Research: Atmospheres 2013, 118. [Google Scholar] [CrossRef]
- Slingo, J.M. The Development and Verification of A Cloud Prediction Scheme For the Ecmwf Model. Quarterly Journal of the Royal Meteorological Society 1987, 113, 899–927. [Google Scholar] [CrossRef]
- Sicard, M.; Pérez, C.; Rocadenbosch, F.; Baldasano, J.M.; García-Vizcaino, D. Mixed-Layer Depth Determination in the Barcelona Coastal Area From Regular Lidar Measurements: Methods, Results and Limitations. Boundary-Layer Meteorology 2006, 119, 135–157. [Google Scholar] [CrossRef]
- Sen; Kumar, P. Estimates of the Regression Coefficient Based on Kendall's Tau. Publications of the American Statistical Association 1968, 63, 1379–1389. [Google Scholar] [CrossRef]
- Mann, H.B. Nonparametric test against trend. Econometrica 1945, 13, 245–259. [Google Scholar] [CrossRef]
- Kendall, M.G. Rank Correlation Methods. British Journal of Psychology 1990, 25, 86–91. [Google Scholar] [CrossRef]
- Meng, X.; Guo, J.; Han, Y. Preliminarily assessment of ERA5 reanalysis data. Journal of Marine Meteorology 2018, 38, 91–99. [Google Scholar]
- K. , N. Mixing depth model using hourly surface observations. Report 7053 1973. [Google Scholar]
- Gu, J.; Zhang, Y.; Yang, N.; Wang, R. Diurnal variability of the planetary boundary layer height estimated from radiosonde data. Earth and Planetary Physics 2020, 4, 1–14. [Google Scholar] [CrossRef]
- Guo, J.; Deng, M.; Lee, S.S.; Wang, F.; Li, Z.; Zhai, P.; Liu, H.; Lv, W.; Yao, W.; Li, X. Delaying precipitation and lightning by air pollution over the Pearl River Delta. Part I: Observational analyses. Journal of Geophysical Research, D. Atmospheres: JGR 2016, 121, 6472–6488. [Google Scholar] [CrossRef]
- Guo; Jianping; Tianning; Zhanqing; Miao; Yucong; Jing; Liu; Huan; Hui. Declining frequency of summertime local-scale precipitation over eastern China from 1970 to 2010 and its potential link to aerosols. Geophysical Research Letters 2017, 44, 5700–5708. [Google Scholar] [CrossRef]
- Zhao, X.; Zhang, X.; Xu, X.; Jing, X.; Meng, W.; Pu, W. Seasonal and diurnal variations of ambient PM2.5 concentration in urban and rural environments in Beijing. Atmospheric Environment 2009, 43, 2893–2900. [Google Scholar] [CrossRef]
- Hua, G.; Ming, X.; Qi, H. Changes in near-surface wind speed in China: 1969–2005. International Journal of Climatology 2011, 31. [Google Scholar]
- Mao, W.; Wang, Q.; Ma, H.; Shen, G. Temporal and Spatial Variation Characteristics of Meiyu in Jianghuai. Transactions of Atmospheric Sciences 2008, 31, 116–122. [Google Scholar]
- Xia, L.; Zhang, Q.; Sun, N.; Wang, S. Climate change characteristics in Jiangsu Province, 1960-2012. Journal of Glaciology and Geocryology 2015, 37, 1480–1489. [Google Scholar]
- Chen, J.; Zhong, S. Analysis of characteristics of relative humidity changes in Nanjing from 1961 to 2010. ENERGY RESEARCH & UTILIZATION, 2018, 45-49.
Figure 1.
Spatial distribution of the radiosonde sounding sites (red dots) and ERA5 grids (black dots) over Jiangsu, China.
Figure 1.
Spatial distribution of the radiosonde sounding sites (red dots) and ERA5 grids (black dots) over Jiangsu, China.
Figure 2.
Scatter plots of planetary boundary layer height (PBLH) derived from ERA5 and sounding datasets at (top) 08:00 LST and (bottom) 20:00 LST at the (a, d) Xuzhou, (b, e) Sheyang and (c, f) Nanjing sites. R is the correlation coefficient, the star superscript indicates the value that is statistically significant (p < 0.01), RMSE is the root mean square error, the black solid line is the fitting line, and the red solid line is the reference line.
Figure 2.
Scatter plots of planetary boundary layer height (PBLH) derived from ERA5 and sounding datasets at (top) 08:00 LST and (bottom) 20:00 LST at the (a, d) Xuzhou, (b, e) Sheyang and (c, f) Nanjing sites. R is the correlation coefficient, the star superscript indicates the value that is statistically significant (p < 0.01), RMSE is the root mean square error, the black solid line is the fitting line, and the red solid line is the reference line.
Figure 3.
Spatial distribution of the average PBLH over Jiangsu, China, from 1979 to 2019: (a) 08:00 LST, (b) 14:00 LST, (c) 20:00 LST, and (d) 02:00 LST.
Figure 3.
Spatial distribution of the average PBLH over Jiangsu, China, from 1979 to 2019: (a) 08:00 LST, (b) 14:00 LST, (c) 20:00 LST, and (d) 02:00 LST.
Figure 4.
Box-and-whisker plot of ERA5 PBLH during each season at 08:00, 14:00, 20:00 and 02:00 LST during the period from 1979 to 2019 over Jiangsu, China.
Figure 4.
Box-and-whisker plot of ERA5 PBLH during each season at 08:00, 14:00, 20:00 and 02:00 LST during the period from 1979 to 2019 over Jiangsu, China.
Figure 5.
Spatial distributions of the seasonal mean PBLH at 08:00 LST (a1-a4), 14:00 LST (b1-b4), 20:00 LST (c1-c4) and 02:00 LST (d1-d4) over Jiangsu, China.
Figure 5.
Spatial distributions of the seasonal mean PBLH at 08:00 LST (a1-a4), 14:00 LST (b1-b4), 20:00 LST (c1-c4) and 02:00 LST (d1-d4) over Jiangsu, China.
Figure 6.
Seasonal diurnal cycles of the mean PBLH throughout Jiangsu, China. The orange, red, blue and black lines indicate spring, summer, fall, and winter, respectively.
Figure 6.
Seasonal diurnal cycles of the mean PBLH throughout Jiangsu, China. The orange, red, blue and black lines indicate spring, summer, fall, and winter, respectively.
Figure 7.
(left) Amplitude and (right) phase of the diurnal cycle of the PBLH over Jiangsu, China, for each season. Amplitudes are mean values of daily maximum minus daily minimum, and phases are the hour of maximum PBLH.
Figure 7.
(left) Amplitude and (right) phase of the diurnal cycle of the PBLH over Jiangsu, China, for each season. Amplitudes are mean values of daily maximum minus daily minimum, and phases are the hour of maximum PBLH.
Figure 8.
Diurnal cycle of the average (a) 2 m temperature (T2m), (b) 10 m wind speed (WS), (c) relative humidity (RH) and (d) lower tropospheric stability (LTS) over Jiangsu, China.
Figure 8.
Diurnal cycle of the average (a) 2 m temperature (T2m), (b) 10 m wind speed (WS), (c) relative humidity (RH) and (d) lower tropospheric stability (LTS) over Jiangsu, China.
Figure 9.
(a) Time series of annual mean BLH (black line) and 5-year moving average BLH (blue line) throughout Jiangsu, China, at 14:00 LST for the 1979-2019 period. The red line denotes the linear fitting trend. (b) Spatial distribution of the trend of annual mean BLH over Jiangsu, China, at 14:00 LST during 1979-2019. Black dots indicate trends that are statistically significant (p < 0.05).
Figure 9.
(a) Time series of annual mean BLH (black line) and 5-year moving average BLH (blue line) throughout Jiangsu, China, at 14:00 LST for the 1979-2019 period. The red line denotes the linear fitting trend. (b) Spatial distribution of the trend of annual mean BLH over Jiangsu, China, at 14:00 LST during 1979-2019. Black dots indicate trends that are statistically significant (p < 0.05).
Figure 10.
Same as in
Figure 9(a), but for the time series of (a) T2m, (b) RH, and (c) LTS. R is the correlation coefficient between the annual mean PBLH and meteorological variables.
Figure 10.
Same as in
Figure 9(a), but for the time series of (a) T2m, (b) RH, and (c) LTS. R is the correlation coefficient between the annual mean PBLH and meteorological variables.
Figure 11.
Same as in
Figure 9(b) but for the spatial distribution of the trend of (a) T2m, (b) RH, and (c) LTS.
Figure 11.
Same as in
Figure 9(b) but for the spatial distribution of the trend of (a) T2m, (b) RH, and (c) LTS.
Figure 12.
Spatial distribution of correlation coefficients (R) between annual mean BLH and (a) T2m, (b) RH and (c) LTS during 1979-2019.
Figure 12.
Spatial distribution of correlation coefficients (R) between annual mean BLH and (a) T2m, (b) RH and (c) LTS during 1979-2019.
Table 1.
Correlation coefficients between the PBLH and meteorological variables (2 m temperature (T2m), 10 m wind speed (WS), relative humidity (RH) and lower tropospheric stability (LTS)) based on the diurnal cycle in the four seasons over Jiangsu, China. Bold figures are significant at the 95% confidence level or greater.
Table 1.
Correlation coefficients between the PBLH and meteorological variables (2 m temperature (T2m), 10 m wind speed (WS), relative humidity (RH) and lower tropospheric stability (LTS)) based on the diurnal cycle in the four seasons over Jiangsu, China. Bold figures are significant at the 95% confidence level or greater.
|
T2m |
WS |
RH |
LTS |
Spring |
0.92 |
0.91 |
-0.83 |
-0.93 |
Summer |
0.94 |
0.93 |
-0.92 |
-0.95 |
Fall |
0.92 |
0.92 |
-0.70 |
-0.93 |
Winter |
0.89 |
0.88 |
-0.57 |
-0.92 |
|
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. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).