Vertical profile observations of water vapor deuterium excess in the lower troposphere

We use airborne measurements of water vapor (H2Ov) stable isotopologues and complementary meteorological observations to examine how boundary layer dynamics, cloud processing, and atmospheric mixing influence the vertical structure of δD, δO, and deuterium-excess (d-excess = δD – 8×δO) in the boundary layer, inversion layer, and lower free 15 troposphere. Flights were conducted around two continental U.S. cities in February – March 2016 and included vertical profiles extending from near the surface to ≤2 km. We examine observations from three unique case study flights in detail. One case study shows observations that are consistent with Rayleigh isotopic distillation theory coinciding with clear skies, dry adiabatic lapse rates within the boundary layer, and relatively constant vertical profiles of wind speed and wind direction. This suggests that the air mass retained the isotopic fingerprint of dehydration during moist adiabatic processes upwind of the study area. Also, 20 observed d-excess values in the free troposphere were sometimes larger than Rayleigh theory predicts, which may indicate mixing of extremely dehydrated air from higher altitudes. The two remaining case studies show isotopic anomalies in the dexcess signature relative to Rayleigh theory, and indicate cloud processes and complex boundary layer development. The most notable case study with stratocumulus clouds present had extremely low (negative) d-excess values at the interface of the inversion layer and the free troposphere, which is possibly indicative of cloud or rain droplet evaporation. We discuss how in situ 25 H2Ov stable isotope measurements, and d-excess in particular, could be useful for improving our understanding of moisture processing and transport mixing occurring between the boundary layer, inversion layer, and free troposphere.

Free troposphere (FT) entrainment and cloud evaporation influence the maintenance of the cloud layer, which in turn influences radiative forcing (Gerber et al., 2013;Yamaguchi and Feingold, 2013).The nature of H 2 O v as a climate feedback agent adds further complexity to our understanding of H 2 O v 's role in weather.Anthropogenic greenhouse gas emissions have resulted in increasing global temperatures, enhanced evaporation from soil and the oceans, and higher atmospheric concentrations of H 2 O v , the dominant absorber of infrared radiation (Held and Soden, 2006;Hurley and Galewsky, 2010;Willett et al., 2007).Warmer temperatures and more humid conditions have caused a shift towards less frequent but more intense precipitation events, increasing the risk of both floods and droughts (Roque-Malo and Kumar, 2017;Trenberth, 2011).H 2 O v also modulates production of the dominant atmospheric oxidant, the hydroxyl radical (Thompson, 1992).Thus, accurately representing H 2 O v in mesoscale processes is of great importance in a warming world.
H 2 O v stable isotopologue measurements are a potential tool for informing our understanding of the distribution and dynamics of H 2 O v in the lower troposphere (see review by Galewsky et al., 2016).H 2 O v stable isotopologue ratios, i.e., the ratio of heavy (HDO or H 18  2 O) to light (H 16 2 O) molecules, capture the origin and mixing of moisture sources as well as the condensation and cloud processes that modified that moisture.The δ notation indicates the sample's heavy-tolight isotope ratio reported relative to an international standard (δ = R sample /R standard − 1), where δ is multiplied by 1000 to report in per mil (‰).Both equilibrium fractionation, which separates heavy and light isotopologues based on their unique vapor-pressure differences, and kinetic fractionation, which is controlled by differences in diffusion rates of the isotopologues, modulate the isotopic signature of atmospheric H 2 O v (Dansgaard, 1964).Isotopic fractionation processes act to enrich and/or deplete both HDO and H 18  2 O relative to H 16  2 O in atmospheric waters, resulting in co-varying δD and δ 18 O signatures.The Rayleigh distillation model describes the equilibrium fractionation between vapor and condensation that dominates the variability of H 2 O v isotopes in the troposphere.This model can be used to calculate the degree to which condensate is removed from an air parcel as it cools, such as when it ascends from the surface to higher altitudes.Perturbations in δD and δ 18 O signatures often co-vary by a factor of approximately 8 because the ratio of the δD and δ 18 O equilibrium fractionation factors is approximately 8 : 1 at typical surface temperatures (Dansgaard, 1964).The second-order isotope parameter deuterium excess (d excess = δD-8 × δ 18 O) can be used to help identify the non-equilibrium kinetic fractionation processes.For example, remoistening of the atmosphere by raindrop evaporation imparts a kinetic isotopic fingerprint on H 2 O v d excess (Field et al., 2010;Lawrence et al., 2004;Risi et al., 2008;Worden et al., 2007;Wright et al., 2009).
Tropospheric vapor d-excess measurements may provide important information about cloud microphysics, convection processes, precipitation recycling, and FT entrainment.Local-scale studies have examined the role of moisture source, condensation history, and land surface evapotranspiration on d-excess signatures (Benetti et al., 2014;Delattre et al., 2015;Griffis et al., 2016;Kelsey et al., 2018;Lai and Ehleringer, 2011;Uemura et al., 2008;Welp et al., 2012).Measurements of d excess in precipitation have been used to investigate below-cloud evaporation (Aemisegger et al., 2015;Froehlich et al., 2008;Wang et al., 2016) and mixing between the boundary layer (BL) and FT from stationary platforms near the surface or at high-altitude mountain sites (Bailey et al., 2015;Benetti et al., 2015;2018;Froehlich et al., 2008;Galewsky, 2015;Lowenthal et al., 2016;Samuels-Crow et al., 2014).Kelsey et al. (2018) conducted mobile traverses along leeward and windward slope roads of a mountain in an instrumented vehicle to report vertical profiles of d excess.While high-elevation surface monitoring has the advantage of sampling BL and FT air over a diurnal cycle, it does not provide a complete picture of the H 2 O v isotope vertical profile (VP) at a discrete time.(Herman et al., 2014;Worden et al., 2012).Airborne platforms are capable of δD, δ 18 O, and d-excess VP measurements at higher spatiotemporal resolution and have been conducted since the 1960s, extending from the lower troposphere to the stratosphere for investigating a variety of scientific questions (overview in Sodemann et al., 2017).Due to either study objectives or limitations of the instrumentation (Dyroff et al., 2015;Herman et al., 2014), however, only one airborne H 2 O v isotope study, to our knowledge, has reported d-excess measurements in the lower troposphere to provide high-vertical-resolution snapshots at discrete time points in the boundary later (Sodemann et al., 2017).Additional airborne d-excess observations may improve our understanding of the role that cloud processes and convection have on determining the moisture distribution of the troposphere (Noone, 2012;Bolot et al., 2013).
In this study we present stable H 2 O v isotope measurements over two continental sites.These measurements provide a unique data set for understanding variability in the lower-troposphere d-excess profile and what it reveals about lower-troposphere moisture processing on relatively small regional scales.H 2 O v stable isotope VP measurements were conducted in the lower troposphere during four flights around the Washington, D.C., and Baltimore, MD, metropolitan area in February 2016 and during five flights around the Indianapolis, IN, metropolitan area in March 2016.We compare and contrast observations of the unique vertical structure of δD, δ 18 O, and d excess from 3 representative case study days.The case studies provide information about meteorological conditions that produce H 2 O v isotopic VP profiles indicative of past Rayleigh dehydration of the air mass and those where other processes must explain the observations.The case study observations reveal d-excess features unique to cloud droplet or rain droplet evaporation and show the influence of synoptic weather patterns and urban versus rural differences on BL development.Interpretations of case study VPs are supported with observations from the remaining flight days in the Washington, D.C., and Baltimore, MD, and Indianapolis areas.

