Ground-based GNSS for climate research: review and perspectives

In climate research, the role of water vapour can hardly be overestimated. Water vapour is the most important natural greenhouse gas and is responsible for the largest known feedback mechanism for amplifying climate change. It also strongly influences atmospheric dynamics and the hydrologic cycle through surface evaporation, latent heat transport and diabatic heating, and is, in particular, a source of clouds and precipitation.Atmospheric water vapour is highly variable, both in space and in time. Therefore, measuring it remains a demanding and challenging task. The Zenith Total Delay (ZTD) estimated from GNSS observations, provided at a temporal resolution of minutes and under all weather conditions, can be converted to Integrated Water Vapour (IWV), if additional meteorological variables are available. Inconsistencies introduced into long-term time series from improved GNSS processing algorithms, instrumental, and environmental changes at GNSS stations make climate trend analyses challenging. Ongoing re-processing efforts using state-of-the-art models aim at providing consistent time series of tropospheric data, using 24+ years of GNSS observations from global and regional networks. GNSS is reaching the “maturity age” of 30 years when climate normal of ZTD/IWV (and horizontal gradients) can be derived. Being not assimilated in numerical weather prediction model reanalyses, GNSS products can also be used as independent datasets to validate climate model outputs (ZTD/IWV). However, what is the actual use of GNSS ZTDs in climate monitoring? What are the advantages of using GNSS ZTDs for climate monitoring? In addition, what would be the best ZTD time series to serve the climate community?The presentation will provide a review of the progress made in and the status of using GNSS tropospheric datasets for climate research, highlighting the challenges and pitfalls, and outlining the major remaining steps ahead. We will show examples demonstrating the benefits for climate monitoring brought by using GNSS ZTD and/or IWV datasets in complement to other observations.This contribution is related to the activities of JWG C.2: Quality control methods for climate applications of geodetic tropospheric parameters, https://iccc.iag-aig.org/joint-work-groups/216, of the IAG Inter-Commission Committee on "Geodesy for Climate Research" (ICCC)

An enhanced integrated water vapour dataset from more than 10 000 global ground-based GPS stations in 2020

We developed a high-quality global integrated water vapour (IWV) dataset from 12 552 ground-based global positioning system (GPS) stations in 2020. It consists of 5 min GPS IWV estimates with a total number of 1 093 591 492 data points. The completeness rates of the IWV estimates are higher than 95 % at 7253 (58 %) stations. The dataset is an enhanced version of the existing operational GPS IWV dataset provided by the Nevada Geodetic Laboratory (NGL). The enhancement is reached by employing accurate meteorological information from the fifth generation of European ReAnalysis (ERA5) for the GPS IWV retrieval with a significantly higher spatiotemporal resolution. A dedicated data screening algorithm is also implemented. The GPS IWV dataset has a good agreement with in situ radiosonde observations at 182 collocated stations worldwide. The IWV biases are within ±3.0 kg m−2 with a mean absolute bias (MAB) value of 0.69 kg m−2. The standard deviations (SD) of IWV differences are no larger than 3.4 kg m−2. In addition, the enhanced IWV product shows substantial improvements compared to NGL\u27s operational version, and it is thus recommended for high-accuracy applications, such as research of extreme weather events and diurnal variations of IWV and intercomparisons with other IWV retrieval techniques. Taking the radiosonde-derived IWV as reference, the MAB and SD of IWV differences are reduced by 19.5 % and 6.2 % on average, respectively. The number of unrealistic negative GPS IWV estimates is also substantially reduced by 92.4 % owing to the accurate zenith hydrostatic delay (ZHD) derived by ERA5. The dataset is available at https://doi.org/10.5281/zenodo.6973528 (Yuan et al., 2022)

Homogenizing GPS integrated water vapour time series: methodology and benchmarking the algorithms on synthetic datasets

We would like to thank the COST Action ES1206 GNSS4SWEC for financial support

