First data set of H2O/HDO columns from the Tropospheric Monitoring Instrument (TROPOMI)

Global measurements of atmospheric water vapour isotopologues aid to better understand the hydrological cycle and improve global circulation models. This paper presents a new data set of vertical column densities of H2O and HDO retrieved from short-wave infrared (2.3 μm) reflectance measurements by the Tropospheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor satellite. TROPOMI features daily global coverage with a spatial resolution of up to 7km× 7km. The retrieval utilises a profile-scaling approach. The forward model neglects scattering, and strict cloud filtering is therefore necessary. For validation, recent ground-based water vapour isotopologue measurements by the Total Carbon Column Observing Network (TCCON) are employed. A comparison of TCCON δD with ground-based measurements by the Multi-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water (MUSICA) project for data prior to 2014 (where MUSICA data are available) shows a bias in TCCON δD estimates. As TCCON HDO is currently not validated, an overall correction of recent TCCON HDO data is derived based on this finding. The agreement between the corrected TCCON measurements and co-located TROPOMI observations is good with an average bias of (−0.2± 3)× 1021 molec cm−2 ((1.1± 7.2)%) in H2O and (−2±7)×1017 molec cm−2 ((−1.1±7.3)%) in HDO, which corresponds to a mean bias of (−14± 17)‰ in a posteriori δD. The bias is lower at lowand mid-latitude stations and higher at high-latitude stations. The use of the data set is demonstrated with a case study of a blocking anticyclone in northwestern Europe in July 2018 using single-overpass data.

Abstract. Global measurements of atmospheric water vapour isotopologues aid to better understand the hydrological cycle and improve global circulation models. This paper presents a new data set of vertical column densities of H 2 O and HDO retrieved from short-wave infrared (2.3 µm) reflectance measurements by the Tropospheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor satellite. TROPOMI features daily global coverage with a spatial resolution of up to 7 km × 7 km. The retrieval utilises a profile-scaling approach. The forward model neglects scattering, and strict cloud filtering is therefore necessary. For validation, recent ground-based water vapour isotopologue measurements by the Total Carbon Column Observing Network (TCCON) are employed. A comparison of TCCON δD with ground-based measurements by the Multi-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water (MUSICA) project for data prior to 2014 (where MUSICA data are available) shows a bias in TCCON δD estimates. As TCCON HDO is currently not validated, an overall correction of recent TCCON HDO data is derived based on this finding. The agreement between the corrected TCCON measurements and co-located TROPOMI observations is good with an average bias of (−0.2 ± 3) × 10 21 molec cm −2 ((1.1 ± 7.2) %) in H 2 O and (−2±7)×10 17 molec cm −2 ((−1.1±7.3) %) in HDO, which corresponds to a mean bias of (−14 ± 17) ‰ in a posteri-ori δD. The bias is lower at low-and mid-latitude stations and higher at high-latitude stations. The use of the data set is demonstrated with a case study of a blocking anticyclone in northwestern Europe in July 2018 using single-overpass data.