Study sites
Flights were conducted around the Washington, D.C. to Baltimore, MD, area in February 2016 and around Indianapolis in March 2016.Washington, D.C., and Baltimore, MD, is a metropolitan area of 9.8 million residents that includes the District of Columbia and encompasses parts of Maryland, Virginia, West Virginia, and Pennsylvania (US Census Bureau, 2018).The Appalachian Mountains lie to the west of Washington, D.C., and Baltimore, MD, and the Chesapeake Bay and the Atlantic Ocean lie on the eastern side of Washington, D.C., and Baltimore.By contrast, Indianapolis has a population of 2.0 million and is relatively isolated from other metropolitan areas by agricultural fields (US Census Bureau, 2018).The closest large body of water to Indianapolis is Lake Michigan, over 200 km to the north.

Airborne Laboratory for Atmospheric Research (ALAR)
The Purdue Airborne Laboratory for Atmospheric Research (ALAR) is a modified twin-engine Beechcraft Duchess aircraft.The ALAR's two rear passenger seats have been removed to make room for scientific instrumentation.Ambient air at the nose of the aircraft is pulled through a forwardfacing unheated 5 cm diameter PFA Teflon inlet called the "main manifold" at a flow rate of 1840 L min −1 using a blower installed at the rear of the aircraft.Residence time in the main manifold is ≤ 0.1 s.Instruments sample from the main manifold with individual Swagelok "T " connections and Teflon sampling lines.The Purdue ALAR is equipped with a global positioning and internal navigation system (GPS/INS) for 50 Hz geopositional measurements and a Best Air Turbulence (BAT) probe for 50 Hz three-dimensional winds and pressure measurements (Crawford and Dobosy, 1992;Garman et al., 2006Garman et al., , 2008)).Temperature measurements are made with a microbead thermistor installed in the center pressure port of the BAT probe (Garman, 2009 S1 in the Supplement).The slower turnover of air in the LGR compared to the Picarro results in LGR H 2 O v measurements that are effectively low-pass filtered when compared to the Picarro measurements.
The LGR TWVIA isotopic measurements were calibrated in the lab for H 2 O v concentration dependence before and after the field campaign using an LGR Water Vapor Isotope Standard Source (WVISS; model: 908-0004-9003) with five standards ranging in isotope enrichment from −39.9 ‰ There was no need for an additional correction to normalize to the VSMOW-SLAP (Vienna Standard Mean Ocean Water -Standard Light Antarctic Precipitation) scale (discussed in Sect.S2; Fig. S2.3).Discussion of the instrument precision and calibration uncertainties is provided in Sect.S3.Total uncertainties, the quadrature sums of instrument precision and calibration uncertainties, increase as the H 2 O v mole fraction decreases below 4000 ppmv (Fig. 1).Flight measurements of δD, δ 18 O, and d excess reported here are smoothed using a 20 s moving average which corresponds to the time required for the TWVIA-reported δ values to stabilize after a change in the sample's H 2 O v mole fraction or isotopic signature (Sect.S3).

Flight design
Nine daytime research flights were conducted around Washington, D.C., and Baltimore, MD, and Indianapolis in February and March 2016 (Table 1; Fig. 2).Flight paths were de-signed to maximize the number of vertical profiles (VPs) conducted while also characterizing upwind-downwind gradients in the H 2 O v isotopic signature.VPs were sometimes conducted in a spiral pattern to limit the horizontal spatial coverage of the measurements, while other VPs were conducted in a sawtooth pattern ("porpoising"; Gerber et al., 2013) between the BL and FT.The aircraft flew up to ∼ 1600 m a.s.l. on average during the VPs.Only data collected on the descents of the VPs, when sampled air transitions from relatively dry to relatively humid are presented here to minimize the potential influence of memory effects.However, similar features were observed on the ascents and descents.The number of VPs (Table 1) conducted on each flight was limited by air traffic and restricted air space (which was worse for the Washington, D.C., and Baltimore, MD, study site), cloud cover, and available flight time.The research aircraft typically does not fly through clouds during experimental flights, as the BAT probe has electronics exposed to air.Flights included other maneuvers, such as transects conducted upwind, intersecting, and downwind of the urban centers, the interpretation of which is beyond the scope of this paper (Fig. 2).