Performance of various homogenization tools on a synthetic benchmark dataset of GPS and ERA-interim IWV differences

Presentación realizada en: IAG-IASPEI 39th Joint Scientific Assembly celebrada en Kobe, Japón, del 30 de julio al 4 de agosto de 2017

Study on homogenization of synthetic GNSS-Retrieved IWV time series and its impact on trend estimates with autoregressive noise

Póster presentado en: EGU General Assembly celebrada del 23 al 28 de abril de 2017 en Viena, Austria.A synthetic benchmark dataset of Integrated Water Vapour (IWV) was created within the activity of “Data homogenisation” of sub-working group WG3 of COST ES1206 Action. The benchmark dataset was created basing on the analysis of IWV differences retrieved by Global Positioning System (GPS) International GNSS Service (IGS) stations using European Centre for Medium-Range Weather Forecats (ECMWF) reanalysis data (ERA-Interim). Having analysed a set of 120 series of IWV differences (ERAI-GPS) derived for IGS stations, we delivered parameters of a number of gaps and breaks for every certain station. Moreover, we estimated values of trends, significant seasonalities and character of residuals when deterministic model was removed. We tested five different noise models and found that a combination of white and autoregressive processes of first order describes the stochastic part with a good accuracy. Basing on this analysis, we performed Monte Carlo simulations of 25 years long data with two different types of noise: white as well as combination of white and autoregressive processes. We also added few strictly defined offsets, creating three variants of synthetic dataset: easy, less complicated and fully complicated. The synthetic dataset we present was used as a benchmark to test various statistical tools in terms of homogenisation task. In this research, we assess the impact of the noise model, trend and gaps on the performance of statistical methods to detect simulated change points

Homogenization of tropospheric data: evaluating the algorithms under the presence of autoregressive process

Presentación realizada en: IX Hotine-Marussi Symposium celebrado en Roma del 18 al 22 de junio de 2018.This research was supported by the Polish National Science Centre, grant No. UMO-2016/21/B/ST10/02353

Advanced exploitation of Ground-Based measurements for Atmospheric Chemistry and Climate Applications "AGACC"