Introduction
Atmospheric water vapour represents the strongest natural greenhouse gas and transports a large amount of energy via latent heat; thus, it plays a fundamental role in shaping weather and climate (Kiehl and Trenberth, 1997;Harries, 1997). However, uncertainties in the quantification of the two abovementioned effects are still large and represent one of the key uncertainties in current climate prediction (Stevens and Bony, 2013). Improvement upon current climate prediction requires new observations on a global scale and with a long-term perspective. To this end, satellite observations from space are considered to be the most promising approach (Rast et al., 2014).
Constraints for the hydrological cycle are offered by observations of isotopologues of water vapour. Different equilibrium vapour pressures and diffusion constants of different isotopologues lead to isotopic fractionation whenever a phase change occurs. Isotopic fractionation occurs at the point of A. Schneider et al.: H 2 O/HDO from TROPOMI phase change, partitioning the heavier and lighter isotopologues, depending on the thermodynamic conditions of the environment. The relative abundance of a heavy isotopologue with respect to the light isotopologue in an air parcel is therefore dependent on the source region's temperature and relative humidity, the source water's isotopic composition as well as the entire transport history of the air parcel, including all evaporation, condensation and mixing events (e.g. Dansgaard, 1964;Craig and Gordon, 1965). This makes measurements of water vapour isotopologues a unique diagnostic of the hydrological cycle (Dansgaard, 1964) and a valuable benchmark for the evaluation and further development of global and regional circulation models (e.g. Joussaume et al., 1984;Hoffmann et al., 1998;Yoshimura et al., 2008;Risi et al., 2010;. The usual notation to describe the isotopological abundance variations is the relative difference of the ratio of the heavy and the light isotopologues, here HDO and H 2 O, (Coplen, 2011). The commonly used standard ratio is Vienna Standard Mean Ocean Water (VSMOW), R D,std = 3.1152 × 10 −4 . Measurements of atmospheric water vapour isotopologues are not very common. In situ observations are performed from aircrafts and balloons (e.g. Rinsland et al., 1984;Dyroff et al., 2010Dyroff et al., , 2015Herman et al., 2014;Sodemann et al., 2017) and on the ground (e.g. Wen et al., 2010;Aemisegger et al., 2012;Bastrikov et al., 2014) using laser spectrometers or cryogenic trapping techniques. Remote sensing instruments exist on the ground and on space-or balloon-based platforms. The former are usually Fourier transform infrared (FTIR) spectrometers. Ground stations are often organised in networks. The largest networks are the Total Carbon Column Observing Network (TCCON, Wunch et al., 2011) and the Network for the Detection of Atmospheric Composition Change (NDACC, De Mazière et al., 2018). The data product of the former includes H 2 O and HDO, whereas the latter involves water vapour isotopologue measurements retrieved by the Multi-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water (MUSICA) project . With respect to satellites, H 2 O and HDO were first retrieved by Zakharov et al. (2004) using thermal infrared measurements from the Interferometric Monitor for Greenhouse gases (IMG) sensor onboard the Advanced Earth Observing Satellite (ADEOS). Later, this was followed by the Tropospheric Emission Spectrometer (TES) on the Earth Observing System (EOS) Aura satellite (Worden et al., 2006), the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) onboard the European Space Agency (ESA)'s environmental satellite (ENVISAT) (Steinwagner et al., 2007;Payne et al., 2007), the SCanning Imaging Absorption spectroMeter for Atmo-spheric CHartographY (SCIAMACHY) instrument on EN-VISAT (Frankenberg et al., 2009;Scheepmaker et al., 2015;Schneider et al., 2018), the Infrared Atmospheric Sounding Interferometer (IASI) onboard the MetOP satellites (Herbin et al., 2009;Schneider and Hase, 2011;Schneider et al., 2016;Lacour et al., 2012), the Greenhouse Gases Observing Satellite (GOSAT) (Frankenberg et al., 2013;Boesch et al., 2013) and the Atmospheric Infrared Sounder (AIRS) onboard the NASA Aqua satellite (Worden et al., 2019). The sensitivity of instruments observing in the thermal infrared (IMG, TES, MIPAS, IASI and AIRS) is very different from that of instruments measuring in the short-wave infrared, such as SCIAMACHY and GOSAT. While the former are mainly sensitive in the stratosphere and free troposphere, the latter have good sensitivity in the lower troposphere, including the boundary layer. On 13 October 2017, the Tropospheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor (S5P) satellite (Veefkind et al., 2012) was launched. It has a short-wave infrared band in heritage of SCIAMACHY with a spectral range of 2305-2385 nm and a spectral resolution of 0.25 nm, although its signal-to-noise ratio is much better than SCIAMACHY and it has an unprecedented spatial resolution of 7 km×7 km (in the centre of the swath). This work presents a new H 2 O and HDO column data set from TROPOMI observations starting at first light of the instrument on 9 November 2017. Section 2 introduces the retrieval method. Section 3 presents a ground-based data set to validate the satellite observations against, and the comparison between both data sets is shown in Sect. 4. Section 5 provides a first insight into the data set's use with respect to studying synoptic-scale variability in the atmospheric branch of the water cycle. Finally, the summary of the results and the conclusions are given in Sect. 6.

