A cloud clustering and classification algorithm is developed for a ground-based Ka-band radar system in the vertically pointing mode. Cloud profiles are grouped based on the combination of a time–height clustering me...A cloud clustering and classification algorithm is developed for a ground-based Ka-band radar system in the vertically pointing mode. Cloud profiles are grouped based on the combination of a time–height clustering method and the k-means clustering method. The cloud classification algorithm, developed using a fuzzy logic method, uses nine physical parameters to classify clouds into nine types: cirrostratus, cirrocumulus, altocumulus, altostratus, stratus, stratocumulus, nimbostratus,cumulus or cumulonimbus. The performance of the clustering and classification algorithm is presented by comparison with all-sky images taken from January to June 2014. Overall, 92% of the cloud profiles are clustered successfully and the agreement in classification between the radar system and the all-sky imager is 87%. The distribution of cloud types in Beijing from January 2014 to December 2017 is studied based on the clustering and classification algorithm. The statistics show that cirrostratus clouds have the highest occurrence frequency(24%) among the nine cloud types. High-level clouds have the maximum occurrence frequency and low-level clouds the minimum occurrence frequency.展开更多
The automatic all-sky imager developed by the Institute of Atmospheric Physics,Chinese Academy of Sciences,provides all-sky visible images in the red,green,and blue channels.This paper presents three major cali-bratio...The automatic all-sky imager developed by the Institute of Atmospheric Physics,Chinese Academy of Sciences,provides all-sky visible images in the red,green,and blue channels.This paper presents three major cali-bration experiments of the all-sky imager,geometric an-gular calibration,optical calibration,and radiometric calibration,and then infers an algorithm to retrieve rela-tive radiance from the all-sky images.Field experiments show that the related coefficient between retrieved radi-ance and measured radiance is about 0.91.It is feasible to use the algorithm to retrieve radiance from images.The paper sets up a relationship between radiance and the im-age,which is useful for using the all-sky image in nu-merical-simulations that predict more meteorological pa-rameters.展开更多
In the present reported study, the vertical distributions of local atmospheric refractivity were retrieved from ground- based GPS observations at low elevation angles. An improved optimization method was implemented a...In the present reported study, the vertical distributions of local atmospheric refractivity were retrieved from ground- based GPS observations at low elevation angles. An improved optimization method was implemented at altitudes of 0-10 km to search for a best-fit refractivity profile that resulted in atmospheric delays most similar to the delays calculated from the observations. A ray-tracing model was used to simulate neutral atmospheric delays corresponding to a given refractivity profile. We initially performed a ＂theoretical retrieval＂, in which no observation data were involved, to verify the optimization method. A statistical relative error of this ＂theoretical retrieval＂ （-2% to 2%） indicated that such a retrieval is effective. In a practical retrieval, observations were obtained using a dual-frequency GPS receiver, and its initial value was provided by CIRA86aQ_UoG data. The statistical relative errors of the practical retrieval range from -3% to 5% were compared with co-located radiosonde measurements, Results clearly revealed diurnal variations in local refractivity prc,files, The results also suggest that the general vertical distribution of refractivity can be derived with a high temporal resolution. However, further study is needed to describe the vertical refractivity gradient clearly.展开更多
The 10-year climatological features related to midlatitude double tropopause events （DTs） are examined using ERA- Interim data from 2003 to 2012. The analysis is based on tropopauses defined by lapse rate. Results s...The 10-year climatological features related to midlatitude double tropopause events （DTs） are examined using ERA- Interim data from 2003 to 2012. The analysis is based on tropopauses defined by lapse rate. Results show that DTs are permanent or semi-permanent in the midlatitudes, and high DT frequency bands move poleward in winter and equatorward in summer, which is consistent with the seasonal movement of the subtropical jet. Based on our statistics, the second tropopause is found at about 100 hPa in the subtropics and at slightly lower altitudes in sub-polar regions. The thickness between the first and second tropopause is smaller in the subtropics and increases with latitude. Next, the origin of air sandwiched between the first and second tropopause of DTs is studied with a revised version of the UK Universities Global Atmospheric Modelling Programme Offline Trajectory Code （Version 3） diabatic trajectory model. The results show that, in the lower or middle troposphere, air is transported into the DTs from lower latitudes, mainly in the tropics. The dominant source regions are mainly areas of deep convection and steep orography, e.g., the western Pacific and Himalayan Mountains, and they show strong seasonality following the seasonal shift of these strong upwelling regions.展开更多
Previous studies have shown that year-to-year incremental prediction （YIP） can obtain considerable skill in seasonal forecasts. This study analyzes the mathematical deRnition of YiP and derives its formula in the no...Previous studies have shown that year-to-year incremental prediction （YIP） can obtain considerable skill in seasonal forecasts. This study analyzes the mathematical deRnition of YiP and derives its formula in the nonlinear time series prediction （NP） method, it is shown that the two methods are equivalent when the prediction time series is embedded in one-dimensional phase space. Compared to previous NP models, the new one introduces multiple external forcings in the form of year-to-year increments. The year-to-year increments have physical meaning, which is better than the NP model with empirically chosen parameters. The summer rainfall over the middle to lower reaches of the Yangtze River is analyzed to examine the prediction skill of the NP models. Results show that the NP model with year-to-year increments can reach a similar skill as the YiP model. When the embedded number of dimensions is increased to two, more accurate prediction can be obtained. Besides similar results, the NP method has more dynamical meaning, as it is based on the classical reconstruction theory. Moreover, by choosing different embedded dimensions, the NP model can reconstruct the dynamical curve into phase space with more than one dimension, which is an advantage of the NP model. The present study suggests that YIP has a robust dynamical foundation, besides its physical mechanism, and the modified NP model has the potential to increase the operationaJ skill in short- term climate prediction.展开更多
The relationship between the radiance ratio （radiance at wavelength 450 nm to 650 nm） and aerosol optical depth （AOD） is analyzed in this paper by numerical simulation and a ＂LUT＂ （look-up table） approach is t...The relationship between the radiance ratio （radiance at wavelength 450 nm to 650 nm） and aerosol optical depth （AOD） is analyzed in this paper by numerical simulation and a ＂LUT＂ （look-up table） approach is then presented for the retrieval of AOD from the radiance ratio. In this LUT approach, the AOD retrieval error depends mainly on the assumption of aerosol types. From the preliminary simulation, a typical error of 15%–20% in AOD obtained from the radiance ratio is estimated, due to neglecting changes in the ground albedo and background atmosphere. At its worst, the AOD error reached a maximum of around 50%, which will be refined in the future. In the latter part of the paper, comparisons are made between AOD from the imager and from the CE-318 sun photometer, both located at Xianghe observatory in Hebei Province （39.75° N, 116.96° E）. This field experiment shows that AOD from the imager is highly correlated with AOD from the sun photometer, with a correlation coefficient of 0.95 and an average retrieval error of around 7%. A contrast experiment confirms the feasibility of retrieving AOD from all-sky images, but more analysis and future research are required to improve the accuracy.展开更多
In order to provide a long time-series,high spatial resolution,and high accuracy dataset of land surface temperature（LST） for climatic change research,a modified Becker and Li＇s split-window approach is proposed in...In order to provide a long time-series,high spatial resolution,and high accuracy dataset of land surface temperature（LST） for climatic change research,a modified Becker and Li＇s split-window approach is proposed in this paper to retrieve LST from the measurements of Advanced Very High Resolution Radiometer（AVHRR） onboard National Oceanic and Atmospheric Administration（NOAA）-7 to-18 and the Visible and InfraRed Radiometer（VIRR） onboard FY-3A.For this purpose,the Moderate Resolution Transmittance Model（MODTRAN） 4.1 was first employed to compute the spectral radiance at the top of atmosphere（TOA） under a variety of surface and atmosphere conditions.Then,a temperature dataset consists of boundary temperature T s（which is one of the input parameters to MODTRAN）,and channels 4 and 5 brightness temperatures（T 4 and T 5） were constructed.Note that channels 4 and 5 brightness temperatures were simulated from the MODTRAN output spectral radiance by convolving them with the spectral response functions（SRFs） of channels 4 and 5 of AVHRRs and VIRR.The coefficients of modified Becker and Li＇s split-window approach for various AVHRRs and VIRR were subsequently regressed based on this temperature dataset using the least square method.As an example of validation,one AVHRR satellite image over Beijing acquired at 0312 UTC 27 April 2008 by AVHRR onboard NOAA-17 was selected to retrieve the LST image using the modified Becker and Li＇s approach.The comparison between this LST image and that from the MODIS level-2 LST product provided by the University of Tokyo in Japan indicates that the correlation coefficient is 0.88,the bias is 0.6 K,and the root mean square deviation（RMSD） is 2.1 K.Furthermore,about 70% and 37% pixels in the LST difference image,which is the result of retrieved LST image from AVHRR minus the corresponding MODIS LST image,have the values within ± 2 and ± 1 K,respectively.展开更多
It has been several years since the Greenhouse Gases Observing Satellite （GOSAT） began to observe the distribution of CO2 and CH4 over the globe from space. Results from Thermal and Near-infrared Sensor for Carbon O...It has been several years since the Greenhouse Gases Observing Satellite （GOSAT） began to observe the distribution of CO2 and CH4 over the globe from space. Results from Thermal and Near-infrared Sensor for Carbon Observation-Cloud and Aerosol Imager （TANSO-CAI） cloud screening are necessary for the retrieval of CO2 and CH4 gas concentrations for GOSAT TANSO-Fourier Transform Spectrometer （FTS） observations. In this study, TANSO-CAI cloud flag data were compared with ground-based cloud data collected by an all-sky imager （ASI） over Beijing from June 2009 to May 2012 to examine the data quality. The results showed that the CAI has an obvious cloudy tendency bias over Beijing, especially in winter. The main reason might be that heavy aerosols in the sky are incorrectly determined as cloudy pixels by the CAI algorithm. Results also showed that the CAI algorithm sometimes neglects some high thin cirrus cloud over this area.展开更多
The radiances scattered or emitted by clouds demonstrate diverse features at different wavelengths due to different cloud physical structures. This paper presents a method(the smallest-radiance-distance method, SRaDM)...The radiances scattered or emitted by clouds demonstrate diverse features at different wavelengths due to different cloud physical structures. This paper presents a method(the smallest-radiance-distance method, SRaDM) of revealing the physical structures of clouds. The method is based on multi-spectral radiances measured by the Moderate Resolution Imaging Spectroradiometer(MODIS)onboard Aqua. The principle and methodology of SRaDM is deduced and provided in this paper. Correlation analysis based on data from MODIS and Cloud Profiling Radar(onboard CloudSat), collected from January 2007 to December 2010 over an ocean area(15°N–45°N,145°E–165°E), led to selection of radiances at 13 wavebands of MODIS that demonstrated high sensitivity to cloud physical structures;radiances at the selected wavebands were subjected to SRaDM. The Standardized Euclidean distance is introduced to quantify the degree of changes in multi-spectral radiances(termed D_(rd)) and in physical structures(termed D_(st)) between cloud profiles. Statistics based on numerous cloud profiles show that D_(rd) decreases monotonically with a decrease in D_(st), which implies that small D_(rd) always accompanies small D_(st). According to the law of D_(rd) and D_(st), the new method, SRaDM, for revealing physical structures of clouds from the collocation of cloud profiles of similar multi-spectral radiances, is presented. Then, two successful experiments are presented in which cloud physical structures are captured using multi-spectral radiances. SRaDM provides a way to obtain knowledge of the physical structures of clouds over relatively larger areas, and is a new approach to obtaining 3D cloud fields.展开更多
In this study,the authors focus on the cut-off low pressure systems（COLs）lingering over East Asia in late spring and early summer and quantify the two-way stratosphere–troposphere exchange（STE）by 3D trajectory in...In this study,the authors focus on the cut-off low pressure systems（COLs）lingering over East Asia in late spring and early summer and quantify the two-way stratosphere–troposphere exchange（STE）by 3D trajectory integrations,achieved using a revised version of the UK Universities Global Atmospheric Modelling Programme Offline Trajectory Code（Version 3）.By selecting 10 typical COLs and calculating the cross-tropopause air mass fluxes,it is found that stratosphere-to-troposphere transport（STT）fluxes exist in the center of COLs;and in the periphery of the COL center,troposphereto-stratosphere transport（TST）fluxes and STT fluxes are distributed alternately.Net transport fluxes in COLs are from stratosphere to troposphere,and the magnitude is about 10-4 kg m-2 s-1.The ratio between the area-averaged STT and TST fluxes increases with increasing strength of the COLs.By adopting appropriate residence time,the spurious transports are effectively excluded.Finally,the authors compare the results with previous studies,and find that the cross-tropopause fluxes（CTFs）induced by COLs are about one to two orders of magnitude larger than global CTFs.COLs play a significant role in local,rapid air mass exchanges,although they may only be responsible for a fraction of the total STE.展开更多
基金supported by the National Natural Science Foundation of China (Grant Nos. 41775032 and 41275040)
文摘A cloud clustering and classification algorithm is developed for a ground-based Ka-band radar system in the vertically pointing mode. Cloud profiles are grouped based on the combination of a time–height clustering method and the k-means clustering method. The cloud classification algorithm, developed using a fuzzy logic method, uses nine physical parameters to classify clouds into nine types: cirrostratus, cirrocumulus, altocumulus, altostratus, stratus, stratocumulus, nimbostratus,cumulus or cumulonimbus. The performance of the clustering and classification algorithm is presented by comparison with all-sky images taken from January to June 2014. Overall, 92% of the cloud profiles are clustered successfully and the agreement in classification between the radar system and the all-sky imager is 87%. The distribution of cloud types in Beijing from January 2014 to December 2017 is studied based on the clustering and classification algorithm. The statistics show that cirrostratus clouds have the highest occurrence frequency(24%) among the nine cloud types. High-level clouds have the maximum occurrence frequency and low-level clouds the minimum occurrence frequency.
基金supported by the national natural science foundation of China (Grant Nos.40505006 and 40775026)
文摘The automatic all-sky imager developed by the Institute of Atmospheric Physics,Chinese Academy of Sciences,provides all-sky visible images in the red,green,and blue channels.This paper presents three major cali-bration experiments of the all-sky imager,geometric an-gular calibration,optical calibration,and radiometric calibration,and then infers an algorithm to retrieve rela-tive radiance from the all-sky images.Field experiments show that the related coefficient between retrieved radi-ance and measured radiance is about 0.91.It is feasible to use the algorithm to retrieve radiance from images.The paper sets up a relationship between radiance and the im-age,which is useful for using the all-sky image in nu-merical-simulations that predict more meteorological pa-rameters.
基金supported by the National Basic Research Program of China (973 Program) under Grant No. 2010CB428601the Open Fund of the State Key Laboratory of Satellite Ocean Environment Dynamics under Contract No. SOED0705the China Postdoctoral Science Foundation
文摘In the present reported study, the vertical distributions of local atmospheric refractivity were retrieved from ground- based GPS observations at low elevation angles. An improved optimization method was implemented at altitudes of 0-10 km to search for a best-fit refractivity profile that resulted in atmospheric delays most similar to the delays calculated from the observations. A ray-tracing model was used to simulate neutral atmospheric delays corresponding to a given refractivity profile. We initially performed a ＂theoretical retrieval＂, in which no observation data were involved, to verify the optimization method. A statistical relative error of this ＂theoretical retrieval＂ （-2% to 2%） indicated that such a retrieval is effective. In a practical retrieval, observations were obtained using a dual-frequency GPS receiver, and its initial value was provided by CIRA86aQ_UoG data. The statistical relative errors of the practical retrieval range from -3% to 5% were compared with co-located radiosonde measurements, Results clearly revealed diurnal variations in local refractivity prc,files, The results also suggest that the general vertical distribution of refractivity can be derived with a high temporal resolution. However, further study is needed to describe the vertical refractivity gradient clearly.
基金supported by the Special Fund for Strategic Pilot Technology of the Chinese Academy of Sciences(Grant No.XDA05040300)
文摘The 10-year climatological features related to midlatitude double tropopause events （DTs） are examined using ERA- Interim data from 2003 to 2012. The analysis is based on tropopauses defined by lapse rate. Results show that DTs are permanent or semi-permanent in the midlatitudes, and high DT frequency bands move poleward in winter and equatorward in summer, which is consistent with the seasonal movement of the subtropical jet. Based on our statistics, the second tropopause is found at about 100 hPa in the subtropics and at slightly lower altitudes in sub-polar regions. The thickness between the first and second tropopause is smaller in the subtropics and increases with latitude. Next, the origin of air sandwiched between the first and second tropopause of DTs is studied with a revised version of the UK Universities Global Atmospheric Modelling Programme Offline Trajectory Code （Version 3） diabatic trajectory model. The results show that, in the lower or middle troposphere, air is transported into the DTs from lower latitudes, mainly in the tropics. The dominant source regions are mainly areas of deep convection and steep orography, e.g., the western Pacific and Himalayan Mountains, and they show strong seasonality following the seasonal shift of these strong upwelling regions.
基金supported by the National Natural Sciences Foundation of China[41375112],[41530426],[41575058]the Key Technology Talent Program of the Chinese Academy of Sciencesthe Public Science and Technology Research Funds Projects of Ocean[201505013]
文摘Previous studies have shown that year-to-year incremental prediction （YIP） can obtain considerable skill in seasonal forecasts. This study analyzes the mathematical deRnition of YiP and derives its formula in the nonlinear time series prediction （NP） method, it is shown that the two methods are equivalent when the prediction time series is embedded in one-dimensional phase space. Compared to previous NP models, the new one introduces multiple external forcings in the form of year-to-year increments. The year-to-year increments have physical meaning, which is better than the NP model with empirically chosen parameters. The summer rainfall over the middle to lower reaches of the Yangtze River is analyzed to examine the prediction skill of the NP models. Results show that the NP model with year-to-year increments can reach a similar skill as the YiP model. When the embedded number of dimensions is increased to two, more accurate prediction can be obtained. Besides similar results, the NP method has more dynamical meaning, as it is based on the classical reconstruction theory. Moreover, by choosing different embedded dimensions, the NP model can reconstruct the dynamical curve into phase space with more than one dimension, which is an advantage of the NP model. The present study suggests that YIP has a robust dynamical foundation, besides its physical mechanism, and the modified NP model has the potential to increase the operationaJ skill in short- term climate prediction.
基金supported bythe National Natural Science Foundation of China (GrantNo. 40500506)
文摘The relationship between the radiance ratio （radiance at wavelength 450 nm to 650 nm） and aerosol optical depth （AOD） is analyzed in this paper by numerical simulation and a ＂LUT＂ （look-up table） approach is then presented for the retrieval of AOD from the radiance ratio. In this LUT approach, the AOD retrieval error depends mainly on the assumption of aerosol types. From the preliminary simulation, a typical error of 15%–20% in AOD obtained from the radiance ratio is estimated, due to neglecting changes in the ground albedo and background atmosphere. At its worst, the AOD error reached a maximum of around 50%, which will be refined in the future. In the latter part of the paper, comparisons are made between AOD from the imager and from the CE-318 sun photometer, both located at Xianghe observatory in Hebei Province （39.75° N, 116.96° E）. This field experiment shows that AOD from the imager is highly correlated with AOD from the sun photometer, with a correlation coefficient of 0.95 and an average retrieval error of around 7%. A contrast experiment confirms the feasibility of retrieving AOD from all-sky images, but more analysis and future research are required to improve the accuracy.
基金Supported by the National Science and Technology Special Funds for Infrastructure Work Projects of China (2006DAK31700)the GF Verification Program of the National Satellite Meteorological Center of China (220043001011003-1)
文摘In order to provide a long time-series,high spatial resolution,and high accuracy dataset of land surface temperature（LST） for climatic change research,a modified Becker and Li＇s split-window approach is proposed in this paper to retrieve LST from the measurements of Advanced Very High Resolution Radiometer（AVHRR） onboard National Oceanic and Atmospheric Administration（NOAA）-7 to-18 and the Visible and InfraRed Radiometer（VIRR） onboard FY-3A.For this purpose,the Moderate Resolution Transmittance Model（MODTRAN） 4.1 was first employed to compute the spectral radiance at the top of atmosphere（TOA） under a variety of surface and atmosphere conditions.Then,a temperature dataset consists of boundary temperature T s（which is one of the input parameters to MODTRAN）,and channels 4 and 5 brightness temperatures（T 4 and T 5） were constructed.Note that channels 4 and 5 brightness temperatures were simulated from the MODTRAN output spectral radiance by convolving them with the spectral response functions（SRFs） of channels 4 and 5 of AVHRRs and VIRR.The coefficients of modified Becker and Li＇s split-window approach for various AVHRRs and VIRR were subsequently regressed based on this temperature dataset using the least square method.As an example of validation,one AVHRR satellite image over Beijing acquired at 0312 UTC 27 April 2008 by AVHRR onboard NOAA-17 was selected to retrieve the LST image using the modified Becker and Li＇s approach.The comparison between this LST image and that from the MODIS level-2 LST product provided by the University of Tokyo in Japan indicates that the correlation coefficient is 0.88,the bias is 0.6 K,and the root mean square deviation（RMSD） is 2.1 K.Furthermore,about 70% and 37% pixels in the LST difference image,which is the result of retrieved LST image from AVHRR minus the corresponding MODIS LST image,have the values within ± 2 and ± 1 K,respectively.
基金support from the Strategic Pilot Science and Technology project of the Chinese Academy of Sciences(Grant No.XDA05040200)the National Natural Science Foundation of China(Grant No.41275040)
文摘It has been several years since the Greenhouse Gases Observing Satellite （GOSAT） began to observe the distribution of CO2 and CH4 over the globe from space. Results from Thermal and Near-infrared Sensor for Carbon Observation-Cloud and Aerosol Imager （TANSO-CAI） cloud screening are necessary for the retrieval of CO2 and CH4 gas concentrations for GOSAT TANSO-Fourier Transform Spectrometer （FTS） observations. In this study, TANSO-CAI cloud flag data were compared with ground-based cloud data collected by an all-sky imager （ASI） over Beijing from June 2009 to May 2012 to examine the data quality. The results showed that the CAI has an obvious cloudy tendency bias over Beijing, especially in winter. The main reason might be that heavy aerosols in the sky are incorrectly determined as cloudy pixels by the CAI algorithm. Results also showed that the CAI algorithm sometimes neglects some high thin cirrus cloud over this area.
基金supported by the National Natural Science Foundation of China (Grant 41775032)Key Research Program of Frontier Sciences, Chinese Academy of Sciences (Grant QYZDY-SSW-DQC027)
文摘The radiances scattered or emitted by clouds demonstrate diverse features at different wavelengths due to different cloud physical structures. This paper presents a method(the smallest-radiance-distance method, SRaDM) of revealing the physical structures of clouds. The method is based on multi-spectral radiances measured by the Moderate Resolution Imaging Spectroradiometer(MODIS)onboard Aqua. The principle and methodology of SRaDM is deduced and provided in this paper. Correlation analysis based on data from MODIS and Cloud Profiling Radar(onboard CloudSat), collected from January 2007 to December 2010 over an ocean area(15°N–45°N,145°E–165°E), led to selection of radiances at 13 wavebands of MODIS that demonstrated high sensitivity to cloud physical structures;radiances at the selected wavebands were subjected to SRaDM. The Standardized Euclidean distance is introduced to quantify the degree of changes in multi-spectral radiances(termed D_(rd)) and in physical structures(termed D_(st)) between cloud profiles. Statistics based on numerous cloud profiles show that D_(rd) decreases monotonically with a decrease in D_(st), which implies that small D_(rd) always accompanies small D_(st). According to the law of D_(rd) and D_(st), the new method, SRaDM, for revealing physical structures of clouds from the collocation of cloud profiles of similar multi-spectral radiances, is presented. Then, two successful experiments are presented in which cloud physical structures are captured using multi-spectral radiances. SRaDM provides a way to obtain knowledge of the physical structures of clouds over relatively larger areas, and is a new approach to obtaining 3D cloud fields.
基金supported by the Special Fund for Strategic Pilot Technology,Chinese Academy of Sciences[grant number XDA05040300]
文摘In this study,the authors focus on the cut-off low pressure systems（COLs）lingering over East Asia in late spring and early summer and quantify the two-way stratosphere–troposphere exchange（STE）by 3D trajectory integrations,achieved using a revised version of the UK Universities Global Atmospheric Modelling Programme Offline Trajectory Code（Version 3）.By selecting 10 typical COLs and calculating the cross-tropopause air mass fluxes,it is found that stratosphere-to-troposphere transport（STT）fluxes exist in the center of COLs;and in the periphery of the COL center,troposphereto-stratosphere transport（TST）fluxes and STT fluxes are distributed alternately.Net transport fluxes in COLs are from stratosphere to troposphere,and the magnitude is about 10-4 kg m-2 s-1.The ratio between the area-averaged STT and TST fluxes increases with increasing strength of the COLs.By adopting appropriate residence time,the spurious transports are effectively excluded.Finally,the authors compare the results with previous studies,and find that the cross-tropopause fluxes（CTFs）induced by COLs are about one to two orders of magnitude larger than global CTFs.COLs play a significant role in local,rapid air mass exchanges,although they may only be responsible for a fraction of the total STE.
