Miniaturized laser heterodyne radiometer for measurements of CO 2 in the atmospheric column

We have developed a low-cost, miniaturized laser heterodyne radiometer for highly sensitive measurements of carbon dioxide (CO 2 ) in the atmospheric column. In this passive design, sunlight that has undergone absorption by CO 2 in the atmosphere is collected and mixed with continuous wave laser light that is step-scanned across the absorption feature centered at 1,573.6 nm. The resulting radio frequency beat signal is collected as a function of laser wavelength, from which the total column mole fraction can be de-convolved. We are expanding this technique to include methane (CH 4 ) and carbon monoxide (CO), and with minor modiﬁcations, this technique can be expanded to include species such as water vapor (H 2 O) and nitrous oxide (N 2 O).


Introduction
Surface-observing networks have made major contributions to understand the carbon cycle as well as atmospheric distributions of key gases such as CO 2 , CH 4 , and CO. However, these observations have not been sufficient to constrain regional CO 2 fluxes, and many parameters related to the emission of these gases still vary widely. Significant uncertainty remains in the source and sink distributions of these gases as well as the processes that govern them and their interannual variability [1,2]. Fluxes inferred from surface measurements alone are subject to errors in atmospheric transport models [3,4]. However, vertical column measurements are less affected by model transport errors, making them particularly useful in constraining surface fluxes [5]. Rayner and O'Brien [6] demonstrated that monthly averaged column CO 2 from satellites with precisions of 2.5 ppmv on a 8 9 10°scale could provide a constraint on surface fluxes comparable to the existent surface network. A more recent study by Chevallier et al. [7] used data from 14 TCCON [8] stations to estimate fluxes and compared the results with estimates based on surface and aircraft observations. Their results indicated that TCCON-based flux estimates were consistent with northern hemisphere seasonal cycles and regional carbon budgets based on flux estimates using the in situ data sources. Though the TCCON results were less precise, in part due to the very sparse nature of the network, this demonstrated for the first time the utility of CO 2 column measurements in constraining carbon fluxes. The Chevallier et al. results demonstrate the scientific value of a high-precision, groundbased, column observation network.
Here, we present a miniaturized laser heterodyne radiometer (mini-LHR) that measures CO 2 in the atmospheric column and offers a low-cost solution for filling in gaps in the current observational network. Miniaturization was possible through recent developments in the telecommunications industry that resulted in low-cost, commercially available, distributive feedback (DFB) lasers that produce stable near infrared light. The mini-LHR has been designed to operate in tandem with an AERONET sun tracker [9][10][11]. With more than 450 instrument sites worldwide, AERONET offers an established framework for future expansion of the mini-LHR and includes key arctic regions (not covered by GOSAT or the planned OCO-2 mission) [12] where release of CO 2 and CH 4 from thawing tundra and permafrost is a growing concern. An additional benefit to tandem operation with AERONET is the quantification of cloud and aerosol light scattering through simultaneous aerosol measurements. In their studies of sensitivity of a space-based measurement, Mao and Kawa [13] find that aerosol and cloud data are key to reduce errors in the CO 2 column because scattering by cloud and aerosol particles changes the pathlength of sunlight through the column and consequently changes the total column CO 2 absorption. A network of continuous, colocated measurements of CO 2 and aerosol abundance solves this pathlength issue and provides the opportunity to better understand the processes controlling carbon flux. Aerosols induce a radiative effect that is an important modulator of regional carbon cycles. Changes in the diffuse radiative flux fraction (DRF) due to aerosol loading have the potential to alter the terrestrial carbon exchange [6,14].

Background of laser heterodyne radiometry
Since the 1970s, laser heterodyne radiometry has become a well established receiver technique [15][16][17][18][19] and has been used to measure atmospheric gases such as ozone (O 3 ) [20,21], water vapor (H 2 O) [22], methane (CH 4 ) [22], ammonia (NH 3 ) [23], and chlorine monoxide (ClO-an atmospheric free radical) [24], with efforts for additional species such as nitrous oxide (N 2 O) underway. We present a passive variation of the laser heterodyne radiometer that uses sunlight as the light source for absorption of CO 2 in the infrared. Sunlight is superimposed with laser light in a single mode fiber coupler. The signals are mixed in a fast photoreceiver, and the RF beat signal is extracted. Changes in concentration of the trace gas are realized through analyzing changes in the beat frequency amplitude. Miniaturization is possible through the use of smaller DFB lasers and related fiber optic components that have recently become commercially available and inexpensive through progress in the telecommunications industry [9].