Retrieval method
The retrievals are performed with SICOR (short-wave infrared CO retrieval algorithm), which utilises a profilescaling approach and is described in detail by Scheepmaker et al. (2016), Landgraf et al. (2016) and Borsdorff et al. (2014). In the following, the most important features are summarised and the specific setup is given.
Using the spectral window from 2354.0 to 2380.5 nm , the algorithm fits the total columns of H 2 O, HDO, CH 4 and CO as well as a Lambertian surface albedo in the form of a Legendre polynomial of order 1. The isotopologue H 2 18 O is included in the forward model but not fitted. A priori profiles of water vapour are adapted from the European Centre for Medium-Range Weather Forecasts (ECMWF) analysis product. As the ECMWF data product does not distinguish between individual isotopologues, H 2 O, HDO and H 2 18 O profiles are obtained from the water vapour profile by scaling it with the respective average relative natural abundances. That implicitly corresponds to a Corresponding residuals (defined as measured minus modelled radiances, in blue) and its root mean square (rms, in cyan), precision of the radiance (in red) and its rms (in purple). (c) Simulated absorption by H 2 O (red), HDO (green) and CH 4 (yellow).
prior of δD of 0 ‰. A priori profiles of CH 4 and CO are taken from TM5 simulations (Krol et al., 2005). Scattering crosssections are taken from HITRAN 2016 (Gordon et al., 2017). The forward model ignores scattering, so that strict filtering for clear-sky scenes is necessary. To this end, co-located measurements from the Visible Infrared Imaging Radiometer Suite (VIIRS) instrument onboard the Suomi National Polar-orbiting Partnership (S-NPP) satellite, which flies in formation with S5P, are used (Siddans, 2016). The cloud cover threshold is 1 % for both the inner field of view and the outer field of view. Moreover, soundings with a high aerosol load are filtered out by a two-band filter as introduced by Scheepmaker et al. (2016) and Hu et al. (2018), which in the present configuration requires that the ratio of retrieved methane in bands with weak and strong absorption (2310-2315 and 2363-2373 nm respectively) is between 0.94 and 1.06. Furthermore, scenes with a solar zenith angle greater than 75 • are discarded because they are prone to errors due to more scattering and diffraction effects, which are not cov-ered well by the forward model, and due to typically low radiances, meaning low signal-to-noise ratios.
An exemplary spectral fit and the resulting residuals (which are defined as measured minus modelled radiances) are shown in Fig. 1. The root-mean-square (rms) residual (cyan horizontal line in Fig. 1b) is in the order of the rms uncertainty of the radiance (purple horizontal line in Fig. 1b).
The sensitivity of a retrieved column to changes in a given altitudinal region is described by the column averaging kernel (Rodgers, 2000). The ideal averaging kernel is unity at all altitudes, but in practice the sensitivity changes with height. Figure 2 depicts examples of column averaging kernels for different solar zenith angles. The sensitivity for the two isotopologues are significantly different. For H 2 O, the highest sensitivity is in the lowest layer (where most water vapour typically resides) and decreases with increasing altitude. The sensitivity in the stratosphere is small; however, the amount of water vapour in this altitudinal region is very small and contributes little to the total column. The sensitivity of HDO does not deviate as much from unity as that of H 2 O. In the lower troposphere it increases slightly with increasing altitude before reaching a maximum depending on the solar zenith angle, above which it decreases. The differences in the column averaging kernel are due to the different absorption strengths of the two isotopologues and mean that a posteriori δD is sensitive to the profile shapes, particularly of the main isotopologue H 2 O; this is due to the fact that the averaging kernel for H 2 O deviates considerably from unity at higher altitudes.