Atmospheric layer identification
We classify regions of the atmosphere into the boundary layer (BL), INV, and FT to compare and contrast observed isotopic features during the research flights.The altitude at the base of the INV (z INV ) is defined as the lowest altitude at which the change in potential temperature (θ ) exceeds 0.5 K for a 10 m change in altitude (dθ/dz>0.05K m −1 ).Rates of dθ/dz>0.05K m −1 were commonly observed within the INV during the research flights.The altitude at the base of the FT (z FT ) is defined as the altitude above z INV at which dθ/dz transitions to < 0.05 K m −1 .A recent evaluation of methods for determining boundary layer height from aircraft measurements indicates that the potential temperature gradient approach is most reliable (Dai et al., 2014).However, if layers are not definable using the dθ/dz > or < 0.05 K m −1 criterion, the secondary criterion of |d(H 2 O v )|/dz>20 and |d(H 2 O v )|/dz<20 ppmv m −1 is used to define z INV and z FT , respectively.These threshold values are appropriate for our wintertime, midlatitude observations but may not be universally appropriate in different locations or seasons.If neither criterion is met, the profiles of θ , dθ/dz, H 2 O v , and d(H 2 O v )/dz are collectively considered, and z INV is visually defined as the point at which H 2 O v and θ begin decreasing and increasing, respectively.Similarly, z FT is visually determined as the altitude at which the rate of change of H 2 O v and θ with altitude begins to decrease.

Isotope theory
There are many processes that can influence the isotopic composition of atmospheric H 2 O v ; here we examine  two common ways: condensation and mixing.We use the Rayleigh distillation model to represent condensation of an ascending air parcel and a simple two-member mixing model to represent atmospheric mixing.We calculate the influence of these processes using R notation but present the results by converting to delta notation (δ = R sample /R standard − 1) using VSMOW as the standard.
The Rayleigh distillation model describes the effects of equilibrium fractionation on the isotopic composition of a dehydrating air parcel, for example during ascent (Dansgaard, 1964).Condensate that is formed as an ascending air parcel expands and cools is enriched in the heavier isotope relative to the vapor, and in the open form of the Rayleigh model, is assumed to be immediately removed from the system (Eq.1): Here R o and R Ray are the heavy-to-light isotopologue ra- ) of the parcel prior to the ascent and at any point throughout the ascent, respectively.The remaining fraction of H 2 O v left in the ascending parcel relative to initial conditions is given by H 2 O vRay H 2 O vo .We determined the initial R o and H 2 O vo input values for each day from the average BL values measured along the VP descents.The temperaturedependent equilibrium fractionation factor, α e , for each isotopologue is calculated using the temperature of the air parcel's lifting condensation level (LCL) altitude using findings of Horita and Wesolowski (1994) for LCL temperatures greater than 0 • C and those of Ellehøj et al. (2013) for LCL temperatures less than 0 • C. The LCL is the height at which an air parcel would become saturated along a dry adiabatic ascent and is often used as an estimate of cloud base height (Romps, 2017).The VP observations show that ambient temperatures vary with altitude along the vertical profiles.However, Rayleigh distillation curves calculated with α e values defined by the varying ambient temperatures measured along the vertical profiles are nearly identical to Rayleigh curves calculated with a single LCL-defined α e value (Fig. S4).
The mixing of two air parcels (A and B) results in a heavyto-light isotopologue ratio of an air parcel, R mix , given by Eq. (2) using HDO v as an example.R mix is the ratio of the weighted average of the heavy isotopologue to the weighted average of the light isotopologue.The fractions of air parcel A, f A , and air parcel B, f B , sum to unity.The mixture's H 2 O v mole fraction is simply the weighted average of each parcels' individual H 2 O v mole fraction.H 18 2 O v can replace HDO v in Eq. ( 2): We consider mixing processes in our observations by choosing two regions of the VPs as potential end members, for example, the BL and the FT, and use the observed H 2 O v and isotopic ratios to define end-member values.
3 Airborne campaign observations of H 2 O v d excess in different layers of the lower troposphere The campaign-wide observations show that values of δD and δ 18 O decrease as H 2 O v mole fractions approach zero (Fig. 3).This relationship results from preferential removal of the heavier isotopologues during condensation processes, which becomes more pronounced at colder temperatures, and is consistent with previous airborne and high-altitude measurements of H 2 O v stable isotopologues (Bailey et al., 2013;Galewsky et al., 2007;He and Smith, 1999;Noone et al., 2013;Samuels-Crow et al., 2014;Sodemann et al., 2017).Later we will return to examine individual profiles of δD and δ 18 O on case study days, but here we focus on identifying common patterns in d-excess signatures observed in the BL, INV, and FT that can be used to understand processes controlling moisture in the lower troposphere.The Washington, D.C., and Baltimore, MD, and Indianapolis VP d-excess observations within the FT, INV, and BL are presented as a function of the H 2 O v mole fraction in Fig. 4a-c, respectively.Generally, the air became drier and the d-excess signature exhibited greater variability with increasing altitude (Fig. 4).The greatest variability was observed in the INV (Fig. 4b), where the d-excess signature deviated both positively and negatively from the global average precipitation d-excess value of 10 ‰ (which is provided in Fig. 4 for reference only).The FT showed low H 2 O v mole fractions as well as large positive d-excess values (Fig. 4a), which are predicted by Rayleigh distillation theory for very low H 2 O v mole fractions (Bony et al., 2008).High FT dexcess signatures have been reported by other studies, which hypothesize that extremely dry, depleted air masses in the mid-to-upper troposphere with large positive d-excess signatures mix downward towards flight-level altitudes in the  lower FT (Bony et al., 2008;Samuels-Crow et al., 2014;Sodemann et al., 2017).These FT air masses likely originated from another source region and possibly underwent multiple condensation cycles to achieve such isotopic depletion prior to mixing with more humid air across the INV.Thus, FT d-excess values likely act as a record of the condensation and isotopic depletion history of a transported air parcel.
The 7 March 2016 flight in Indianapolis (MAR07) reveals an unusual set of meteorological conditions because, unlike all other flight days of the campaign, the H 2 O v mole fraction increased with altitude in the lower troposphere (Fig. 4a-c; gold trace), and it was the most humid day of the entire campaign (Fig. 3; gold trace).A warm, southerly front moved into the Indianapolis study area on this day, and rain preceded the flight measurements.The relatively high H 2 O v mole fractions in the INV and FT likely reflect residual humidity from the storm.Overall, the MAR07 VP observations do not exhibit distinctive d-excess features, i.e., not deviating much from 10 ‰ in the INV and FT (Fig. 4).