We live in an era in which human activities are causing significant changes to the atmospheric environment which result in local to global consequences on the ecosystems. Changes in the atmospheric composition impact our climate via chemical and dynamical feedback mechanisms; in many instances they also affect air quality, and the health of the biosphere. Monitoring and understanding those changes and their consequences is fundamental to establish adequate actions for adaptation to and mitigation of the environmental changes. Furthermore, after implementation of regulatory measures like the Montreal Protocol, it is necessary to verify whether the measures are effective. This can only be achieved if we have adequate detection methods and a reliable long record of a series of key geophysical parameters. Thus the AGACC project contributes to the provision of basic new knowledge regarding the atmospheric composition and its changes, based on advanced groundbased monitoring, in combination with satellite and numerical modelling data. Its results are integrated in ongoing international research programmes. The general objective of AGACC has been to improve and extend the groundbased detection capabilities for a number of climate-related target species and, based hereupon, analyse past and present observations to derive new information about the atmospheric composition, its variability and long-term changes. Despite the advent of a growing and more performant fleet of Earth Observation satellites, ground-based observations are still indispensable to (1) guarantee long-term continuity, homogeneity and high quality of the data, and (2) to underpin the satellite data for calibration and (long-term) validation. A first target gas is atmospheric water vapour. It is the key trace gas controlling weather and climate. It is also the most important greenhouse gas in the Earth’s atmosphere. Its amount and vertical distribution are changing, but how and why? Especially in the upper troposphere - lower stratosphere, the radiative effects of changes in the water vapour are significant and should be quantified. The measurement of water vapour is a hot topic since several years. It is a challenge, because water vapour exhibits a large gradient in its concentration when going from the ground to the stratosphere, and because it is highly variable in time and space. For example, we have found that the time scale of the variations of the total water vapour amount at Jungfraujoch is in the order of minutes. In AGACC, we have therefore investigated various experimental techniques to measure the concentration of water vapour in the atmosphere, focusing on the total column as well as on the vertical distribution in the troposphere up to the lower stratosphere. The retrieval of water vapour vertical profiles and total columns from ground-based FTIR data has been initiated at three very different stations where correlative data for verification are available, namely Ukkel (± sea level, mid-latitude), Ile de La Réunion (± sea level, tropical) and Jungfraujoch (high altitude, mid-latitude), with promising results. In particular, at Jungfraujoch, it has been demonstrated that the precision of the FTIR integrated water vapour (IWV) measurements is of order 2%. The capability to retrieve individual isotopologues of water vapour, and to monitor their daily and diurnal variations, has also been demonstrated. This could open new ways to study in the future the role of water vapour in the radiative balance, the global circulation, precipitation etc. We also started joint exploitation of ground-based FTIR and satellite IASI data for water vapour and its isotopologues, in order to exploit fully the potential of the existing instrumentation. A correction method for the radiosoundings at Ukkel has been successfully implemented, resulting in a homogeneous and reliable time series from 1990 to 2008 from which trends in upper troposphere humidity (UTH) and tropopause characteristics have been derived. One observes a rising UTH until September 2001, followed by a decline, accompanied by a descent and heating of the tropopause up to the turning point and an ascent and cooling afterwards. The changes after September 2001 in the upper troposphere can be explained by surface heating and convective uplift. At Jungfraujoch, one does not observe any significant trend in the total water vapour abundance above the station over the 1988-2010 time period, although significant positive summer and negative winter trends have been detected. We have made a quantitative statistical comparison between ground-based FTIR, CIMEL, GPS and integrated (corrected) radio sounding measurements of the IWV at Ukkel. This work is important to better characterize the different sensors in order to exploit together different observations made by different instruments. A second target species is atmospheric aerosol. There is a very large variety of aerosol both from natural or anthropogenic origin. One of the reasons why they are so important is that they affect the optical properties of the atmosphere. In particular, it has been demonstrated in previous studies that the aerosols have a large impact on the quantity of harmful UV-B radiation received at the Earth’s surface. The latest IPCC Report also stressed that the radiative forcing caused by atmospheric aerosols is one of the largest uncertainties in determining the total radiative forcing in the atmosphere. Better monitoring capabilities of aerosol properties can therefore improve our understanding and forecasting of the atmospheric processes and evolution, and in particular of UV-B and climate changes. Several measurement techniques are now operational in the AGACC consortium for the ground-based monitoring of aerosol properties. These are the Brewer spectrometer and CIMEL observations at Ukkel, the latter contributing also to the AERONET network since July 2006, and the newly developed MAXDOAS observations. Unlike CIMEL and Brewer measurements, that provide the total Aerosol Optical Depth, it has been demonstrated that the MAXDOAS measurements also provide additional information about the vertical distribution of the aerosol extinction in the lowest kilometres of the troposphere. A better understanding of the ultimate capabilities of MAXDOAS aerosol remote sensing has been gained through participation to the international CINDI campaign (Cabauw Intercomparison Campaign of Nitrogen Dioxide measuring Instruments ) in summer 2009. The combination of Brewer, CIMEL and MAXDOAS instruments gives us a remote-sensing dataset that will enable a more comprehensive characterization of the tropospheric aerosol optical properties. The usefulness of these aerosol observations has already been demonstrated in the improvement of the UVindex predictions for the general public. Another application is their use as input data in the retrieval of vertical profiles of tropospheric pollutants from MAXDOAS measurements, like tropospheric NO2 and formaldehyde. Third we have focused on a few climate-related trace gases. Changing greenhouse gas and aerosol concentrations directly affect the radiative budget of the atmosphere, and therefore climate. But many species known as pollutants like carbon monoxide (CO), nitrogen oxides (NOx) and hydrocarbons, - often related to fossil fuel or biomass burning -, also affect climate through their role in chemical reactions that produce tropospheric ozone, which is a well-known greenhouse gas, or that modify the lifetime of gases like methane, or the oxidation capacity of the atmosphere. Therefore in AGACC, we have focused on the measurement of a number of trace gases that are subject to changing concentrations, that directly or indirectly affect climate, and that are either difficult to monitor or that have not yet been measured from the ground. We have included attempts to observe distinctly some isotopologues, because the isotopic ratios observed in an airmass provide information on its history, and because the FTIR solar absorption measurements provide a rather unique capability hereto. The investigated species are the isotopologues of CH4 and CO, and hydrogen cyanide (HCN), as examples of biomass burning tracers, some hydrocarbons like formaldehyde (HCHO), ethylene (C2H4) and acetylene (C2H2), and HCFC-142b, a replacement product for CFCs and a greenhouse gas. In many cases, retrieval strategies had to be adapted when going from one site to another with different atmospheric conditions, especially when the local humidity and abundances are very different as is the case between Jungfraujoch (dry, high altitude, mid-latitude) and Ile de La Réunion (humid, low altitude, low latitude). Still we have been able to show the feasibility of retrieving particular trace gas information even under difficult conditions. Many of our results have been compared to correlative data, to validate the approach and to gain complementary information. It is also important to note that the retrieval strategies developed in AGACC have regularly been presented to the global Network for the Detection of Atmospheric Composition Change (NDACC) UV-Vis and Infrared communities and have often been adopted by others or even proposed for adoption as a standard in the community (e.g., for hydrogen cyanide (HCN)). In particular: We have been able to study the seasonal variations of HCN at the Jungfraujoch and at Ile de La Réunion, and to show the dominant impact of biomass burning. Formaldehyde was studied in much detail at Ukkel, Jungfraujoch and Ile de la Réunion. The challenge for detection at Jungfraujoch is the small abundance (about 10 times smaller than at Ukkel and Ile de La Réunion); a particular observation strategy was developed successfully, resulting in a time series that already shows the day-to-day and seasonal variations. At Ile de La Réunion, comparisons of FTIR, MAXDOAS, satellite and model data have (1) shown the good agreement between the various data sets, but also, (2), the variability of HCHO (diurnal, seasonal, day-to-day), and (3), thanks to the complementarities of the various data sets, they have enabled us to learn more about the long-range transport of Non-methane Volatile Organic Compounds (NMVOCS, precursors of HCHO) and deficiencies in the models. It was shown that fast, direct transport of NMVOCS from Madagascar has a significant impact on the HCHO abundance and its variability at Ile de La Réunion, and that this is underestimated in the model. Significant progress was made as to the detection of 13CH4 and CH3D from ground-based FTIR observations, both at Jungfraujoch and Ile de La Réunion. To our knowledge, it is the first time that a d13C data set is derived from ground-based FTIR observations. More work is needed to improve the CH3D retrieval at Ile de La Réunion, and to interpret the results, in combination with models. Also for the first time, 12CO and 13CO have been retrieved individually at Jungfraujoch. The d13C time series shows significant seasonal and interannual changes. As to the hydrocarbon ethylene, it is shown that it can be detected at Jungfraujoch only in spectra at low solar elevation, given its small atmospheric abundance. Regarding acetylene, the observed time series at Jungfraujoch and Ile de La Réunion show clear seasonal variations and enhancements due to the impact of biomass burning events, correlated with enhancements in CO, C2H6 and HCN. It is not clear yet whether we can reliably retrieve the concentration of HCFC- 142b, a replacement product that is increasing strongly in the troposphere. New line parameters for the interfering species HFC-134a are required to confirm/infirm the preliminary results. This highlights again the importance of the laboratory work for providing such parameters. Improved line parameters have been obtained for water vapour and its isotopologues, ethylene and formic acid. These AGACC results have been integrated in the international spectroscopic databases. We also showed that line intensities available around 2096 cm–1 for the 13C16O isotopologue of carbon monoxide in the HITRAN database seem to be accurate to 2%. We failed to improve line intensities for the 13.6 μm region of acetylene. The new data sets that have been derived in AGACC from FTIR and MAXDOAS observations have been archived in the NDACC data centre, where they are available for users (generally modelers and satellite teams). In addition, they are stored locally and are available to users upon request. AGACC results have been reported to the international scientific community, via the literature, via integration in geophysical or spectroscopic databases, and via participation to international research initiatives like the Atmospheric Water Vapour in the Climate System (WAVACS) Cost Action, the International Space Science Institute (ISSI) Working Group on Atmospheric Water Vapour, the International Union of Pure and Applied Chemistry (IUPAC) project, the International CINDI campaign, etc. The results have already found important scientific applications. A few examples are worth mentioning: the re-evaluation of methane emissions in the tropics from SCIAMACHY based on the new H2O spectroscopy, and the improved retrievals of HCOOH from the satellite experiments ACE-FTS and IASI, and from the ground. In the longer-term, the AGACC results will no doubt benefit the research in atmospheric sciences –in particular in the monitoring of its composition changes–, which is the fundamental basis of environmental assessment reports for supporting policy makers.Advanced exploitation of ground-based measurements for atmospheric chemistry and climate applications "AGACC