Ground-based FTIR data sets
To validate the TROPOMI retrievals, ground-based Fourier transform infrared (FTIR) measurements are used. HDO is a product of NDACC-MUSICA ( Barthlott et al., 2017) and TCCON . NDACC-MUSICA provides two products: type 1 is the direct retrieval output, and type 2 contains a posteriori processed output that reports the optimal estimation of (H 2 O, δD) pairs; here the type 2 product is used because it is recommended for isotopologue analyses (  co-locating with a maximal time difference of 1 h. An example for Wollongong is plotted in Fig. 3. A comparison between MUSICA and TCCON was also performed by , who compared the MUSICA type 1 product with TC-CON and found a bias in δD of 40 ‰ on average. MUSICA is explicitly created for isotopologue studies, and δD profiles have been validated against aircraft measurements in an altitudinal range between 2 and 7 km during a dedicated campaign in summer 2013 Dyroff et al., 2015). However, data are only available until 2014; thus, there is no temporal overlap with TROPOMI which was launched in October 2017. TCCON H 2 O total columns are calibrated with in situ measurements (mainly radiosondes); a so-called aircraft correction factor of 1.0183 is applied to match the reference . However, TCCON HDO is currently not verified; thus, no correction factor is applied to it. Therefore, it is assumed that TCCON HDO has to be corrected.
In order to correct for the discrepancy, the idea is to scale TCCON HDO to match MUSICA δD. Scaling HDO by a factor a, i.e. c HDO −→ a c HDO , is equivalent to the linear transformation in δD. Figure 4a depicts a correlation histogram of TCCON δD vs. MUSICA δD for the Wollongong station. Here, the relation between MUSICA and TCCON is described to a large degree by a simple scaling of the column. The result of a fit of Eq.
(2) to the data is plotted as a blue line, giving the scaling factor for the TCCON HDO column. To demonstrate that this approach does not involve intercept issues, a linear fit of slope and intercept (red line) as well as the confidence interval computed using the bootstrap method (i.e. by fitting a randomly reduced data set 10 000 times, shown using red shading) has also been plotted in the figure. Both the slope and the offset are similar to the approach using Eq. (2), and the latter lies within the confidence band of the former. The bar chart in Fig. 4c visualises fit results for all stations in both networks. It shows that the correction factor does not change much between stations. The large difference in fit error is mostly due to the large difference in the amount of data (Fig. 4b). Thus, it is meaningful to scale HDO at all TC-CON stations by the error-weighted average correction factor a = 1.0778 in order to correct TCCON's bias in HDO and, in turn, δD.