Case studies
We focus on three particular Indianapolis flights as representative case studies based on their distinct features and vertical profiles over large altitude ranges (Table 1).H 2 O v isotope measurements conducted within and above the cloudless, well-mixed boundary layer on 6 March (CLR) represent the meteorologically simplest observations of the entire campaign, as no clouds, precipitation, or shifting synoptic conditions were observed.Isotopic observations on 4 March (STC) may reflect moisture processing in a stratocumulus topped BL, and the 18 March (developing boundary layer -DBL) observations may reveal differences in urban versus rural BL development and the influence of changing synoptic conditions.

Clear-sky observations of a well-mixed boundary layer (CLR)
Four VPs were conducted on 6 March 2016 ("CLR" case study) around Indianapolis during clear-sky conditions (Fig. 5a; weather map is presented in Fig. S5.1).Cloud-top height estimated from the Terra MODIS satellite retrievals (https://worldview.earthdata.nasa.gov/,last access: 5 February 2018) indicates that the sparse cloud cover shown in Fig. 5a corresponds to higher-altitude (> 4800 m) clouds.The CLR measurements were made below 1400 m a.s.l.(Fig. 5b) and, as a result, were likely not impacted by higher-cloud processes.In terms of meteorology, the CLR case study is the simplest flight day of the airborne campaign.It is a useful case study for examining isotopic signatures across the BL, INV, and FT without the influence of complex atmospheric circulations or vapor-condensate interactions from clouds or precipitation.
The VP measurements made during CLR indicate a wellmixed BL and FT, and wind speed and wind direction were relatively constant from the BL to FT, indicating that BL, INV, and FT air parcels shared recent advection histories.δ 18 O, δD, and d excess are relatively constant with altitude within the BL (surface to z INV ), varying by 1.2 ‰, 15.3 ‰, and 10.9 ‰, respectively (Fig. 6a).VP2 observations are presented as a representative example of CLR because VP2 was conducted approximately midway through the flight (Fig. 5b), it covers the largest vertical range, and it was conducted in a spiral formation to minimize the horizon- tal spatial extent over which the measurements were made (Fig. 5a).The ambient temperature profile approximately follows the dry adiabatic lapse rate to the top of the BL (Fig. 6a Observations of δD, δ 18 O, and d excess measured along the CLR VP descents are plotted as a function of the H 2 O v mole fraction in Fig. 7 along with predictions from Rayleigh distillation theory (Fig. 7a-c) and different BL-FT mixing scenarios (Fig. 7d-f; Sect.2.5).The measured δ values more closely match Rayleigh curves than the mixing lines from the BL through the INV (Fig. 7).However, positive deviations relative to Rayleigh d excess exist in the upper INV and lower FT (Fig. 7c).We hypothesize that dry, isotopically depleted FT air parcels carrying large, positive d-excess values can mix downward into more humid air parcels of smaller d-excess values near the top of the INV (Sodemann et al., 2017).As the H 2 O v mole fraction approaches zero, Rayleigh-predicted d excess approaches 7000 ‰ (Bony et al., 2008) 2015), who report lower-troposphere δD observations over the Atlantic Ocean and explain the vertical structure of δD at lower altitudes using Rayleigh theory, while mixing scenarios dictate the δD profile based on higher-altitude observations.Despite apparent similarities between observed and Rayleigh d excess in the BL and INV, the meteorological conditions within the BL are contrary to assumptions of the Rayleigh distillation model.Rayleigh distillation theory describes a two-phase system at saturation; however, no clouds (i.e., condensate) were observed at the flight altitudes where the CLR measurements were made.Furthermore, the BL was dry adiabatic (Fig. 6a).The agreement between the CLR observations and the Rayleigh-predicted d excess, however, could be described by upwind condensation consistent with Rayleigh distillation theory followed by subsequent advection to the Indianapolis study site.Given that wind direction and wind speed were observed to be relatively constant throughout the lower troposphere (Fig. 6a), air parcels in the BL, INV, and FT likely share similar trajectories.Thus, the Rayleigh-consistent CLR VP observations at lower flight altitudes likely result from upwind condensation that imprinted an isotopic signature on the air parcels that was maintained during transport, while lower FT observations point to mixing between subsiding FT and INV air.Past studies have also reported on the usefulness of H 2 O v stable isotopes as "imprints" of condensation, mixing, and atmospheric transport (Bailey et al., 2013;Brown et al., 2008;Galewsky et al.,  2007 ;Gedzelman, 1988;He and Smith, 1999;Samuel-Crow et al., 2014;Taylor, 1984).

Stratocumulus-topped boundary layer observations (STC)
The central and eastern portions of the Indianapolis study area were covered by stratocumulus clouds during part of the 4 March 2016 flight ("STC" case study) on which five VPs were conducted (Fig. 8).The cloud cover map in Fig. 8a is provided to show the cloud type and extent during the afternoon of STC; however, it does not necessarily represent the cloud cover conditions throughout the 2.5 h flight.Satellite images (not shown) captured during the early afternoon of STC show that a thick cloud cover was sustained from the beginning into the middle of the flight, particularly over the city of Indianapolis, but eventually transitioned to scattered cover throughout the afternoon.This is consistent with visual observations made by the pilot and mission scientist that thick clouds persisted from the beginning of the flight until approximately 15:30 LT (local time).Vertical profile temperature and RH measurements collected on this day support the presence of a cloud layer (Fig. 6b).BL air was nearly saturated at 788 m (z CB for "cloud base") along VP2 on STC (Fig. 6b).The ambient temperature lapse rate is 8.8 K km −1 (close to the dry adiabatic lapse rate of 9.8 K km −1 ) near the surface until an altitude of z CB , where the lapse rate transitions to 2.8 K km −1 .These observations are indicative of a  6).Contrary to the CLR INV d-excess profile (Fig. 6a), d excess on STC first increases with altitude within the lower INV before decreasing to a minimum at z FT (Fig. 6b).Obvious anomalies in the STC d-excess signature relative to CLR are also apparent when plotted as a function of the H 2 O v mole fraction (Figs.9c and  7c).
While the atmosphere on STC was unsaturated at most flight altitudes (except near the cloud base; Fig. 6b), the comparison to open-system Rayleigh curves is a useful exercise, as we showed above in the discussion of the CLR case study.STC VP δ 18 O and δD observations deviated negatively from the Rayleigh curve in drier portions of the INV and FT and were more pronounced for δD than δ 18 O (Fig. 9a-b).Most of the VPs' δD values are more negative relative to Rayleigh in the INV and plateau in the FT, which is unusual in water vapor isotope observations (Fig. 9b).Mixing processes can produce δ 18 O and δD values that plot very nearly on or above the open-system Rayleigh curve depending on the mixing end members (e.g., Fig. 9d, e, and f), whereas processes involving non-equilibrium liquid-vapor interactions or closed-system Rayleigh processes can plot below the Rayleigh curve.The FT-BL mixing scenario is a poor match to observations, especially for d excess, whereas the FT-zFT and zFT-BL mixing scenarios show remarkable agreement with the VP dexcess observations on STC (Fig. 9d, e, and f).However, the mixing scenarios do not identify the third end-member source responsible for the minimum in d excess at z FT .
The d-excess measurements along VP2 through VP5 (Fig. 9c) reveal two anomalies that yield some insight: (1) the slight increase in d excess in the middle of the inversion layer (particularly for VP2 and VP5) and (2) the d-excess minimum at the INV-FT interface (z FT ).From this minimum at z FT , the FT d-excess signature becomes more positive with increasing altitude (as the air becomes drier) and eventually transitions to being more positive than the Rayleigh curve (Fig. 9c).
The d-excess values observed during VP1 represent the only VP measurements on STC that do not exhibit a minimum at the INV-FT interface (Fig. 9c).Only VP1 was conducted before the research aircraft encountered the thick stratocumulus cloud layer over Indianapolis (Fig. 8a).Unlike VP, VP1 d excess tracks the Rayleigh curve at the INV-FT interface.Slightly above z FT , the VP1 H 2 O v mole fraction increases and d excess approaches the mixing line (Fig. 9c  and f).We believe that this VP represents conditions prior to cloud process influences.
The presence of the stratocumulus cloud layer is a defining meteorological characteristic of the STC case study day; therefore, we evaluate the potential for cloud and rain processes to cause the observed d-excess anomalies.Sodemann et al. (2017)

