Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/ doi:10.5194/amt-9-2845-2016 Author(s) 2016. CC Attribution 3.0 License.
Accomplishments of the MUSICA project to provide accurate, long-term, global and high-resolution observations of tropospheric {H2O, D} pairs a review
Matthias Schneider1, Andreas Wiegele1, Sabine Barthlott1, Yenny Gonzlez2,3,a, Emanuel Christner1, Christoph Dyroff1,b, Omaira E. Garca2, Frank Hase1, Thomas Blumenstock1, Eliezer Seplveda2, Gizaw Mengistu Tsidu4,5, Samuel Takele Kenea4,c, Sergio Rodrguez2, and Javier Andrey6,d
1Institute of Meteorology and Climate Research (IMK-ASF), Karlsruhe Institute of Technology, Karlsruhe, Germany
2Izaa Atmospheric Research Center, Agencia Estatal de Meteorologa (AEMET), Santa Cruz de Tenerife, Spain
3Sieltec Canarias, S.L., Hbitat 2, 38204, San Cristbal de La Laguna, Santa Cruz de Tenerife, Spain
4Department of Physics, Addis Ababa University, P.O. Box 1176, Addis Ababa, Ethiopia
5Botswana International University of Technology and Science (BIUST) Priv. Bag 16, Palapye, Botswana
6Area de Investigacin e Instrumentacin Atmosfrica, INTA, Torrejn de Ardoz, Spain
anow at: Dept. of Earth, Atmospheric, and Planetary Sciences, Massachusetts Institute of Technology,77 Massachusetts Avenue, Cambridge, MA 02139-4307, USA
bnow at: Aerodyne Research Inc., 45 Manning Road, Billerica, MA 01821, USA
cnow at: Department of Physics, Samara University, P.O. Box 132, Samara, Ethiopia
dnow at: CNRM-GAME, Mto France and CNRS, Toulouse, France
Correspondence to: Matthias Schneider ([email protected])
Received: 29 October 2015 Published in Atmos. Meas. Tech. Discuss.: 18 January 2016 Revised: 31 May 2016 Accepted: 4 June 2016 Published: 7 July 2016
Abstract. In the lower/middle troposphere, {H2O, D} pairs are good proxies for moisture pathways; however, their observation, in particular when using remote sensing techniques, is challenging. The project MUSICA (MUlti-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water) addresses this challenge by integrating the remote sensing with in situ measurement techniques. The aim is to retrieve calibrated tropospheric {H2O, D} pairs from the middle infrared spectra measured from ground by FTIR (Fourier transform infrared) spectrometers of the NDACC (Network for the Detection of Atmospheric Composition Change) and the thermal nadir spectra measured by IASI (Infrared Atmospheric Sounding Interferometer) aboard the MetOp satellites. In this paper, we present the nal MUSICA products, and discuss the characteristics and potential of the NDACC/FTIR and MetOp/IASI {H2O, D} data pairs.
First, we briey resume the particularities of an {H2O, D} pair retrieval. Second, we show that the remote sensing data
of the nal product version are absolutely calibrated with respect to H2O and D in situ prole references measured in the subtropics, between 0 and 7 km. Third, we reveal that the {H2O, D} pair distributions obtained from the different remote sensors are consistent and allow distinct lower/middle tropospheric moisture pathways to be identied in agreement with multi-year in situ references. Fourth, we document the possibilities of the NDACC/FTIR instruments for climatological studies (due to long-term monitoring) and of the MetOp/IASI sensors for observing diurnal signals on a quasi-global scale and with high horizontal resolution. Fifth, we discuss the risk of misinterpreting {H2O, D} pair distributions due to incomplete processing of the remote sensing products.
Published by Copernicus Publications on behalf of the European Geosciences Union.
2846 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
1 Introduction
Atmospheric moisture (condensed water and vapour) strongly interacts with solar as well as thermal radiances and distributes energy in the form of latent heat. In consequence, it has a wide impact on the atmospheric energy budget and strongly affects circulation on regional and global scales. The insufcient understanding of tropospheric moisture pathways and their coupling to atmospheric circulation is seen as a major challenge for climate system modelling (e.g.Stevens and Bony, 2013), and a focus on three specic research areas is recommended (http://www.wcrp-climate.org/grand-challenges
Web End =http://www.wcrp-climate.org/ http://www.wcrp-climate.org/grand-challenges
Web End =grand-challenges http://www.wcrp-climate.org/grand-challenges
Web End = ): rst, the low cloud feedback and the responsible physical processes (modelled equilibrium climate sensitivity is strongly linked to the low cloud feedback, e.g.Sherwood et al., 2014); second, the climate response of large-scale tropospheric circulation systems and precipitation patterns (whereby palaeoclimate archives offer valuable possibilities, e.g. Ortega et al., 2015); third, the coupling between small-scale processes and large-scale dynamics (e.g. cloud processes that take place on diurnal timescale can signicantly affect large-scale circulation, Marsham et al., 2013).
Simultaneous observations of different tropospheric water isotopologues can aid advancements in these research areas. The ratio between the different water isotopologues provides information on processes related to moisture up-take, exchange, clouds and atmospheric transportation up-wind of the detected air mass (e.g. Dansgaard, 1964; Gat, 2000; Yoshimura et al., 2004), thereby offering potential for investigating the coupling between atmospheric circulation and moisture pathways.
Furthermore, the water isotopologue ratios archived in ice cores, speleothems, lake sediments or tree rings contain information about past climate conditions (e.g. Ortega et al., 2015), whose reconstruction, however, relies on a comprehensive understanding of the linkages between the archived isotopologues, on the one hand, and the tropospheric isotopologues, temperature and circulation, on the other hand.
The isotopologue ratios are typically expressed in the notation, which relates the observed ratio to the standard ratio VSMOW (Vienna Standard Mean Ocean Water). For instance, the HDO / H2O ratio is typically expressed as D =
HDO/H2O VSMOW
1. Here and in the following we use H2O and
HDO as equivalent to H162O and HD16O, respectively. The consideration of other isotopologues will be specied explicitly (e.g H182O or H172O).
During the last years there has been great progress in observations of the tropospheric water vapour isotopologues, whereby remote sensing observations are particularly interesting since they can be performed continuously (for cloud-free conditions). Ground-based remote sensing can offer long-term data records. There are the ground-based water vapour isotopologue remote sensing retrievals using the NDACC (Network for the Detection of Atmospheric Com-
position Change) middle infrared spectra (Schneider et al., 2012, and references therein) and retrievals that use the TCCON (Total Carbon Column Observing Network) near-infrared spectra (e.g. Rokotyan et al., 2014). Tropospheric water vapour isotopologue data sets have also been presented using space-based sensors. There are different research groups using short-wave infrared (SWIR) spectra measured by the satellite sensors SCIAMACHY (Frankenberg et al., 2009; Scheepmaker et al., 2015) or GOSAT (Boesch et al., 2013; Frankenberg et al., 2013) as well as the thermal nadir spectra of TES (Thermal Emission Spectrometer; Worden et al., 2006, 2012) or IASI (Infrared Atmospheric Sounding Interferometer; Schneider and Hase, 2011;Lacour et al., 2012; Wiegele et al., 2014).
While in the dry upper troposphere and stratosphere D observations alone allow signicant conclusions on the moisture pathways from the troposphere to the stratosphere and on stratospheric circulation (e.g. Kuang et al., 2003; Stein-wagner et al., 2010), the situation is different in the lower and middle troposphere. There, humidity is much more variable and the moisture pathways can be best investigated by analysing the distribution of the {H2O, D} pairs (e.g.
Galewsky et al., 2005; Noone, 2012; Gonzlez et al., 2016).
Recently, there have been a variety of publications that use remote sensing observations of tropospheric {H2O, D} pairs for tropospheric moisture pathway studies: for instance, for estimating the importance of rain recycling (e.g. Worden et al., 2007), for investigating the dynamics of the Madden-Julian oscillation (e.g. Berkelhammer et al., 2012; Tuinenburg et al., 2015) or for drawing conclusions on vertical mixing processes (e.g. Noone, 2012; Risi et al., 2012a; Sutanto et al., 2015). However, there are very few studies so far where attempts have been made to empirically validate these {H2O, D} pairs (Schneider et al., 2015; Lacour et al., 2015). Further and more detailed validation efforts for the {H2O, D} pairs are urgently needed because an estimation of {H2O, D} pairs is complex and not the same as an individual optimal estimation of H2O and HDO (nor of H2O and D). The reason is that the sensitivities for H2O and HDO (and also for H2O and D) are generally different (Schneider et al., 2006; Worden et al., 2006; Schneider et al., 2012).
A further pending detail with the isotopologue remote sensing data is the unclear bias in D, which can signicantly compromise their scientic usefulness (e.g. Risi et al., 2012b; Field et al., 2014; Yoshimura et al., 2014). First empirical bias assessment studies were presented by Schneider and Hase (2011) and Worden et al. (2011). For a reliable bias documentation, we need vertical isotopologue reference proles measured by well-calibrated in situ instrumentation (Herman et al., 2014; Dyroff et al., 2015). Currently we are only aware of one campaign (the summer 2013 MUSICA (MUlti-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water) campaign, Dyroff et al., 2015), where such proles are measured in coincidence with ground- and space-based remote sensing observations
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2847
and over the wide altitude range where the remote sensors are sensitive.
Removing the shortcomings in tropospheric water vapour isotopologue remote sensing data has been a focus of the project MUSICA (http://www.imk-asf.kit.edu/english/musica.php
Web End =http://www.imk-asf.kit.edu/english/ http://www.imk-asf.kit.edu/english/musica.php
Web End =musica.php ). In this paper, we summarize the nal results of the MUSICA project, whose most relevant papers are collected in Appendix A. We demonstrate that calibrated, long-term, global and high-resolution remote sensing data of tropospheric {H2O, D} pairs can be produced. In Sect. 2, we give a brief review on the theory of {H2O, D} pair estimations and of the MUSICA NDACC/FTIR (Fourier transform infrared) and MetOp/IASI retrievals, in particular. Section 3 documents that the nal MUSICA data are well-calibrated with respect to in situ H2O and D references. Section 4 shows an empirical validation of the {H2O, D} pair distributions by intercomparing the {H2O, D} signals ob-served by in situ instruments, NDACC/FTIR, MetOp/IASI-A and MetOp/IASI-B for different atmospheric situations. Sections 5 and 6 document the possibilities offered by the unique long-term characteristics of the NDACC/FTIR data and by the unique spatial and temporal coverage of the MetOp/IASI observations. In Sect. 7 we discuss risks for defective interpretations of the {H2O, D} remote sensing data pairs. Section 8 gives a summary.
2 Remote sensing of water vapour isotopologues and their ratios
In this section, we give a very brief overview on the challenges of tropospheric water vapour isotopologue remote sensing. Then we resume the retrieval approaches as developed in preparation for and continuously improved during the project MUSICA. In addition to the MUSICA products, there are other ground- and space-based (non-MUSICA) tropospheric water vapour isotopologue remote sensing products. A brief overview and a short discussion of the differences to the MUSICA products is given in Appendix B.
2.1 The challenge
In situ instruments analyse a clearly dened air mass and from the H2O and HDO measurements the ratio ( D) can be directly calculated. This is more complicated for remote sensing observations. There, the amount of H2O and HDO along the line of sight is retrieved from the measured spectra.The sensitivity of the retrieval depends on the noise in the spectra and on the shape and strength of the different spectral lines. Typically, the sensitivity for H2O is different to the sensitivity for HDO and the retrieved H2O value represents a different altitude range than the retrieved HDO value (H2O and HDO have different averaging kernels). In consequence, the D values calculated from individual H2O and
HDO retrievals can be rather misleading. Instead of individual retrievals, a combined H2O and HDO retrieval is needed.
A logarithmic-scale retrieval of H2O and HDO together with a constraint of ln[HDO] ln[H2O] has been proposed
by Schneider et al. (2006) and Worden et al. (2006) for generating D values by remote sensing techniques. This approach means actually an optimal estimation of (ln[H2O] +
ln[HDO])/2 and ln[HDO] ln[H2O], which are good prox
ies for H2O and D (Schneider et al., 2012); i.e. it is a quasi-optimal estimation of H2O and D. However, it is an individual optimal estimation of H2O and D, and thus such a product is still not comparable to the {H2O, D} pairs obtained from the in situ measurements. The reason is that the remote sensing system is more sensitive to atmospheric H2O than to atmospheric D and in addition there is a slight cross-dependency of the retrieved D on atmospheric H2O (e.g.
Schneider et al., 2012).
2.2 The MUSICA {H2O, D} pair product
The aforementioned problems can be overcome by an a posteriori processing of the retrieval output. The result of the a posteriori processing is an estimation of {H2O, D} pairs that are comparable to the {H2O, D} pairs obtained from in situ measurements. In the following, we give a very brief explanation of the necessity and functionality of the a posteriori processing.
The remote sensing retrievals produce state vectors x, averaging kernels A, error covariance matrices S, etc., for the {ln[H2O],ln[HDO]} basis system (or for the
{ln[H162O],ln[H182O],ln[HD16O]} basis system, if isotopo
logues with different oxygen atoms are distinguished). However, besides H2O (99.7 % in form of H162O), we are actually interested in D, and eventually deuterium excess (d = D
8 18, with 18 = H
18
2 O/H162O
VSMOW
1). These parameters, their variations and uncertainties can be captured and characterized in an elegant manner by using the {H2O, D} (or {H2O, D,d})
proxy basis system: the basis {(ln[H2O] + ln[HDO])/2} (or
{(ln[H162O] + ln[H182O] + ln[HD16O])/3}, if we distinguish
between the isotopologues with different oxygen atoms) well captures the parallel variations of the different isotopologues and is thus a good proxy for the dominant H2O variations.
Variations in D are 1 order of magnitude smaller, for which the basis {ln[HDO]ln[H2O]} (or {ln[HD16O]ln[H162O]})
is a good proxy. Variations in deuterium excess are a further order of magnitude smaller, for which we can use the proxy basis {7ln[H162O]8ln[H182O]+ln[HD16O]} (Barthlott et al.,
2016a).
The transformation of this proxy basis system can easily be realized by a transformation operator P, and the transformed state vector, averaging kernel and any covariance matrix (x[prime], A[prime] and S[prime]) are then given by
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2848 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
x[prime] = PxA[prime] = PAP1
S[prime] = PSPT . (1)
In the proxy basis system, the different sensitivities with respect to H2O and D and the cross-dependencies become clearly visible. We can remove these decits by the a posteriori processing, which consists of simple matrix multiplications (using the a posteriori correction operator C). The a posteriori corrected state vector, averaging kernel and any covariance matrix (x , A and S ) can be calculated as
x = C(x[prime] x[prime]a) + x[prime]a A = CA[prime]
S = CS[prime]CT . (2)
In the MUSICA papers, this a posteriori processed product is also often called the Type 2 product, and since it is most useful for isotopologue studies (H2O and D have almost the same averaging kernels), it will be the product that is generally used in this work. In addition, MUSICA offers a Type 1 product which is not a posteriori processed and it is the kind of product that is generally distributed by other data producers. It consists of optimal H2O data, but has limited possibilities for isotopologue studies (see discussion in Sect. 7).
Readers that are interested in more details about the proxy state method, the a posteriori correction and the operators P and C as used in Eqs. (1) and (2) are recommended to study Schneider et al. (2012), Wiegele et al. (2014) and Barthlott et al. (2016a). However, we would like to note that the MUSICA remote sensing data users do not have to be concerned about details of the a posteriori processing. The processing does not have to be performed by the data users because the data are provided as an a posteriori processed product (Type 2) and also in the form of the direct retrieval output (H2O and HDO or H2O and D products that have not been a posteriori processed, Type 1 product).
2.3 The MUSICA ground-based products
The MUSICA ground-based remote sensing retrieval uses the PROFFIT retrieval code (Hase et al., 2004) and the middle-infrared spectra recorded within the NDACC (http://www.acom.ucar.edu/irwg/
Web End =www. http://www.acom.ucar.edu/irwg/
Web End =acom.ucar.edu/irwg/ ). The NDACC spectra are of very high spectral resolution (0.005 cm1) and offer H162O, H182O and
HD16O absorption lines of similar strength.
For the nal MUSICA retrieval version (v2015), we perform an optimal estimation of (ln[H162O] + ln[H182O] + ln[HD16O])/3, ln[HD16O] ln[H162O] and 7ln[H162O]
8ln[H182O] + ln[HD16], which are good proxies for H2O, D
and d (deuterium excess). In addition to the cross-constrained
t of the water vapour isotopologues H162O, H182O and HD16O, we perform simultaneous but individual ts (no cross-constraints) for proles of the water vapour isotopologue H172O, the temperature and the interfering species CO2,
O3, N2O, CH4 and HCl. The code PROFFIT uses the root mean square value of the residual (difference between simulated and measured radiances) as noise level for constraining the retrievals. We use HITRAN 2012 parameters (Rothman et al., 2013) optimized for speed-dependent Voigt line parameterization. Details on the general retrieval set-up are given in Schneider et al. (2012) and the modications made for v2015 are summarized in Barthlott et al. (2016a).
The Type 1 MUSICA NDACC/FTIR product provides water vapour proles for the lower, middle and upper tropo-sphere (DOFS, degree of freedom of signal, of almost 3). The Type 2 product (the a posteriori corrected product) offers consistent {H2O, D} pairs, which are sensitive to the lower and the middle troposphere, whereby it is possible to reasonably separate both altitude regions (DOFS of about 1.7). This is illustrated in Fig. 1, which shows the rows of the averaging kernel matrix (A ) in the {H2O, D,d} proxy basis system and after applying the a posteriori correction (see Eq. 2). The full averaging kernels matrix consists of nine blocks, each of which is a {nol [notdef] nol} matrix, whereby {nol} is the number
of the vertical atmospheric grid points used for the retrieval. In total, A has the dimension {(nol [notdef] 3) [notdef] (nol [notdef] 3)}:
A = 0
@
A 11 A 12 A 13 A 21 A 22 A 23
A 31 A 32 A 33
1
A. (3)
The three blocks along the diagonal describe the direct responses; i.e. they represent the averaging kernels for the H2O, D and d proxies. Figure 1 demonstrates that the sensitivity with respect to H2O and D is very similar (compare plots showing the entries of A 11 and A 22), which is achieved by the a posteriori processing. The outer diagonal blocks describe the cross-responses, whereby in Fig. 1 the respective x axes are scaled, thereby accounting for the different magnitudes of the H2O, D and d variations. Concerning the cross-responses between H2O and D, we have to consider that ln[H2O] variations are 1 order of magnitude larger than D
variations. This means that the entries in the A 12 block must be 10 times larger than entries in the A 11 block in order to be of similar importance. Vice versa, entries in the A 21 block can be 1 order of magnitude smaller than entries in the A 22 block and still have a similar importance. The blocks A 31,
A 32 and A 33 describe how the retrieved d proxy responds to real atmospheric variations. Although there is sensitivity with respect to real atmospheric deuterium excess (see entries of A 33), the retrieved d signals are mainly caused by variations in H2O and D (signicant entries in blocks A 31 and A 32). The situation can be improved by an alternative a posteriori processing being dedicated to {H2O, D,d} triplets (see Sect. 4.4.2 in Barthlott et al., 2016a), which is, however,
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2849
2.4 The MUSICA space-based products
For the MUSICA MetOp/IASI retrievals, a nadir version of
PROFFIT is used. The basic retrieval set-up is presented and analysed in detail in Schneider and Hase (2011) and Wiegele et al. (2014). It has been developed in consistency with the NDACC/FTIR retrieval. It uses a broad spectral window (11901400 cm1), and performs an optimal estimation of the humidity and D proxies ((ln[H2O]+ln[HDO])/2
and ln[HDO] ln[H2O]). Simultaneously, we retrieve sur
face skin temperature and atmospheric temperature, as well as the interfering species CH4, N2O, CO2 and HNO3. Our
MetOp/IASI retrieval does not distinguish between H162O and H182O because the spectroscopic signatures of the latter are very weak in the IASI spectra. So, we treat all water molecules with two hydrogen atoms as a single molecule (in the following referred to as H2O), whose variability is mainly the one of H162O. For our atmospheric temperature t, we constrain strongly towards the EUMETSAT L2 atmospheric temperature. The root mean square value of the residual (dif-
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
Figure 1. Row entries of the nine blocks of the NDACC/FTIR full averaging kernel matrix in the {H2O, D,d} proxy basis system and after a posteriori correction (A as given in Eq. 3). The kernel is for the retrieval of a ground-based FTIR spectrum measured on 09:41 UT on 24
Jul 2013 (the retrieval result of this measurement is one of the data points shown in Fig. 3).
not subject of this paper. For this paper and in order to understand the {H2O, D} pair signals, it is sufcient to work with the kernel blocks A 11, A 12, A 21 and A 22.
Error estimations are discussed in detail in Schneider et al. (2012). The random error is about 2 % for H2O and 25 for D, whereby the leading error source is uncertainty in the atmospheric temperature proles and artefacts in the spectral baseline (like channelling or offset). Systematic errors are dominated by uncertainties in the spectroscopic parameters. Already for small uncertainties of 1 and 2 % for line intensity and pressure broadening parameters, the systematic errors in H2O and D proles can reach 10 % and 150 , respectively. An empirical study indicates that for the v2015 data, the bias is actually much smaller (see Sect. 3). In addition, for an atmosphere with ne vertical structures, the D cross-dependency on H2O can reach 15 (see Fig. 9 in Schneider et al., 2012).
2850 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
ference between simulated and measured radiances) is used as noise level for constraining the retrievals.
The previous version of the MUSICA MetOp/IASI retrieval (Schneider and Hase, 2011; Wiegele et al., 2014) worked with HITRAN 2008 spectroscopic water vapour line parameters (Rothman et al., 2009). For the nal MUSICA retrieval version (v2015) we use the HITRAN 2012 parameters (Rothman et al., 2013) and modied the line intensities (S) for the HDO absorption signatures by +10 %. We also
tested modications of S for the H2O signatures and changes of air (pressure broadening parameter), but nally found that a modication of S of HDO works most effectively for correcting biases in D. In this context, we would like to remark that the bias correction as suggested for TES (Worden et al., 2011; Herman et al., 2014) is also consistent with a positive change of S of HDO.
The retrieval provides a Type 1 product of H2O proles that is sensitive to variations between the surface and about 15 km altitude (DOFS of about 4). The Type 2 product (consistent {H2O, D} pairs) has a typical DOFS of 0.71.0 (after ltering according to Sect. 6.2), whereby the sensitivity is mainly limited to the middle troposphere. This is illustrated in Fig. 2, which plots the rows of the averaging kernel matrix (A ) in the {H2O, D} proxy basis system and after applying the a posteriori correction. For the MetOp/IASI retrieval we have no deuterium excess basis and the dimension of A is {(nol [notdef] 2) [notdef] (nol [notdef] 2)}
A = [parenleftbigg]
A 11 A 12 A 21 A 22
. (4)
As for the NDACC/FTIR kernels, A 11 and A 22 describe the sensitivities with respect to H2O and D, and A 12 and
A 21 their cross-responses, respectively.
The random error is about 5 % for H2O and 20 for D, whereby atmospheric temperature uncertainties, thin elevated clouds and noise in the spectra are the leading uncertainty sources. The systematic error is dominated by uncertainties in the spectroscopic parameters. It can easily reach 5 % and 50 for H2O and D, respectively (in case of a 5 % uncertainty in the spectroscopic line intensity parameters). However, an empirical study as shown in Sect. 3 reveals that for v2015 the bias is actually much smaller. In addition, for an atmosphere with ne vertical structures, the D cross-dependency on H2O can occasionally be larger than 40 (see Fig. 5 in Wiegele et al., 2014).
2.5 Uniform a priori settings
For v2015, we work with the same globally constant water vapour isotopologue a priori data for all ground-based retrievals (different globally distributed NDACC/FTIR stations) and for all MetOp/IASI retrievals (for the whole globe). Thereby, we assure that observations at different locations are not affected by the use of different a priori data, and therefore an interpretation of regional differences (e.g. lati-
Figure 2. Row entries of the four blocks of the MetOp/IASI full averaging kernel matrix in the {H2O, D} proxy basis system and after a posteriori correction (A as given in Eq. 4). It is for a retrieval of an IASI spectrum measured on 11:07 UT on 24 Jul 2013 close to the southern coast of Tenerife (the retrieval result of this measurement is one of the data points shown in Fig. 4).
tudinal gradients) becomes rather straightforward. This is a further development of the previous retrieval version, where we used different a priori data for the (sub)tropics, the mid-latitudes and the polar regions.
The a priori proles used are mean values of LMDZiso calculations (Risi et al., 2012b) and are as depicted in Fig. 2 of Lacour et al. (2012). As H2O a priori variability, we assume 75 % in the boundary layer, 150 % in the middle and upper troposphere and 30 % in the stratosphere. For D, the respective variability values are 60 , 120 and 50 . The a priori covariances are then calculated by assuming a correlation length of 2 km in the boundary layer, 4 km in the middle and upper troposphere and 8 km in the stratosphere.
3 Calibrated remote sensing products
The summer 2013 MUSICA campaign generated reliable H2O and HDO (and D) in situ reference proles on 6 days between the sea surface and 67 km altitude (Dyroff et al., 2015; Schneider et al., 2015). From the D proles we calculate 18 and then H182O proles by assuming D = 8 [notdef] 18
(H182O is needed for validating the v2015 NDACC/FTIR product). The in situ proles have been measured in coincidence with high-resolution ground-based FTIR observations and with IASI observations and are unique for documenting the bias in the remote sensing data because they cover the
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2851
Figure 3. Correlation between the reference proles (in situ data measured by the ISOWAT instrument for 07 km and climatology above, then smoothed with FTIR kernels) and FTIR data for FTIR measurements made in the morning (08:1509:45 UT), i.e. before the aircraft ights, but reasonably representative of the free troposphere. Left panels are for H2O and right panels are for D. Upper panels for retrievals at 2.4 km and bottom panels for retrievals at 5 km. The black line is the 1 : 1 diagonal. The error bars represent the uncertainty estimations
for the reference and FTIR data. Note that for the comparison at 5 km a large part of the uncertainty in the reference data is due to the fact that there are no ISOWAT measurements above the aircrafts ceiling altitude (Schneider et al., 2015, and discussion in Appendix C).
whole altitude range from the surface up to 67 km. In this section, we show that the MUSICA v2015 remote sensing products are well calibrated with respect to these reference data.
3.1 NDACC/FTIR
In a rst study with data from the previous MUSICA
NDACC/FTIR retrieval version, we found a bias of 25 70 for D with respect to the prole references (see right panels in Figs. 9 and 10 of Schneider et al., 2015). That study was done with a limited data set (only one exemplary remote sensing observation per day). For the v2015 product we perform a comprehensive empirical bias assessment and compare the ground-based FTIR data obtained for all optimal coincidences (about 10 observations each day between 08:15 and 09:45 UT) with the reference data (one prole per day:21 Jul 2013, 22 Jul 2013, 24 Jul 2013, 25 Jul 2013, 30 Jul 2013 and 31 Jul 2013).
Figure 3 shows the plots for the correlations between the reference and the FTIR data. The upper panel depicts the comparison for the lower free troposphere and the bottom panel for the middle free troposphere. The references are constructed from the in situ prole measurements (surface up to ceiling altitude of 67 km) and a climatology for higher altitudes by convolution with the FTIR averaging kernels. The climatology above the ceiling altitude is the same as the a priori data as discussed in Sect. 2.5. The technical details for this comparison, like the need for applying the averaging kernels to the reference proles, are the same as for the exemplary study of Schneider et al. (2015). The error bars on the reference data are largely due to the unknown humidity and D values above the aircrafts ceiling altitude. This can easily be understood from the averaging row kernels of the matrix blocks A 11 and A 22 as shown in Fig. 1. They reveal that the atmospheric state above 67 km does affect the retrieval for lower altitudes. Appendix C provides a brief discussion on the importance of reaching high ceiling altitudes.
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2852 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
Table 1. Empirical validation and bias assessment of the NDACC/FTIR H2O and D products in the lower troposphere.
Sensor Altitude Number of Number of R2/slope Mean bias range remote sensing reference from t [notdef] condence
observations observation (NR) H2O D H2O (%) D ()
NDACC/FTIR 2.44 km 65 6 70 %/0.86 68 %/0.86 +2.1 [notdef] 12.4 12.1 [notdef] 16.6
Table 2. Empirical validation and bias assessment of the NDACC/FTIR and MetOp/IASI H2O and D products in the middle troposphere.
Sensor Altitude Number of Number of R2/slope Mean bias range remote sensing reference from t [notdef] condence
observations observation (NR) H2O D H2O (%) D ()
NDACC/FTIR 37 km 65 6 86 %/0.96 91 %/0.95 0.8 [notdef] 8.2 2.7 [notdef] 7.4
MetOp/IASI 28 km 10 4 97 %/0.89 88 %/0.74 +0.6 [notdef] 3.7 +8.5 [notdef] 7.9
We observe a reasonable correlation and no signi-cant bias between the reference data and the MUSICA NDACC/FTIR v2015 data. We also see variations in the FTIR data measured between 08:15 and 09:45 UT on a single day, which cannot be attributed to changes in the averaging kernels because there is no similar variation seen in the smoothed reference data (there is only a single reference prole per day). This variation is seen in the lower free tropo-sphere and in the middle free troposphere and is very likely a true variation in the free tropospheric humidity and D elds.
The relatively high variability of the atmospheric state is a problem when comparing the different measurements because we do not know whether the different measurements detect air masses with the same atmospheric characteristics well. Our experience is that in the surroundings of Tenerife, the middle tropospheric humidity elds can be reasonably compared if they are made within 2 h (Schneider et al., 2010a, b) and within a horizontal distance of about 100 km (Wiegele et al., 2014). For these coincidence criteria, the mismatch between lower and middle tropospheric H2O and D can be assumed to be within 10 % and 10 , respectively.However, the comparison of aircraft measurements (free troposphere) with measurements made at Izaa (on a mountain ridge) is more difcult due to the local thermal circulation that starts on the island during the morning hours (Schneider et al., 2010b; Gonzlez et al., 2016). For that reason, we compare the in situ aircraft references with FTIR data measured in the early morning hours (between 08:15 and 09:45 UT) and dene these measurements the optimal coincidences with the free tropospheric aircraft measurements made between 10:30 and 13:30 UT. As a consequence the mismatch uncertainty for the comparison between aircraft and FTIR data will be larger (within 30 % for H2O and very likely within 30 for D). For a detailed discussion on optimal coincidences between the aircraft-based in situ and the remote sensing measurements, please refer to Appendix C.
The Tables 1 and 2 resume the results of the empirical validation and bias assessment. The reference and the NDACC/FTIR products detect the same variations to a large extent (R2 values of about 70 % for the lower and almost 90 % for the middle troposphere). The calculated mean NDACC/FTIR biases are not signicant if we take the condence range of these assessments into account. The condence range is calculated as the standard deviation of the bias divided by pNR 1 (with NR being the number of indepen
dent reference observations). For the lower troposphere, we are able to determine the bias with a condence of 12.4 % and 16.6 for H2O and D, respectively. In the middle troposphere, the assessment is even more reliable. There, the condence ranges are 8.2 % and 7.4 for H2O and D, respectively. These condence ranges are in good agreement with the estimated mismatch uncertainties. In summary, the lower tropospheric bias is very likely somewhere between
10 and +15 % for H2O and between 30 and +5 for
D. The middle tropospheric bias is very likely between 10
and +8 % for H2O and between 10 and +5 for D.
3.2 MetOp/IASI
In an exemplary study with the previous MUSICA MetOp/IASI retrieval version, Schneider et al. (2015) reported a bias between the IASI D product and the D reference of about 60 (see Fig. 11 therein). For the comprehensive bias assessment of the v2015 data, we follow the procedure as described in the context of that exemplary study.
Figure 4 shows the correlation plots between the reference and the MetOp/IASI v2015 data. The data points plotted by smaller symbols and by black error bars are for non-ideal coincidences, and the rest of the data points are for good coincidences (Appendix C gives more details on the coincidences and the error bars). For this comparison, we smooth the reference prole data (in situ measurements and climatological data above ceiling altitude, as described in Sect. 3.1) with the respective IASI averaging kernels. Since there is
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2853
Figure 4. Same as Fig. 3 but for the correlation between reference data (in situ data measured by the ISOWAT instrument for 07 km and climatology above, then smoothed with IASI kernels) and IASI data for 5 km altitude. The colour code is for the different days. The data points for non-optimal coincidences as discussed in Appendix C can be identied by the smaller symbols and the black error bars.
only one reference prole per day, the small variations in the smoothed reference data on an individual day must be due to varying averaging kernels (for instance, the small variability of the green dots for 24 Jul 2013 in parallel to the x axis). The respective variations of the IASI data (variations in parallel to the y axis) are due to variations in the sensitivity of IASI (variation in the kernels) and due to variations in the real atmospheric state encountered at the different observational pixels. The in situ data and the remote sensing data observed for good coincidences are well correlated. We estimate a mismatch uncertainty of about 10 % and 10 for the H2O and D comparisons, respectively (for more details see discussions in Sect. 3.1 and Appendix C).
We are able to assess the bias for the middle tropospheric
IASI data with a condence of 3.7 % and 7.9 for H2O and D, respectively (see Table 2). The condence ranges are in agreement with the expected mismatch uncertainties. The obtained mean bias values lie within these condence ranges (for H2O) or are only very slightly outside this range (for D), meaning that the actual biases are very likely between
2.5 and +5 % for H2O and between 0 and +15 for D.
4 Validation of {H2O, D} pair distributions
An individual validation of H2O and D (as shown in the previous section) is important for documenting that the data are calibrated to the reference scales for H2O and D. However, it is not sufcient. It is the tropospheric {H2O, D} pair distribution that gives insight into tropospheric moisture pathways and therefore, it is the distribution of these pairs that has to be validated. Wiegele et al. (2014) and Schneider et al. (2015) presented approaches for such kinds of validation exercises. Here, we present a further rened {H2O, D} pair validation and compare in situ, MUSICA v2015 NDACC/FTIR, MUSICA v2015 MetOp/IASI-A and MUSICA v2015 MetOp/IASI-B data sets.
4.1 Moisture pathways to the North Atlantic subtropical free troposphere
In the surroundings of Tenerife there are three distinct moisture transport pathways that control free tropospheric humidity. Gonzlez et al. (2016) showed that these three pathways have a distinct {H2O, D} pair distribution, which offers a unique opportunity for validating the different {H2O, D} pair data sets.
Generally, the free troposphere in the subtropics receives air that has been transported from higher latitudes and altitudes and subsides to the subtropics (e.g. Galewsky et al., 2005). In the following, we call this pathway ATL, desc. However, the summertime free troposphere close to West Africa is often affected by the Saharan air layer (SAL). The SAL is a well-mixed planetary boundary layer that can expand up to 67 km and has its origin in the strong vertical mixing (dry convection) over the summertime Sahara. This dry convection process mixes boundary layer air with free tropospheric air. The SAL is then often advected westward over the Atlantic, where it can be identied by high dust concentrations (Rodrguez et al., 2011) and increased humidity levels. The free troposphere above Tenerife is also particularly humid when the air has been transported from lower altitudes over the tropical/subtropical Atlantic (Gonzlez et al., 2016). This mainly occurs in the late summer and early autumn. We call this pathway ATL, asc in the following. The three distinct moisture transport pathways, their identication methods and the prevailing occurrences are summarized in Table 3 and discussed in great detail in Gonzlez et al. (2016).
For our validation purpose we sort the {H2O, D} observations according to the moisture pathway that has been prevailing during the observation. The sorting is done according to the identication method as given in Table 3 and for all the different data sets independently: the Picarro
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2854 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
Table 3. The dominant moisture pathways to the free troposphere in the surroundings of Tenerife.
Pathway Description Identication method Season
ATL, desc. Atlantic air mass low aerosol load (measurement) and low temperature at point of NovemberMay descending from high latitudes last condensation (trajectories), Gonzlez et al. (2016)
ATL, asc. Atlantic air mass low aerosol load (measurement) and high temperature at point of JulyOctober ascending from lower latitudes last condensation (trajectories), Gonzlez et al. (2016)
SAL Saharan air layer high aerosol load (measurement), July + August
advected over the Atlantic Rodrguez et al. (2011); Gonzlez et al. (2016)
Figure 5. Free tropospheric {H2O, D} pair distributions as obtained by different measurement techniques in the surroundings of Tenerife (for orientation see also map as shown in Fig. C1). The contour lines indicate areas with the highest densities of the {H2O, D} pairs: red denotes air masses that are clearly affected by dry convection over the African continent; blue and green denote Atlantic air masses with different pathways (see Table 3). The thin dashed and thick solid lines mark the areas that include 95 % and 66 % of all data, respectively. Left: two Picarro in situ instruments (L2120-i) measuring during nighttime at 2390 and 3550 m a.s.l. (Izaa Observatory and Teide). Middle: ground-based NDACC/FTIR located at Izaa. Right: space-based MetOp/IASI-A and MetOp/IASI-B observing in a 2 [notdef] 2 area south of the
island. In addition, the panels show three simulated curves (thin black lines): a Rayleigh curve for initialization with T = 25 C, RH = 80%
and D = 80 and two mixing curves (rst line for mixing between H2O[1] = 25 000 ppmv; D[1] = 80 and H2O[2] = 900 ppmv;
D[2] = 430 and second line for mixing between H2O[1] = 25 000 ppmv; D[1] = 80 and H2O[2] = 200 ppmv; D[2] = 610 ).
The yellow star marks the a priori value used for the remote sensing retrievals at 4.9 km.
surface-based in situ observations (at 2390 and 3550 m a.s.l.), the NDACC/FTIR measurements made from Izaa and the MetOp/IASI observations made close to Tenerife. Then we calculate the {H2O, D} pair density distribution for each observational technique and each moisture pathway ensemble. The results are the {H2O, D} pair contour plots as depicted in Fig. 5 (from the left to the right for the different data sets as described in the panels). We present these plots on a logarithmic scale for H2O and maintain the scale for D (the D scale is in a rst approximation the same as a ln[HDO] ln[H2O]
scale). These are the scales on which the optimal estimation of the remote sensing products is performed. This largely facilitates the interpretation of the remote sensing data because then the Type 2 kernels are very similar on the x and y scales. Furthermore, on these scales a Rayleigh process will become visible by a straight line.
For air descending from the Northern Atlantic (ATL, desc, blue contours), the data points are well distributed between a typical Rayleigh line (gradual dehydration due to condensation) and mixing lines (mixing of two end members with {H2O, D} on the exemplary Rayleigh line). These water masses have gone through different condensation and mixing processes.
For air ascending from the tropical/subtropical Atlantic (ATL, asc, green contours), the air is more humid and the {H2O, D} pairs generally group around the Rayleigh line with some data points lying signicantly below the Rayleigh line. These water masses are strongly depleted in HDO, which indicates rain re-evaporation or gradual dehydration (Rayleigh distillation) after evaporation over a warm ocean.
For SAL conditions (red contours), the air is also humid but HDO is signicantly enriched when compared to the typical Rayleigh distribution. This can be well explained by
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2855
Figure 6. {H2O, D} pair plots for 4.9 km altitude for coincident FTIR and IASI measurements for the three locations: Tenerife, Karlsruhe, and Kiruna. The FTIR data, the IASI data and FTIR data smoothed with the IASI averaging kernels are plotted (from the top to the bottom).
The colour code displays the upper 10 % and lower 10 % of D values as identied in the FTIR data. The black lines and the yellow stars are the same as in Fig. 5 (Rayleigh line, mixing lines and a priori value for 5 km altitude, respectively). For more details on this kind of validation approach, please refer to Wiegele et al. (2014).
the mixing of planetary boundary layer humidity with middle/upper tropospheric humidity.
We observe that the three different data sets reveal very similar {H2O, D} pair distributions depending on the history of the detected air mass. The study is statistically very robust, since it uses several hundred individual observations (number given by parameter N in the legends of Fig. 5) made on many different days (number given by parameter D). This conrms that the v2015 MUSICA remote sensing data are reasonably well bias corrected and proves that they are capable of tracking different moisture pathways. The smaller variability in the remote sensing data is due to the fact that they represent averages for layers of several kilometres (see averaging kernel plots of Figs. 1 and 2). For the MetOp/IASI data the variability is particularly small for dry air (blue contours) because the drier the atmosphere, the lower IASIs sensitivity for middle tropospheric {H2O, D} pairs.
4.2 {H2O, D} extremes on global scale
Figure 5 shows validations of {H2O, D} pairs for a subtropical site. In order to perform a similar study for other sites, we would need respective middle tropospheric in situ references, which are not available. The ISOWAT prole and surface-based Izaa and Teide Picarro in situ references as observed in the surroundings of Tenerife are unique and a
generation of similar data sets for middle or high latitudes would be expensive (it would require a large number of aircraft campaigns).
Here, we show a comparison between NDACC/FTIR,
MetOp/IASI-A and MetOp/IASI-B, which is of global validity. Our argument is that a global agreement between the different remote sensing data sets would suggest that the in situ validations made for Tenerife are of global validity.
For this kind of validation, we work with {H2O, D} extremes. For this purpose, we identify anomalous or extreme {H2O, D} distributions in one remote sensing data set and document to what extent these extremes are seen in another remote sensing data set (by comparing coincident observations). The validation approach with {H2O, D} extremes was rst presented by Wiegele et al. (2014), which should be consulted for more details.
4.2.1 NDACC/FTIR vs. MetOp/IASI
We compare the FTIR and IASI data for three rather different sites: Tenerife (subtropical Atlantic), Karlsruhe (central Europe) and Kiruna (northern Scandinavia). At these sites, we have ground-based FTIR observations of NDACC that contribute to MUSICA and we performed continuous IASI retrievals around the FTIR locations. Figure 6 shows the {H2O, D} distributions as retrieved at the three sites from
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2856 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
Figure 7. Similar to Fig. 6 but for coincidences of IASI-A and IASI-B measurements. The products retrieved at 4.9 km for all coincidences within 1 h and 0.25 [notdef] 0.25 are presented for 16 August 2014. The left panels show the morning overpasses and right panels the evening
overpasses. The anomalies are identied in the IASI-A observations (upper panels) and then checked in the IASI-B observations (bottom
panels).
coincident FTIR and IASI measurements (IASI observations made in an 110 [notdef] 110 km area south of the FTIR instru
ments). As temporal coincidence criterium we required that the two measurements were made within 1 h. The left column of plots shows data for Tenerife (coincidences between 2007 and 2013), the central column data for Karlsruhe (coincidences between 2010 and 2013) and the right column data for Kiruna (coincidences between 2007 and 2012). The rst row of plots depicts the FTIR data, the second row of plots shows the IASI data and the third row of plots the FTIR data smoothed by the IASI averaging kernels. In all plots we show retrievals for 4.9 km altitude. The grey data points represent all data. The FTIR observations that show unusual low or strong HDO depletion (high or low D values) are marked in red and green, respectively. These anomalies or extremes have been identied by a second-order least squares t to the {ln[H2O], D} distribution. The 10 % of the data points that
have the largest positive/negative D difference to the regression curve are dened as the extreme values.
First, comparing the {H2O, D} distribution relative to the unique a priori point, we see a good agreement between both data sets. In both data sets and from Tenerife via Karlsruhe to Kiruna, the water masses get generally more and more depleted in HDO and the {ln[H2O], D} slopes become more
and more shallow. Second, both data sets reveal very similar anomalies. If the FTIR observes an anomalously weak depletion, IASI also does (red dots in the IASI plots are situated at the upper end of the D distribution). The same is true for the anomalies with strong depletion (green dots).
In summary, the {H2O, D} distribution patterns (latitudinal gradients, {ln[H2O], D} slopes, anomalous distribu
tions) as observed in the MUSICA NDACC/FTIR and MU
SICA MetOp/IASI data are in good agreement. This nding suggests that the {H2O, D} pair validation as shown in the previous sections for the surroundings of Tenerife is valid for very different geophysical locations.
4.2.2 IASI-A vs. IASI-B
Since 2013 two IASI instruments (A and B) on two different satellites (MetOp-A and MetOp-B) have provided operational spectra. Their respective overpasses take place typically within 30 min, which offers very good conditions for cross-validating the IASI-A and IASI-B products.
Figure 7 depicts {H2O, D} distributions considering all valid observations on 16 August 2014 (left columns for the morning overpass and right columns for the evening over-pass). The colour code is the same as in Fig. 6. The grey data points show all data, the red data points mark the observations that have been identied in the IASI-A data as a positive D extreme and the green data points mark the observations that correspond to a negative IASI-A D extreme.
The top panels show the IASI-A data. These data are used for identifying the extremes and the red and green data points are of course separated. The bottom panels show the IASI-B data and green and red mark the IASI-B observations that are conducted in coincidence with the extreme IASI-A observations. The coincidence criteria were measurements within 1 h
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2857
Table 4. List of current MUSICA NDACC/FTIR sites (ordered from north to south) and available MUSICA data record of quality-ltered data. DOFS (degrees Of freedom of signal) for the optimal estimation of {H2O, D} pairs (trace of A 11 A 22; an example of these matrix
blocks is plotted in Fig. 1). This table is adopted from Barthlott et al. (2016a).
Site Location Altitude Data record No. of meas. (N) Meas. days (D) DOFS
Eureka, Canada 80.1 N, 86.4 W 610 m a.s.l. 20062014 1890 398 1.7 Ny lesund, Norway 78.9 N, 11.9 E 21 m a.s.l. 20052014 730 251 1.6
Kiruna, Sweden 67.8 N, 20.4 E 419 m a.s.l. 19962014 1981 969 1.6 Bremen, Germany 53.1 N, 8.9 E 27 m a.s.l. 20042014 582 316 1.6
Karlsruhe, Germany 49.1 N, 8.4 E 110 m a.s.l. 20102014 1756 425 1.6 Jungfraujoch, Switzerland 46.6 N, 8.0 E 3580 m a.s.l. 19962014 1884 1175 1.6
Izaa/Tenerife, Spain 28.3 N, 16.5 W 2367 m a.s.l. 20012014 9350 1210 1.7 Altzomoni, Mexico 19.1 N, 98.7 W 3985 m a.s.l. 20122014 1489 234 1.7
Addis Ababa, Ethiopia 9.0 N, 38.8 E 2443 m a.s.l. 20092013 528 154 1.6 Wollongong, Australia 34.5 S, 150.9 E 30 m a.s.l. 20072014 5834 927 1.6
Lauder, New Zealand 45.1 S, 169.7 E 370 m a.s.l. 19972014 3533 1653 1.6 Arrival Heights, Antarctica 77.8 S, 166.7 E 250 m a.s.l. 20022014 374 287 1.4
and within an area of 0.25 [notdef] 0.25 , whereby the compared
IASI-A and IASI-B data are often measured with rather different swath angles. We nd that IASI-B detects very similar D extremes as IASI-A, which demonstrates the good global consistency of the IASI-A and IASI-B {H2O, D} pairs, the robustness of the retrieval (compared are observations with different swath angles) and that our coincidence criteria used in the context of Sect. 3 are reasonable. In summary, we can use the IASI-A and IASI-B products as a uniform and consistent {H2O, D} pair data set.
5 Consistent long-term observation with NDACC/FTIR
Ground-based FTIR high-resolution solar absorption spectra have been measured within the NDACC for many years and can be used for generating long-term data sets of tropospheric {H2O, D} pairs (Schneider et al., 2012). For MU
SICA, the principal investigators of the individual FTIR stations send the spectra to the MUSICA retrieval team where they are centrally evaluated. This strategy assures highest consistency between the retrieval products for the different stations.
The NDACC/FTIR activities complement the surface-based in situ isotopologue monitoring activities. While the data obtained from the latter represent near-surface small-scale variations, which are often difcult for models to capture, the MUSICA NDACC/FTIR isotopologue data are representative of different altitudes and for processes that take place over deeper layers (the data represent vertical layers averaged over 25 km; see typical averaging kernels in Fig. 1). Due to their long-term data characteristics, the NDACC/FTIR data are particularly interesting for climatological studies.
5.1 Contributing stations and currently available data volume
The number of stations contributing to the MUSICA activities is gradually increasing and the data sets have been updated. Table 4 gives an overview of the 12 NDACC/FTIR sites that currently contribute to the MUSICA activities and the available data records. The stations are well distributed from the Arctic to the Antarctic and in some occasions have provided data since the late 1990s (a plot of time series until 2012 is given in Fig. 12 of Schneider et al., 2012). A further extension of the data set to other sites or for some stations to measurements made in the beginning of the 1990 is feasible, but has not been possible with the funds available for the MUSICA project.
5.2 Data ltering
All the data pass through a quality lter with different criteria. First, we require total DOFS for the three water vapour isotopologues (H162O, HD16O and H182O) of at least 4.0. Second, we analyse the position of solar lines with respect to terrestrial lines and require a [Delta1] / within 5 [notdef] 106 ([Delta1] is
the difference in the line shift of solar and terrestrial lines and the position of the solar line, both given in cm1).
This method allows observations made with incorrect pointing of the solar tracker to be excluded (Gisi et al., 2011).Third, we require that the tted phase error of the instrumental line shape does not change by more than 0.02 rad during the period of the data record. Thereby, we exclude data that have been measured under anomalous instrumental characteristics. Fourth, we perform XCO2 retrievals using the same spectra as for the water vapour isotopologue retrievals. We compare the retrieved XCO2 values with XCO2 as obtained from a multi-parameter model for XCO2 (Barthlott et al., 2015). We require that the measured and modelled XCO2 data agree within 2 %.
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2858 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
Figure 8. Example of seasonal {H2O, D} pair distribution climatologies obtained from FTIR retrievals for 4.9 km at Izaa/Tenerife (top panels) and Addis Ababa (bottom panels). The left panels show the distribution density plots, whereby the thin dashed and thick solid lines mark the areas that include 95 and 66 % of all data, respectively. The different colours correspond to the seasons as explained in the legends. The yellow star and black lines are the a priori and the simulated lines as in Figs. 5 and 6. The right panels show corresponding typical row kernels for 4.9 km for the H2O and D proxies as well as for their cross-responses (for more general details about the FTIR kernels, see Fig. 1 and corresponding discussions).
5.3 Seasonal {H2O, D} climatologies
Figure 8 gives an example of the seasonal cycles in the {H2O, D} distributions (around 5 km altitude) obtained from the observations made on Tenerife in the subtropical North Atlantic (14 years: 20012014) and in Addis Ababa in East Africa (4 years: 20092013).
For the subtropical North Atlantic (upper panels), the seasonal cycle can be explained by the seasonality of the prevailing moisture pathways as summarized in Table 3, so the {H2O, D} pair distribution plot shown here is very similar to the one shown in Fig. 5 (central panel). The right panels show typical row kernels for the three different seasons and
for the altitude of 4.9 km. Due to the a posteriori processing, the sensitivity with respect to H2O and D is almost identical (compare row kernels of A 11 and A 22). However, the sensitivities slightly depend on the season. It seems that during winter and late summer the 4.9 km retrieval is more sensitive to the lower troposphere than to the middle troposphere, whereas the situation is vice versa for retrievals of midsummer observations.
Over East Africa (bottom panel) the air is generally more humid and less depleted in HDO than over the subtropical North Atlantic. Furthermore, we observe different {H2O, D}
distributions for the different seasons. Directly after the rain season (mid-OctoberDecember, green contours) the vapour
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2859
tween 0 (nadir) and 48.3 , whereby the ground pixel at nadir has a diameter of 12 km. MetOp-A with IASI-A was launched in October 2006 and MetOp-B with IASI-B in September 2012. Currently, both IASI instruments are operative.
6.1 Spatial and temporal resolution, coverage and currently available data volume
For one morning or evening overpass the IASI-B swaths typically complement the area left out by the IASI-A swaths, and vice versa. Since the MUSICA IASI-A and IASI-B data are very consistent (see Fig. 7), we can treat them as a uniform data set and create extremely dense global data point maps for each daily morning and evening overpass. Figure 9 depicts typical maps for single day morning and evening over-passes during boreal winter and summer, respectively. The areas with missing data are cloudy areas or correspond to scenarios where the retrieval has rather low sensitivity.
Currently, MUSICA MetOp/IASI-A and MetOp/IASI-B data (Type 1 and Type 2 products) are available on a global scale for 6 days in February 2014 and for the whole months of August 2014. In addition, we have retrievals results for longer time periods and for 2 [notdef] 2 areas around the three
ground-based FTIR sites, Izaa (Tenerife), Karlsruhe and
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
Figure 9. Daily morning and evening coverage for the MUSICA IASI-A and IASI-B Type 2 retrieval product that passed the sensitivity lter (Eq. 5) for 4.9 km altitude and the retrieval quality lter. Left panels show the morning overpass and the right panels show the corresponding evening overpass. Top panels are for a typical situation during northern winter (example 16 February 2014) and bottom panels for northern summer (example 16 August 2014). The marked areas (blueish and reddish) are discussed in detail in the context of Fig. 10.
is most depleted in HDO. In January and February (blue contours) the air remains similarly dry, but becomes more enriched in HDO. Then, before the rain season (MarchJune, purple contours), the air gets more humid, but D gets only slightly enriched in HDO. The typical row kernels depicted in the right panels, reveal very similar sensitivities for all the different seasons.
The main intention of this gure is to briey demonstrate the potential of the NDACC/FTIR data for climatological {H2O, D} pair distribution analyses. A deeper scientic discussion of the climatological signals would need model calculations and is beyond the scope of this paper.
6 Quasi-global and high-resolution observations with
MetOp/IASI
IASI sensors are aboard the MetOp satellites, which is a series of three satellites (MetOp-A, MetOp-B and MetOp-C) for covering the time period from 2006 to the beginning of the 2020s. MetOp has 14 orbits per day at about 817 km altitude, which, together with the swath width of about 2200 km of the IASI instruments, leads to a quasi-global coverage of morning overpasses (at about 10:00 LT) as well as evening overpasses (at about 22:00 LT). The swath angles are be-
2860 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
-
A
A
Figure 10. Example of regional-scale signals (top panels), and seasonal- and diurnal-scale signals (bottom panels) in the {H2O, D} pair distribution as obtained from IASI-A and IASI-B spectra. The left panels show the distribution density plots for retrieval results of 4.9 km, whereby the thin dashed and thick solid lines mark the areas that include 95 and 66 % of all data, respectively. The different colours correspond to the regions, seasons and daytime as marked in Fig. 9 and as explained in the legends. The yellow star and black lines are the a priori and the simulated lines as in Figs. 5, 6 and 8. The right panels show the corresponding typical row kernels for 4.9 km for the H2O and D proxies as well as for their cross-responses (for more general details about the IASI kernels, see Fig. 2 and corresponding discussions).
Kiruna. More MUSICA MetOp/IASI retrievals are planned, but depend on future funding.
6.2 Sensitivity and retrieval quality lter
The height region around 5 km altitude (about 500 hPa) is generally most sensitive with respect to the {H2O, D} pairs.
However, occasionally (e.g. for an extremely dry or humid troposphere) the sensitivity peaks at other altitudes.
To lter out data with low sensitivity at a certain altitude, we set up a matrix in representation of the atmospheric covariances (the matrixs elements represent the different altitude levels). This matrix (S c) has unity values on the diago-
nal, and the outer diagonal elements are obtained by assuming an inter-level correlation length of 5 km. Then we calculate the error covariance in the retrieved data as
S [epsilon1] = (A I)S c(A I)T . (5)
Here, A is the averaging kernel matrix (see Eq. 4 and Fig. 2 for an example) and I the identity matrix. If we are interested in data that represent the atmosphere at altitude X, we do not consider retrievals, for which the respective diagonal element of S [epsilon1] is larger than 0.52, thereby requiring that at least 50 %
of the atmospheric variation at altitude X is seen in the retrieved data. In Fig. 9 we only plot data points for which this sensitivity criterion is fullled for the altitude of 4.9 km (i.e.
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
Ke
-
-
1217 Feb. 2014:
1217 Aug. 2014:
Ke
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2861
it only shows data that reect variations of a relatively deep layer at about 5 km altitude).
We lter out poor quality observational data by only considering retrievals for which the root mean square value of the residuals (difference between measured and simulated radiances) relative to the maximum value of the radiances is smaller than 0.0065.
The cloud, retrieval quality and sensitivity lter for 4.9 km altitude leaves us with about 120 000 and 110 000 valid data points for each single morning and evening overpass in August, respectively. In February there are typically 100 000 and 95 000 valid morning and evening observations for each day, respectively. Each of these data points represents the middle tropospheric situation of a small area (12 km diameter at nadir).
6.3 Regional-scale signals
In order to demonstrate the high potential of IASI for a daily detection of regional-scale moisture transport pathways, we analyse the {H2O, D} pair distribution measured on 16 August 2014 in different regions during the morning overpass.We investigate observations over land (two distinct locations: Alaska and South Africa) and sea (two distinct locations: subtropical North Atlantic and Gulf of Persia). The analysed regions are marked by a bluish colour in Fig. 9 and the respective {H2O, D} pair density distributions are plotted in the upper left panel of Fig. 10.
IASI detects very distinct {H2O, D} pair distributions for the different regions. The air over Alaska and South Africa is similarly dry; however, the D values differ systematically by more than 150 . A look on the row kernels (right panels) reveals that in both regions, IASI has a very similar sensitivity, which peaks between 2 and 5 km altitude. An apparent explication is that in Alaska, the drying happens by condensation (via Rayleigh distillation), while in South Africa, the drying is due to mixing with very dry air (subsidence from the upper troposphere).
The air over the North Atlantic and the Gulf of Persia is similarly humid, but there is a clear difference in D. A reason for the difference in D might be that over the North Atlantic, rain re-evaporation is important, whereas the air over the Gulf of Persia is strongly affected by dry convection processes (vertical mixing) over the Arabian peninsula. The row kernels for both scenes are not the same however, both show peak sensitivities for altitudes between 4 and 8 km, meaning that for both regions, IASI should be able to consistently capture variations in the {H2O, D} pair distributions that take place within a deep middle tropospheric layer.
6.4 Seasonal- and diurnal-scale signals
The bottom panels of Fig. 10 demonstrate IASIs potential for detecting seasonal- and diurnal-scale signals in the {H2O, D} pair distribution. This is done by analysing the
{H2O, D} pair distributions over the Sahara (22.5 to 32.5 N and 10 W to 30 E, region marked by a reddish colour in Fig. 9) for six consecutive winter and summer days and for morning and evening overpasses. Morning overpasses of the IASI instruments are at about 10:00 LT and evening over-passes at about 22:00 LT.
In the context of Fig. 5 we discussed the SAL events that can be observed during July and August over the Atlantic Ocean. Actually, the dry convection process that is responsible for the distinct {H2O, D} pair distribution under SAL conditions in the surroundings of Tenerife takes place over the Sahara. The strong heating of the Earths surface during the day in summer is the main driver of these processes, which should be manifested by a pronounced diurnal cycle over the Sahara in August. Indeed, for August we ob-serve such diurnal signals in the MUSICA {H2O, D} pair data. For the morning overpasses we observe a fractionation that is similar to the situation in winter ( D values between
300 and 200 ), whereas for the evening overpass the D
values veer away from the Rayleigh line and group around a line that simulates mixing between the planetary boundary layer air and middle free tropospheric air ( D values between
250 and 140 ). This evening distribution is very similar
to the distribution that IASI typically observes over the Atlantic under SAL conditions (red contours in the right panel of Fig. 5). For February, the surface heating is much weaker than in summer and dry convection processes are unlikely.As a result, we observe no signicant difference between the morning and evening {H2O, D} pair distribution.
Regarding the morning data, the moistening from winter to summer does not happen perfectly parallel to a Rayleigh line.This is probably because in summer the evaporation takes place over a warmer ocean than in winter and we have to consider different Rayleigh lines for winter and summer.
The main intention of Fig. 10 is to give examples of the large potential of IASI {H2O, D} products. A more profound scientic study of these {H2O, D} pair distribution patterns can only be done with model calculations and is out of the scope of this paper.
7 Defective interpretation of H2O and D remote sensing data
The previous sections have demonstrated the feasibility of the remote sensing of {H2O, D} pair distributions. In this context, it is important to recall that so far we have only presented examples of {H2O, D} pair distributions using a posteriori processed remote sensing products (thereby assuring the same sensitivities for H2O and D; see right panels in Figs. 8 and 10). In this section, we will briey discuss the difculty of correctly interpreting the distribution of the {H2O, D} pairs obtained from remote sensing data that are not a posteriori processed, which is the type of data that are typically provided by other data producers.
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2862 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
Figure 11. Same as Fig. 8, but instead of the {H2O, D} pair product (a posteriori corrected product), the plots are for the direct H2O and D retrieval output (Type 1, no a posteriori correction). Please note the different scale on the x axis of the row kernel plots compared to Fig. 8.
7.1 Discussion of example cases
Figure 11 shows exactly the same as Fig. 8, but for data that have not undergone the a posteriori processing. In MUSICA we call them the Type 1 products and they represent individual optimal estimations of H2O and D. For the Type 1 product, H2O and D are not sensitive for the same air mass and the retrieval response is much more sensitive to atmospheric H2O than to D variations. This is revealed by the typical row kernels as depicted on the right panels. There is clearly a higher sensitivity for H2O than for D (compare the row kernels of the matrix blocks A 11 and A 22).
The different sensitivities affect the slopes in the {H2O, D} distribution plots. Concretely, the NDACC/FTIR
Izaa Type 1 data for winter (blue contour lines in upper panel of Fig. 11) show a {H2O, D} pair distribution
where dry air can occasionally be weakly depleted (significant number of data points with H2O concentrations below 102.7 ppmv 500ppmv and D above 400 ). It is very
likely that these data points only appear there because of the
H2O sensitivity being higher than the D sensitivity. Actually, when accounting for the different H2O and D sensitivities, this observation is not made (see Fig. 8). For midsummer, the Type 1 data show a {H2O, D} pair distribution that is reasonably parallel to a Rayleigh line (red contour lines in upper left panel of Fig. 11), whereas the distribution when H2O and D have almost identical sensitivities is more along a mixing line (red contour lines in upper left panel of Fig. 8).
Concerning East Africa (lower panel of Fig. 11), the
Type 1 {H2O, D} pairs are generally distributed around a line with a slope being less steep than the slope of a Rayleigh line. However, the conclusion that mixing with dry air is the
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2863
plotted in the right panels, and the very different entries in the matrix blocks A 11 and A 22 can be clearly observed.
The cyan contour lines in the upper panel of Fig. 12 are reasonably close to the exemplary Rayleigh line. This means that the moisture above the Gulf of Persia might be defectively interpreted as being partly dried by Rayleigh distillation processes. This is in contrast to what is indicated by the more reliable {H2O, D} pair distribution pattern from
Fig. 10. There the cyan contour lines are close to a mixing line, revealing that mixing with dry air is actually much more important than dehydration by condensation. Vice versa, for Alaska a D value of 300 for H2O concentrations below
102.8 ppmv 630ppmv (see purple contour lines in the up
per panel of Fig. 12) indicates drying by mixing with dry air; however, when considering the different sensitivities in H2O
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
Figure 12. Same as Fig. 10, but instead of the {H2O, D} pair product (a posteriori corrected product), the plots are for the direct H2O and D retrieval output (Type 1, no a posteriori correction). Please note the different scale on the x axis of the row kernel plots compared to
Fig. 10.
dominating drying process might be wrong because when analysing the a posteriori processed {H2O, D} pair distribution, the slope is steeper and in parallel to a Rayleigh line (see Fig. 8).
The risk of such defective interpretations is larger the more pronounced the difference between the H2O and D sensitivities is. For the MetOp/IASI products the sensitivity difference is especially important. Figure 12 shows exactly the same as Fig. 10, but for data that have not undergone the a posteriori processing. We observe completely changed {H2O, D} pair distribution patterns. These changed patterns reveal the risk of a defective interpretation of the real atmospheric situation when using H2O and D data that have signicantly different averaging kernels. The row kernels are
remote sensing data with the model because the modelled {H2O, D} pair distributions can be strongly camouaged by averaging kernels properties.
The different issues with the averaging kernels are less important when using a posteriori processed data because the respective kernels are less complex and have less of an effect on the {H2O, D} pair distributions (e.g. the slopes of the {H2O, D} pair distributions are not signicantly affected by the kernels). Then the application of the averaging kernel to the model data is less critical and a rst-order comparison between modelled and measured {H2O, D} pair distributions is even possible without applying the averaging kernels to the model (if variations of deep layers are of interest). This possibility becomes evident from Fig. 5. There, the in situ data as depicted in the left panel correspond to point measurements made at 2390 and 3550 m a.s.l., whereas the NDACC/FTIR and MetOp/IASI data reect the middle troposphere according to their averaging kernels (A 11 and A 22 in Figs. 1 and 2).
The {H2O, D} pair distributions of the a posteriori processed NDACC/FTIR and MetOp/IASI data compare well with the in situ data that are not affected by averaging kernels, suggesting that a posteriori processed remote sensing data can be used as a direct reference for a rst-order validation of modelled {H2O, D} pair distributions (there is no need for applying the kernels to the model as long as the dominating {H2O, D} pair signal of a deep layer is of interest).
A further improved integration of model and remote sensing {H2O, D} pairs might be achieved by the development of a {H2O, D} pair remote sensing retrieval simulator (the principle idea of a retrieval simulator is presented in Field et al., 2012).
8 Summary and conclusions
The MUSICA ground- and space-based water vapour isotopologue remote sensing products (generated from NDACC/FTIR and MetOp/IASI spectra with the nal MUSICA retrieval version, v2015) are calibrated and their quality is empirically documented. Compared to previous versions, v2015 improves the consistency between the different locations (uniform a priori for all retrievals and improved spectral windows for the NDACC/FTIR retrievals) and the consistency with respect to in situ references (data are calibrated using the aircraft prole references measured between the surface and 7 km altitude during the summer 2013 MUSICA campaign). The remaining bias in the v2015 data is very likely smaller than 15 % and 25 for the NDACC/FTIR lower tropospheric H2O and D products, respectively. In the middle troposphere, the remaining bias for the NDACC/FTIR and MetOp/IASI products is estimated to be within 10 % (for H2O) and 15 (for D).
Tropospheric D values are most interesting for science if provided together with H2O, i.e. it is important to validate the {H2O, D} pair distributions that are obtained
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
2864 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
and D it seems that drying by mixing is rather unlikely. Instead, dehydration by condensation is suggested (see corresponding contour line in Fig. 10).
A further example for defective interpretations is the seasonal cycle over the Sahara. The {H2O, D} pair distributions as depicted in the bottom panel of Fig. 12 suggest strong differences between the summer and winter humidity levels, but only small differences in D. Whereas during summer the {H2O, D} pairs are reasonably close to the exemplary Rayleigh line, they are far away from this line in winter. When removing the inconsistencies between the H2O and D sensitivities, this behaviour is very different. Then, the seasonal variation happens much closer to the exemplary Rayleigh line (see bottom panel of Fig. 10).
Already the few examples discussed here clearly demonstrate that the interpretation of {H2O, D} pair remote sensing data has to be done with great care. It is important to have a look at the sensitivities of H2O and D (as well as at their cross-sensitivities). The transformation of the {H2O, D} proxy basis system helps to gain insight into the complex characteristics of these remote sensing data products (see the discussions in Sect. 2 and the references cited therein). In order to reduce the risk of defective data interpretation, we strongly recommend the a posteriori data processing, which provides {H2O, D} pair distributions that can be interpreted in a straightforward manner (see Figs. 8 and 10).
7.2 Combination of remote sensing data and models
A profound interpretation of the {H2O, D} pair distributions is only possible in combination with isotopologue models. In this context and if model data are available, it might be argued that the a posteriori processing is actually not needed because one can compare Type 1 data with model outputs that have been convolved with the Type 1 averaging kernels, thereby simulating the effect of the different H2O and D kernels. However, model data that have been convolved with Type 1 averaging kernels will be strongly different from the original model data and they will, to a large extent, reect averaging kernel properties instead of real atmospheric signals (e.g. the slopes of the {H2O, D} pair distributions will strongly depend on the differences between the Type 1 H2O and D kernels). It will be difcult to understand what {H2O, D} pair signals are introduced by the averaging kernels and what {H2O, D} pair signals are actually modelled.Furthermore, the differences between the Type 1 H2O and D kernels depend on the atmospheric and geophysical conditions (surface and atmospheric temperatures, atmospheric humidity concentrations, etc.; see, for instance, the example row kernels of the A 11 and A 22 matrix blocks in Fig. 12), and since modelled and measured atmospheric state often differ signicantly (high small-scale variability of tropospheric humidity), there is a high risk when using inadequate Type 1 H2O and D kernels when convolving the model data. Under such circumstances it will be very difcult to compare the
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2865
using the remote sensing techniques. MUSICAs surface-based in situ measurements made on Tenerife at 2390 and 3550 m a.s.l. (subtropical North Atlantic) provide a continuous free tropospheric in situ reference record of {H2O, D}
pairs for validating the remote sensing data. We nd that the in situ and the calibrated and a posteriori processed remote sensing products reveal similar {H2O, D} pair distributions and consistently capture the three principle moisture pathways to the subtropical free troposphere: transport from the upper troposphere of the extratropics, transport from the lower troposphere over the subtropical/tropical Ocean and uplift via dry convection over the Sahara followed by advection over the Atlantic. We show that the space- and ground-based MUSICA v2015 data are rather consistent on global scale. First, there is no signicant bias between both data sets and second, the space- and ground-based products consistently detect extremes in the {H2O, D} distribution at different globally distributed locations. This suggests that the calibrations with respect to H2O and D reference scales and the validations with respect to the reference {H2O, D} pair distributions have global validity.
We present examples of seasonal cycles in the NDACC/FTIR {H2O, D} pair distribution and briey discuss possible links between the seasonality in this distribution and the seasonality in moisture sources and transport processes. Since NDACC/FTIR provides long-term data records, there are good possibilities for studying the {H2O, D} pair distributions in a climatological context. In the framework of follow-up studies, it would be interesting to examine whether there is a signicant long-term change in the seasonality of the {H2O, D} pair distributions at any
MUSICA NDACC/FTIR site.
MetOp/IASI offers high horizontal resolution, on a quasi-global scale, and morning as well as evening observations.We present and briey discuss examples of {H2O, D} pair distribution patterns on regional scales as well as on seasonal and diurnal timescales. The regional {H2O, D} patterns give insight into the horizontal distribution of humidity control mechanisms; for instance, they suggest regions where humidity is determined mainly by mixing and regions where drying by condensation or moistening by rain re-evaporation is dominant. The diurnal timescale patterns allow conclusions about the mechanisms that drive the diurnal cycle of atmospheric moisture; for instance in the summertime, Sahara dry convection seems to be very important. As soon as more MUSICA MetOp/IASI {H2O, D} data are produced, seasonal cycle analyses on global scales could be performed, and specic horizontal {H2O, D} patterns for different atmospheric and climate modes (NAO, ENSO, etc.) could be investigated.
The MUSICA {H2O, D} pair remote sensing data discussed here are produced from H2O and D optimal estimation retrievals by an a posteriori correction method. The a posteriori correction is needed for generating {H2O, D}
pairs, i.e. H2O and D products with the same sensitivities.
MUSICA data are available as a posteriori processed data, and we strongly recommend their usage whenever {H2O, D}
pair distributions are of interest. At the same time, we advise against using the H2O and D optimal estimation products (original retrieval output, not a posteriori processed) for {H2O, D} pair distribution studies. The reason is the different sensitivities of H2O and D, which imply a signicant risk for defective interpretations of {H2O, D} pair distributions.
This paper shows that reliable and carefully characterized {H2O, D} pair remote sensing observations can be made available by retrievals of spectra measured by NDACC/FTIR instruments and by MetOp/IASI sensors. The data have long-term characteristics and offer global coverage and high resolution (in space and time), thereby opening up new opportunities for addressing the focus research areas that are briey described in the Introduction section. In the next step, concrete research opportunities should be identied and evaluated. For this purpose, atmospheric models (like COSMO-iso, Pfahl et al., 2012) will be essential. Sensitivity tests with the models are needed in order to reveal the links between different moisture processes and distinct {H2O, D} pair distributions and in order to estimate whether the characteristics of the remote sensing data (e.g. their limited sensitivity and vertical resolution) allow the expected {H2O, D} pair signals to be detected.
9 Data availability
The MUSICA NDACC/FTIR data are available via the database of NDACC (ftp://ftp.cpc.ncep.noaa.gov/ndacc/MUSICA/
Web End =ftp://ftp.cpc.ncep.noaa.gov/ndacc/ ftp://ftp.cpc.ncep.noaa.gov/ndacc/MUSICA/
Web End =MUSICA/ ) and via ZENODO (Barthlott et al., 2016b). Details are given in Barthlott et al. (2016a).
In the long run, we also plan to disseminate the MUSICA MetOp/IASI data via a database and in a standard data format. At this stage the data are available as ascii data les and can be requested from the MUSICA team (by email to the leading author of this paper). Correct data usage will be assured by direct contact between the data users and the MUSICA team.
The MUSICA aircraft-based and ground-based in situ water vapour isotopologue observations used as the references in this study are discussed in detail in Dyroff et al. (2015) and Gonzlez et al. (2016), respectively. The data can be requested from the MUSICA team (by email to the leading author of this paper). It is foreseen to provide the data via a dedicated international database at Laboratoire des Science du Climat et de lEnvironnement, which is currently in development (see status at https://waterisotopes.lsce.ipsl.fr/
Web End =https://waterisotopes.lsce.ipsl.fr/ ).
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2866 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
Appendix A: Milestones of the MUSICA project
MUSICA ends in July 2016. During the last few years, methods for a theoretical characterization and empirical validation of {H2O, D} pair remote sensing have been developed.
The most important milestones and the corresponding publications are collected in Table A1.
Table A1. Important developments/milestones in the context of the MUSICA activities.
Development/milestone References
Optimal estimation of H2O, HDO and D Schneider et al. (2006)
Improving spectroscopic line parameterization using atmospheric spectra Schneider and Hase (2009); Schneider et al. (2011)
The MUSICA FTIR/NDACC retrieval, H2O and D proxies Schneider et al. (2012); Barthlott et al. (2016a) A posteriori processing for generating optimal {H2O, D} pairs Schneider et al. (2012)
The MUSICA MetOp/IASI retrieval Schneider and Hase (2011); Wiegele et al. (2014) Validation of {H2O, D} pairs Wiegele et al. (2014); Schneider et al. (2015)
Using XCO2 for quality ltering of MUSICA NDACC/FTIR Barthlott et al. (2015)
In situ prole references (ISOWAT aircraft campaign, 07 km) Dyroff et al. (2015); Schneider et al. (2015)
Continuous in situ reference for the free troposphere Gonzlez et al. (2016)
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2867
Appendix B: MUSICA in the context of other isotopologue ratio remote sensing data sets
We would like to remark that the results as shown in this paper are only valid for the MUSICA products. This appendix gives a brief overview on other (non-MUSICA) tropospheric water vapour isotopologue remote sensing products and briey discusses their differences to the MUSICA products.
B1 The TCCON XH2O and XHDO data
The ground-based FTIR water vapour isotopologue products that are made available via TCCON (http://www.tccon.caltech.edu/
Web End =www.tccon.caltech. http://www.tccon.caltech.edu/
Web End =edu/ ) are fundamentally different from the MUSICA ground-based FTIR isotopologue products.
A TCCON-like product is discussed in Rokotyan et al. (2014). It relies on near-infrared absorption lines (where HD16O is a rather weak absorber) and the isotopologue ratios are calculated a posteriori from independently retrieved H162O and HD16O column amounts. Such a posteriori calculated ratios are affected by the different sensitivities of the individual H162O and HD16O retrievals. At the moment, the
TCCON kernels do not give information about the cross-correlations between the H162O and HD16O product, and it is not possible to calculate kernels for humidity and D proxies. This means that no Type 2 product can be calculated. The TCCON retrievals use NCEP (National Centers for Environmental Prediction) humidity analyses as H2O a priori and construct the D a priori proles by assuming a xed linear correlation between ln[H2O] and D ( D =
0.0695 [notdef] ln[H2O] + 0.28).
B2 Satellite-based tropospheric water vapour isotopologue data
A brief overview of available products of tropospheric water vapour isotopologues and the respective satellite sensors is given in the Tables B1 and B2.
The thermal nadir sensors TES and IASI have the best sensitivity with respect to the water vapour isotopologues in the middle troposphere (about 28 km altitude). In addition to the MUSICA research team, a group at the University of Brussels (ULB) is working on IASI water vapour isotopologue retrievals. The ULB IASI retrieval uses two small spectral microwindows (11931223 and 12511253 cm1) and ts the proxies for humidity, D and CH4 below 10 km altitude as well as ground temperature (Lacour et al., 2012).It uses the EUMETSAT Level 2 temperature output for the whole atmosphere and the EUMETSAT Level 2 humidity output for altitudes above 10 km (no t). The ULB group uses the same globally uniform a priori data as the MUSICA group. For the ULB retrieval, Type 2 products (a posteriori processed, H2O and D have almost identical sensitivities)
and Type 1 products (individual optimal estimation of H2O and D) can be made available.
Another tropospheric isotopologue product is generated from AURA/TES spectra. It has rst been presented by Worden et al. (2006), whereby small spectral microwindows have been tted (similar to the ULB IASI retrieval). The TES version 5 product used nowadays is discussed in Worden et al. (2012), and the respective retrieval set-up is rather similar to the MUSICA IASI retrieval set-up: a broad spectral window, simultaneous t of proxies for humidity and D as well as of temperature and the interfering gases throughout the atmosphere. However, for the TES retrieval the H2O a priori assumption comes from the NCEP humidity analyses, and the D a priori has a latitudinal dependency. The generation of a Type 2 product is theoretically possible, but it is currently not provided. TES measures limb and thermal nadir spectra (the isotopologue data are generated from the nadir spectra). It has a similar spectral coverage to IASI, but a signicantly higher spectral resolution and, on the other hand, much sparser daily horizontal coverage. Lacour et al. (2015) showed a cross-validation of the TES version 5 products and the IASI products generated by the ULB group.
Space-based sensors measuring solar short-wave infrared spectra (SWIR) reected on the Earths surface theoretically have better sensitivity in the lower troposphere than the thermal nadir sensors. Retrievals using the sensors SCIAMACHY and GOSAT have been presented and assessed by Frankenberg et al. (2009), Frankenberg et al. (2013), Boesch et al. (2013) and Scheepmaker et al. (2015). All use humidity analyses (NCEP or ECMWF) as humidity a priori, but a globally uniform D a priori. The respective retrievals work independently for the different isotopologues, and the isotopologue ratio is calculated after the retrieval process. This is an important difference to the thermal nadir retrievals, which optimally estimate the proxies of H2O and HDO / H2O. The near-infrared retrievals are thus affected by the different sensitivities for the different isotopologues. A further difference is that in the near-infrared, the absorption signatures of the secondary isotopologue (HD16O) is signicantly smaller than in the thermal infrared. The daily horizontal coverage of these sensors is much sparser than for IASI. For the current SWIR retrievals, no humidity and D proxy kernels are available, and it is not possible to assess the difference between the H2O and D kernels and to correct for it; i.e. it is not possible to perform an a posteriori correction and generate a Type 2 product.
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2868 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
Table B1. Overview of tropospheric water vapour isotopologue retrievals using space-based observations.
Research group/sensor Spectral window Fitted parameter References and remarks
MUSICA (IMK-ASF)/IASI 11901400 cm1 (ln[H2O] + ln[HDO])/2 Schneider and Hase (2011), Wiegele et al. (2014),
ln[HDO] ln[H2O] atm. temp. constrained to EUMETSAT L2,
CH4, N2O, CO2 and HNO3 only clear sky retrievals atmospheric temperatureground temperature
ULB (U. Brussels)/IASI 11931223 cm1 (ln[H2O] + ln[HDO])/2 Lacour et al. (2012), Lacour et al. (2015)
12511253 cm1 ln[HDO] ln[H2O] only t for 010 km,
CH4 atm. temp. from EUMETSAT L2, ground temperature only clear sky retrievals
NASA/TES 11701320 cm1 (ln[H2O] + ln[HDO])/2 Worden et al. (2006), Worden et al. (2012)
ln[HDO] ln[H2O]
CH4 and N2O atmospheric temperature ground temperature cloud ( and pressure)
SRON/SCIAMACHY 42124248 cm1 H16
2 O, H
18
2 O and HD16O Frankenberg et al. (2009),
CH4 and CO Scheepmaker et al. (2015)
NASA/GOSAT 63116441 cm1 H2O and HDO Frankenberg et al. (2013)
U. Leicester/GOSAT 64396464 cm1 H2O and HDO Boesch et al. (2013),uses CH4 obtained from extra CH4 retrieval
Table B2. Space-based sensors with available tropospheric water vapour isotopologue retrieval products.
Sensor Meas. geo. Pixel size Meas. per day Temporal coverage Spectral res.
IASI thermal nadir 12 km diameter 1.3 million IASI-A: since 2007 0.5 cm
(at nadir) IASI-B: since 2013
IASI-C: scheduled for 2018
TES thermal nadir 5 [notdef] 8 km 2100 since 2002 0.1 cm
(at nadir) (since 2010 temporarily)
SCIAMACHY SWIR 120 [notdef]30 km 32 000 20032012 0.45 cm
(after 2007 increased detector degradation)
GOSAT SWIR 10 km diameter 10 000 since 2009 0.4 cm
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
1
1
1
1
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2869
Appendix C: Reference proles
C1 Coincidences in space and time
During July 2013 we performed an aircraft campaign in the surroundings of Tenerife. We operated the ISOWAT instrument (Dyroff et al., 2010) aboard the aircraft and measured highly resolved vertical proles of H2O and D from sea surface up to almost 7 km altitude on 6 individual days (21 Jul 2013, 22 Jul 2013, 24 Jul 2013, 25 Jul 2013, 30 Jul 2013 and31 Jul 2013; Dyroff et al., 2015; Schneider et al., 2015). The aircrafts ascent and descent took place between 10:30 and 13:30 UT. Figure C1 shows a site map indicating the horizontal ight track of the aircraft (grey line) as well as the location of the Tenerife FTIR instrument (green star) and the IASI observation pixels (coloured squares and diamonds).
A main characteristic of tropospheric humidity is the high short-term and small-scale variability, which make inter-comparison studies difcult. Concerning temporal variability, Steinke et al. (2015) reported good correlation of integrated water vapour (IWV) observations and simulations made within 3 h (see their Fig. 7). Vogelmann et al. (2015) use DIAL (differential absorption lidar) for estimating temporal mismatches in vertical proles. Their Fig. 6 reports middle tropospheric (46 km) H2O variabilities within 2 h of about 20 % and within 1 h of about 15 %. This variability increases strongly in the upper troposphere where they found about 60 % within 2 h. When interpreting that study, we have to take into account that the DIAL detects vertically high-resolution proles, whereas the FTIR and IASI remote sensing data represent deep layers (where variability largely cancels out). Spatial variability can be estimated by space-based observations or by models. Steinke et al. (2015) used ICON simulations and estimated the scatter for humidity encountered at a distance of 10 km to be about 4 % (right column in their Fig. 4 using an IWV of 12 kg m2 from their
Fig. 2). In Wiegele et al. (2014), we estimated a variability of about 19 % and 17 for middle tropospheric H2O and D measured within an 110 km [notdef] 110 km area around Tenerife,
whereby at least half of this variability is due to the random errors in the IASI data.
The temporal and spatial variability of humidity elds in the surroundings of Tenerife has been estimated in different studies (Schneider et al., 2010a, b; Wiegele et al., 2014;Gonzlez et al., 2016) using NDACC/FTIR, MetOp/IASI, ground-based in situ, radiosondes, GPS, radiometers and sun photometer data. Observations made for the middle tropo-sphere within 2 h have a scatter of typically 10 % (for H2O)
and 10 (for D) that cannot be explained by uncertainties of the instruments. However, close to the surface at the Izaa Observatory, the variability is larger and shows a strong diurnal cycle (Schneider et al., 2010b; Gonzlez et al., 2016), which is caused by local circulation (thermal upslope ow reaches the islands mountains in the late morning hours).These local circulations have to be considered when compar-
Figure C1. Site map indicating the location of the different instruments and ground pixels during the aircraft campaign on 6 days in July 2013. Green star: Izaa Observatory (location of the rst Picarro and the FTIR; the green dashed lines indicate the line of sight of the FTIR between 08:15 and 13:30 UT); blue star: Teide Observatory (location of the second Picarro); grey lines: aircraft ight track during ISOWAT measurements; black squares and diamonds: cloud-free ground pixels of IASI-A and IASI-B, respectively, during the six aircraft ights; red lled squares and diamonds: pixels that full our coincidence criteria for IASI-A and IASI-B, respectively, whereby the different lling colour corresponds to the 6 different days as in Figs. 3 and 4.
ing measurements made in Izaa with measurements made at the free tropospheric location of the aircraft. Since the H2O and D signals as measured in Izaa in the later morning are already strongly affected by the upslope airow that is developing during the late morning, early morning data are better representative of the free troposphere. For this reason, we compare FTIR data measured between 08:15 UT and 09:45 UT (for July this means solar elevation between 25 and 45 , which is about 2 and 3.5 h after sunrise) with the free tropospheric aircraft measurements made between 10:30 and 13:30 UT (see Fig. 3 in Sect. 3). This means that there is a temporal mismatch between the FTIR and the ISOWAT observations between 45 min and 5 h for the measurement pairs we dene as optimal coincidences. This increases the mismatch uncertainty by about a factor of 3, i.e. from 10 % and 10 to 30 % and 30 (this estimatiuon is based on Fig. 12 of Schneider et al., 2010b). Nevertheless, this is what we dene as optimal coincidences.
Figure C2 shows the same as Fig. 3, but for FTIR data that have been measured during the time of the aircraft prole measurements, i.e. in the late morning hours and during midday. The comparison plots show relatively large variability in FTIR data that represent the atmospheric layer just above the island, revealing the large impact of the local diurnal upslope ow on the FTIR observations. Although for this comparison the temporal mismatches are rather low, it does not represent
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2870 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
Figure C2. Same as Fig. 3, but for FTIR measurements corresponding to best temporal coincidences (FTIR observations made during the 3 h of the aircraft ights, typically 10:3013:30 UT).
optimal coincidences because the strong local upslope ow on the island means that the FTIR and the ISOWAT instrument detect different air masses.
The coloured squares and diamonds in Fig. C1 indicate the locations of the IASI observation pixels, whereby the different colours correspond to the different days (see legend in Fig. 4). We require, as coincidence criteria, that the observation pixel is not further away than 50 km from the aircrafts track (coloured squares and diamonds group around the track, which is indicated as the grey line). On the 2 days21 Jul 2013 and 22 Jul 2013 there are no IASI pixels within 50 km of the aircrafts track and for those days we also include observation pixels that are located more than 100 km away from the ight track. These are the grey and red pixels marked by a black edge, indicating that they correspond to non-optimal coincidences. In addition, there are three magenta-coloured pixels (representing day 31 Jul 2013) that are marked by black edges. These are also representative of non-optimal coincidences, since on day 31 July 2013 there
was a very sharp gradient from the southeast of the ight track (air mass with strong SAL conditions) to the northwest of the ight track (air mass with weaker SAL conditions), and the aircrafts ISOWAT measurements and IASI detect air masses of different characteristics (for a more detailed discussion of this day, 31 Jul 2013, please see Appendix A of Schneider et al., 2015).
C2 Ceiling altitude and uncertainties
The reference proles are constructed from the ISOWAT measurements (surface up to almost 7 km) and climatological values assumed above the ceiling altitude (see Sect. 2.5). This prole (ISOWAT + climatology) is then smoothed by
the averaging kernel of the remote sensor. The uncertainty of ISOWAT is 4 % for H2O and typically better than 10 for D (it can only reach high values for rather dry air, e.g. 35 when the H2O concentration is below 500 ppm,
Dyroff et al., 2015). However, above the ceiling altitude we have no measurements and the uncertainty is signicantly
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2871
larger. We assume 100 % and 80 for H2O and D, respectively. These large uncertainties for the atmosphere above the ceiling altitude propagate to lower altitudes due to the smoothing with the averaging kernels (see typical row kernels in Figs. 1 and 2). In fact, the error bars on the smoothed reference data for 4.9 km as well as for 2.4 km as depicted in Figs. 3, 4 and C2 are dominated by unavailable reference H2O and D values above 67 km (the D error bars are only dominated by the uncertainty in the D measurements up to 7 km for very dry conditions).
For the references used for the FTIR data validation at2.4 km altitude (top panels of Figs. 3 and C2), the uncertainty introduced from missing data above the ceiling altitude can reach 15 % for H2O and 12 for D. For the FTIR data validation at 4.9 km altitude (bottom panels of Figs. 3 and C2), the respective uncertainty can reach 25 % for H2O and 20 for D. For the IASI data validation at 4.9 km, the respective uncertainties are 256 % for H2O and 205 for D. The values depend on both the averaging kernels as well as the ceiling altitude (e.g. on 21 Jul 2013 we only reached6.0 km, leading to higher uncertainties than for days when we reached 6.8 km). Please consider that the kernels in Figs. 1 and 2 are on a logarithmic scale; i.e. a value of 0.1 of the4.9 km row kernel at 6.5 km altitude together with a prole uncertainty of 80 (no measurement at 6.5 km) means an uncertainty in the smoothed prole value at 4.9 km of 8 .The error sums up to 15 by adding up the contributions from all altitudes above 6.5 km.
For the validation/calibration of the remote sensing data, it is essential to have reference measurements that cover the troposphere from the surface up to high altitudes. During the ISOWAT campaign we almost reached 7 km during most of the ights and we are not aware of another D reference prole data set with similarly good altitude coverage.
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2872 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
Author contributions. A. Wiegele performed the IASI retrievals and S. Barthlott the FTIR retrievals. Y. Gonzlez, E. Christner andC. Dyroff produced the in situ reference data and helped with their interpretation. F. Hase developed the PROFFIT and PROFFIT-nadir retrieval codes. O. E. Garca, E. Seplveda, S. Barthlott, F. Hase,T. Blumenstock, G. Mengitsu Tsidu and S. Takele Kenea performed and maintained the NDACC/FTIR measurements at Izaa and Addis Ababa. S. Rodrguez and Y. Gonzlez helped in interpreting the aerosol data and the SAL events. J. Andrey coordinated the data exchange between KIT and INTA in the context of the July 2013 aircraft campaign. M. Schneider coordinated and designed the MUSICA project and prepared the manuscript with contributions from all co-authors.
Acknowledgements. This study has been conducted in the framework of the project MUSICA which is funded by the European Research Council under the European Communitys Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement number 256961.E. Seplveda is supported by the Ministerio de Economa and Competitividad of Spain for the project NOVIA (CGL2012-37505).
The aircraft campaign has been co-funded by the project MUSICA and the Spanish national project AMISOC (CGL2011-24891).
We are grateful to INTA Aerial Platforms, a branch of the Spanish ICTS program, and the Spanish Air Force for their efforts in maintaining and operating the aircraft.
The AERONET sun photometer at Izaa (principal investigator: Emilio Cuevas) has been calibrated within AERONET EUROPE TNA supported by the European Community Research Infrastructure Action under the FP7 Capacities program for Integrating Activities, ACTRIS grant agreement number 262254.
The Izaa aerosol in situ measurements are part of the project POLLINDUST (CGL2011-26259) funded by the Minister of Economy and Competitiveness of Spain.
We thank all the personnel from the Izaa Atmospheric Research Center (IARC) of the Agencia Estatal de Meteorologa (AEMET).Our study has strongly benetted from this great support and important measurements have been made in IARC research facilities.
We acknowledge the support by the Deutsche Forschungsgemeinschaft and the Open Access Publishing Fund of the Karlsruhe Institute of Technology.
The article processing charges for this open-access publication were covered by a ResearchCentre of the Helmholtz Association.
Edited by: D. Feist
The AERONET sun photometer at Izaa (PI: Dr. Emilio Cuevas) has been calibrated within AERONET EUROPE European Community Research Infrastructure Action under the FP7 Capacities program for Integrating Activities, number 262254.
The Izaa aerosol in-situ measurements are part of the project POLLINDUST (CGL2011-26259) funded by the Competitiveness of Spain.
References
Barthlott, S., Schneider, M., Hase, F., Wiegele, A., Christner,E., Gonzlez, Y., Blumenstock, T., Dohe, S., Garca, O. E., Seplveda, E., Strong, K., Mendonca, J., Weaver, D., Palm, M., Deutscher, N. M., Warneke, T., Notholt, J., Lejeune, B., Mahieu,E., Jones, N., Grifth, D. W. T., Velazco, V. A., Smale, D., Robinson, J., Kivi, R., Heikkinen, P., and Raffalski, U.: Using XCO2 retrievals for assessing the long-term consistency of
NDACC/FTIR data sets, Atmos. Meas. Tech., 8, 15551573, doi:http://dx.doi.org/10.5194/amt-8-1555-2015
Web End =10.5194/amt-8-1555-2015 http://dx.doi.org/10.5194/amt-8-1555-2015
Web End = , 2015.
Barthlott, S., Schneider, M., Hase, F., Blumenstock, T., Kiel, M., Dubravica, D., Garca, O. E., Seplveda, E., Mengistu Tsidu, G., Takele Kenea, S., Grutter, M., Plaza, E. F., Stremme, W., Strong,K., Weaver, D., Palm, M., Warneke, T., Notholt, J., Mahieu, E., Servais, C., Jones, N., Grifth, D. W. T., Smale, D., and Robinson, J.: Tropospheric water vapour isotopologue data (H16
2 O,
H18
2 O and HD16O) as obtained from NDACC/FTIR solar absorption spectra, Earth Syst. Sci. Data Discuss., doi:http://dx.doi.org/10.5194/essd-2016-9
Web End =10.5194/essd- http://dx.doi.org/10.5194/essd-2016-9
Web End =2016-9 , in review, 2016a.
Barthlott, S., Schneider, M., Hase, F., Blumenstock, T., Mengistu Tsidu, G., Grutter de la Mora, M., Strong, K., Notholt,J., Mahieu, E., Jones, N., and Smale, D.: The ground-based MUSICA dataset: Tropospheric water vapour isotopologues (H16
2 O, H
182 O and HD16O) as obtained from NDACC/FTIR solar absorption spectra, doi:http://dx.doi.org/10.5281/zenodo.48902
Web End =10.5281/zenodo.48902 http://dx.doi.org/10.5281/zenodo.48902
Web End = , 2016b.
Berkelhammer, M., Risi, C., Kurita, N., and Noone, D. C.: The moisture source sequence for the Madden-Julian Oscillation as derived from satellite retrievals of HDO and H2O, J. Geophys.
Res.-Atmos., 117, D03106, doi:http://dx.doi.org/10.1029/2011JD016803
Web End =10.1029/2011JD016803 http://dx.doi.org/10.1029/2011JD016803
Web End = , 2012. Boesch, H., Deutscher, N. M., Warneke, T., Byckling, K., Cogan,A. J., Grifth, D. W. T., Notholt, J., Parker, R. J., and Wang, Z.: HDO/H2O ratio retrievals from GOSAT, Atmos. Meas. Tech., 6, 599612, doi:http://dx.doi.org/10.5194/amt-6-599-2013
Web End =10.5194/amt-6-599-2013 http://dx.doi.org/10.5194/amt-6-599-2013
Web End = , 2013.
Dansgaard, W.: Stable isotopes in precipitation, Tellus, 16, 436468, doi:http://dx.doi.org/10.1111/j.2153-3490.1964.tb00181.x
Web End =10.1111/j.2153-3490.1964.tb00181.x http://dx.doi.org/10.1111/j.2153-3490.1964.tb00181.x
Web End = , 1964.
Dyroff, C., Ftterer, D., and Zahn, A.: Compact diode-laser spectrometer ISOWAT for highly sensitive airborne measurements of water-isotope ratios, Appl. Phys. B, 98, 537548, doi:http://dx.doi.org/10.1007/s00340-009-3775-6
Web End =10.1007/s00340-009-3775-6 http://dx.doi.org/10.1007/s00340-009-3775-6
Web End = , 2010.
Dyroff, C., Sanati, S., Christner, E., Zahn, A., Balzer, M., Bouquet,H., McManus, J. B., Gonzlez-Ramos, Y., and Schneider, M.: Airborne in situ vertical proling of HDO / H16
2 O in the subtropical troposphere during the MUSICA remote sensing validation campaign, Atmos. Meas. Tech., 8, 20372049, doi:http://dx.doi.org/10.5194/amt-8-2037-2015
Web End =10.5194/amt- http://dx.doi.org/10.5194/amt-8-2037-2015
Web End =8-2037-2015 , 2015.
Field, R. D., Risi, C., Schmidt, G. A., Worden, J., Voulgarakis,A., LeGrande, A. N., Sobel, A. H., and Healy, R. J.: A Tropospheric Emission Spectrometer HDO/H2O retrieval simulator for climate models, Atmos. Chem. Phys., 12, 1048510504, doi:http://dx.doi.org/10.5194/acp-12-10485-2012
Web End =10.5194/acp-12-10485-2012 http://dx.doi.org/10.5194/acp-12-10485-2012
Web End = , 2012.
Field, R. D., Kim, D., LeGrande, A. N., Worden, J., Kelley,M., and Schmidt, G. A.: Evaluating climate model performance in the tropics with retrievals of water isotopic composition from Aura TES, Geophys. Res. Lett., 41, 60306036, doi:http://dx.doi.org/10.1002/2014GL060572
Web End =10.1002/2014GL060572 http://dx.doi.org/10.1002/2014GL060572
Web End = , 2014.
5
We thank all the personal from the Izaa Atmospheric Research Center (IARC) of the Agencia Estatal destudy has strongly benetted from this great support and important measurements have been made in IARC research
We acknowledge the support by the Deutsche Forschungsgemeinschaft and the Open Access Publishing Fund of of Technology.10
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2873
Frankenberg, C., Yoshimura, K., Warneke, T., Aben, I., Butz, A., Deutscher, N., Grifth, D., Hase, F., Notholt, J., Schneider, M., Schrejver, H., and Rckmann, T.: Dynamic processes governing lower-tropospheric HDO/H2O ratios as observed from space and ground, Science, 325, 13741377, doi:http://dx.doi.org/10.1126/science.1173791
Web End =10.1126/science.1173791 http://dx.doi.org/10.1126/science.1173791
Web End = , 2009.
Frankenberg, C., Wunch, D., Toon, G., Risi, C., Scheepmaker, R., Lee, J.-E., Wennberg, P., and Worden, J.: Water vapor isotopologue retrievals from high-resolution GOSAT shortwave infrared spectra, Atmos. Meas. Tech., 6, 263274, doi:http://dx.doi.org/10.5194/amt-6-263-2013
Web End =10.5194/amt-6- http://dx.doi.org/10.5194/amt-6-263-2013
Web End =263-2013 , 2013.
Galewsky, J., Sobel, A., and Held, I.: Diagnosis of Subtropical Humidity Dynamics Using Tracers of Last Saturation, J. Atmos.Sci., 62, 33533367, 2005.
Gat, J. R.: Atmospheric water balance the isotopic perspective, Hydrol. Process., 14, 13571369, doi:http://dx.doi.org/10.1002/1099-1085(20000615)14:8<1357::AID-HYP986>3.0.CO;2-7
Web End =10.1002/1099- http://dx.doi.org/10.1002/1099-1085(20000615)14:8<1357::AID-HYP986>3.0.CO;2-7
Web End =1085(20000615)14:8<1357::AID-HYP986>3.0.CO;2-7 , 2000.Gisi, M., Hase, F., Dohe, S., and Blumenstock, T.: Camtracker: a new camera controlled high precision solar tracker system for FTIR-spectrometers, Atmos. Meas. Tech., 4, 4754, doi:http://dx.doi.org/10.5194/amt-4-47-2011
Web End =10.5194/amt-4-47-2011 http://dx.doi.org/10.5194/amt-4-47-2011
Web End = , 2011.
Gonzlez, Y., Schneider, M., Dyroff, C., Rodrguez, S., Christner,E., Garca, O. E., Cuevas, E., Bustos, J. J., Ramos, R., Guirado-Fuentes, C., Barthlott, S., Wiegele, A., and Seplveda, E.: Detecting moisture transport pathways to the subtropical North Atlantic free troposphere using paired H2O- D in situ measurements, Atmos. Chem. Phys., 16, 42514269, doi:http://dx.doi.org/10.5194/acp-16-4251-2016
Web End =10.5194/acp- http://dx.doi.org/10.5194/acp-16-4251-2016
Web End =16-4251-2016 , 2016.
Hase, F., Hannigan, J. W., Coffey, M. T., Goldman, A., Hpfner, M., Jones, N. B., Rinsland, C. P., and Wood, S.: Intercomparison of retrieval codes used for the analysis of high-resolution, J. Quant.Spectrosc. Ra., 87, 2552, 2004.
Herman, R. L., Cherry, J. E., Young, J., Welker, J. M., Noone, D., Kulawik, S. S., and Worden, J.: Aircraft validation of Aura Tropospheric Emission Spectrometer retrievals of HDO / H2O, Atmos. Meas. Tech., 7, 31273138, doi:http://dx.doi.org/10.5194/amt-7-3127-2014
Web End =10.5194/amt-7-3127-2014 http://dx.doi.org/10.5194/amt-7-3127-2014
Web End = , 2014.
International in situ water isotope database: available at: https://waterisotopes.lsce.ipsl.fr
Web End =https: https://waterisotopes.lsce.ipsl.fr
Web End =//waterisotopes.lsce.ipsl.fr , 2013.
Kuang, Z., Toon, G., Wennberg, P., and Yung, Y.: Measured HDO/H2O ratios across the tropical tropopause, Geophys. Res.
Lett., 30, 251254, doi:http://dx.doi.org/10.1029/2003GL017023
Web End =10.1029/2003GL017023 http://dx.doi.org/10.1029/2003GL017023
Web End = , 2003.
Lacour, J.-L., Risi, C., Clarisse, L., Bony, S., Hurtmans, D., Clerbaux, C., and Coheur, P.-F.: Mid-tropospheric D observations from IASI/MetOp at high spatial and temporal resolution, Atmos. Chem. Phys., 12, 1081710832, doi:http://dx.doi.org/10.5194/acp-12-10817-2012
Web End =10.5194/acp- http://dx.doi.org/10.5194/acp-12-10817-2012
Web End =12-10817-2012 , 2012.
Lacour, J.-L., Clarisse, L., Worden, J., Schneider, M., Barthlott, S., Hase, F., Risi, C., Clerbaux, C., Hurtmans, D., and Coheur, P.-F.: Cross-validation of IASI/MetOp derived tropospheric D with TES and ground-based FTIR observations, Atmos. Meas. Tech., 8, 14471466, doi:http://dx.doi.org/10.5194/amt-8-1447-2015
Web End =10.5194/amt-8-1447-2015 http://dx.doi.org/10.5194/amt-8-1447-2015
Web End = , 2015.
Marsham, J. H., Dixon, N. S., Garcia-Carreras, L., Lister, G. M. S.,
Parker, D. J., Knippertz, P., and Birch, C. E.: The role of moist convection in the West African monsoon system: Insights from continental-scale convection-permitting simulations, Geophys.Res. Lett., 40, 18431849, doi:http://dx.doi.org/10.1002/grl.50347
Web End =10.1002/grl.50347 http://dx.doi.org/10.1002/grl.50347
Web End = , 2013.Noone, D.: Pairing Measurements of the Water Vapor Isotope Ratio with Humidity to Deduce Atmospheric Moistening and Dehydra-
tion in the Tropical Midtroposphere, J. Climate, 25, 44764494, doi:http://dx.doi.org/10.1175/JCLI-D-11-00582.1
Web End =10.1175/JCLI-D-11-00582.1 http://dx.doi.org/10.1175/JCLI-D-11-00582.1
Web End = , 2012.
NDACC database: available at: ftp://ftp.cpc.ncep.noaa.gov/ndacc/MUSICA/
Web End =ftp://ftp.cpc.ncep.noaa.gov/ndacc/ ftp://ftp.cpc.ncep.noaa.gov/ndacc/MUSICA/
Web End =MUSICA/ , 2016.
Ortega, P., Lehner, F., Swingedouw, D., Masson-Delmotte, V., Raible, C. C., Casado, M., and Yiou, P.: A model-tested North Atlantic Oscillation reconstruction for the last millennium, Nature, 523, 7174, doi:http://dx.doi.org/10.1038/nature14518
Web End =10.1038/nature14518 http://dx.doi.org/10.1038/nature14518
Web End = , 2015.
Pfahl, S., Wernli, H., and Yoshimura, K.: The isotopic composition of precipitation from a winter storm a case study with the limited-area model COSMOiso, Atmos. Chem. Phys., 12, 1629 1648, doi:http://dx.doi.org/10.5194/acp-12-1629-2012
Web End =10.5194/acp-12-1629-2012 http://dx.doi.org/10.5194/acp-12-1629-2012
Web End = , 2012.
Risi, C., Noone, D., Worden, J., Frankenberg, C., Stiller, G., Kiefer,M., Funke, B., Walker, K., Bernath, P., Schneider, M., Bony, S., Lee, J., Brown, D., and Sturm, C.: Process-evaluation of tropospheric humidity simulated by general circulation models using water vapor isotopic observations. Part 2: an isotopic diagnostic to understand the mid and upper tropospheric moist bias in the tropics and subtropics, J. Geophys. Res., 117, D05304, doi:http://dx.doi.org/10.1029/2011JD016623
Web End =10.1029/2011JD016623 http://dx.doi.org/10.1029/2011JD016623
Web End = , 2012a.
Risi, C., Noone, D., Worden, J., Frankenberg, C., Stiller, G., Kiefer,M., Funke, B., Walker, K., Bernath, P., Schneider, M., Wunch, D., Sherlock, V., Deutscher, N., Grifth, D., Wennberg, P., Strong,K., Barthlott, S., Hase, F., G. O., Smale, D., Mahieu, E., Sayres,D., Bony, S., Lee, J., Brown, D., Uemura, R., and Sturm, C.: Process-evaluation of tropospheric humidity simulated by general circulation models using water vapor isotopic observations. Part 1: comparison between models and datasets, J. Geophys. Res., 117, D05303, doi:http://dx.doi.org/10.1029/2011JD016621
Web End =10.1029/2011JD016621 http://dx.doi.org/10.1029/2011JD016621
Web End = , 2012b. Rodrguez, S., Alastuey, A., Alonso-Prez, S., Querol, X., Cuevas,E., Abreu-Afonso, J., Viana, M., Prez, N., Pandol, M., and de la Rosa, J.: Transport of desert dust mixed with North African industrial pollutants in the subtropical Saharan Air Layer, Atmos. Chem. Phys., 11, 66636685, doi:http://dx.doi.org/10.5194/acp-11-6663-2011
Web End =10.5194/acp-11-6663- http://dx.doi.org/10.5194/acp-11-6663-2011
Web End =2011 , 2011.
Rokotyan, N. V., Zakharov, V. I., Gribanov, K. G., Schneider, M., Bron, F.-M., Jouzel, J., Imasu, R., Werner, M., Butzin, M., Petri, C., Warneke, T., and Notholt, J.: A posteriori calculation of 18O and D in atmospheric water vapour from ground-based near-infrared FTIR retrievals of H16
2 O, H
182 O, and HD16O, Atmos. Meas. Tech., 7, 25672580, doi:http://dx.doi.org/10.5194/amt-7-2567-2014
Web End =10.5194/amt-7-2567-2014 http://dx.doi.org/10.5194/amt-7-2567-2014
Web End = , 2014.
Rothman, L. S., Gordon, I. E., Barbe, A., Chris Benner, D.,
Bernath, P. F., Birk, M., Boudon, V., Brown, L. R., Campargue, A., Champion, J.-P., Chance, K., Coudert, L. H., Dana,V., Devi, V. M., Fally, S., Flaud, J.-M., Gamache, R. R., Goldman, A., Jacquemart, D., Kleiner, I., Lacome, N., Lafferty, W. J., Mandin, J.-Y., Massie, S. T., Mikhailenko, S. N., Miller, C. E., Moazzen-Ahmadi, N., Naumenko, O. V., Nikitin, A. V., Orphal, J., Perevalov, V. I., Perrin, A., Predoi-Cross, A., Rins-land, C. P., Rotger, M., Simeckov, M., Smith, M. A. H., Sung,K., Tashkun, S. A., Tennyson, J., Toth, R. A., Vandaele, A. C., and Vander-Auwera, J.: The HITRAN 2008 molecular spectroscopic database, J. Quant. Spectrosc. Ra., 110, 533572, doi:http://dx.doi.org/10.1016/j.jqsrt.2009.02.013
Web End =10.1016/j.jqsrt.2009.02.013 http://dx.doi.org/10.1016/j.jqsrt.2009.02.013
Web End = , 2009.
Rothman, L., Gordon, I., Babikov, Y., Barbe, A., Benner, D. C.,
Bernath, P., Birk, M., Bizzocchi, L., Boudon, V., Brown, L.,
Campargue, A., Chance, K., Cohen, E., Coudert, L., Devi, V.,
Drouin, B., Fayt, A., Flaud, J.-M., Gamache, R., Harrison, J.,
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
2874 M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs
Hartmann, J.-M., Hill, C., Hodges, J., Jacquemart, D., Jolly,A., Lamouroux, J., Roy, R. L., Li, G., Long, D., Lyulin, O., Mackie, C., Massie, S., Mikhailenko, S., Mller, H., Naumenko,O., Nikitin, A., Orphal, J., Perevalov, V., Perrin, A., Polovtseva,E., Richard, C., Smith, M., Starikova, E., Sung, K., Tashkun, S., Tennyson, J., Toon, G., Tyuterev, V., and Wagner, G.: The HITRAN2012 molecular spectroscopic database, J. Quant. Spec-trosc. Ra., 130, 450, doi:http://dx.doi.org/10.1016/j.jqsrt.2013.07.002
Web End =10.1016/j.jqsrt.2013.07.002 http://dx.doi.org/10.1016/j.jqsrt.2013.07.002
Web End = , 2013.Scheepmaker, R. A., Frankenberg, C., Deutscher, N. M., Schneider,M., Barthlott, S., Blumenstock, T., Garcia, O. E., Hase, F., Jones,N., Mahieu, E., Notholt, J., Velazco, V., Landgraf, J., and Aben,I.: Validation of SCIAMACHY HDO/H2O measurements using the TCCON and NDACC-MUSICA networks, Atmos. Meas.
Tech., 8, 17991818, doi:http://dx.doi.org/10.5194/amt-8-1799-2015
Web End =10.5194/amt-8-1799-2015 http://dx.doi.org/10.5194/amt-8-1799-2015
Web End = , 2015.
Schneider, M. and Hase, F.: Improving spectroscopic line parameters by means of atmospheric spectra: Theory and example for water vapour and solar absorption spectra, J. Quant. Spectrosc.Ra., 110, 18251839, doi:http://dx.doi.org/10.1016/j.jqsrt.2009.04.011
Web End =10.1016/j.jqsrt.2009.04.011 http://dx.doi.org/10.1016/j.jqsrt.2009.04.011
Web End = , 2009.Schneider, M. and Hase, F.: Optimal estimation of tropospheric
H2O and D with IASI/METOP, Atmos. Chem. Phys., 11, 1120711220, doi:http://dx.doi.org/10.5194/acp-11-11207-2011
Web End =10.5194/acp-11-11207-2011 http://dx.doi.org/10.5194/acp-11-11207-2011
Web End = , 2011.
Schneider, M., Hase, F., Blavier, J.-F., Toon, G. C., and Leblanc,T.: An empirical study on the importance of a speed-dependent Voigt line shape model for tropospheric water vapor prole remote sensing, J. Quant. Spectrosc. Ra., 112, 465474, doi:http://dx.doi.org/10.1016/j.jqsrt.2010.09.008
Web End =10.1016/j.jqsrt.2010.09.008 http://dx.doi.org/10.1016/j.jqsrt.2010.09.008
Web End = , 2011.
Schneider, M., Hase, F., and Blumenstock, T.: Ground-based remote sensing of HDO/H2O ratio proles: introduction and validation of an innovative retrieval approach, Atmos. Chem. Phys., 6, 47054722, doi:http://dx.doi.org/10.5194/acp-6-4705-2006
Web End =10.5194/acp-6-4705-2006 http://dx.doi.org/10.5194/acp-6-4705-2006
Web End = , 2006.Schneider, M., Romero, P. M., Hase, F., Blumenstock, T., Cuevas,E., and Ramos, R.: Continuous quality assessment of atmospheric water vapour measurement techniques: FTIR, Cimel, MFRSR, GPS, and Vaisala RS92, Atmos. Meas. Tech., 3, 323 338, doi:http://dx.doi.org/10.5194/amt-3-323-2010
Web End =10.5194/amt-3-323-2010 http://dx.doi.org/10.5194/amt-3-323-2010
Web End = , 2010a.
Schneider, M., Seplveda, E., Garca, O., Hase, F., and Blumen-stock, T.: Remote sensing of water vapour proles in the framework of the Total Carbon Column Observing Network (TCCON), Atmos. Meas. Tech., 3, 17851795, doi:http://dx.doi.org/10.5194/amt-3-1785-2010
Web End =10.5194/amt-3- http://dx.doi.org/10.5194/amt-3-1785-2010
Web End =1785-2010 , 2010b.
Schneider, M., Barthlott, S., Hase, F., Gonzlez, Y., Yoshimura,K., Garca, O. E., Seplveda, E., Gomez-Pelaez, A., Gisi, M., Kohlhepp, R., Dohe, S., Blumenstock, T., Wiegele, A., Christner,E., Strong, K., Weaver, D., Palm, M., Deutscher, N. M., Warneke,T., Notholt, J., Lejeune, B., Demoulin, P., Jones, N., Grifth, D.W. T., Smale, D., and Robinson, J.: Ground-based remote sensing of tropospheric water vapour isotopologues within the project MUSICA, Atmos. Meas. Tech., 5, 30073027, doi:http://dx.doi.org/10.5194/amt-5-3007-2012
Web End =10.5194/amt- http://dx.doi.org/10.5194/amt-5-3007-2012
Web End =5-3007-2012 , 2012.
Schneider, M., Gonzlez, Y., Dyroff, C., Christner, E., Wiegele, A., Barthlott, S., Garca, O. E., Seplveda, E., Hase, F., Andrey, J., Blumenstock, T., Guirado, C., Ramos, R., and Rodrguez, S.: Empirical validation and proof of added value of MUSICAs tropospheric D remote sensing products, Atmos. Meas. Tech., 8, 483503, doi:http://dx.doi.org/10.5194/amt-8-483-2015
Web End =10.5194/amt-8-483-2015 http://dx.doi.org/10.5194/amt-8-483-2015
Web End = , 2015.
Sherwood, S. C., Bony, S., and Dufresne, J.-L.: Spread in model climate sensitivity traced to atmospheric convective mixing, Nature, 505, 3742, doi:http://dx.doi.org/10.1038/nature12829
Web End =10.1038/nature12829 http://dx.doi.org/10.1038/nature12829
Web End = , 2014.
Steinke, S., Eikenberg, S., Lhnert, U., Dick, G., Klocke, D., Di
Girolamo, P., and Crewell, S.: Assessment of small-scale integrated water vapour variability during HOPE, Atmos. Chem. Phys., 15, 26752692, doi:http://dx.doi.org/10.5194/acp-15-2675-2015
Web End =10.5194/acp-15-2675-2015 http://dx.doi.org/10.5194/acp-15-2675-2015
Web End = , 2015. Steinwagner, J., Fueglistaler, S., Stiller, G. P., von Clarmann, T.,
Kiefer, M., Borsboom, P.-P., van Delden, A., and Rckmann, T.: Tropical dehydration processes constrained by the seasonality of stratospheric deuterated water, Nature Geosciene, 3, 262266, doi:http://dx.doi.org/10.1038/ngeo822
Web End =10.1038/ngeo822 http://dx.doi.org/10.1038/ngeo822
Web End = , 2010.
Stevens, B. and Bony, S.: What Are Climate Models Missing?, Science, 340, 10531054, doi:http://dx.doi.org/10.1126/science.1237554
Web End =10.1126/science.1237554 http://dx.doi.org/10.1126/science.1237554
Web End = , 2013. Sutanto, S. J., Hoffmann, G., Scheepmaker, R. A., Worden, J.,
Houweling, S., Yoshimura, K., Aben, I., and Rckmann, T.: Global-scale remote sensing of water isotopologues in the troposphere: representation of rst-order isotope effects, Atmos. Meas. Tech., 8, 9991019, doi:http://dx.doi.org/10.5194/amt-8-999-2015
Web End =10.5194/amt-8-999-2015 http://dx.doi.org/10.5194/amt-8-999-2015
Web End = , 2015. Tuinenburg, O. A., Risi, C., Lacour, J. L., Schneider, M., Wiegele,A., Worden, J., Kurita, N., Duvel, J. P., Deutscher, N., Bony, S., Coheur, P. F., and Clerbaux, C.: Moist processes during MJO events as diagnosed from water isotopic measurements from the IASI satellite, J. Geophys. Res.-Atmos., 120, 1061910636, doi:http://dx.doi.org/10.1002/2015JD023461
Web End =10.1002/2015JD023461 http://dx.doi.org/10.1002/2015JD023461
Web End = , 2015.
Vogelmann, H., Sussmann, R., Trickl, T., and Reichert, A.: Spatiotemporal variability of water vapor investigated using lidar and FTIR vertical soundings above the Zugspitze, Atmos. Chem. Phys., 15, 31353148, doi:http://dx.doi.org/10.5194/acp-15-3135-2015
Web End =10.5194/acp-15-3135-2015 http://dx.doi.org/10.5194/acp-15-3135-2015
Web End = , 2015. Wiegele, A., Schneider, M., Hase, F., Barthlott, S., Garca, O.E., Seplveda, E., Gonzlez, Y., Blumenstock, T., Raffalski,U., Gisi, M., and Kohlhepp, R.: The MUSICA MetOp/IASI H2O and D products: characterisation and long-term comparison to NDACC/FTIR data, Atmos. Meas. Tech., 7, 27192732, doi:http://dx.doi.org/10.5194/amt-7-2719-2014
Web End =10.5194/amt-7-2719-2014 http://dx.doi.org/10.5194/amt-7-2719-2014
Web End = , 2014.
Worden, J. R., Bowman, K., Noone, D., Beer, R., Clough, S., Eldering, A., Fisher, B., Goldman, A., Gunson, M., Herman, R., Kulawik, S. S., Lampel, M., Luo, M., Osterman, G., Rinsland,C., Rodgers, C., Sander, S., Shephard, M., and Worden, H.: TES observations of the tropospheric HDO/H2O ratio: retrieval approach and characterization, J. Geophys. Res., 11, D16309, doi:http://dx.doi.org/10.1029/2005JD006606
Web End =10.1029/2005JD006606 http://dx.doi.org/10.1029/2005JD006606
Web End = , 2006.
Worden, J., Noone, D., Bowman, K., Beer, R., Eldering, A., Fisher,B., Gunson, M., Goldman, A., Herman, R., Kulawik, S. S., Lampel, M., Osterman, G., Rinsland, C., Rodgers, C., Sander, S., Shephard, M., Webster, R., and Worden, H.: Importance of rain evaporation and continental convection in the tropical water cycle, Nature, 445, 528532, doi:http://dx.doi.org/10.1038/nature05508
Web End =10.1038/nature05508 http://dx.doi.org/10.1038/nature05508
Web End = , 2007. Worden, J., Noone, D., Galewsky, J., Bailey, A., Bowman, K.,
Brown, D., Hurley, J., Kulawik, S., Lee, J., and Strong, M.: Estimate of bias in Aura TES HDO/H2O proles from comparison of TES and in situ HDO/H2O measurements at the
Mauna Loa observatory, Atmos. Chem. Phys., 11, 44914503, doi:http://dx.doi.org/10.5194/acp-11-4491-2011
Web End =10.5194/acp-11-4491-2011 http://dx.doi.org/10.5194/acp-11-4491-2011
Web End = , 2011.
Worden, J., Kulawik, S., Frankenberg, C., Payne, V., Bowman, K., Cady-Peirara, K., Wecht, K., Lee, J.-E., and Noone, D.: Proles of CH4, HDO, H2O, and N2O with improved lower tropospheric vertical resolution from Aura TES radiances, Atmos.
Meas. Tech., 5, 397411, doi:http://dx.doi.org/10.5194/amt-5-397-2012
Web End =10.5194/amt-5-397-2012 http://dx.doi.org/10.5194/amt-5-397-2012
Web End = , 2012.
Atmos. Meas. Tech., 9, 28452875, 2016 www.atmos-meas-tech.net/9/2845/2016/
M. Schneider et al.: MUSICA review on the observation of tropospheric {H2O, D} pairs 2875
Yoshimura, K., Oki, T., and Ichiyanagi, K.: Evaluation of two-dimensional atmospheric water circulation elds in reanalyses by using precipitation isotopes databases, J. Geophys. Res.-Atmos., 109, D20109, doi:http://dx.doi.org/10.1029/2004JD004764
Web End =10.1029/2004JD004764 http://dx.doi.org/10.1029/2004JD004764
Web End = , 2004.
Yoshimura, K., Miyoshi, T., and Kanamitsu, M.: Observation system simulation experiments using water vapor isotope information, J. Geophys. Res.-Atmos., 119, 78427862, doi:http://dx.doi.org/10.1002/2014JD021662
Web End =10.1002/2014JD021662 http://dx.doi.org/10.1002/2014JD021662
Web End = , 2014.
www.atmos-meas-tech.net/9/2845/2016/ Atmos. Meas. Tech., 9, 28452875, 2016
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright Copernicus GmbH 2016
Abstract
In the lower/middle troposphere, {H<sub>2</sub>O,<i>δ</i>D} pairs are good proxies for moisture pathways; however, their observation, in particular when using remote sensing techniques, is challenging. The project MUSICA (MUlti-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water) addresses this challenge by integrating the remote sensing with in situ measurement techniques. The aim is to retrieve calibrated tropospheric {H<sub>2</sub>O,<i>δ</i>D} pairs from the middle infrared spectra measured from ground by FTIR (Fourier transform infrared) spectrometers of the NDACC (Network for the Detection of Atmospheric Composition Change) and the thermal nadir spectra measured by IASI (Infrared Atmospheric Sounding Interferometer) aboard the MetOp satellites. In this paper, we present the final MUSICA products, and discuss the characteristics and potential of the NDACC/FTIR and MetOp/IASI {H<sub>2</sub>O,<i>δ</i>D} data pairs. <br><br> First, we briefly resume the particularities of an {H<sub>2</sub>O,<i>δ</i>D} pair retrieval. Second, we show that the remote sensing data of the final product version are absolutely calibrated with respect to H<sub>2</sub>O and <i>δ</i>D in situ profile references measured in the subtropics, between 0 and 7 km. Third, we reveal that the {H<sub>2</sub>O,<i>δ</i>D} pair distributions obtained from the different remote sensors are consistent and allow distinct lower/middle tropospheric moisture pathways to be identified in agreement with multi-year in situ references. Fourth, we document the possibilities of the NDACC/FTIR instruments for climatological studies (due to long-term monitoring) and of the MetOp/IASI sensors for observing diurnal signals on a quasi-global scale and with high horizontal resolution. Fifth, we discuss the risk of misinterpreting {H<sub>2</sub>O,<i>δ</i>D} pair distributions due to incomplete processing of the remote sensing products.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer