Isotope analysis of water by means of near-infrared dual-wavelength diode laser spectroscopy

A novel diode laser spectrometer was developed using dualwavelength multiplexing, ensuring ideal conditions for high-precision and simultaneous measurements of the H/H, O/O, and O/O isotope ratios in water. A 1.4-μm diode laser probed a HOH / HOH line pair near 7198 cm, while a similar laser observed HOH, HOH, and HOH rovibrational lines around 7183 cm, or a HOH / HOH line pair near 7200 cm. The 1-σ standard deviation is 0.2‰ for O/O, and 0.5‰ for the H/H and O/O isotope ratios. Preliminary experiments with repeated injections of a natural abundance sample point to an accuracy of about 1‰ for all three isotope ratios in natural samples. 2003 Optical Society of America OCIS codes: (300.1030) Spectroscopy: Absorption; (300.6190) Spectroscopy: Spectrometers; (300.6260) Spectroscopy: Spectroscopy, diode lasers References and links 1. E. R. Th. Kerstel: Isotope Ratio Infrared Spectrometry, in Handbook of Stable Isotope Analytical Techniques, ed. P.A. de Groot, Elsevier, accepted for publication. 2. G. Gagliardi, R. Restieri, G. De Biasio, P. De Natale, F. Cotrufo, and L. Gianfrani, “Quantitative diode laser absorption spectroscopy near 2 μm with high precision measurements of CO2 concentration,” Rev. Sci. Instr. 72, 4228-4233 (2001). 3. M. P. Arroyo and R. K. Hanson, “Absorption measurements of water-vapor concentration, temperature, and line-shape parameters using tunable InGaAsP diode laser,” Appl. Opt. 32, 6104-6116 (1993). 4. I. Linnerud, P. Kaspersen, and T. Jaeger, “Gas monitoring in the process industry using diode laser spectroscopy,” Appl. Phys. B 67, 297-305 (1998). 5. D. S. Baer, R. K. Hanson, M. E. Newfield, and N. K. J. M. Gopaul, “Multiplexed diode-laser sensor system for simultaneous H2O, O2, and temperature measurements”, Opt. Lett. 19, 1900-1902 (1994). 6. L. Gianfrani, P. De Natale, and G. De Natale, “Remote sensing of volcanic gases with a DFB-laser-based fiber spectrometer,” Appl. Phys. B 70, 467-470 (2000). 7. M. P. Arroyo, T. P. Birbeck, D. S. Baer, and R. K. Hanson, "Dual diode-Laser fiber-optic diagnostic for water vapor measurements," Opt. Lett. 19, 1091-1093 (1994). 8. G. Totschnig, D. S. Baer, J. Wang, F. Winter, H. Hofbauer, and R. K. Hanson, "Multiplexed continuouswave diode-laser cavity ringdown measurements of multiple species," Appl. Opt. 39, 2009-2016, (2000). 9. D. B. Oh, M. E. Paige, and D. S. Bomse, “Frequency modulation multiplexing for simultaneous detection of multiple gases by use of wavelength modulation spectroscopy with diode laers,” Appl. Opt. 37, 24992501 (1998). 10. E. R. Th. Kerstel, R. van Trigt, N. Dam, J. Reuss, and H. A. J. Meijer, “Simultaneous determination of the H/H, O/O, and O/O isotope abundance ratios in water by means of laser spectrometry,” Anal. Chem. 71, 5297-5303 (1999). 11. R. van Trigt, H. A. J. Meijer, A. E. Sveinbjornsdottir, S. J. Johnsen, and E. R. Th. Kerstel, ” Measuring stable isotopes of hydrogen and oxygen in ice: the Bølling transition in the Dye–3 (south Greenland) ice core,” Ann. Glaciol. 35, 125-130 (2002). 12. R. van Trigt, E. R. Th. Kerstel, R. E. M. Neubert, H. A. J. Meijer, M. McLean, and G. H. Visser, “Validation of the doubly labeled water method in japanese quail at different water fluxes,” J. Appl. Physiol. 93, 2147-2154 (2002). (C) 2003 OSA 30 June 2003 / Vol. 11, No. 13 / OPTICS EXPRESS 1566 #2463 $15.00 US Received May 07, 2003; Revised June 21, 2003 13. E. R. Th. Kerstel, G. Gagliardi, L. Gianfrani, H. A. J. Meijer, R. van Trigt, and R. Ramaker, “Determination of the H/H, O/O, and O/O isotope ratios in water by means of tunable diode laser spectroscopy at 1.39 μm,” Spectrochim. Acta 58, 2389-2396 (2002). 14. R. A. Toth, “Extensive measurements of H2 O line frequencies and strengths: 5750 to 7965 cm,” Appl. Opt. 33, 4851-4867 (1994). 15. R. A. Toth, “Transition frequencies and strengths of H2 O and H2 O: 6600 to 7640 cm,” Appl. Opt. 33, 4868-4879 (1994). 16. R. A. Toth: private communication (2000).


Introduction
Optical techniques for stable isotope analysis have attracted a growing interest in recent years, because of the need of valid alternatives to isotope ratio mass spectrometry (IRMS) [1]. In fact, although commercially available with a potentially very high precision and throughput, IRMS instrumentation is expensive, bulky, and usually requires a skilled technician for its operation. IRMS is also not compatible with condensable gases, such as water, which require chemical conversion of the sample to more easily analyzed gases. These chemical preparatory steps limit the overall throughput and often also the final precision and accuracy. In addition, they complicate the on-line, continuous analysis of samples. Water is certainly one of the unluckiest cases. Reduction of water to H 2 or equilibration with CO 2 is necessary in order to determine 2 H and 18 O isotope ratios, respectively. Isotope ratio measurements of 17 O are virtually impossible, as the 17 OCO and 13 CO 2 molecules appear in the same mass channel.
On the other hand, diode laser spectroscopy has been widely recognized as a powerful tool for trace-gas diagnostics, fast and accurate concentration measurements being possible in a large variety of environments [2][3][4]. Very interesting advances have been made in the field of multiple-gas sensing. The use of multiple independently operated diode lasers in conjunction with optical fiber components has enabled the development of multiple-gas analyzers capable of measuring temperature, pressure and concentrations in a single optical path configuration [5]. Compact and robust systems have been implemented for field measurements, often in hostile conditions like those of a volcanic site [6]. Multiple-gas sensing can be performed using different approaches. In the multiplexed configuration, one laser can be activated at a time providing a sequential operation at a given repetition frequency [7,8]. Recently, continuous and simultaneous probing of different species in a two-laser set-up has been reported, using an original scheme based on wavelength modulation spectroscopy [9]. In this case, the two lasers were modulated at incommensurate frequencies and detected by the same photodiode, whose amplitude-modulated output signal was processed by two different lock-in amplifiers. In this paper, we demonstrate that a similar spectroscopic design can be successfully applied to the precise and simultaneous determination of the 2 H/ 1 H, 17 O/ 16 O, and 18 O/ 16 O isotope abundance ratios in water. This is not an obvious issue, as the precision requirements of isotope ratio measurements are such that cross-talk between the two wavelength channels may not exceed about 1 in 10 4 . In particular for tunable diode laser isotope ratio measurements, wavelength multiplexing provides a major advantage, as most such experiments are severely hampered by the requirement of finding all istopomers of interest absorbing within the limited tuning range of the laser, at near equal intensities, and without interfering spectral overlap or unfavorable temperature dependencies.
At the Groningen Center for Isotope Research, we have already shown the capability of infrared laser spectroscopy to measure isotope ratios in small water samples in the gas phase, exciting the fundamental OH-stretching vibration at 2.7 µm [10]. The isotope ratios were determined from simultaneously recorded, direct absorption spectra in parallel sample and reference gas cells, the latter being filled with a local reference material, back-traceable to international standards such as Standard Light Antarctic Precipitation (SLAP) and Vienna Standard Mean Ocean Water (VSMOW). The instrument used a commercial color center laser as its infrared light source. The accuracy for routine measurements was of the order of 0.5‰ for the oxygen isotope ratios and of about 1‰ for the deuterium isotope ratio for natural ( samples, sufficient for successful application of the instrument to studies in paleoclimatology [11] and biomedicine [12]. More recently, we have also shown that a distributed feedback diode laser, operating near 1.39 µm, can replace the color center laser with a large number of technical advantages, including a higher spectral purity, lower frequency and amplitude noise, and much better performance in terms of frequency control and modulation, in addition to the practical advantages of low cost, small size and weight, and low power consumption. In spite of the relatively small line intensities associated with the overtone or combination vibrations at 1.4 µm, promising results were obtained for the 2 H/ 1 H, 17 O/ 16 O, and 18 O/ 16 O isotope ratios using wavelength modulation spectroscopy [13]. As in the case of the color center laser setup, the simultaneous measurements of these three ratios was possible thanks to the fortuitous occurrence of suitable ro-vibrational lines of the H 16 16 OH / H 18 OH line pair near 7200 cm -1 ), to measure all three isotope ratios (or just the two most frequently desired) with a precision of better than 1‰. Figure 1 shows the experimental setup, consisting of the diode laser system, two multiplereflection cells and the detection system. The diode lasers are of the InGaAsP/InP type, with a distributed feedback (DFB) structure and were bought individually at different times from Sensors Unlimited. Their specifications differ significantly, most prominently in beam quality and output power. The first (DL 1 ) emits about 1 mW in a single mode around 1.392 µm, while the second laser (DL 2 ) can deliver up to 10 mW at a wavelength of 1.389 µm. Both lasers operate at or near room temperature. Homebuilt current generators and (thermo-electric) temperature controllers were used to drive the diode laser sources. In the wavelength multiplexed experiments DL 1 was wavelength tuned over vib-rotational lines of H 16 OH, H 17 OH, and H 18 OH isotopomers, with a periodic modulation of the injection current providing a 0.5 cm -1 wide scan at a repetition rate of 0.125 Hz. At the same time, DL 2 was tuned into coincidence with H 16 OH and HO 2 H lines by means of a frequency scan as large as 0.6 cm -1 , which was performed at the same repetition rate and in phase with that of DL 1 . For this purpose, a single triangular wave was used with an individually adjustable gain on each modulation channel.