Validation of TROPOMI retrievals
For validation, TROPOMI observations are co-located with TCCON measurements with a radius of 30 km, a maximal altitude difference of 500 m, a field of view of 45 • in the FTIR viewing direction and a maximal time difference of 2 h. Here, the TCCON HDO data are corrected according to the approach presented in the previous section. Table 1 gives an overview of all stations used. Other stations have too few (less than 5 d) co-located measurements and have therefore not been included in the validation study. No altitude correction is applied here. The mentioned co-location criterion for altitude is used to ensure that no bias due to the large height difference between the station and satellite ground pixel is introduced (cf. Schneider et al., 2018). For each station, daily averages are computed over all co-located measurements. Figure 5 shows an exemplary time series for Edwards station. The co-located observations of H 2 O and HDO agree very well, and the agreement in δD is also good, with more scatter and a small bias. Corresponding correlation plots are depicted in Fig. 6. Figure 6a and b confirm the excellent agreement in H 2 O and HDO with Pearson correlation coefficients of 0.98 and 0.99 respectively, and a corresponding correlation coefficient of 0.96 for δD. The average difference between TROPOMI and TCCON defines the bias. In δD, a small bias is plain in the correlation plot and amounts to −22 ‰. Figure 7 depicts the validation statistics for all TCCON stations. The correlation in H 2 O and HDO is high for all stations. In δD, the correlation is high except for a low correlation of 0.37 at Saga station, where the amount of data is very small, the variability in δD is small, but H 2 O and HDO vary considerably. Figure 8 shows the biases. At low-and mid-latitude (< 54 • ) stations the bias is as low as (−0.3 ± 3) × 10 21 molec cm −2 (corresponding to a relative bias of (0.6±5.7) %) in H 2 O and (−2±8)×10 17 molec cm −2 ((−0.6 ± 6.3) %) in HDO, which corresponds to (−9 ± 11) ‰ (4.7 ± 6.8 %) in a posteriori δD. At these stations the bias in δD ranges between about −30 ‰ and +15 ‰. At high-latitude stations it can be as high as −45 ‰ to −60 ‰. Possible reasons for these high biases are higher relative biases in H 2 O and/or HDO at these relatively dry locations. At high-latitude, retrievals are generally challenging due to high solar zenith angles and low albedos which lead to low signal-to-noise ratios. The average bias over all stations is (−0.2 ± 3) × 10 21 molec cm −2 or (1.1 ± 7.2) % in H 2 O, (−2±7)×10 17 molec cm −2 or (−1.1±7.3) % in HDO, and (−14 ± 17) ‰ or (5.6 ± 6.7) % in a posteriori δD. This is good considering that δD is very sensitive to small errors in H 2 O or HDO. Figure 9 shows how differences in TROPOMI observations and corrected TCCON measurements depend on the H 2 O column. For H 2 O and HDO there is no such dependence, and the Pearson correlations are very low: 0.06 and 0.08 respectively. In a posteriori δD there is no dependence for low-and mid-latitude stations (correlation coefficient of 0.28). At the high-latitude stations Eureka, Sodankylä and East Trout Lake (marked using red in Fig. 9), a large range of H 2 O columns show no dependence, but at H 2 O columns below ∼ 2×10 21 molec cm −2 (i.e. in dry conditions), the differences between TROPOMI and TCCON increase and do depend on the H 2 O column. This is the main reason for the high biases at these stations, as discussed in the previous paragraph. As mentioned above, retrievals at low albedos and high solar zenith angles are generally challenging. Whether the dependence of the difference in δD on the H 2 O column in dry conditions at high latitudes is due to co-location errors or due to the retrieval is still unclear and has to be examined in future research.