Sodemann et al., 2017). HDO molecules preferentially evaporate relative to H 18
2 O molecules (Dansgaard, 1964).The result is a relatively positive vapor d excess, while the d-excess signature of the residual droplet becomes progressively more negative as it evaporates (Aemisegger et al., 2015).Therefore, if liquid droplet evaporation occurred in separate atmospheric layers from start to finish, such as from the bottom of the INV to the top of the INV, a positive d-excess anomaly could be transferred to the surrounding vapor as the droplet begins evaporating starting in the lower INV.As the liquid droplet is subsequently transported to the upper INV, a negative d-excess anomaly could be transferred to the surrounding vapor as the evaporation of the liquid droplet nears completion.Evidence exists that stratocumulus cloud droplets can evaporate at different altitudes, specifically within the cloud layer and within the inversion layer (Kollias and Albrecht, 2000;de Lozar and Melledo, 2015).Furthermore, inversion layer depths are not homogeneous above a stratocumulus cloud layer.Observations show inversion layers to be thicker above downdrafts and thinner near updrafts (Kollias and Albrecht, 2000).Differing INV depths, and associated thermodynamic properties, above the stratocumulus cloud layer on STC may in part explain why unique d-excess anomalies were observed near the middle of the INV and at z FT (Fig. 6b and 9).
When evaluating the effect that cloud evaporation could have on VPs of H 2 O v isotopologues, we must consider the altitudes at which clouds form and evaporate and the speed of the relevant processes.Stratocumulus cloud tops are typically directly below z INV (Wood, 2012).The top of the INV (z FT ) is approximately the upper limit of BL mixing (Wood, 2012).Lofting of cloud droplets into the INV would cause droplet evaporation, as the INV was not saturated (Fig. 6b).This hypothesis requires that a negative d-excess anomaly is retained by an evaporating droplet while it is transported to the upper INV.This requires us to consider two important timescales: first, the timescale at which a liquid droplet isotopically equilibrates with its surrounding vapor, and second, the speed at which droplets move through the INV by vertical winds.Using the method described by Bolot et al. ( 2013), which is based on prior work by Jouzel (1986), we estimate that the e-folding time required for a cloud droplet with a 15 µm radius to isotopically equilibrate with the surrounding vapor is approximately 2 s under the conditions observed in the middle of the INV or about 5 s at the top of the INV.The time for a droplet to move from the bottom to the top of the INV is 19 s based on vertical wind speed observations.These calculations indicate that a cloud droplet located near the bottom of the INV would approach isotopic equilibrium with surrounding vapor before reaching the top of the INV under the observed flight conditions.Thus, these calculations do not strongly support the hypothesis that evaporation of stratocumulus cloud droplets produced the observed minimum in d excess at z FT .
The question of what process or processes are responsible for the minimum in d excess at the INV-FT interface (z FT ) remains.A key characteristic of stratocumulus cloud layers is the longwave cooling at the cloud tops which maintains and enhances in-cloud turbulence (Wood, 2012).If previous INV conditions on STC were colder, had faster vertical wind speeds, or the liquid droplet radii were larger than 50 µm (e.g., rain drizzle), the timescales associated with isotopic equilibrium and transport across the INV could converge.The main difference between the cloud and rain droplet evaporation processes is the size of the liquid droplets: < 50 µm for cloud versus > 50 µm for rain droplets (Kollias and Albrecht, 2000).It is possible that such conditions were present prior to when the STC measurements were made and that the minimum in d excess is the isotopic imprint of rain droplet evaporation under cold, unsaturated, or turbulent conditions.We find discussion of these potential d-excess anomalies in the literature in reference to rain droplets evaporating below the cloud layer (Aemisegger et al., 2015;Gat, 1996), but we believe that the same isotopic fingerprint on vapor could occur as liquid cloud or rain drizzle droplets lofted in drier environments finish evaporating at the top of the INV, or at the top of the BL, as hypothesized by Sodemann et al. (2017).
The negative H 2 O v d-excess observations at z FT could also have resulted from kinetic fractionation of vapor during deposition on ice crystals or snow (i.e., in ice-supersaturated conditions; Bolot et al., 2013;Casado et al., 2016;Galewsky, 2015;Lowenthal et al., 2016;Moore et al., 2016;Samuels-Crow et al., 2014;Schmidt et al., 2005).Low d-excess values (relative to Rayleigh curves assuming that RH = 100%) were sometimes observed within the INV layers during the Washington, D.C., and Baltimore, MD, flights (Fig. 4b).Ambient temperatures observed in flight in Washington, D.C., and Baltimore, MD, were sometimes less than 0 • C (Table 1); thus vapor deposition on ice crystals could be possible for those scenarios.It is unlikely that vapor deposition on ice occurred during the STC case study flight because temperatures were greater than 0 • C (Table 1).However, as an example, Fig. S6 shows the theoretical d-excess values of STC vapor under ice-supersaturated conditions (Sect.S6).We reiterate that ice supersaturation is an unlikely explanation for the STC z FT d-excess minimum because flight altitudes were less than 2 km, and ice (cirrus) clouds are typically present at ∼ 6 km.It is unlikely that ice hydrometeors falling from higher altitudes could be sustained at the top of the inversion and contribute to the low d-excess signal observed on VP2-VP5 through vapor deposition given the > 0 • C temperatures.There was, however, a region of high humidity upwind (northwest) of the study site at altitudes between 3 and 5.5 km, where ice or mixed-phase condensate could have been present (Fig. S5.2).It is possible that condensation under ice-supersaturated conditions occurred prior to the STC flight and that the resulting isotopic imprint was maintained during transport to Indianapolis and subsequently mixed downward via subsiding FT air (Fig. 9f).Both of these explanations for the minimum in d excess at z FT on STC require an advected signal of a prior process (complete cloud droplet evaporation or vapor deposition on ice).This process must have happened relatively close in time to the STC flight, since the minimum in d excess at z FT was not observed on VP1, but the anomaly was observed 2 h later when VP2-VP5 were conducted.