Experimental setup
The two laser beams were combined through a pair of half-wave plates and a polarizing cube, thus producing a two-color beam with crossed polarization components. This beam was split into two parts of equal intensities, which were subsequently sent to the two gas cells. The optical path-length inside the multiple-pass cells was 20.5 m in a 48-pass configuration. The use of a pair of cells enables one to measure the isotope ratio in a sample gas with respect to that of a reference gas and, hence, to determine the δ-values in a manner as described previously [10]. A vacuum system, based on an oil-free molecular-drag pump, ensured high purity in the gas cells. The liquid water samples were injected into the gas cell through a silicon membrane using a 10 µl syringe. Simultaneous operation of the two lasers was possible using frequency modulation multiplexing. The diode lasers were modulated at incommensurate frequencies of approximately 18 and 27 kHz. The modulation frequencies were chosen with particular care in order to optimize the signal-to-noise ratio and to avoid potential cross-talk effects among the absorption signals at the two wavelengths. This was accomplished using a Fourier transform spectrum analyzer, which enabled us to compare the amplitude noise spectra in the range 0-100 kHz associated with the transmission of one of the gas cells with one laser activated at a time and, subsequently, with both lasers operating simultaneously. A Ge photodetector was used to monitor the transmitted power for each cell. Two pairs of lock-in amplifiers were used to process the signals carrying information on the ( Wavelength modulation spectroscopy with first harmonic (1-f) detection was implemented with the modulation index set to about 2. First-harmonic detection was selected over secondharmonic detection, based on preliminary experiments. It was concluded that the reduced amplitude modulation in a 2-f scheme does not compensate for the loss of signal strength. The lock-in time constants were equal to 10 ms with digital filtering of 24 dB/oct for the digital lock-ins, and 30 ms with 6 dB/oct, for the analog models. Signal averaging over 10 scans was performed to further reduce the equivalent noise detection bandwidth (ENDB) down to about 1 Hz in each case.
As shown schematically in Fig. 1, the digital lock-in amplifiers were interfaced to a personal computer through a GPIB interface, while auxiliary analog-to-digital converters of the digital lock-ins were used to acquire the output signals of the analog lock-in amplifiers with 16-bit resolution. In this way, four data vectors (2048-record length) were initially stored in the buffers of the two digital lock-in amplifiers before being transferred to the personal computer, thus enabling the fully automated acquisition of spectra for many hours without operator intervention.
The alignment of the multiple-reflection cells was performed with utmost care in order to minimize interference effects, which could introduce spurious amplitude modulated signals, thus disturbing the absorption signals. No periodic modulation of the baseline, due to "fringes", was observable in the present setup. In addition, great care was taken to insure mechanical and thermal stability. Particularly, the entire spectrometer was enclosed in a sealed plexi-glas box, which was shown to be already quite effective in reducing thermal variations between the gas cells. A second box of thermally insulating foam provided further passive thermal stabilization of the gas cells. This is a crucial point, as a temperature difference of 1 °C between the sample and reference cells may easily introduce a systematic error in the relative isotope ratio values of the order of 1‰, depending on the temperature coefficients of the ro-vibrational lines. The temperature coefficient is defined as the relative variation of the line strength corresponding to a temperature increase of 1 K at 296 K, and may be expressed in ‰/K. Table 1 gives the temperature coefficients of the ro-vibrational transitions used in this study, together with their ro-vibrational assignments based on the spectroscopic data reported by R. A. Toth [14][15][16]. The temperature coefficients were calculated from the ground state energies E" using the following formula for the temperature dependence of the line strength S at constant volume: This formula neglects the effects of stimulated emission (a very good approximation in the infrared), and was derived assuming a classical expression for the rotational partition function, while ignoring the temperature dependence of the vibrational partition function (see, e.g., formulas 5 and 6 in Ref. [1]).
The temperature sensitivity of the isotope ratio measurement (δ-value; see below) is in very good approximation equal to the difference in temperature coefficients of the rovibrational transitions involved (i.e., of rare and abundant isotopomer).

Results and discussion
An example of acquired spectra is shown in Fig. 2.  We may note the excellent signal-to-noise ratio, which was measured to be about 10 4 for the H 16 OH line at 7197.5505 cm -1 and 3 10 3 for the H 18 OH line at 7183.5871 cm -1 , with an ENBW of 1.2 Hz and 0.8 Hz, respectively. Using the so-called δ-notation, the isotopic composition of a sample with respect to that of a reference material is given by, for instance in the case of 18 O in water: n 18 and n 16 being the H 18 OH and H 16 OH number densities. Since the δ-value is generally very small, it is usually expressed as per mil (‰). The usual international standard reference material for water is VSMOW, although in practice a local reference material is used as the working standard, which is well-known with respect to VSMOW. The water vapor pressure was about 10 Torr, which lead to a maximum fractional absorption at line center of less than 5% for the strongest lines, due to H 18 OH and H 16 OH. At this value of absorption, the deviation from linearity of the Lambert-Beer law is about 2.6%. Since for all practical purposes the abundance of H 16 OH does not change, we need only worry about underestimating the absorption of the rare isotopomer. However, only at relatively large isotopic abundance differences between sample and reference material, will this lead to a sizable error in the determined isotope δ-value of Eq. (2). For example, if one measures a sample with δ 18 Ο = -50‰ (this would be about the most depleted natural water to be found on Earth on Antarctica), the systematic error on this account would be about 1.3‰. In our case of samples of identical or very similar isotopic composition, the maximum absorption of 5% is sufficiently low to ensure that we may safely ignore the effect of Lambert-Beer non-linearity. The δ-values are retrieved by means of a line-by-line non-linear least squares fit of the sample spectrum to the sum of the reference spectrum and a linear or quadratic baseline, multiplied by a factor representing the ratio of the isotopomer concentrations, in a manner as described in Ref. [10]. This procedure requires a careful determination of the line centers and is performed, for each isotopomer, over a frequency range as wide as 3 times the full-width-at-halfmaximum, symmetric with respect to the line center. The analytical capabilities of the dual-wavelength spectrometer were carefully tested and optimized through a large number of measurements performed with initially one laser activated at a time and, subsequently, with both lasers operating simultaneously. For this purpose, both gas cells were filled with 10 µl of identical de-mineralized water. We thus expect to measure a δ-value of zero. The most interesting performances were obtained with the δ 2 H determination, using DL 2 . In Fig. 3, the results of precision tests, over short (upper part) and long (lower part) time spans, are shown.
As each measurement took 160 s, the total acquisition times were equal to 80 and 400 minutes for series of 30 and 150 measurements. Of these, measurements that were removed more than 3.5 times the absolute deviation from the median were excluded from the figure and further statistical analysis. In this way, in a typical experiment between 5% and 10% of the measurements were considered to be outliers. For the measurements of Fig. 3 we found a 1-σ precision of 0.5‰ and 0.8‰ for the 80 and 400 minute time-span series, respectively. It is interesting to note that the δ-values appear to change periodically between -0.5‰ and +2.5‰. This may be attributed to a small temperature difference between the two cells, the temperature coefficients of the HO 2 H and H 16  We also checked the reproducibility over several days through different series of measurements, each time filling both cells with 10 µl of the same de-mineralized water, followed by evacuation, flushing, and refilling, before starting a new series. Over a time interval of six days, we found mean values ranging from -0.73‰ to +1.68‰, with an overall average value of 0.48‰ and a standard deviation (1-σ) of 1‰. This test also provides preliminary information on the level of accuracy. Through the measurement of two or more reference materials of well-known isotopic composition, the reproducibility may be transferred into a measurement accuracy over the range spanned by the reference materials. Although a non-zero average δ-value indicates the presence of a systematic deviation, the average δ 2 H-value was well within the long-term precision level. The results of the long-term reproducibility test are presented in Table 2.  Similar tests were performed for the detection of H 16 OH, H 17 OH, and H 18 OH isotopomers by means of DL 1 . Recording several series of 30 spectra we found values of the short-term reproducibility in the range 0.6‰ to 1.1‰ and 1.4‰ to 1.8‰, for the δ 18 O and δ 17 O determinations, respectively. Hence, in spite of the much lower available power with respect to DL 2 , we were able to obtain only slightly inferior results in terms of short-term precision. Nonetheless, the major problem with δ 17 O and δ 18 O measurements was related to the accuracy, as systematic errors emerged in each series. In Table 3, we report the results of δ 17 O and δ 18 O determinations in different sets of measurements performed over two days, each time using a new sample of de-mineralized water. We may note different apparent δ-shifts, with the δ 17 O-and δ 18 O-values always exhibiting opposite signs. This point deserves a more detailed discussion. We believe that these apparent δ-shifts are due to residual temperature and pressure differences between the two cells, as well as to the partial overlapping of the H 16 OH and H 17 OH lines at 7183 cm -1 . Particularly, using the temperature coefficients reported in Table 1, we may calculate shifts in the δ 17 O and δ 18 O values of about +9 and -1 ‰/K, due to the differential temperature effect alone. In other words, the observed fluctuations in the δ 17 O value would suggest a temperature difference between the two cells of 0.7 K. But this temperature difference should lead to much smaller variations in the δ 18 O-value than those observed experimentally, and reported in Table 3. However, the spectral overlapping of the H 16 OH and H 17 OH lines, clearly visible in Fig. 2, may enhance the effect of a differential temperature on the δ 18 O-value because of the opposite sign of the H 16 OH and H 17 OH temperature coefficients. Indeed, as a consequence of the partial line overlapping, the effective temperature coefficient for the H 16 OH signal would be lower than that given in Table 1, thus increasing the temperature dependence of the δ 18 O-value. Furthermore, a differential pressure may limit the accuracy in the fit of the sample spectrum to the reference one, for each isotopomer individually, because of the different line shapes and widths. Line overlapping may be responsible for pressure-dependent asymmetries in the 1-f profiles, further limiting the accuracy in the spectra analysis. It should be noted that these effects are a direct consequence of the data analysis procedure used, in which a section around each spectral feature in the sample spectrum is fit to the same region in the reference spectrum (with the isotopomer ratio as a multiplicative free parameter). An analysis procedure based on the simultaneous fit of multiple (frequency modulated) absorption profiles to all spectral features present in the recorded spectrum would not suffer from such a differential pressure (line width and line shape) effect. This is because in this case the absorption intensities may be calculated from the line height and line width parameters of the individual line profiles. In the calculation of the δ-value according to Eq. (2), the pressure induced line shape change of the spectral feature corresponding to the rare isotopomer would then be canceled by the same change in the abundant isotopomer line shape. A disadvantage of a simultaneous fit to all spectral features present in the recorded spectrum is the requirement of a proper calibration of the frequency scale. We have observed and extensively studied the above mentioned effects with our color center laser spectrometer at 2.7 µm [10]. We concluded that in most cases the, line-by-line, sample to reference fitting procedure is superior to the multiple line (Voigt) profile procedure, provided that well-calibrated (thus reproducible) corrections are made for the differential pressure between the sample and reference gas cells (as measured by the line-widths). In fact, we were even able to reproduce the experimentally determined correction factors in a simulation of the data analysis procedure with synthetic spectra. These aspects of the present apparatus are currently being studied in more detail. The analytical performance of our spectrometer was also tested with both lasers operating simultaneously in the wavelength multiplexed mode. We found the simultaneous operation to be effective, as the precision and reproducibility (upon repeated injection of a fresh sample) levels achieved with each laser separately were not measurably affected by the presence of the other laser.
We also used the better, more powerful diode laser DL 2 , instead of DL 1 , to determine δ 18 O and δ 17 O. For this purpose, DL 2 was temperature tuned to 7183 cm -1 to probe the H 16 Fig. 4(a) and Fig. 4(b). This improvement is attributed to the better technical features of DL 2 with respect to DL 1 , in terms of emitted power and spatial beam quality. Also, it must be noted that the gain of the pre-amplified detectors was optimized on the transmission of DL 2 , thus resulting in a less than optimal detection of the DL 1 beam. Also in this configuration systematic errors similar to those discussed above for the experiment with DL 1 were clearly present.
Subsequently, we tuned DL 2 in coincidence with a pair of H 16 OH and H 18 OH lines at 7200 cm -1 , which exhibits the important features of having a low temperature dependence and being well resolved and isolated from other lines. Here, we found an interesting result, which is shown in Fig. 4(c). The mean δ 18 O value was 0.1‰, with a standard deviation of 0.2‰. The much better performance with respect to the oxygen isotope ratio measurements at 7183 cm -1 , using the same laser, confirms our hypothesis of line overlap as a major source of systematic errors in the δ 17 O and δ 18 O measurements at 7183 cm -1 .

Conclusions
We have demonstrated that dual-wavelength diode laser spectroscopy enables the simultaneous determination of 2  Precision levels are better that those achieved in previous experiments with a color center laser at 2.7 µm. Herewith, laser spectrometry becomes a viable alternative to conventional IRMS, particularly for deuterium, and considering the laser spectrometry advantages of fast and direct measurements on water samples, no sample pretreatment, and a small and relatively cheap instrument. Preliminary tests of the reproducibility upon injection of fresh samples, performed in different spectral regions, revealed that δ 18 O and δ 2 H measurements can be carried out with systematic errors smaller than 1‰, in the spectral windows around 7200 and 7198 cm -1 , respectively. On the other hand, δ 17 O and δ 18 O determinations performed at 7183 cm -1 require a correction of the apparent δ-values based on the measured line width difference between sample and reference gas, or an altogether different data analysis procedure, in order to arrive at the same level of accuracy, due to the partial overlapping of H 16 OH and H 17 OH lines (which enhances differential temperature and pressure effects in our current data analysis strategy). More measurements on samples, that cover a range of isotopic abundances, including standard materials back-traceable to international standards, are required in order to calibrate the instrument and truly assess its accuracy.