Atmospheric measurement approach
With the exception of satellite measurements such as GOSAT [25] and the upcoming OCO-2 mission [12,[26][27][28], there are two general approaches that are in regular use for measuring carbon cycle gases in the atmosphere: (1) in situ instruments that sample gases from a specific location and then analyze samples with an analytical instrument such as an infrared gas analyzer and (2) suntracking instruments that measure gases in the atmospheric column using sunlight as the light source for absorption measurements. FLUXNET [29] and NOAA's ESRL/GMD Surface flask sampling network [30] support the first type of measurement. These approaches offer global coverage, but because sampling occurs near the surface (and not in the entire atmospheric column), these measurements are typically not adequate for validation and calibration of satellite data. Other examples of the first type are NOAA's tall tower [31] and CCGG aircraft measurements [32]. These retrieve some altitude information-but are limited by the availability of tall tower sites and flight sampling locations. The second type of instrument measures gases in the atmospheric column and consequently provides a better comparison with satellite retrievals where trace gases are also measured in the column. The main network of these ground instruments, total carbon column-observing network (TCCON) [8], is comprised of 16 operational global instrument sites. TCCON data are already being used for validation of GOSAT [25] data and in the future, will be used for OCO-2 [33] validation. These Fouriertransform spectrometers (FTS) can measure the largest range of trace gases, but expansion of the network is limited due to cost.
The mini-LHR is in the second category of instruments and offers a low-cost solution for filling in gaps in the TCCON network with targeted measurements of key carbon cycle gases. The passive mini-LHR has been designed to operate in tandem with an AERONET sun photometer. With more than 450 instrument sites worldwide, AER-ONET offers an established framework for future expansion of the mini-LHR. While only CO 2 measurements are presented here, the mini-LHR can be expanded to include methane (CH 4 ) and carbon monoxide (CO) measurements. Also, because the mini-LHR is operated in tandem with AERONET, it will be possible to identify interferences from cloud and aerosol scattering through simultaneous aerosol measurements.

Sun tracking and collimation
The mini-LHR uses the sun tracker technology of an AERONET sun photometer [9-11] for collection of sunlight in the atmospheric column. Collimation optics for the mini-LHR is noninvasively ''piggy-backed'' to the sun tracker through a custom-made lightweight band of aluminum (shown in Figs. 1, 3). While the mini-LHR collimation optics can be connected to a variety of different sun trackers, the AERONET sun tracker has been found to be reliable and the link to the AERONET network provides an existing global framework for future deployments of the mini-LHR. AERONET sun-tracking instruments are commercially produced by CIMEL Electronique (model CE 318, CIMEL.FR). A typical AERONET installation consists of three components: an optical head, an electronics box, and a robot. The optical head consists of two collimators (with and without lenses). The electronics box runs the motion control as well as the data collection through separate processors. The robot has two stepper motors that control the zenith and azimuth of the tracker which are used in combination with a 4-quadrant detector for a tracking accuracy of better than 0.1°.
For AERONET data collection, these sun trackers operate every 15 min: returning to a ''rest'' position between measurements. For our experiments, we operated the sun tracker continuously with an update interval of 30 s. A single mode fiber optic cable brings light from the collimator into the mobile instrument. The field of view (FOV) for each of the collimators was *0.2°( compared to the 1.2°FOV of AERONET and the *0.5°F OV of the sun). The tracking accuracy of AERONET is 0.1°and results in a minor oscillation that can be seen in the baseline of the absorption scans. To mitigate this, a second LHR collimator serves as a reference channel to monitor the signal from the sun in this wavelength region. This channel is sensitive to the accuracy of the tracker as well as drops in the signal due to cloud cover or aerosols.

Miniaturized laser heterodyne radiometer (mini-LHR)
In this passive instrument, sunlight that has undergone absorption by CO 2 is collected and mixed with laser light that scans across a CO 2 absorption feature at 1,573.6 nm. This line was selected at a wavelength region that contains minimal interferences from other species such as water vapor, but also that was within the C-L bands for commercially available DFB lasers. A schematic of the instrument set-up is shown in Fig. 1. Incoming light is modulated with an optical chopper and introduced into the Fig. 1 Schematic of the mini-LHR. Sunlight that has undergone absorption by CO 2 is collected with fiber-coupled collimation optics that is connected to an AERONET sun tracker. This incoming sunlight is chopped and then superimposed with light from the DFB laser (local oscillator) in a single mode fiber coupler. Superimposed light is mixed in an InGaAs detector (fast photoreceiver) to produce the RF beat signal. The RF receiver amplifies and detects the beat signal that is then measured through a data acquisition device that is synchronized with the chopper frequency. To measure the CO 2 absorption feature, the laser is scanned over the appropriate IR wavelength, and changes in the resulting RF beat signal are monitored Fig. 2 The compact RF receiver amplifies and detects the beat signal that is produced in the fast photoreceiver. Components from upper right include 1 bias tee with 50 ohm resistor, 2 high-gain, low-noise amplifiers, 3 square-law detector, and 4 video amplifier mini-LHR through a single mode optical fiber. Light from the local oscillator (a DFB laser) is superimposed on this incoming sunlight in a 50/50 single mode fiber coupler. Superimposed light is mixed in a fast photoreceiver (a 5 GHz InGaAs detector) to produce an RF beat signal. The beat signal enters the RF receiver/amplification stage (shown in Figs. 1, 2) through a bias tee with a 50 ohm resistor (the required InGaAs detector DC impedance) to separate RF and DC outputs of the InGaAs detector. The RF signal passes through a gain stage to amplify and set the bandwidth of the measurement (twice the RF bandwidth due to the double-sideband property of the optical-RF mixing). RF detection is performed with a square-law detector that outputs a voltage that is proportional to the square of the input voltage, allowing the measurement of the power of the RF signal. Output from the RF detector is amplified and low-pass filtered with a video amplifier circuit to produce a DC signal equal to the averaged power in the RF circuit over a chosen time interval. The final output voltage from the video amplifier is a chopped signal with the power of the input signal. When the optical chopper is blocking the incoming light, the RF subsystem output is proportional to the system noise. The amplitude of this chopped signal is monitored as the laser is scanned through the absorption feature and provides an output proportional to the incoming light around the laser wavelength with a bandwidth equal to twice the RF bandwidth. The instrument bandwidth/resolution was measured by mixing light from duplicate DFB lasers and measuring the emergence of the beat signal. The first laser was set to a constant wavelength and the second laser was scanned across the first wavelength while monitoring the beat signal after the RF subsystem. We found the bandwidth of the system to be *1.3 GHz (*0.011 nm at 1,573 nm).

Data acquisition
The small changes in the modulated signal amplitude due to perturbations in the CO 2 content within the planetary boundary layer are easily captured with a commercially available data acquisition cards, DAQ, of modest speed but high digitization fidelity. We use a National Instruments NI-USB-6255 16 bit digitizer operating at 20Ks/s to capture the variation in signal voltage from the two wide-band RF detection channels of the mini-LHR. The measured signals are a modulated, near square-wave voltage of a few 10-100 s of mV at 200 Hz with a 50/50 duty cycle. The DAQ is electronically triggered from the mechanical chopper, an SR540 chopper with negligible phase drift. The captured signals are process in real time using custom software written in LabVIEW running on a Windows laptop. The custom software essentially performs a boxcar operation to measure the amplitude of the modulation. The amplitude of the signal and reference are then recorded as a function of wavelength for later analyses.

Measurements of CO 2 in the atmospheric column
The mini-LHR is being tested at NASA's Goddard Space Flight Center (38°59 0 33 00 N, 76°50 0 23 00 W). The LHR is located on the roof of a building, allowing for a clear view of the sky throughout the entirety of the day.
The LHR instrument is comprised of an AERONET sun tracker and a weatherproof module (shown in Fig. 3) that houses the laser, the chopper, the InGaAs detectors, and the electronic components of the instrument. The module was not originally thermally controlled with the exception of some ventilation. Because summer roof temperatures exceed the upper operating range of the laser and detectors (operating range for both is 0-40°C), the module was moved to a climate-controlled shipping container that is held at 25°C. Future versions of the module will include both thermal control and internal temperature monitoring to minimize temperature-related wavelength uncertainty. Two collimators are attached by a mount fitted to the AERONET sun tracker (one for the collection of the heterodyne signal and other as a reference). The mounts have 4°of adjustment in both pitch and yaw axes. The collimators are aligned to collect the maximum amount of sunlight (generally 2-3 lW) in the morning, and they Fig. 3 The mini-LHR (shown housed in a weather-proof case) collects sunlight with fiber-coupled collimation optics connected to an AERONET sun tracker (shown at top of photo) maintain good alignment throughout the day as the AER-ONET device tracks the sun. Sunlight is brought to the laser/detector module via FC/APC single mode fibers. Figure 4 shows a typical measurement of CO 2 in the atmospheric column. The laser is scanned in 0.005 nm increments across the wavelength region of the absorption feature. Two-hundred data points are collected and averaged for an integration time of one second per increment. Raw data was fit with a simulated atmospheric CO 2 absorption feature, the details of which will be described in the following section.

Data handling and analysis
Raw data that have been collected from the LHR sensor are analyzed with a retrieval algorithm to (1) remove data outliers due to cloud and aerosol interferences, (2) remove data collected when the sun tracker was not able to lock onto the sun (due to weather or tracker maintenance issues), (3) fit the RF scans of the trace gas absorption features through a spectral simulation program (described below), and (4) extract mole fractions of the trace gases in the atmospheric column as a function of time of day.
For more than 20 years, a spectral simulation program, SpecSyn, has been under development at GWU [34][35][36]. For this project, SpecSyn is coupled with well known algorithms developed for atmospheric radiation transmission. For the latter, we borrow heavily from coding developed under the LOWTRAN and MODTRAN programs by AFOSR (and others) [37]. The interface between these two codes is detailed below.
In the case of light transmission through the atmosphere, it is necessary to account for variation of pressure, temperature, composition, and refractive index through the atmosphere that are all functions of latitude, longitude, time of day, altitude, etc. The calculation of air mass is a well-studied problem, and sophisticated algorithms are available that can follow a light path through the atmosphere.
In MODTRAN (and most other air-mass calculations), the atmosphere is modeled as a series of spherically symmetric shells with boundaries specified at defined altitudes. Species mixing ratios are defined at the boundaries and are assumed to be invariant with each layer. Temperature and pressure are allowed to vary continuously with altitude as described below. In the case of the LHR instrument, the light path between an observer at the earth's surface and the sun passes through these layers and may be slightly curved due to variation in refractive index (See Fig. 5). The observed spectrum at the LHR instrument is the convolution of the contributions along this light path.
To quantify the composition for different layers, we implement a multilinear regression algorithm. For any absorption measurement the signal at a particular spectral frequency is a linear combination of spectral line contributions from several species. For each species that might absorb in a spectral region, we have precalculated its contribution as a function of temperature and pressure. The integrated path absorption spectrum can then by calculated using the initial sun angle (from location, date, and time) and assumptions about pressure and temperature profiles from an atmospheric model (as detailed in the MODTRAN documentation). Initial guesses for compositions of analytes and confounders at several atmospheric layer boundaries will be made from one of the six model atmospheres provided by MODTRAN. We will iterate the modeled spectrum to match the experimental observation using standard multilinear regression techniques. In addition to the layer concentrations, the numerical technique also provides uncertainty estimates for these quantities as well as dependencies on assumptions inherent in the atmospheric models.
For this initial proof-of-principle work, several simplifications are made in the light path definitions. First, refraction is ignored and the paths are assumed to be linear. This assumption is reasonable as most of the measurements reported here were made mid day (mostly in the 10:00 a.m.-2:00 p.m. window) in summer months, and thus with relatively small zenith angles (refraction is most important at larger zenith angles [37]. For this work, estimates of the variation of temperature and pressure with altitude were extracted from the 1976 Standard US atmosphere (referenced as Model 6 in the Modtran2/3 Report [37]). The following empirical results were used in the data fits: ln pðtorrÞ ð Þ¼À1:49889 Â 10 À1 ZðkmÞ þ 6:73791 ð1Þ For concentration determinations, 4 layers were considered in the atmosphere, defined at lower altitude boundaries of 0, 1, 5, and 10 km.
The resolution of the spectra collected using the mini-LHR is determined by the width of the RF filter spectral width, which here is less than line widths of the spectral features in the atmosphere. Thus, the variation of absorption with wavelength the Beer-Bougert-Lambert Law gives the magnitude of the absorption: where S is the line strength of the absorption, g(m) is a line shape function, q is the gas density, v j is the mole fraction of the absorber, l is the pathlength though the absorber, and a is the extinction coefficient. In this case, transmitted near infrared light intensity is proportional to the measured RF power. The intensity of light in the absence of absorption can be determined by directly measuring solar intensity in a spectral absorption line, or by fitting a baseline from adjacent, nonabsorbing wavelengths. In either case, measured absorbance, defined in Eq. 3 is compared to that calculated using SpecSyn. Regression on the concentrations of analytes is then used to minimize differences between measured and modeled signals. A Monte Carlo simulation was performed to assess sensitivity of reported concentrations to random noise in absorbance data. The fitting program provides an option for fitting one or multiple concentrations of CO 2 in the atmospheric column. Not surprisingly, when multiple CO 2 concentrations were fit, a given noise in the spectrum produced a greater uncertainty in CO 2 concentrations than if a uniform atmospheric column concentrations is assumed.

Sensitivity estimate
Signal to noise ratio for the instrument was calculated from using the following equation, where T BB is the temperature of the black body, g e is the effective quantum efficiency of the InGaAs photodiode, B is the bandwidth of the intermediate frequency (or IF), m is the frequency, h is Planck's constant, j is Boltzmann's constant, s is integration time and T 0 is the total transmission factor [18]. Total transmission factor was estimated by combining expected optical loss. The optical loss in the system comes from coupling efficiency of the collimators, insertion loss of connecting fibers, chopping of the signal, phase misalignment, and polarization mismatch between the source and the local oscillator. Combining the sensitivity of the photodiode and the total transmission gave a total degeneracy factor of *70 (1/g e T 0 ). For the instrument the SNR expanded with an increasing integration time according to Fig. 6. For an integration time of one second, an estimated SNR of 477 was calculated. To measure the actual SNR, the data Fig. 5 Illustration of data treatment algorithm. For each laser wavelength, a spectrum is collected at a sensor location and sun angle, time, and date are recorded. A standard atmospheric model is used as input to a SpecSyn calculation that calculates a response matrix with each species at each layer as the columns and spectral frequencies as the rows. Multilinear regression is then used to determine concentrations (e.g., [CO 2 ] in layer 0, [CO 2 ] in layer 1, etc.) from August 16, 2012 was utilized. This day was chosen because it was a particularly clear day in which consistently uninterrupted data was acquired. The baseline signal and uncertainty for each spectral sweep was calculated by fitting a Lorentzian line to the data. These values were then weighted together to give a final mean of baseline signal and uncertainty. The ratio of this output resulted in a measured SNR of 342. The difference in the calculated SNR and the measured SNR is most likely due to an under estimation of the effects of the beam filling factor, mixerpreamplifier noise, and lineshape distribution of signal.
Long-term stability of the components of the LHR instrument is a minor factor in the source of error. The local oscillator used in this setup is rated by the manufacture at having a power stability of 0.01 dB, and a wavelength stability of 0.002 nm over a period of 24 h continuous use. The battery in the photodiode is rated for reliable operation for 190 mA h [38]. At the average power of the local oscillator this rating returns *20 h of reliable battery operation. This time is taken into account and batteries are frequently replaced for the photodiode.
The instrument sensitivity was estimated through Monte Carlo simulations that calculated the CO 2 concentration given varying levels of noise. Figure 7 shows the standard deviation in the calculated CO 2 as a function of absorbance noise. In this simulation, the concentration of CO 2 in the atmosphere was treated as a constant at 390 ppmv, and the zenith angle was assumed to be 25°. The noise intervals tested were 0.001, 0.002, 0.005, 0.01, 0.02, 0.05, and 0.1 which correspond to 0.1 % absorption, 0.2 % absorption, 0.5 % absorption, etc. The data was fitted with a power regression line. As shown in the fit, to achieve a sensitivity of 1 ppmv, noise in the measurement cannot exceed 1 %. Currently the measured instrument noise is approximately 10 % for 1 s of averaging of each data point corresponding to a sensitivity of *8 ppmv. Future improvements to the data acquisition speed and transitioning from laser step scanning to a modulated sweep scanning are anticipated to reduce this noise to the \1 % requirement.

Conclusions
We have presented a passive mini-LHR and preliminary measurements of carbon dioxide in the atmospheric column at 1,573.6 nm. The mini-LHR has been designed to operate in tandem with the passive aerosol sensor currently used in AERONET. Through leveraging of AERONET's existing global aerosol monitoring network, the mini-LHR has a clear pathway to future deployment in key locations that reduce flux uncertainty. Tandem operation will offer correction of cloud and aerosol light scattering through the simultaneous aerosol measurements. We report a precision of *8 ppmv for 1 s of averaging for column CO 2 measurements in cloud-free conditions with the current instrument configuration. Future changes in the data acquisition speed and laser scanning method are expected to improve this to a better than 1 ppmv. Fig. 6 Calculated signal to noise ratio for the mini-LHR for integration times ranging from 0 to 100 s Fig. 7 Simulations indicate that to achieve a 1 ppmv sensitivity, the noise in the measurement can not exceed 1 %