Developing boundary layer observations (DBL)
The final airborne case study, conducted on 18 March 2016 in Indianapolis (Fig. 10), is referred to as DBL because measurements on this day reveal considerable spatiotemporal variability in the vertical structure of the observed meteorological and isotopic variables.The boundary layer height increased over the course of the flight and may reflect a combination of a residual layer from the previous day, urban vs. rural differences in BL development, and the effects of a frontal pattern moving across the Indianapolis study area.
A defining characteristic of the DBL case study is the variability in both the meteorological and isotopic variables between each of the four VPs (Fig. 11).Observations along VP1 show an INV layer, marked by a characteristic increase in θ and a corresponding decrease in H 2 O v with altitude (Fig. 11a), separating the BL and FT.However, there appears to be two distinct atmospheric layers between the BL and FT in VP2 (Fig. 11b).The layer directly below the FT in VP2 is consistent with a residual layer (RL) from the previous day's boundary layer (Fig. 11b).We define the base of the RL using the same approach described in Sect.2.4 (dθ/dz and |d(H 2 O v )|/dz threshold values) for determining the base of the INV (z INV ).Both the RL and the INV (directly below the RL) show characteristic decreases in H 2 O v , δ 18 O, and δD values with altitude.The presence of multiple layers is supported by the increase in the variance of the vertical wind speed (W σ 2 ), indicating wind shear, at the interface of atmospheric layers (Fig. 11b).Similarities in the vertical structure of H 2 O v and d excess between VP2 and VP3 give indications that the RL persisted for an additional hour after VP2 was conducted (Fig. 11c).Despite a temporary failure of the aircraft's wind measurement system halfway through the VP3 descent, available measurements show an increase in wind speed and vertical wind variance (W σ 2 ) as well as sharp temperature changes at the base and top of the RL, supporting VP3's d-excess indications of the persisting RL (Fig. 11c).
This residual layer hypothesis is supported by ambient temperature profiles from the Rapid Refresh model-data assimilation system (Iowa State University, https://mesonet.agron.iastate.edu/archive/,last access: 7 May 2018) for the Indianapolis International Airport (KIND) from the day prior to DBL, 17 March 2016, which follow a nearly dry adiabatic lapse rate all the way to an altitude of 3 km (Fig. S7; archived at https://mesonet.agron.iastate.edu/archive/,last access: 7 May 2018).This was a relatively warm, turbulent day.A cold front moved into the Indianapolis study area on the DBL (18 March 2016).The ambient temperature profile on the DBL shows the previous day's residual layer persisted into the early afternoon, between approximately 1 and 3 km (Fig. S7), before being incorporated into the BL.The distinction between the RL and BL blurs as surface heating progresses throughout the day and the RL is incorporated into the BL.It would also be expected that the downwind edge of the Indianapolis city boundaries (winds were from the northwest; Fig. 11a) would have a more well-mixed BL due to stronger turbulent mixing from the urban heat island and increased surface roughness (Grimmond et al., 2010;Stull, 1988).Support for this is given by measurements made along VP3, which were conducted on the downwind edge of the Indianapolis city boundaries and reveal a considerably more homogenous structure in δD and δ 18 O relative to VP1 and VP2 (Fig. 11c).The ambient temperatures measured along VP3 in the FT and RL are warmer relative to ambient temperatures along the other three VPs (Fig. 11), demonstrating the influence of the urban heat island.Although the VP3 H 2 O v , δD, and δ 18 O values are relatively more homogenous in the vertical dimension, the d-excess signatures still maintain indications of the RL, as the vertical structure of d excess is similar to VP2 d-excess observations.
Measurements of δ 18 O and δD values along DBL VP3 are the only observations during our airborne campaign measurements of enriched, and relatively constant, δ values that extend from the BL through to the FT (Fig. 11c).FT H 2 O v mole fractions along VP3 only decreased to ∼ 2500 ppmv, whereas H 2 O v mole fractions of 1700 ppmv and less were observed in the FT of .Relatively humid conditions were also observed at the highest altitudes flown in the FT along VP4 (Fig. 11d).Enriched, vertically unvarying δ 18 O and δD values are unique to DBL VP3.VP4 δ 18 O and δD values decrease across the INV (Fig. 11d) and track a δ-H 2 O v path similar to VP1 and VP2 from the BL up to the INV-FT interface, which are similar in shape to both Rayleigh predictions or BL-FT mixing scenarios (Fig. 12).Interestingly, VP4 FT δ 18 O and δD values become enriched at higher altitudes, corresponding to an increase in H 2 O v mole fractions (Fig. 11d), and appear to track a mixing line with the VP3 z INV end-member δ values (red trace in Fig. 12d-e).Despite differences in δ-value features along the four DBL VPs, the relationship between d excess and the H 2 O v mole fraction appears to be remarkably consistent throughout the day (Fig. 12c).
The relatively humid, isotopically enriched FT air observed along VP3 and VP4 could have been caused by the shortwave trough in the mid-troposphere (3-5 km) which carried moist air into the Indianapolis study air in the late afternoon on this day (Fig. S5.3).The dew-point profile in Fig. S7 shows this relatively moist mid-tropospheric air subsiding over the course of the afternoon, reaching flight altitudes by the time VP3 and VP4 were flown.The vertically unvarying, isotopically enriched VP3 observations likely reflect a combination of enhanced turbulence from the urban area and humid air from the shortwave trough mixing down- ward into the lower FT.In contrast, VP4 was conducted over a rural area north of Indianapolis and likely did not experience the same degree of vertical mixing as the atmosphere downwind of Indianapolis (where VP3 was flown; Fig. 10a).Potentially due to weaker turbulence in the rural area, a sharp decrease in δD and δ 18 O was observed at z FT on VP4 before increasing with altitude, reaching enriched δ values observed in the VP4 BL and throughout VP3.
The DBL case study shows how isotopic water vapor measurements can be used as tracers to track the development of different atmospheric structures and circulations, including residual layers, urban heat island impacts, and passing fronts.δD, δ 18 O, and particularly d excess can support meteorological measurements in identifying distinct atmospheric layers that are difficult to identify solely based on meteorological data.

Conclusions
The aim of this study is to observe and interpret the vertical structure of H 2 O v stable isotopic composition, specifically d excess, in the continental lower troposphere.Previous studies presenting d-excess observations are heavily focused on ocean evaporation at coastal or island surface sites (Benetti et al., 2014(Benetti et al., , 2015(Benetti et al., , 2018;;Delattre et al., 2015;Steen-Larsen et al., 2014;Uemura et al., 2008).Few reported observations of d excess in the INV and FT exist (Galewsky et al., 2015;Lowenthal et al., 2016;Kelsey et al., 2018;Samuels-Crow et al., 2014;Schmidt et al., 2005;Sodemann et al., 2017), and, to our knowledge, only one study has used airborne measurements to provide high-vertical-resolution snapshots of the lower-troposphere d-excess profile at discrete time points (Sodemann et al., 2017).Our stable H 2 O v isotope measurements over two continental sites are a starting point in filling the field's gap in understanding variability in the lowertroposphere d-excess profile and what it reveals about lowertroposphere moisture processing on relatively small regional scales.
Our observations reaffirm the dominant role that Rayleigh distillation processes have on H 2 O v isotopic signatures and that these isotopic signatures can be retained by air masses as they are advected from previous condensation events.Our measurements indicate that H 2 O v isotopologues, and d excess in particular, can act as fingerprints of earlier processes.This is true for the CLR case study, in which Rayleigh distil-lation curves represent vertical profile observations well from near the surface up through the atmospheric inversion layer despite the temperature profile within the BL following the dry rather than the (Rayleigh-assumed) moist adiabat.Similarly, the evidence of upwind processes was also retained in the H 2 O v isotopic signature on STC.The STC measurements give clues that the observed air mass experienced prior evaporation of cloud or rain droplets.Measurements made during the rapidly changing atmospheric conditions on the DBL case study also show that the H 2 O v isotopic signature can reflect near-instantaneous changes in the atmosphere.
These new results highlight the potential for H 2 O v isotope ratios, especially d excess, to identify complex processes across the atmospheric inversion layer, including cloud condensation, evaporation, and mixing or entrainment of free tropospheric air into the boundary layer.On each case study day, the combination of all three isotope metrics provided additional information that may have been missed from either δD or δ 18 O alone.For example, our results show that d excess clearly distinguishes Rayleigh from mixing processes in the FT on the CLR case study day.The potential role of cloud or raindrop evaporation is most obvious in the d-excess observations followed by some indication in the δD vertical profile as well on the STC case study day.Our observations of the d-excess profile in a stratocumulus cloud-topped boundary layer and the d-excess observations reported by Sodemann et al. (2017) near marine cumulus clouds represent an opportunity to investigate the sensitivity of the dexcess signature to different classes of clouds and their associated unique cloud processes.Finally, on the DBL case study day, both d excess and the individual isotope ratios were able to distinguish unique air masses mixing at different altitudes evolving throughout the afternoon.These types of measurements may become increasingly valuable as we seek to understand the physical processes that sustain cloud layers and spatiotemporally variable boundary layer mixing depths.
Our interpretation of the d-excess VPs could be further evaluated by isotope-enabled circulation and weather models (Aemisegger et al., 2015;Pfahl et al., 2012;Schmidt et al., 2005).However, the simulation of convective boundary layer processes with isotope-enabled models is complex (Bolot et al., 2013;Benetti et al., 2018).The measurements reported here could help further develop current and forthcoming isotope enabled models, particularly for simulating wintertime, continental lower-troposphere processes or stratocumulus evaporation.
Data availability.Geopositional, meteorological, greenhouse gas, and water vapor isotope measurements are available for the Washington, D.C., and Baltimore, MD, and Indianapolis flight days (Table 1) are available by request or through the Yale Stable Water Vapor Isotopes Database: https://vapor-isotope.yale.edu/(Salmon and Welp, 2019).The data have also been archived at the Purdue Uni-versity Research Repository under https://doi.org/10.4231/1MZN-1C18(Salmon et al., 2019).The authors request that they be notified if the data are to be used in publication.

Figure 1 .
Figure 1.Total uncertainty of δD, δ 18 O, and d excess over the range of H 2 O v mole fractions observed in flight.

Figure 2 .
Figure 2. Flight paths conducted around the (a) Washington, D.C. and Baltimore, MD, and (b) Indianapolis study sites for the research flights listed in Table 1.Case study flight paths are indicated with solid lines, while all other flight paths are indicated with dotted lines.

Figure 3 .
Figure 3. Measurements of (a) δD and (b) δ 18 O versus H 2 O v mole fraction made during vertical profile (VP) descents on the nine research flights (Table1).

Figure 4 .
Figure 4. Observations of d excess versus H 2 O v mole fraction in the (a) free troposphere (FT), (b) inversion (INV), and (c) boundary layer (BL) during all the VP descents.Colors indicate each research flight.BL, INV, and FT observations are defined by the rate of change of atmospheric variables.The dashed lines correspond to the global average precipitation d-excess value of 10 ‰ for reference only.

Figure 5 .
Figure 5. (a) CLR flight path overlaying the study site's cloud cover (https://worldview.earthdata.nasa.gov/,last access: 5 February 2018) captured at 12:30 LT.The teal outline indicates the Indianapolis city boundaries.The West Lafayette, IN, Purdue airport is indicated by the airplane marker.(b) CLR case study altitude time series.Colored segments indicate vertical profile (VP) descents.
), above which temperature increases within the INV layer.The H 2 O v mole fraction decreases rapidly in the INV between z INV and z FT before becoming relatively constant in the FT.The δ 18 O and δD values track the H 2 O v profile within the INV, decreasing by 30.8 ‰ and 193.2 ‰, respectively.Observed d-excess values in the INV initially decrease slightly with altitude and then increase, varying overall by 66.6 ‰.Above ∼ 1100 m in the FT, the H 2 O v , δ 18 O, δD, and dexcess signatures are relatively constant with altitude.

Figure 6 .
Figure 6.Observations of meteorological and isotope variables along the second VP (VP2) conducted on (a) CLR (MAR06) and (b) STC (MAR04).Shading around the δ 18 O, δD, and d-excess VP measurements defines measurement uncertainty.Altitudes of the boundary layer (BL), inversion layer (INV), and free troposphere (FT) are indicated for reference.The inversion layers, which are bound at the bottom and top by z INV and z FT , respectively, are identified by blue horizontal bands.The dashed blue line in (b) corresponds to stratocumulus cloud base (z CB ).

Figure 7 .
Figure 7.Comparison of vertical profile δ 18 O (a, d), δD (b, e), and d-excess (c, f) measurements to Rayleigh theory (a, b, c) and mixing (d, e, f) curves for CLR (MAR06).Individual VP descents are indicated by the differently colored points.The bounds of the inversion layer (INV), indicated by grey shading, are defined by the average H 2 O v mole fractions observed at z INV and z FT .

Figure 8 .
Figure 8.(a) STC flight path overlaying the study site's cloud cover (https://worldview.earthdata.nasa.gov/,last access: 5 February 2018) captured at 12:30 LT.The teal outline indicates the Indianapolis city boundaries.The West Lafayette, IN, Purdue airport is indicated by the airplane marker.(b) STC case study altitude time series.Colored segments indicate vertical profile (VP) descents.
also describe VP observations of a minimum in d excess (negative values) at the BL top and hypothesize that the negative d-excess signal results from rain droplet evaporation directly below cumulus clouds at the top of the boundary layer.Rehydration processes, like cloud and rain droplet evaporation, have been proposed as mechanisms that could produce negative d-excess anomalies(Bolot et al., 2013;

Figure 9 .
Figure 9.Comparison of vertical profile δ 18 O (a, d), δD (b, e), and d excess (c, f) measurements to Rayleigh theory (a, b, c) and mixing (d, e, f) curves for STC (MAR04).Individual VP descents are indicated by the differently colored points.The bounds of the inversion layer (INV), indicated by grey shading, are defined by the average H 2 O v mole fractions observed at z INV and z FT .

Figure 10
Figure 10.(a) DBL flight path overlaying the study site's cloud cover (https://worldview.earthdata.nasa.gov/,last access: 5 February 2018) captured at approximately 12:30 LT.The teal outline indicates the Indianapolis city boundaries.The West Lafayette, IN, Purdue airport is indicated by the airplane marker.(b) DBL case study altitude time series.Colored segments indicate vertical profile (VP) descents.

Figure 12 .
Figure 12.Comparison of vertical profile δ 18 O (a, d), δD (b, e), and d-excess (c, f) measurements to Rayleigh theory (a, b, c) and mixing (d, e, f) curves for DBL (MAR18).Individual VP descents are indicated by the differently colored points.The bounds of the inversion layer (INV), indicated by grey shading, are defined by the average H 2 O v mole fractions observed at z INV and z FT .
Satellites can provide discrete VP measurements of HDO v /H 2 O v but have low vertical resolution and currently cannot give information about d excess, as they cannot retrieve H 18 2 O v /H 2 O v

Table 1 .
Flight date, research flight code, flight time (local time, LT), number of vertical profiles conducted, and the observed range of potential temperature (θ ) and ambient temperature (T ) during the flights.