Demonstration of applications of the data set
An illustration of the TROPOMI retrievals on the global and monthly scale is depicted in Fig. 10 for September 2018. There are no data over the oceans because water is too dark in the short-wave infrared and glint measurements are not taken into account. The data gaps in tropical regions are due to persistent clouds. The data quality in terms of noise is significantly better than for a multiyear average of SCIAMACHY observations (cf. Schneider et al., 2018, Fig . 7). In the spatial distribution shown in Fig. 10 the major isotopic effects formulated by Dansgaard (1964) can be recognised. The general latitudinal gradient due to the temperature dependence of the fractionation effects and progressive rain out of heavy isotopologues, the so-called latitudinal effect, is clearly visible. The continental effect of depletion due to the rain out of the heavy isotopologue is visible on all continents, including Australia. The altitude effect, which describes depletion above high ground due to lower temperature and increasing rain out, can be seen, for example, over the Andes and the Himalayas.
To demonstrate the quality and the possibilities of the new data set of water vapour isotopologues from TROPOMI, a case study using single-overpass results over Europe on 30 July 2018 is presented in Fig. 11. The summer 2018 was one of the hottest and driest in central and northern Europe (Copernicus Climate Service, 2018;Gubler et al., 2018) with forest fires in Scandinavia, dry fields and low river stages all over the central and northern parts of the continent. The reason for this exceptionally hot and dry summer was the presence of a high-pressure system over northern Europe that blocked the otherwise predominant westerly moist flow from the North Atlantic. Synoptic-scale atmospheric blocking situations can lead to hot temperature ex- tremes due to adiabatic warming of the descending air in the core of the anticyclone . The descending vertical motion favours clear-sky conditions and, thus, further contributes to surface warming via radiative effects in the centre of the anticyclone (Trigo et al., 2004). In particular, the end of July 2018 was characterised by a stationary blocking anticyclone extending over the entire troposphere over northwestern Russia and Scandinavia. This blocking led to large-scale descent and to a divergent flow near the surface in its core, resulting in clear-sky conditions over northwestern Russia and Finland (see Fig. 11c). The iso-topic signature of the blocking anticyclone in Fig. 11b reflects this synoptic flow configuration with low δD signals of between −250 ‰ and −200 ‰ in the centre of the anticyclone. The depleted total column vapour in this region is due to the large-scale subsidence transporting depleted (Fig. 11b) and dry (Fig. 11d) upper tropospheric air towards lower levels. The near-surface divergent wind exports more enriched freshly evaporated moisture that is taken up near the surface towards the edges of the blocking. The anticyclone area is characterised by clear skies (Fig. 11c) with low specific humidity (1-3 g kg −1 at 700 hPa, Fig. 11d), low relative hu- Figure 11. TROPOMI single-overpass results for H 2 O column (a) and δD (b) over Europe on 30 July 2018; VIIRS cloud fraction on the same day (c); specific humidity (d), relative humidity (e), and potential temperature (f) at 700 hPa from the ECMWF analysis product over Europe at 12:00 UTC on 30 July 2018. The 700 hPa level is chosen for the thermodynamic variables because it reflects the large-scale conditions in the lower troposphere above the continental boundary layer. The overlaying contours in all panels show mean sea-level pressure from ECMWF at 12:00 UTC with a contour line distance of 2 hPa. midity (10 %-30 % at 700 hPa, Fig. 11e) and high potential temperature associated with the dry subsiding (adiabatically warming) air masses (Fig. 11f). The dry low-level outflow encounters moister and warmer air at the edge of the surface anticyclone, leading to a very strong horizontal gradient of specific and relative humidity (Fig. 11d, e) in the lower troposphere. As a consequence, the warm moist air is forced to rise, localised instabilities occur, and isolated convective cells develop leading to condensation and the formation of a ring of clouds around the blocking anticyclone. A distinct arc-like feature of enriched total column water vapour at the edge of the anticyclone can be distinguished and is slightly displaced from the first clouds in the northwest (Fig. 11b). Turbulent mixing and convection that inject more enriched, freshly evaporated moisture advected with the large-scale flow from marine environments (Barents Sea, North Sea and Black Sea) could be the reason for this interesting enriched ring-like water vapour isotopologue pattern. A very depleted cloud-free area south of the Ob River with δD values below −250 ‰ (Fig. 11b) might be connected to anomalously strong subsidence of northerly continental air masses.
The large range of measured δD and H 2 O mixing ratios during the northeastern European blocking on 30 July 2018 (Fig. 11a, b) becomes apparent in Fig. 12, where a twodimensional histogram and the cumulative density of the TROPOMI data in the 50-70 • N, 20-60 • E region is shown along with two different types of idealised air mass transformation scenarios (coloured lines). These simple idealised scenarios are frequently used in the literature to guide the interpretation of stable isotope measurements in (H 2 O, δD) diagrams (e.g. Rozanski and Sonntag, 1982;Worden et al., 2007;Noone, 2012). The first scenario, illustrated by the dashed green line in Fig. 12, is that an air parcel with a humidity of 15 000 ppm and δD = −80 ‰ (typical for the continental boundary layer in this northerly continental region, Bastrikov et al., 2014) experienced moist adiabatic ascent with condensation following a Rayleigh process (dashed green line; Rayleigh, 1902;Dansgaard, 1964). The progressive decrease of δD with a decreasing water vapour mixing ratio would thus be due to preferential condensation of HDO compared with H 2 O and the subsequent removal of hydrometeors by precipitation. The dashed green Rayleigh curve in Fig. 12 shows a behaviour that is different from the TROPOMI data points. Given the clear-sky conditions and the subsidising movement of the air masses within the block-ing anticyclone, the assumptions needed for a Rayleigh distillation process are hardly fulfilled. The second scenario, illustrated by the solid green, blue, and cyan lines in Fig. 12, is that two air parcels with distinct humidity and δD are mixed due to turbulent and convective mixing to yield different blends that follow the so-called mixing lines in the (H 2 O, δD) space. The highest density of observed blocking anticyclone points retrieved by TROPOMI is located in the region spanned by the green and the cyan mixing lines in the (H 2 O, δD) space. The 25 %, 50 % and 75 % contours of the cumulative density of points are aligned with the blue mixing line. This suggests that a two end-member mixing process describes the data much better than an idealised Rayleigh process (dashed green line in Fig. 12). In this particular synoptic situation, this corresponds to the moistening of a subsidising air mass from the mid troposphere.
In future work, the nature and occurrence of these features should be analysed in more detail, including a catalogue of different continental blocking events with observations from TROPOMI.
Apart from investigations on the water cycle dynamics associated with continental blockings, many other dynamically interesting contexts exist where TROPOMI could present an important added value for further investigations. These comprise, among others, the region of the heat low over the Sahara (e.g. Schneider et al., 2015;González et al., 2016;Lacour et al., 2017) or continental regions upstream of cold air surges leading to events of strong ocean evaporation along the warm ocean western boundary currents (Aemisegger and Papritz, 2018;Aemisegger and Sjolte, 2018).