Updated trends of the stratospheric ozone vertical distribution in the 60° S–60° N latitude range based on the LOTUS regression model

This study presents an updated evaluation of stratospheric ozone profile trends in the 60° S–60° N latitude range over the 2000–2020 period using an updated version of the Long-term Ozone Trends and Uncertainties in the Stratosphere (LOTUS) regression model that was used to evaluate such trends up to 2016 for the last WMO Ozone Assessment (2018). In addition to the derivation of detailed trends as a function of latitude and vertical coordinates, the regressions are performed with the datasets averaged over broad latitude bands, i.e. 60–35° S, 20° S–20° N and 35–60° N. The same methodology as in the last assessment is applied to combine trends in these broad latitude bands in order to compare the results with the previous studies. Longitudinally resolved merged satellite records are also considered in order to provide a better comparison with trends retrieved from ground-based records, e.g. lidar, ozonesondes, Umkehr, microwave and Fourier transform infrared (FTIR) spectrometers at selected stations where long-term time series are available. The study includes a comparison with trends derived from the REF-C2 simulations of the Chemistry Climate Model Initiative (CCMI-1). This work confirms past results showing an ozone increase in the upper stratosphere, which is now significant in the three broad latitude bands. The increase is largest in the Northern and Southern Hemisphere midlatitudes, with ∼2.2 ± 0.7 % per decade at ∼2.1 hPa and ∼2.1 ± 0.6 % per decade at ∼3.2 hPa respectively compared to ∼1.6 ± 0.6 % per decade at ∼2.6 hPa in the tropics. New trend signals have emerged from the records, such as a significant decrease in ozone in the tropics around 35 hPa and a non-significant increase in ozone in the southern midlatitudes at about 20 hPa. Non-significant negative ozone trends are derived in the lowermost stratosphere, with the most pronounced trends in the tropics. While a very good agreement is obtained between trends from merged satellite records and the CCMI-1 REF-C2 simulation in the upper stratosphere, observed negative trends in the lower stratosphere are not reproduced by models at southern and, in particular, at northern midlatitudes, where models report an ozone increase. However, the lower-stratospheric trend uncertainties are quite large, for both measured and modelled trends. Finally, 2000–2020 stratospheric ozone trends derived from the ground-based and longitudinally resolved satellite records are in reasonable agreement over the European Alpine and tropical regions, while at the Lauder station in the Southern Hemisphere midlatitudes they show some differences

Homogenization of the long-term global ozonesonde records

Póster presentado en: WMO Technical Conference on Meteorological and Environmental Instruments and Methods of Observation celebrada del 10 al 13 de octubre de 2022 en París