The“Outline for High-Quality Development of Meteorology(2022-2035)”was published by the State Council of the People’s Republic of China in April 2022(1).This publication highlights the essential role of meteorologi...The“Outline for High-Quality Development of Meteorology(2022-2035)”was published by the State Council of the People’s Republic of China in April 2022(1).This publication highlights the essential role of meteorological disaster warnings,particularly cold spell warnings,in preventing and mitigating disasters.The outline also sets forth stricter expectations for enhancing the overall capacity of society in terms of meteorological disaster preparedness and response.展开更多
Sea surface temperature(SST)is one of the important parameters of global ocean and climate research,which can be retrieved by satellite infrared and passive microwave remote sensing instruments.While satellite infrare...Sea surface temperature(SST)is one of the important parameters of global ocean and climate research,which can be retrieved by satellite infrared and passive microwave remote sensing instruments.While satellite infrared SST offers high spatial resolution,it is limited by cloud cover.On the other hand,passive microwave SST provides all-weather observation but suffers from poor spatial resolution and susceptibility to environmental factors such as rainfall,coastal effects,and high wind speeds.To achieve high-precision,comprehensive,and high-resolution SST data,it is essential to fuse infrared and microwave SST measurements.In this study,data from the Fengyun-3D(FY-3D)medium resolution spectral imager II(MERSI-II)SST and microwave imager(MWRI)SST were fused.Firstly,the accuracy of both MERSIII SST and MWRI SST was verified,and the latter was bilinearly interpolated to match the 5km resolution grid of MERSI SST.After pretreatment and quality control of MERSI SST and MWRI SST,a Piece-Wise Regression method was employed to correct biases in MWRI SST.Subsequently,SST data were selected based on spatial resolution and accuracy within a 3-day window of the analysis date.Finally,an optimal interpolation method was applied to fuse the FY-3D MERSI-II SST and MWRI SST.The results demonstrated a significant improvement in spatial coverage compared to MERSI-II SST and MWRI SST.Furthermore,the fusion SST retained true spatial distribution details and exhibited an accuracy of–0.12±0.74℃compared to OSTIA SST.This study has improved the accuracy of FY satellite fusion SST products in China.展开更多
Since the Beijing 2022 Winter Olympics was the first Winter Olympics in history held in continental winter monsoon climate conditions across complex terrain areas,there is a deficiency of relevant research,operational...Since the Beijing 2022 Winter Olympics was the first Winter Olympics in history held in continental winter monsoon climate conditions across complex terrain areas,there is a deficiency of relevant research,operational techniques,and experience.This made providing meteorological services for this event particularly challenging.The China Meteorological Administration(CMA)Earth System Modeling and Prediction Centre,achieved breakthroughs in research on short-and medium-term deterministic and ensemble numerical predictions.Several key technologies crucial for precise winter weather services during the Winter Olympics were developed.A comprehensive framework,known as the Operational System for High-Precision Weather Forecasting for the Winter Olympics,was established.Some of these advancements represent the highest level of capabilities currently available in China.The meteorological service provided to the Beijing 2022 Games also exceeded previous Winter Olympic Games in both variety and quality.This included achievements such as the“100-meter level,minute level”downscaled spatiotemporal resolution and forecasts spanning 1 to 15 days.Around 30 new technologies and over 60 kinds of products that align with the requirements of the Winter Olympics Organizing Committee were developed,and many of these techniques have since been integrated into the CMA’s operational national forecasting systems.These accomplishments were facilitated by a dedicated weather forecasting and research initiative,in conjunction with the preexisting real-time operational forecasting systems of the CMA.This program represents one of the five subprograms of the WMO’s high-impact weather forecasting demonstration project(SMART2022),and continues to play an important role in their Regional Association(RA)II Research Development Project(Hangzhou RDP).Therefore,the research accomplishments and meteorological service experiences from this program will be carried forward into forthcoming highimpact weather forecasting activities.This article provides an overview and assessment of this program and the operational national forecasting systems.展开更多
The continuous rainy precipitation process from February to March in 2019 was selected to analyze the effect of meteorological service in Tangpu Reservoir basin,so as to sum up service experience and then lay a better...The continuous rainy precipitation process from February to March in 2019 was selected to analyze the effect of meteorological service in Tangpu Reservoir basin,so as to sum up service experience and then lay a better foundation for subsequent services.In response to the rainy weather from December 2018 to early 2019,three rounds of flood discharge were carried out in Tangpu Reservoir.During February-March in 2019,the hit rate of short-term area rainfall forecast for Tangpu Reservoir was 80.0%.Compared with the median of forecast interval,the average absolute error was 7.6 mm,and the relative error was 32.7%.The large deviation in the forecast from March 27 to 28 was deeply analyzed,and it is found that the main reasons were excessive reliance on and trust in a single model,insufficient correction of the actual situation,and insufficient judgment of the nature of precipitation.For the future reservoir meteorological service,three aspects of thinking were put forward,such as further strengthening the sharing of hydrological and meteorological information,improving the forecasting ability,and deepening the research of runoff forecast models.展开更多
It is not only meteorological problems for the medium-range numerical weather prediction (NWP) research to be in operation,but also engineering and technological problems.Here we gener- ally described the results of r...It is not only meteorological problems for the medium-range numerical weather prediction (NWP) research to be in operation,but also engineering and technological problems.Here we gener- ally described the results of research,engineering construction,operation information and testing,in the course of set-up of medium-range NWP operation system in the China National Meteorological Center.展开更多
On January 1,2021,the newly revised Archives Law of the People's Republic of China(hereinafter referred to as the new Archives Law)was officially implemented,providing legal protection for archive management,archi...On January 1,2021,the newly revised Archives Law of the People's Republic of China(hereinafter referred to as the new Archives Law)was officially implemented,providing legal protection for archive management,archive work and archive undertaking in the new era.Starting from the relevant provisions of the new Archives Law,the collection of precious meteorological archives was discussed,and it is believed that only on the basis of clarifying the definition,collection responsibility subject,collection scope,ways and methods of precious meteorological archives can the relevant units be convenient to practice and obtain substantive collection effect.展开更多
CO_(2)is one of the most important greenhouse gases(GHGs)in the earth’s atmosphere.Since the industrial era,anthropogenic activities have emitted excessive quantities of GHGs into the atmosphere,resulting in climate ...CO_(2)is one of the most important greenhouse gases(GHGs)in the earth’s atmosphere.Since the industrial era,anthropogenic activities have emitted excessive quantities of GHGs into the atmosphere,resulting in climate warming since the 1950s and leading to an increased frequency of extreme weather and climate events.In 2020,China committed to striving for carbon neutrality by 2060.This commitment and China’s consequent actions will result in significant changes in global and regional anthropogenic carbon emissions and therefore require timely,comprehensive,and objective monitoring and verification support(MVS)systems.The MVS approach relies on the top-down assimilation and inversion of atmospheric CO_(2)concentrations,as recommended by the Intergovernmental Panel on Climate Change(IPCC)Inventory Guidelines in 2019.However,the regional high-resolution assimilation and inversion method is still in its initial stage of development.Here,we have constructed an inverse system for carbon sources and sinks at the kilometer level by coupling proper orthogonal decomposition(POD)with four-dimensional variational(4DVar)data assimilation based on the weather research and forecasting-greenhouse gas(WRF-GHG)model.Our China Carbon Monito ring and Verification Support at the Regional level(CCMVS-R)system can continuously assimilate information on atmospheric CO_(2)and other related information and realize the inversion of regional and local anthropogenic carbon emissions and natural terrestrial ecosystem carbon exchange.Atmospheric CO_(2)data were collected from six ground-based monito ring sites in Shanxi Province,China to verify the inversion effect of regio nal anthropogenic carbon emissions by setting ideal and real experiments using a two-layer nesting method(at 27 and 9 km).The uncertainty of the simulated atmospheric CO_(2)decreased significantly,with a root-mean-square error of CO_(2)concentration values between the ideal value and the simulated after assimilation was close to 0.The total anthropogenic carbon emissions in Shanxi Province in 2019 from the assimilated inversions were approximately 28.6%(17%-38%)higher than the mean of five emission inventories using the bottomup method,showing that the top-down CCMVS-R system can obtain more comprehensive information on anthropogenic carbon emissions.展开更多
Accurate radar quantitative precipitation estimation(QPE)plays an essential role in disaster prevention and mitigation.In this paper,two deep learning-based QPE networks including a single-parameter network and a mult...Accurate radar quantitative precipitation estimation(QPE)plays an essential role in disaster prevention and mitigation.In this paper,two deep learning-based QPE networks including a single-parameter network and a multi-parameter network are designed.Meanwhile,a self-defined loss function(SLF)is proposed during modeling.The dataset includes Shijiazhuang S-band dual polarimetric radar(CINRAD/SAD)data and rain gauge data within the radar’s 100-km detection range during the flood season of 2021 in North China.Considering that the specific propagation phase shift(KDP)has a roughly linear relationship with the precipitation intensity,KDP is set to 0.5°km^(-1 )as a threshold value to divide all the rain data(AR)into a heavy rain(HR)and light rain(LR)dataset.Subsequently,12 deep learning-based QPE models are trained according to the input radar parameters,the precipitation datasets,and whether an SLF was adopted,respectively.The results suggest that the effects of QPE after distinguishing rainfall intensity are better than those without distinguishing,and the effects of using SLF are better than those that used MSE as a loss function.A Z-R relationship and a ZH-KDP-R synthesis method are compared with deep learning-based QPE.The mean relative errors(MRE)of AR models using SLF are improved by 61.90%,51.21%,and 56.34%compared with the Z-R relational method,and by 38.63%,42.55%,and 47.49%compared with the synthesis method.Finally,the models are further evaluated in three precipitation processes,which manifest that the deep learning-based models have significant advantages over the traditional empirical formula methods.展开更多
Generalized Additive Models(GAMs)are widely employed in ecological research,serving as a powerful tool for ecologists to explore complex nonlinear relationships between a response variable and predictors.Nevertheless,...Generalized Additive Models(GAMs)are widely employed in ecological research,serving as a powerful tool for ecologists to explore complex nonlinear relationships between a response variable and predictors.Nevertheless,evaluating the relative importance of predictors with concurvity(analogous to collinearity)on response variables in GAMs remains a challenge.To address this challenge,we developed an R package named gam.hp.gam.hp calculates individual R^(2) values for predictors,based on the concept of'average shared variance',a method previously introduced for multiple regression and canonical analyses.Through these individual R^(2)s,which add up to the overall R^(2),researchers can evaluate the relative importance of each predictor within GAMs.We illustrate the utility of the gam.hp package by evaluating the relative importance of emission sources and meteorological factors in explaining ozone concentration variability in air quality data from London,UK.We believe that the gam.hp package will improve the interpretation of results obtained from GAMs.展开更多
Understanding the structure of tropical cyclone(TC)hydrometeors is crucial for detecting the changes in the distribution and intensity of precipitation.In this study,the GMI brightness temperature and cloud-dependent ...Understanding the structure of tropical cyclone(TC)hydrometeors is crucial for detecting the changes in the distribution and intensity of precipitation.In this study,the GMI brightness temperature and cloud-dependent 1DVAR algorithm were used to retrieve the hydrometeor profiles and surface rain rate of TC Nanmadol(2022).The Advanced Radiative Transfer Modeling System(ARMS)was used to calculate the Jacobian and degrees of freedom(△DOF)of cloud water,rainwater,and graupel for different channels of GMI in convective conditions.The retrieval results were compared with the Dual-frequency Precipitation Radar(DPR),GMI 2A,and IMERG products.It is shown that from all channels of GMI,rain water has the highest△DOF,at 1.72.According to the radiance Jacobian to atmospheric state variables,cloud water emission dominates its scattering.For rain water,the emission of channels 1–4 dominates scattering.Compared with the GMI 2A precipitation product,the 1DVAR precipitation rate has a higher correlation coefficient(0.713)with the IMERG product and can better reflect the location of TC precipitation.Near the TC eyewall,the highest radar echo top indicates strong convection.Near the melting layer where Ka-band attenuation is strong,the double frequency difference of DPR data reflects the location of the melting.The DPR drop size distribution(DSD)product shows that there is a significant increase in particle size below the melting layer in the spiral rain band.Thus,the particle size may be one of the main reasons for the smaller rain water below the melting layer retrieved from 1DVAR.展开更多
Atmospheric ammonia(NH_(3)) is a chemically active trace gas that plays an important role in the atmospheric environment and climate change. Satellite remote sensing is a powerful technique to monitor NH_(3) concentra...Atmospheric ammonia(NH_(3)) is a chemically active trace gas that plays an important role in the atmospheric environment and climate change. Satellite remote sensing is a powerful technique to monitor NH_(3) concentration based on the absorption lines of NH_(3) in the thermal infrared region. In this study, we establish a retrieval algorithm to derive the NH_(3)column from the Hyperspectral Infrared Atmospheric Sounder(HIRAS) onboard the Chinese Feng Yun(FY)-3D satellite and present the first atmospheric NH_(3) column global map observed by the HIRAS instrument. The HIRAS observations can well capture NH_(3) hotspots around the world, e.g., India, West Africa, and East China, where large NH_(3) emissions exist. The HIRAS NH_(3) columns are also compared to the space-based Infrared Atmospheric Sounding Interferometer(IASI)measurements, and we find that the two instruments observe a consistent NH_(3) global distribution, with correlation coefficient(R) values of 0.28–0.73. Finally, some remaining issues about the HIRAS NH_(3) retrieval are discussed.展开更多
Seasonal snow cover is a key global climate and hydrological system component drawing considerable attention due to glob-al warming conditions.However,the spatiotemporal snow cover patterns are challenging in western ...Seasonal snow cover is a key global climate and hydrological system component drawing considerable attention due to glob-al warming conditions.However,the spatiotemporal snow cover patterns are challenging in western Jilin,China due to natural condi-tions and sparse observation.Hence,this study investigated the spatiotemporal patterns of snow cover using fine-resolution passive mi-crowave(PMW)snow depth(SD)data from 1987 to 2018,and revealed the potential influence of climate factors on SD variations.The results indicated that the interannual range of SD was between 2.90 cm and 9.60 cm during the snowy winter seasons and the annual mean SD showed a slightly increasing trend(P>0.05)at a rate of 0.009 cm/yr.In snowmelt periods,the snow cover contributed to an increase in volumetric soil water,and the change in SD was significantly affected by air temperature.The correlation between SD and air temperature was negative,while the correlation between SD and precipitation was positive during December and March.In March,the correlation coefficient exceeded 0.5 in Zhenlai,Da’an,Qianan,and Qianguo counties.However,the SD and precipitation were neg-atively correlated over western Jilin in October,and several subregions presented a negative correlation between SD and precipitation in November and April.展开更多
A second rain belt sometimes occurs ahead of a frontal rain belt in the warm sector over coastal South China,leading to heavy precipitation.We examined the differences in the mesoscale characteristics and microphysics...A second rain belt sometimes occurs ahead of a frontal rain belt in the warm sector over coastal South China,leading to heavy precipitation.We examined the differences in the mesoscale characteristics and microphysics of the frontal and warm sector rain belts that occurred in South China on May 10–13,2022.The southern rain belt occurred in an environment with favorable mesoscale conditions but weak large-scale forcing.In contrast,the northern rain belt was related to low-level horizontal shear and the surface-level front.The interaction between the enhanced southeasterly winds and the rainfall-induced cold pool promoted the persistent growth of convection along the southern rain belt.The convective cell propagated east over the coastal area,where there was a large temperature gradient.The bow-shaped echo in this region may be closely related to the rear-inflow jet.By contrast,the initial convection of the northern rain belt was triggered along the front and the region of low-level horizontal shear,with mesoscale interactions between the enhanced warm-moist southeasterly airflow and the cold dome associated with the earlier rain.The terrain blocked the movement of the cold pool,resulting in the stagnation of the frontal convective cell at an early stage.Subsequently,a meso-γ-scale vortex formed during the rapid movement of the convective cell,corresponding to an enhancement of precipitation.The representative raindrop spectra for the southern rain belt were characterized by a greater number and higher density of raindrops than the northern rain belt,even though both resulted in comparable hourly rainfalls.These results help us better understand the characteristics of double rain belts over South China.展开更多
Tree radial growth can have significantly differ-ent responses to climate change depending on the environ-ment.To elucidate the effects of climate on radial growth and stable carbon isotope(δ^(13)C)fractionation of Q...Tree radial growth can have significantly differ-ent responses to climate change depending on the environ-ment.To elucidate the effects of climate on radial growth and stable carbon isotope(δ^(13)C)fractionation of Qing-hai spruce(Picea crassifolia),a widely distributed native conifer in northwestern China in different environments,we developed chronologies for tree-ring widths and δ^(13)C in trees on the southern and northern slopes of the Qilian Mountains,and analysed the relationship between these tree-ring variables and major climatic factors.Tree-ring widths were strongly influenced by climatic factors early in the growing season,and the radial growth in trees on the northern slopes was more sensitive to climate than in trees on the southern.Tree-ring δ^(13)C was more sensitive to climate than radial growth.δ^(13)C fractionation was mainly influenced by summer temperature and precipitation early in the growing season.Stomatal conductance more strongly limited stable carbon isotope fractionation in tree rings than photosynthetic rate did.The response between tree rings and climate in mountains gradually weakened as climate warmed.Changes in radial growth and stable carbon isotope fractionation of P.crassifolia in response to climate in the Qilian Mountains may be further complicated by continued climate change.展开更多
This study compares the summer atmospheric water cycle,including moisture sources and consumption,in the upstream,midstream,and downstream regions of the Yarlung Zangbo River Basin in the southern Tibetan Plateau.The ...This study compares the summer atmospheric water cycle,including moisture sources and consumption,in the upstream,midstream,and downstream regions of the Yarlung Zangbo River Basin in the southern Tibetan Plateau.The evolutions of moisture properties under the influence of the westerly and summer southerly monsoon are examined using 5-yr multi-source measurements and ERA5 reanalysis data.Note that moisture consumption in this study is associated with clouds,precipitation,and diabatic heating.Compared to the midstream and downstream regions,the upstream region has less moisture,clouds,and precipitation,where the moisture is brought by the westerly.In early August,the vertical wet advection over this region becomes enhanced and generates more high clouds and precipitation.The midstream region has moisture carried by the westerly in June and by the southerly monsoon from July to August.The higher vertical wet advection maximum here forms more high clouds,with a precipitation peak in early July.The downstream region is mainly affected by the southerly-driven wet advection.The rich moisture and strong vertical wet advection here produce the most clouds and precipitation among the three regions,with a precipitation peak in late June.The height of the maximum moisture condensation is different between the midstream region(325 hPa)and the other two regions(375 hPa),due to the higher upward motion maximum in the midstream region.The diabatic heating structures show that stratiform clouds dominate the upstream region,stratiform clouds and deep convection co-exist in the midstream region,and deep convection systems characterize the downstream region.展开更多
High spatiotemporal resolution radiances from the advanced imagers onboard the new generation of geostationary weather satellites provide a unique opportunity to evaluate the abilities of various reanalysis datasets t...High spatiotemporal resolution radiances from the advanced imagers onboard the new generation of geostationary weather satellites provide a unique opportunity to evaluate the abilities of various reanalysis datasets to depict multilayer tropospheric water vapor(WV),thereby enhancing our understanding of the deficiencies of WV in reanalysis datasets.Based on daily measurements from the Advanced Himawari Imager(AHI)onboard the Himawari-8 satellite in 2016,the bias features of multilayer WV from six reanalysis datasets over East Asia are thoroughly evaluated.The assessments show that wet biases exist in the upper troposphere in all six reanalysis datasets;in particular,these biases are much larger in summer.Overall,we find better depictions of WV in the middle troposphere than in the upper troposphere.The accuracy of WV in the ERA5 dataset is the highest,in terms of the bias magnitude,dispersion,and pattern similarity.The characteristics of the WV bias over the Tibetan Plateau are significantly different from those over other parts of East Asia.In addition,the reanalysis datasets all capture the shift of the subtropical high very well,with ERA5 performing better overall.展开更多
“Diurnal variation of CH4 at the surface from spring to winter.The time units are in local time(+8 h UTC).The error bar is 1σfor all the observed hourly mean data within that season at that local time.”in the capti...“Diurnal variation of CH4 at the surface from spring to winter.The time units are in local time(+8 h UTC).The error bar is 1σfor all the observed hourly mean data within that season at that local time.”in the caption of Fig.8 on Page 604 should be“Diurnal variation of CH4 at the surface from spring to winter.The time units are in UTC.The error bar is 1σfor all the observed hourly mean data within that season at that local time.”展开更多
Based on the lightning observation data from the Fengyun-4A(FY-4A)Lightning Mapping Imager(FY-4A/LMI)and the Lightning Imaging Sensor(LIS)on the International Space Station(ISS),we extract the“event”type data as the...Based on the lightning observation data from the Fengyun-4A(FY-4A)Lightning Mapping Imager(FY-4A/LMI)and the Lightning Imaging Sensor(LIS)on the International Space Station(ISS),we extract the“event”type data as the lightning detection results.These observations are then compared with the cloud-to-ground(CG)lightning observation data from the China Meteorological Administration.This study focuses on the characteristics of lightning activity in Southeast China,primarily in Jiangxi Province and its adjacent areas,from April to September,2017–2022.In addition,with the fifth-generation European Centre for Medium-Range Weather Forecasts reanalysis data,we further delved into the potential factors influencing the distribution and variations in lightning activity and their primary related factors.Our findings indicate that the lightning frequency and density of the FY-4A/LMI,ISS-LIS and CG data are higher in southern and central Jiangxi,central Fujian Province,and western and central Guangdong Province,while they tend to be lower in eastern Hunan Province.In general,the high-value areas of lightning density for the FY-4A/LMI are located in inland mountainous areas.The lower the latitude is,the higher the CG lightning density is.High-value areas of the CG lightning density are more likely to be located in eastern Fujian and southeastern Zhejiang Province.However,the high-value areas of lightning density for the ISS-LIS are more dispersed,with a scattered distribution in inland mountainous areas and along the coast of eastern Fujian.Thus,the mountainous terrain is closely related to the high-value areas of the lightning density.The locations of the high-value areas of the lightning density for the FY-4A/LMI correspond well with those for the CG observations,and the seasonal variations are also consistent.In contrast,the distribution of the high-value areas of the lightning density for the ISS-LIS is more dispersed.The positions of the peak frequency of the FY-4A/LMI lightning and CG lightning contrast with local altitudes,primarily located at lower altitudes or near mountainsides.K-index and convective available potential energy(CAPE)can better reflect the local boundary layer conditions,where the lightning density is higher and lightning seasonal variations are apparent.There are strong correlations in the annual variations between the dew-point temperature(Td)and CG lightning frequency,and the monthly variations of the dew-point temperature and CAPE are also strongly correlated with monthly variations of CG lightning,while they are weakly correlated with the lightning frequency for the FY-4A/LMI and ISS-LIS.This result reflects that the CAPE shows a remarkable effect on the CG lightning frequency during seasonal transitions.展开更多
The Tibet Plateau is one of the regions with the richest solar energy resources in the world.In the process of achieving carbon neutrality in China,the development and utilization of solar energy resources in the regi...The Tibet Plateau is one of the regions with the richest solar energy resources in the world.In the process of achieving carbon neutrality in China,the development and utilization of solar energy resources in the region will play an important role.In this study,the gridded solar resource data with 1km resolution in Tibet were obtained by spatial correction and downscaling of SMARTS model.On this basis,the spatial and temporal distribution characteristics of solar energy resources in the region in the past 30 years(1991–2020)are finely evaluated,and the annual global horizontal radiation resource is calculated.The results show that:1)The average annual global horizontal radiation amount in Tibet is 1816 kWh/m^(2).More than 60%of the area belongs to the“Most abundant”(GHI≥1750 kWh/m^(2))area of China’s solar energy resources category A,and nearly 40%belongs to the“Quite abundant”(1400≤GHI<1750)area of China’s solar energy resource category B.2)In space,the solar energy resources in Tibet increased gradually from north to south and from east to west.Lhasa,Central and Eastern Shigatse,Shannan,and Southwestern Ali are the most abundant cities,with a maximum annual radiation level of 2189 kWh/m2.3)In terms of time,the total horizontal radiation in Tibet was the highest in May and the lowest in December.74%of the total area belongs to the“Very stable”(R_(w)≥0.47)area of solar resource stability category A,and 26%belongs to the“stable”(0.36≤R_(w)<0.47)area of solar resource stability category B.Solar energy resources in the region show the characteristics of both strong and stable.Average solar energy resources in the region have shown a fluctuating downward trend over the past 30 years,with an average decline of about 12.86(kWh/m2)per decade.4)In terms of solar radiation resources reaching the earth’s surface,the theoretical total amount of annual horizontal radiation in Tibet is about 240.07 billion tons of standard coal or 222.91 billion kilowatts on average.展开更多
文摘The“Outline for High-Quality Development of Meteorology(2022-2035)”was published by the State Council of the People’s Republic of China in April 2022(1).This publication highlights the essential role of meteorological disaster warnings,particularly cold spell warnings,in preventing and mitigating disasters.The outline also sets forth stricter expectations for enhancing the overall capacity of society in terms of meteorological disaster preparedness and response.
文摘Sea surface temperature(SST)is one of the important parameters of global ocean and climate research,which can be retrieved by satellite infrared and passive microwave remote sensing instruments.While satellite infrared SST offers high spatial resolution,it is limited by cloud cover.On the other hand,passive microwave SST provides all-weather observation but suffers from poor spatial resolution and susceptibility to environmental factors such as rainfall,coastal effects,and high wind speeds.To achieve high-precision,comprehensive,and high-resolution SST data,it is essential to fuse infrared and microwave SST measurements.In this study,data from the Fengyun-3D(FY-3D)medium resolution spectral imager II(MERSI-II)SST and microwave imager(MWRI)SST were fused.Firstly,the accuracy of both MERSIII SST and MWRI SST was verified,and the latter was bilinearly interpolated to match the 5km resolution grid of MERSI SST.After pretreatment and quality control of MERSI SST and MWRI SST,a Piece-Wise Regression method was employed to correct biases in MWRI SST.Subsequently,SST data were selected based on spatial resolution and accuracy within a 3-day window of the analysis date.Finally,an optimal interpolation method was applied to fuse the FY-3D MERSI-II SST and MWRI SST.The results demonstrated a significant improvement in spatial coverage compared to MERSI-II SST and MWRI SST.Furthermore,the fusion SST retained true spatial distribution details and exhibited an accuracy of–0.12±0.74℃compared to OSTIA SST.This study has improved the accuracy of FY satellite fusion SST products in China.
基金This work was jointly supported by the National Natural Science Foundation of China(Grant Nos.41975137,42175012,and 41475097)the National Key Research and Development Program(Grant No.2018YFF0300103).
文摘Since the Beijing 2022 Winter Olympics was the first Winter Olympics in history held in continental winter monsoon climate conditions across complex terrain areas,there is a deficiency of relevant research,operational techniques,and experience.This made providing meteorological services for this event particularly challenging.The China Meteorological Administration(CMA)Earth System Modeling and Prediction Centre,achieved breakthroughs in research on short-and medium-term deterministic and ensemble numerical predictions.Several key technologies crucial for precise winter weather services during the Winter Olympics were developed.A comprehensive framework,known as the Operational System for High-Precision Weather Forecasting for the Winter Olympics,was established.Some of these advancements represent the highest level of capabilities currently available in China.The meteorological service provided to the Beijing 2022 Games also exceeded previous Winter Olympic Games in both variety and quality.This included achievements such as the“100-meter level,minute level”downscaled spatiotemporal resolution and forecasts spanning 1 to 15 days.Around 30 new technologies and over 60 kinds of products that align with the requirements of the Winter Olympics Organizing Committee were developed,and many of these techniques have since been integrated into the CMA’s operational national forecasting systems.These accomplishments were facilitated by a dedicated weather forecasting and research initiative,in conjunction with the preexisting real-time operational forecasting systems of the CMA.This program represents one of the five subprograms of the WMO’s high-impact weather forecasting demonstration project(SMART2022),and continues to play an important role in their Regional Association(RA)II Research Development Project(Hangzhou RDP).Therefore,the research accomplishments and meteorological service experiences from this program will be carried forward into forthcoming highimpact weather forecasting activities.This article provides an overview and assessment of this program and the operational national forecasting systems.
文摘The continuous rainy precipitation process from February to March in 2019 was selected to analyze the effect of meteorological service in Tangpu Reservoir basin,so as to sum up service experience and then lay a better foundation for subsequent services.In response to the rainy weather from December 2018 to early 2019,three rounds of flood discharge were carried out in Tangpu Reservoir.During February-March in 2019,the hit rate of short-term area rainfall forecast for Tangpu Reservoir was 80.0%.Compared with the median of forecast interval,the average absolute error was 7.6 mm,and the relative error was 32.7%.The large deviation in the forecast from March 27 to 28 was deeply analyzed,and it is found that the main reasons were excessive reliance on and trust in a single model,insufficient correction of the actual situation,and insufficient judgment of the nature of precipitation.For the future reservoir meteorological service,three aspects of thinking were put forward,such as further strengthening the sharing of hydrological and meteorological information,improving the forecasting ability,and deepening the research of runoff forecast models.
文摘It is not only meteorological problems for the medium-range numerical weather prediction (NWP) research to be in operation,but also engineering and technological problems.Here we gener- ally described the results of research,engineering construction,operation information and testing,in the course of set-up of medium-range NWP operation system in the China National Meteorological Center.
基金Supported by the Key Research and Development Planning Project of Anhui Province(2022h11020002).
文摘On January 1,2021,the newly revised Archives Law of the People's Republic of China(hereinafter referred to as the new Archives Law)was officially implemented,providing legal protection for archive management,archive work and archive undertaking in the new era.Starting from the relevant provisions of the new Archives Law,the collection of precious meteorological archives was discussed,and it is believed that only on the basis of clarifying the definition,collection responsibility subject,collection scope,ways and methods of precious meteorological archives can the relevant units be convenient to practice and obtain substantive collection effect.
基金supported by the General Project of Top-Design of Multi-Scale Nature-Social ModelsData Support and Decision Support System for NSFC Carbon Neutrality Major Project(42341202)the Basic Scientific Research Fund of the Chinese Academy of Meteorological Sciences(2021Z014)。
文摘CO_(2)is one of the most important greenhouse gases(GHGs)in the earth’s atmosphere.Since the industrial era,anthropogenic activities have emitted excessive quantities of GHGs into the atmosphere,resulting in climate warming since the 1950s and leading to an increased frequency of extreme weather and climate events.In 2020,China committed to striving for carbon neutrality by 2060.This commitment and China’s consequent actions will result in significant changes in global and regional anthropogenic carbon emissions and therefore require timely,comprehensive,and objective monitoring and verification support(MVS)systems.The MVS approach relies on the top-down assimilation and inversion of atmospheric CO_(2)concentrations,as recommended by the Intergovernmental Panel on Climate Change(IPCC)Inventory Guidelines in 2019.However,the regional high-resolution assimilation and inversion method is still in its initial stage of development.Here,we have constructed an inverse system for carbon sources and sinks at the kilometer level by coupling proper orthogonal decomposition(POD)with four-dimensional variational(4DVar)data assimilation based on the weather research and forecasting-greenhouse gas(WRF-GHG)model.Our China Carbon Monito ring and Verification Support at the Regional level(CCMVS-R)system can continuously assimilate information on atmospheric CO_(2)and other related information and realize the inversion of regional and local anthropogenic carbon emissions and natural terrestrial ecosystem carbon exchange.Atmospheric CO_(2)data were collected from six ground-based monito ring sites in Shanxi Province,China to verify the inversion effect of regio nal anthropogenic carbon emissions by setting ideal and real experiments using a two-layer nesting method(at 27 and 9 km).The uncertainty of the simulated atmospheric CO_(2)decreased significantly,with a root-mean-square error of CO_(2)concentration values between the ideal value and the simulated after assimilation was close to 0.The total anthropogenic carbon emissions in Shanxi Province in 2019 from the assimilated inversions were approximately 28.6%(17%-38%)higher than the mean of five emission inventories using the bottomup method,showing that the top-down CCMVS-R system can obtain more comprehensive information on anthropogenic carbon emissions.
基金supported by National Key R&D Program of China(Grant No.2022YFC3003903)the S&T Program of Hebei(Grant No.19275408D),the Key-Area Research and Development Program of Guangdong Province(Grant No.2020B1111200001)+1 种基金the Key Project of Monitoring,Early Warning and Prevention of Major Natural Disasters of China(Grant No.2019YFC1510304)the Joint Fund of Key Laboratory of Atmosphere Sounding,CMA,and the Research Centre on Meteorological Observation Engineering Technology,CMA(Grant No.U2021Z05).
文摘Accurate radar quantitative precipitation estimation(QPE)plays an essential role in disaster prevention and mitigation.In this paper,two deep learning-based QPE networks including a single-parameter network and a multi-parameter network are designed.Meanwhile,a self-defined loss function(SLF)is proposed during modeling.The dataset includes Shijiazhuang S-band dual polarimetric radar(CINRAD/SAD)data and rain gauge data within the radar’s 100-km detection range during the flood season of 2021 in North China.Considering that the specific propagation phase shift(KDP)has a roughly linear relationship with the precipitation intensity,KDP is set to 0.5°km^(-1 )as a threshold value to divide all the rain data(AR)into a heavy rain(HR)and light rain(LR)dataset.Subsequently,12 deep learning-based QPE models are trained according to the input radar parameters,the precipitation datasets,and whether an SLF was adopted,respectively.The results suggest that the effects of QPE after distinguishing rainfall intensity are better than those without distinguishing,and the effects of using SLF are better than those that used MSE as a loss function.A Z-R relationship and a ZH-KDP-R synthesis method are compared with deep learning-based QPE.The mean relative errors(MRE)of AR models using SLF are improved by 61.90%,51.21%,and 56.34%compared with the Z-R relational method,and by 38.63%,42.55%,and 47.49%compared with the synthesis method.Finally,the models are further evaluated in three precipitation processes,which manifest that the deep learning-based models have significant advantages over the traditional empirical formula methods.
基金supported by the National Natural Science Foundation of China (32271551)National Key Research and Development Program of China (2023YFF0805803)the Metasequoia funding of Nanjing Forestry University。
文摘Generalized Additive Models(GAMs)are widely employed in ecological research,serving as a powerful tool for ecologists to explore complex nonlinear relationships between a response variable and predictors.Nevertheless,evaluating the relative importance of predictors with concurvity(analogous to collinearity)on response variables in GAMs remains a challenge.To address this challenge,we developed an R package named gam.hp.gam.hp calculates individual R^(2) values for predictors,based on the concept of'average shared variance',a method previously introduced for multiple regression and canonical analyses.Through these individual R^(2)s,which add up to the overall R^(2),researchers can evaluate the relative importance of each predictor within GAMs.We illustrate the utility of the gam.hp package by evaluating the relative importance of emission sources and meteorological factors in explaining ozone concentration variability in air quality data from London,UK.We believe that the gam.hp package will improve the interpretation of results obtained from GAMs.
基金funded by the National Key Research and Development Program of China(Grant No.2022YFC3004202)the National Natural Science Foundation of China(Grant Nos.U2142212 and 42105136)。
文摘Understanding the structure of tropical cyclone(TC)hydrometeors is crucial for detecting the changes in the distribution and intensity of precipitation.In this study,the GMI brightness temperature and cloud-dependent 1DVAR algorithm were used to retrieve the hydrometeor profiles and surface rain rate of TC Nanmadol(2022).The Advanced Radiative Transfer Modeling System(ARMS)was used to calculate the Jacobian and degrees of freedom(△DOF)of cloud water,rainwater,and graupel for different channels of GMI in convective conditions.The retrieval results were compared with the Dual-frequency Precipitation Radar(DPR),GMI 2A,and IMERG products.It is shown that from all channels of GMI,rain water has the highest△DOF,at 1.72.According to the radiance Jacobian to atmospheric state variables,cloud water emission dominates its scattering.For rain water,the emission of channels 1–4 dominates scattering.Compared with the GMI 2A precipitation product,the 1DVAR precipitation rate has a higher correlation coefficient(0.713)with the IMERG product and can better reflect the location of TC precipitation.Near the TC eyewall,the highest radar echo top indicates strong convection.Near the melting layer where Ka-band attenuation is strong,the double frequency difference of DPR data reflects the location of the melting.The DPR drop size distribution(DSD)product shows that there is a significant increase in particle size below the melting layer in the spiral rain band.Thus,the particle size may be one of the main reasons for the smaller rain water below the melting layer retrieved from 1DVAR.
基金supported by the Feng Yun Application Pioneering Project (FY-APP-2022.0502)the National Natural Science Foundation of China (Grant No. 42205140)。
文摘Atmospheric ammonia(NH_(3)) is a chemically active trace gas that plays an important role in the atmospheric environment and climate change. Satellite remote sensing is a powerful technique to monitor NH_(3) concentration based on the absorption lines of NH_(3) in the thermal infrared region. In this study, we establish a retrieval algorithm to derive the NH_(3)column from the Hyperspectral Infrared Atmospheric Sounder(HIRAS) onboard the Chinese Feng Yun(FY)-3D satellite and present the first atmospheric NH_(3) column global map observed by the HIRAS instrument. The HIRAS observations can well capture NH_(3) hotspots around the world, e.g., India, West Africa, and East China, where large NH_(3) emissions exist. The HIRAS NH_(3) columns are also compared to the space-based Infrared Atmospheric Sounding Interferometer(IASI)measurements, and we find that the two instruments observe a consistent NH_(3) global distribution, with correlation coefficient(R) values of 0.28–0.73. Finally, some remaining issues about the HIRAS NH_(3) retrieval are discussed.
基金Under the auspices of the Strategic Priority Research Program of the Chinese Academy of Sciences(No.XDA28110502)Science and Technology Development Plan Project of Jilin Province(No.20220202035NC)+1 种基金National Natural Science Foundation of China(No.41871248)Changchun Science and Technology Development Plan Project(No.21ZY12)。
文摘Seasonal snow cover is a key global climate and hydrological system component drawing considerable attention due to glob-al warming conditions.However,the spatiotemporal snow cover patterns are challenging in western Jilin,China due to natural condi-tions and sparse observation.Hence,this study investigated the spatiotemporal patterns of snow cover using fine-resolution passive mi-crowave(PMW)snow depth(SD)data from 1987 to 2018,and revealed the potential influence of climate factors on SD variations.The results indicated that the interannual range of SD was between 2.90 cm and 9.60 cm during the snowy winter seasons and the annual mean SD showed a slightly increasing trend(P>0.05)at a rate of 0.009 cm/yr.In snowmelt periods,the snow cover contributed to an increase in volumetric soil water,and the change in SD was significantly affected by air temperature.The correlation between SD and air temperature was negative,while the correlation between SD and precipitation was positive during December and March.In March,the correlation coefficient exceeded 0.5 in Zhenlai,Da’an,Qianan,and Qianguo counties.However,the SD and precipitation were neg-atively correlated over western Jilin in October,and several subregions presented a negative correlation between SD and precipitation in November and April.
基金National Natural Science Foundation of China(41930972,52239006,41975001)。
文摘A second rain belt sometimes occurs ahead of a frontal rain belt in the warm sector over coastal South China,leading to heavy precipitation.We examined the differences in the mesoscale characteristics and microphysics of the frontal and warm sector rain belts that occurred in South China on May 10–13,2022.The southern rain belt occurred in an environment with favorable mesoscale conditions but weak large-scale forcing.In contrast,the northern rain belt was related to low-level horizontal shear and the surface-level front.The interaction between the enhanced southeasterly winds and the rainfall-induced cold pool promoted the persistent growth of convection along the southern rain belt.The convective cell propagated east over the coastal area,where there was a large temperature gradient.The bow-shaped echo in this region may be closely related to the rear-inflow jet.By contrast,the initial convection of the northern rain belt was triggered along the front and the region of low-level horizontal shear,with mesoscale interactions between the enhanced warm-moist southeasterly airflow and the cold dome associated with the earlier rain.The terrain blocked the movement of the cold pool,resulting in the stagnation of the frontal convective cell at an early stage.Subsequently,a meso-γ-scale vortex formed during the rapid movement of the convective cell,corresponding to an enhancement of precipitation.The representative raindrop spectra for the southern rain belt were characterized by a greater number and higher density of raindrops than the northern rain belt,even though both resulted in comparable hourly rainfalls.These results help us better understand the characteristics of double rain belts over South China.
基金supported by Basic Research Operating Expenses of the Central level Non-profit Research Institutes (IDM2022003)National Natural Science Foundation of China (42375054)+2 种基金Regional collaborative innovation project of Xinjiang (2021E01022,2022E01045)Young Meteorological Talent Program of China Meteorological Administration,Tianshan Talent Program of Xinjiang (2022TSYCCX0003)Youth Innovation Team of China Meteorological Administration (CMA2023QN08).
文摘Tree radial growth can have significantly differ-ent responses to climate change depending on the environ-ment.To elucidate the effects of climate on radial growth and stable carbon isotope(δ^(13)C)fractionation of Qing-hai spruce(Picea crassifolia),a widely distributed native conifer in northwestern China in different environments,we developed chronologies for tree-ring widths and δ^(13)C in trees on the southern and northern slopes of the Qilian Mountains,and analysed the relationship between these tree-ring variables and major climatic factors.Tree-ring widths were strongly influenced by climatic factors early in the growing season,and the radial growth in trees on the northern slopes was more sensitive to climate than in trees on the southern.Tree-ring δ^(13)C was more sensitive to climate than radial growth.δ^(13)C fractionation was mainly influenced by summer temperature and precipitation early in the growing season.Stomatal conductance more strongly limited stable carbon isotope fractionation in tree rings than photosynthetic rate did.The response between tree rings and climate in mountains gradually weakened as climate warmed.Changes in radial growth and stable carbon isotope fractionation of P.crassifolia in response to climate in the Qilian Mountains may be further complicated by continued climate change.
基金supported by The Second Tibetan Plateau Scientific Expedition and Research(STEP)program(2019QZKK0105)the National Natural Science Foundation of China(91437221,91837204).
文摘This study compares the summer atmospheric water cycle,including moisture sources and consumption,in the upstream,midstream,and downstream regions of the Yarlung Zangbo River Basin in the southern Tibetan Plateau.The evolutions of moisture properties under the influence of the westerly and summer southerly monsoon are examined using 5-yr multi-source measurements and ERA5 reanalysis data.Note that moisture consumption in this study is associated with clouds,precipitation,and diabatic heating.Compared to the midstream and downstream regions,the upstream region has less moisture,clouds,and precipitation,where the moisture is brought by the westerly.In early August,the vertical wet advection over this region becomes enhanced and generates more high clouds and precipitation.The midstream region has moisture carried by the westerly in June and by the southerly monsoon from July to August.The higher vertical wet advection maximum here forms more high clouds,with a precipitation peak in early July.The downstream region is mainly affected by the southerly-driven wet advection.The rich moisture and strong vertical wet advection here produce the most clouds and precipitation among the three regions,with a precipitation peak in late June.The height of the maximum moisture condensation is different between the midstream region(325 hPa)and the other two regions(375 hPa),due to the higher upward motion maximum in the midstream region.The diabatic heating structures show that stratiform clouds dominate the upstream region,stratiform clouds and deep convection co-exist in the midstream region,and deep convection systems characterize the downstream region.
基金partly supported by the National Natural Science Foundation of China(Grant Nos.41975020 and 41975031)(Jun LI)。
文摘High spatiotemporal resolution radiances from the advanced imagers onboard the new generation of geostationary weather satellites provide a unique opportunity to evaluate the abilities of various reanalysis datasets to depict multilayer tropospheric water vapor(WV),thereby enhancing our understanding of the deficiencies of WV in reanalysis datasets.Based on daily measurements from the Advanced Himawari Imager(AHI)onboard the Himawari-8 satellite in 2016,the bias features of multilayer WV from six reanalysis datasets over East Asia are thoroughly evaluated.The assessments show that wet biases exist in the upper troposphere in all six reanalysis datasets;in particular,these biases are much larger in summer.Overall,we find better depictions of WV in the middle troposphere than in the upper troposphere.The accuracy of WV in the ERA5 dataset is the highest,in terms of the bias magnitude,dispersion,and pattern similarity.The characteristics of the WV bias over the Tibetan Plateau are significantly different from those over other parts of East Asia.In addition,the reanalysis datasets all capture the shift of the subtropical high very well,with ERA5 performing better overall.
文摘“Diurnal variation of CH4 at the surface from spring to winter.The time units are in local time(+8 h UTC).The error bar is 1σfor all the observed hourly mean data within that season at that local time.”in the caption of Fig.8 on Page 604 should be“Diurnal variation of CH4 at the surface from spring to winter.The time units are in UTC.The error bar is 1σfor all the observed hourly mean data within that season at that local time.”
基金National Natural Science Foundation of China(42175014,42205137)Open Research Fund of Institute of Meteorological Technology Innovation,Nanjing(BJG202202)+3 种基金Joint Research Project of Typhoon Research,Shanghai Typhoon Institute,China Meteorological Administration(TFJJ202209)Innovation Development Project of China Meteorological Administration(CXFZ2023P001)Open Project of KLME&CIC-FEMD(KLME202311)Jiangxi MDIA-ASI Fund。
文摘Based on the lightning observation data from the Fengyun-4A(FY-4A)Lightning Mapping Imager(FY-4A/LMI)and the Lightning Imaging Sensor(LIS)on the International Space Station(ISS),we extract the“event”type data as the lightning detection results.These observations are then compared with the cloud-to-ground(CG)lightning observation data from the China Meteorological Administration.This study focuses on the characteristics of lightning activity in Southeast China,primarily in Jiangxi Province and its adjacent areas,from April to September,2017–2022.In addition,with the fifth-generation European Centre for Medium-Range Weather Forecasts reanalysis data,we further delved into the potential factors influencing the distribution and variations in lightning activity and their primary related factors.Our findings indicate that the lightning frequency and density of the FY-4A/LMI,ISS-LIS and CG data are higher in southern and central Jiangxi,central Fujian Province,and western and central Guangdong Province,while they tend to be lower in eastern Hunan Province.In general,the high-value areas of lightning density for the FY-4A/LMI are located in inland mountainous areas.The lower the latitude is,the higher the CG lightning density is.High-value areas of the CG lightning density are more likely to be located in eastern Fujian and southeastern Zhejiang Province.However,the high-value areas of lightning density for the ISS-LIS are more dispersed,with a scattered distribution in inland mountainous areas and along the coast of eastern Fujian.Thus,the mountainous terrain is closely related to the high-value areas of the lightning density.The locations of the high-value areas of the lightning density for the FY-4A/LMI correspond well with those for the CG observations,and the seasonal variations are also consistent.In contrast,the distribution of the high-value areas of the lightning density for the ISS-LIS is more dispersed.The positions of the peak frequency of the FY-4A/LMI lightning and CG lightning contrast with local altitudes,primarily located at lower altitudes or near mountainsides.K-index and convective available potential energy(CAPE)can better reflect the local boundary layer conditions,where the lightning density is higher and lightning seasonal variations are apparent.There are strong correlations in the annual variations between the dew-point temperature(Td)and CG lightning frequency,and the monthly variations of the dew-point temperature and CAPE are also strongly correlated with monthly variations of CG lightning,while they are weakly correlated with the lightning frequency for the FY-4A/LMI and ISS-LIS.This result reflects that the CAPE shows a remarkable effect on the CG lightning frequency during seasonal transitions.
基金This work was supported by the Major Science and Technology Project of the Science and Technology Department of Tibet under Grant Number XZ202101ZD0015Gthe Second Tibet Plateau Scientific Expedition and Research Program(STEP)under Grant Number 2019QZKK0804.
文摘The Tibet Plateau is one of the regions with the richest solar energy resources in the world.In the process of achieving carbon neutrality in China,the development and utilization of solar energy resources in the region will play an important role.In this study,the gridded solar resource data with 1km resolution in Tibet were obtained by spatial correction and downscaling of SMARTS model.On this basis,the spatial and temporal distribution characteristics of solar energy resources in the region in the past 30 years(1991–2020)are finely evaluated,and the annual global horizontal radiation resource is calculated.The results show that:1)The average annual global horizontal radiation amount in Tibet is 1816 kWh/m^(2).More than 60%of the area belongs to the“Most abundant”(GHI≥1750 kWh/m^(2))area of China’s solar energy resources category A,and nearly 40%belongs to the“Quite abundant”(1400≤GHI<1750)area of China’s solar energy resource category B.2)In space,the solar energy resources in Tibet increased gradually from north to south and from east to west.Lhasa,Central and Eastern Shigatse,Shannan,and Southwestern Ali are the most abundant cities,with a maximum annual radiation level of 2189 kWh/m2.3)In terms of time,the total horizontal radiation in Tibet was the highest in May and the lowest in December.74%of the total area belongs to the“Very stable”(R_(w)≥0.47)area of solar resource stability category A,and 26%belongs to the“stable”(0.36≤R_(w)<0.47)area of solar resource stability category B.Solar energy resources in the region show the characteristics of both strong and stable.Average solar energy resources in the region have shown a fluctuating downward trend over the past 30 years,with an average decline of about 12.86(kWh/m2)per decade.4)In terms of solar radiation resources reaching the earth’s surface,the theoretical total amount of annual horizontal radiation in Tibet is about 240.07 billion tons of standard coal or 222.91 billion kilowatts on average.
基金supported by the Strategic Priority Research Program of Chinese Academy of Sciences[grant numbers XDA23090102]the National Natural Science Foundation of China[grant numbers 42175078 and 42075040]+1 种基金the Health Meteorological Project of Hebei Province[grant number FW202150]the National Key Research and Development Program of China[grant number 2018YFA0606203].