Summary and conclusions
This work presents a new data set of H 2 O and HDO columns retrieved from TROPOMI short-wave infrared observations. Scattering is ignored in the forward model so that a strict cloud filtering is necessary, which is performed with colocated VIIRS measurements. The data quality is such that single overpasses yield meaningful results, which is a huge step forward compared with previous missions like SCIA-MACHY.
For validation of the TROPOMI data product, particular attention must be paid to the reference data sets. At this stage, there are two data products of ground-based observations of the HDO total column available, provided by the TCCON and NDACC-MUSICA networks. Comparing these two data products for stations in both networks reveals a large bias between the ground-based products of 58 ‰ on average in δD. NDACC-MUSICA was decidedly developed for water vapour isotopologue studies and is validated in δD with aircraft measurements; however, data are only available until 2014. TCCON provides recent data with temporal overlap with TROPOMI observations, and its H 2 O total column data product is validated against in situ measurements; however, its HDO data product is not verified. In order to obtain a suitable validation data set, TCCON HDO columns are scaled by a factor of 1.0778 to match the MUSICA δD over the common observation time period.
The use of the new data set is demonstrated in a case study of an atmospheric blocking event with a single TROPOMI overpass over northeastern Europe on 30 July 2018. Depleted air masses are found in the core of the anticyclone due to subsidence transporting upper tropospheric air towards lower levels. At the edge of the anticyclone a ring of enriched air is observed. A climatological study on the water vapour isotopic signature of continental summer blocking events could provide promising insights into the atmospheric water cycling associated with such systems that frequently lead to heat waves and hot temperature extremes. This case study shows the quality of the new data set and the added value for isotopologue studies, enabling studies on a day-by-day basis with high spatial resolution over continental regions.
Due to the restrictive filter for clear-sky scenes, the data coverage is limited. To improve on this, cloudy-sky retrievals over low clouds will be considered in a future study by using a forward model that accounts for scattering. Moreover, a calibration and validation of the TCCON HDO product is necessary. Additionally, it would be beneficial if recent NDACC-MUSICA data became available. Finally, an improvement in the consistency between the networks would be very valuable.