Lightweight mid ‐ infrared methane sensor for unmanned aerial systems

Methane is a strong greenhouse gas and plays a key role in determining the climate impact of sources such as wetlands, livestock, and oil and gas equipment [1–3]. There is a continuing need for technology to better understand and track the sources and sinks of methane at the local scale. Unmanned aerial systems (UAS) are attractive platforms for such measurements due to their flexibility and autonomous nature. In particular, small (<25 kg) UAS have low operation costs, relaxed regulatory requirements, and can operate near the surface, where most emissions originate. As an example, a near-infrared sensor (100 ppbv precision) equipped on a remote-controlled helicopter has been used to estimate the magnitude and temporal variability of methane emissions from a natural gas compressor station [4], and serves as an alternative to mobile laboratory and ground survey techniques which are more time-consuming or may miss elevated leak sources. To leverage the potential of UAS, there is a need for high performance trace gas sensors meeting the payload constraints of typical fixed wing and copter platforms, including for methane [5]. Due in part to this gap, air quality and trace gas monitoring with UAS is still a relatively new area of research. Methane instruments are available that have excellent performance, but are generally designed for aircraft[6] or ground-based applications [7] where size, power consumption and mass are not significantly restricted. A sensor based on integrated cavity output spectroscopy (ICOS) technology from Los Gatos Research Abstract The design and field performance of a compact diode laser-based instrument for measuring methane on unmanned aerial systems (UAS) is described. The system is based on open-path, wavelength modulation spectroscopy with a 3.27 μm GaSb laser. We design two versions of the sensor for a long-endurance fixed wing UAS and a rotary wing hexacopter, with instrument masses of 4.6 and 1.6 kg, respectively. The long-endurance platform was used to measure vertical profiles of methane up to 600 m in altitude and showed repeatability of 13 ppbv between multiple profiles. Additionally, the hexacopter system was used to evaluate the evolution of methane in the nocturnal boundary layer during the ScaleX field campaign in Germany, where measured data is consistent with supporting ground-based methane and meteorological measurements.


Introduction
Methane is a strong greenhouse gas and plays a key role in determining the climate impact of sources such as wetlands, livestock, and oil and gas equipment [1][2][3]. There is a continuing need for technology to better understand and track the sources and sinks of methane at the local scale. Unmanned aerial systems (UAS) are attractive platforms for such measurements due to their flexibility and autonomous nature. In particular, small (<25 kg) UAS have low operation costs, relaxed regulatory requirements, and can operate near the surface, where most emissions originate. As an example, a near-infrared sensor (100 ppbv precision) equipped on a remote-controlled helicopter has been used to estimate the magnitude and temporal variability of methane emissions from a natural gas compressor station [4], and serves as an alternative to mobile laboratory and ground survey techniques which are more time-consuming or may miss elevated leak sources.
To leverage the potential of UAS, there is a need for high performance trace gas sensors meeting the payload constraints of typical fixed wing and copter platforms, including for methane [5]. Due in part to this gap, air quality and trace gas monitoring with UAS is still a relatively new area of research. Methane instruments are available that have excellent performance, but are generally designed for aircraft- [6] or ground-based applications [7] where size, power consumption and mass are not significantly restricted. A sensor based on integrated cavity output spectroscopy (ICOS) technology from Los Gatos Research Abstract The design and field performance of a compact diode laser-based instrument for measuring methane on unmanned aerial systems (UAS) is described. The system is based on open-path, wavelength modulation spectroscopy with a 3.27 µm GaSb laser. We design two versions of the sensor for a long-endurance fixed wing UAS and a rotary wing hexacopter, with instrument masses of 4.6 and 1.6 kg, respectively. The long-endurance platform was used to measure vertical profiles of methane up to 600 m in altitude and showed repeatability of 13 ppbv between multiple profiles. Additionally, the hexacopter system was used to evaluate the evolution of methane in the nocturnal boundary layer during the ScaleX field campaign in Germany, where measured data is consistent with supporting ground-based methane and meteorological measurements. weighing 19.5 kg has been flown on the NASA SIERRA UAS platform [8], which like many long-endurance UAS has a large payload capacity. More compact methane sensors are available that typically either use standoff detection for path integrated measurements [9,10], or have been targeted at high concentrations such as for direct sampling from natural gas leaks. In addition, several lightweight in situ methane sensors have previously been demonstrated for high-altitude balloon-borne or small UAS use [4,11,12]. These employed an open-path design, where air flows directly into the optical cell rather than using sampling equipment allowing for low weight and inherently high time response [7]. However, each used near-infrared laser sources, so a further improvement in performance and weight is expected using a mid-infrared laser probing methane's strong fundamental absorption band in the 3.3 μm region. A notable lightweight, mid-infrared system is the Tunable Laser Spectrometer, which has been used to constrain the possible levels of methane on Mars [13,14].
We demonstrate for the first time a lightweight, midinfrared, in situ methane sensor that was modified to operate on two different UAS platforms. The sensor employs open-path wavelength modulation spectroscopy (WMS) for sensitive detection of methane in the atmosphere. The first platform, a long-endurance fixed-wing UAS, had modest payload constraints and allowed us to test precision and stability on a multi-hour flight. The hexacopter system required additional reductions in mass and operates with much shorter flight times, but was well-suited for a nocturnal boundary layer application [15,16], where methane accumulates in the stable, shallow layer near the surface. The sensor's precision of 5-10 ppbv was adequate to resolve vertical gradients of atmospheric background methane (~2 ppmv) in both day-and night-time conditions, and meets previously estimated criteria of <40 ppbv needed to quantify emissions at a landfill [5], which may produce enhancements of 100-1000 s ppbv above the background level. The combination of both large, long-endurance, and small, flexible platforms allows for a variety of applications for probing the atmospheric boundary layer.

Description of instrument for both configurations
The common attributes of the sensor on both platforms are the core optical cell and most of the electronics components. The optical cell housing was fabricated using 3D printed ABS with spacing and rigidity maintained by three 8 mm hollow carbon fiber rods. A Herriott cell [17] was chosen due to its well-established stability in field environments and compact geometry, with an 11.2 cm base path and 2.7 m optical path length (N = 24 passes) confined by Ø (5.08 cm) gold-coated mirrors. As an open-path design, air moves directly through the optical cell based on atmospheric motions relative to the platform movement rather than using a pump. The laser and detector were attached to opposite sides of the cell and aligned using additional steering mirrors. We employed a GaSb distributed feedback (DFB) laser (Norcada), a LDTC0520 laser temperature and current controller (Wavelength Electronics), and an optically immersed HgCdTe photodetector with an AC coupled preamplifier (Vigo). A single board computer (Advantech) and PCIe-6251 data acquisition board (National Instruments) were dedicated to system control and data processing. Finally, a BMP180 sensor (Bosch) was situated at the edge of the electronics box to measure barometric pressure. While these aspects are the same for both configurations, details of the electronics box and mounting differed due to integration and flight operation needs.

Long-endurance system
The platform that formed the basis of the long-range design was an American Aerospace Technologies, Inc. (AATI) RS-20. The RS-20 has a payload of 25 kg, flight times of up to 16 h, and a maximum speed of 46.3 m s −1 . Specific to this UAS, the electronics and sensor were designed to mount underneath one wing of the RS-20, enabled by its ability to take off with an asymmetric load. DC power and an Ethernet link for real-time communications were connected through the wing to the UAS fuselage. The sensor electronics were enclosed in a 20 cm × 26 cm × 11 cm solid aluminum case and placed behind the optical cell relative to the flow of air during flight (Fig. 1a), with both components mounted on an aluminum plate connecting to a hard point on the wing. The overall mounting scheme was designed to ensure air-readiness and to withstand the pneumatic launch that produces forces of up to 10-g. Altitude was estimated using barometric pressure from the BMP180.

Hexacopter system
The hexacopter platform that formed the basis of this configuration was a Pro X-3 Lite by Stein Elektronik (Fig. 1b). This medium lift hexacopter has a 2.1 kg payload capacity remaining after flight battery and a 5-6 min flight duration with a take-off weight of 4.9 kg. The platform necessitated a further reduction in sensor mass to meet these constraints, which was achieved with a lighter, plastic electronics box and a smaller optical cell housing. The hexacopter was operated using a Pixhawk controller and the ArduPilot Copter 3.3 software to program flight patterns and access telemetry data. Here the sensor was powered by its own Page 3 of 9 170 0.14 kg, 1300 mAh 4-cell LiPo battery. The detector controller was updated from STCC-04 (Vigo) to the newer PTCC-01-OEM, which has a smaller form factor. The overall specifications of the methane sensor as deployed on both platforms are shown (Table 1). For the hexacopter system, an Innovative Sensor Technology HYT 271 was added directly outside of the multipass cell to measure the ambient temperature and humidity. This is based on a proportional to absolute temperature (PTAT) measurement and capacitive polymer relative humidity measurement. Altitude and position were obtained from the hexacopter's flight control, which incorporates both pressure and GPS information.

Spectroscopy and data processing
The laser is tuned over a 1.2 cm −1 region centered at 3057.4 cm −1 containing strong absorption features for both methane and water vapor (Fig. 2). The water vapor feature is composed of one dominant line, while methane is a composite of three lines [S = 2.106 × 10 −19 , 1.26 × 10 −19 , 1.26 × 10 −19 cm −1 /(molecule cm −2 )] [18] that are broadened at atmospheric pressure. A 200 Hz sawtooth ramp was applied with additional 20 kHz sinusoidal modulation to control laser drive current. Increasing the ramp and modulation rates is advantageous for increasing signal-to-noise but is limited by the ability to tune fully over the water and methane features. Data were averaged and logged at 20 Hz for the fixed-wing configuration and 10 Hz for the hexacopter configuration.
The harmonic signals were generated in real-time using a LabVIEW-based lock-in detection scheme [19] with methane and water vapor concentrations retrieved  using a second harmonic (2f) peak-to-trough height calculation normalized by the first harmonic (1f) at the peak centers to account for intensity variations. The 2f/1f values were calibrated either by placing the whole sensor in an enclosed box and flowing calibration standards or with comparison to reference instruments, as described in Sect. 2.3. In both cases, a two-point calibration approach was used with a zero and linear slope.
In post processing, methane and supporting data were time synchronized to the GPS time and location information at 1 Hz. An internal GPS was used for the fixed-wing system, while for the hexacopter the copter GPS was used to avoid the need for an additional antenna. Data with low signal strength (defined as detector voltage <0.2 V) were then excluded to remove incidences where the beam was blocked. To account for pressure effects a correction function, κ(P) is applied [7], which was derived based on simulations using HITRAN 2012 data [18] to account for spectroscopic and density (ideal gas law) effects defined relative to standard conditions (296 K, 1013 hPa). This correction additionally accounts for the fact that a constant modulation amplitude was used, leading to changes in the modulation index and the shape of the 2f signal. For the fixedwing configuration, temperature measurements were not available in-flight and a ground-based, fixed station was used. For the flight altitude envelope of 600 m during midday, the expected change in temperature is expected to be 3-6 K colder than the surface based upon environmental/ adiabatic lapse rates. This temperature decrease results in an increase in the CH 4 concentrations (i.e., absorption) by 1-1.5% from the ideal gas law and a decrease in absorption from the linestrength of 1-1.5%. In combination, these effects nearly cancel out and thus the overall uncertainty in the absolute concentration from ignoring temperature is negligible compared to instrument noise (5-10 ppbv).
An additional limitation was that active line locking (stabilization of the laser tuning) was not used, so the methane concentration was prone to drift if there were changes in the peak position. This was primarily an issue for the hexacopter flights, which were short in duration and had more rapidly varying environmental conditions than with the fixed-wing configuration. This was accounted for using a polynomial to correct for changes to the extent that drift in concentration was correlated with peak position. Individual hexacopter flights were excluded if peak position drifted outside of the preset window used by the single board computer for calculation of peak-to-trough heights.
A final aspect to the approach was inclusion of an inline reference-cell containing the target gas (methane) at reduced pressure. This is a new variation on previous use of in-line reference cells which have targeted a separate gas that is spectrally non-overlapping or distinguishable using higher harmonics [19,20]. Here the in-line reference cell was 2 cm long and contained 0.4% methane in N 2 at 30 Torr, which can be effectively measured using direct absorption spectroscopy (DAS). The signal can be extracted using a time-multiplexed scheme where there are both scans with and without sinusoidal modulation every second [21]. A duty cycle of 19:1 modulated to non-modulated scans was used, so the sensor was predominately aimed at WMS, while the hexacopter system ran in fully WMS mode since the reference cell had lost pressure. The parameters were optimized so that the direct absorption current offset was increased slightly, but not so much as to cause laser instability. Experimental spectra will be discussed in Sect. 3.1.

Long-endurance UAS
Both systems were deployed for field testing and measurements of methane mixing ratios in the lower atmosphere. Fixed-wing UAS flights were conducted on September 9, 2015 and February 24-25, 2016 operated out of the Cape May Airport, a Federal Aviation Administration (FAA) UAS test site in southern New Jersey. The first flight encompassed a 12 × 15 km area over the nearby coastline (Fig. 3a) and is the focus of our analysis here. The sensor was calibrated beforehand using nitrogen as a zero and a span gas of 2.045 ppmv methane (GTS-Welco 2.12 ± 5% ppmv methane in air, cross-calibrated against a NOAA standard to 2.045 ppmv ± 1% using a LI-COR LI-7700). A second flight was flown on February 24 but there were difficulties with wireless communication and the sensor's internal computer restarting. On February 25, a controlled methane release was conducted, with the UAS making downwind passes at the Cape May Airport. Analysis of the data did not see a methane enhancement from the release since wind meteorological conditions were not favorable. Information from a 3D sonic anemometer was used with WindTrax (http://www.thunderbeachscientific.com) to model the release, which also supported that an observable enhancement would not be expected above 50 m.

Hexacopter
Flights with the hexacopter system were conducted in July 2016 as part of the ScaleX campaign (http://scalex. imk-ifu.kit.edu) [22], a "scale-crossing" experiment at the Fendt field site within the TERestrial ENvironmental Observatories (TERENO) network [23] in southern Germany. An innovative approach had been employed in ScaleX 2015 where a UAS was equipped with a Teflon drag-tube connected to a Picarro methane analyzer on the ground, where results have been described elsewhere [22]. Here, we use our new sensor to obtain direct, in situ measurements with efforts focused during an intensive observational period from 21:00 to 06:00 (local time) on July 06-07.
Flight tracks are shown at the Fendt field site (Fig. 3b), along with the location of an instrumented 9 m tower containing air inlets at 1, 3, and 9 m. A gas switching system was used in connection with a Picarro G2508 that measured at these heights sequentially, sampling each height for 2.5 min. A container also located in Fendt (not shown) was equipped with a Los Gatos DLT-100, providing additional long-term information about methane patterns at the site [24].
In total, ten flights were conducted during the UAS intensive observational period, along with five additional daytime flights in Fendt and one at the IMK-IFU on July 18, 2016. One extended (2 h) intercomparison was conducted on July 11 with the sensor collocated at the Picarro tower's 1 m inlet. A second intercomparison (2 h) collocated with the Los Gatos inlet on the roof of a measurement container was conducted on July 17 starting at sunset. Reference gas measurements (CH 4 = 1893 ± 1% ppbv) were made several times during the campaign using a closed chamber, and zero gas measurements were done using nitrogen in a lab at the IMK-IFU campus. A constant linear calibration slope was used throughout campaign, but an offset was obtained by comparison with the Picarro each flight. This was a solution to help account for drifts by comparing to a highly stable ground-based instrument without having to do full re-calibrations in the field.

Fixed-wing based system
A 30-min section of data obtained while flying over the Cape May region at constant altitude was used as a measure of the sensor precision in the in-flight environment. An Allan-Werle deviation plot [25] with this data shows a precision of 5 ppbv at 1 Hz (Fig. 4). Short-term results follow the behavior of white noise (slope ∼ 1/ √ τ ) while there is a relatively short turn-over time of 10-15 s. By comparison, an experiment in the laboratory while flowing a calibration standard (2.045 ppmv, 1% uncertainty) found a precision of 4.5 ppbv at 1 Hz, a decrease of only 12% in the field environment suggesting a robust optomechanical structure and  insensitivity to the vibrations on this platform. Both experiments also likely included effects from uncontrolled temperature variations and in the field-case, real atmospheric variability in methane when traveling over lateral distances of 10's of km.
To test stability and reproducibility on a longer timescale, a sequence of six vertical profiles was also obtained over a 15 min time period during the same flight (Fig. 5). Measurements were conducted while making transects along the tidal wetlands between Cape May and Wildwood. The UAS started and ended the measurement sequence at 430 m and profiled between 120 and 583 m. The atmosphere is generally considered well-mixed (i.e. constant with altitude) in the afternoon boundary layer, but variation is possible at the 10's of ppbv level. Even at the nearby NOAA profiling site located over the Atlantic Ocean ~50 km SE of Cape May, methane mixing ratios often deviate from flat or decreasing with altitude. The UAS profiles we observed showed measurable structure, with minima at 130 and 400 m and maxima at 240 and 540 m, indicating that the sensor can resolve atmospheric variability since results were repeatable across multiple up-down profiles. The deviation between profiles (1σ) had a mean of 13 ppbv and was as low as 8 ppbv near 330 m, which is an additional proxy for the upper limit of in-flight precision.
The multiplexed DAS approach provides additional information about the stability of the sensor during this flight. An example of an in-flight non-modulated scan and subsequent processing sequence is shown (Fig. 6). The spectra are averaged and saved at 15 s interval, but because only 1 in 20 scans are non-modulated, each recorded spectrum represents an integration time of 0.75 s. A fifth order polynomial was used to extract the absorption baseline, calculated from nine fitting points located away from any absorption features. The baseline-subtracted absorption signal of water vapor and methane and the residual between measurement and Voigt line fitting is also shown (Fig. 6). The fit is accomplished using nonlinear least squares to solve for the CH4 ambient , CH4 ref. cell , and H 2 O mixing ratios. Three sharp peaks from the reduced pressure reference cell are clearly visible overlaying the ambient methane feature, which enables it to be retrieved separately from the ambient signal. The standard deviation (σ) of the residual was calculated between 3057.4 and 3057.6 cm −1 , away from the line centers where the residual may be driven by the accuracy of the Voigt line shape function rather than random noise, was approximately 10 −3 . We did not use the ambient methane concentration data directly because of the higher noise compared to WMS, but rather used CH4 ref. cell as a continuous proxy short-term instrumental noise.

Hexacopter-based system
Ten night-time flights were taken during ScaleX, while six passed the quality check criteria and are shown as vertical profiles sequenced by the time of night/morning (Fig. 7), excluding the sixth where pressure and temperature were not saved. The red dots show methane mixing ratios binned by altitude, while the trend is illustrated using Gaussian kernel smoothing. We use 2 m bins since the air is likely mixed below this scale due to the circulation produced by the hexacopter rotors. The downdraft effect on trace gas measurements is an active area of research, but is thought to cause downward mixing of air and may disturb fine-scale features of the profile. The 21:00 profile was well-mixed up to 75 m, consistent with expected results during daytime, and a layer top at 100 m that was also seen in the temperature profile. The profile at midnight saw a decreasing concentration from ground level up to about 125 m. The remaining profiles take a markedly different shape, likely due to local meteorological conditions. Transport of methane from nearby dairy farms is known to cause enhanced methane [24], with the closest farm located 0.5 km southwest of the main UAS launch location. The Fendt site is also influenced by a nocturnal low level jet which tend to mix turbulence downwards from aloft [26]. The predominantly northerly wind was seen to shift to southwesterly around 02:00.
The evolution of methane seen by the Picarro-instrumented tower is shown, along with a comparison to the average measurement from the UAS at the equivalent height (Fig. 8). The ground based measurement resolves some of the complexity, with methane greater at the surface than at 9 m up until 23:00 and after 04:00, but inverted for much of the time in-between. The UAS found consistent results, with small gradients at 21:00 and 0:00 and a larger, inverted gradient at 03:00. Up until 02:30 there was an accumulation seen by the Picarro analyzer. Then the concentration as well as the gradient decreased followed by an again developing strong gradient from 04:00 until 05:15. The gradient at 05:00 from the UAS was consistent in sign but about half the magnitude of the Picarro. This is likely the result of the strong vertical, and likely horizontal, gradient that makes it difficult to directly compare the measurements. Similarly large (>100 ppbv) differences were also sometimes observed between the tower and UAS-drag tube during the ScaleX 2015 campaign [22]. The comparison is also more uncertain where the mixing ratio is changing quickly, since the UAS is measuring over a period of time and the Picarro also only samples at one inlet at a time.

Summary and conclusions
We demonstrated a compact, open-path methane sensor with measurements conducted on both small and large unmanned aerial systems and which weighed 4.6 and 1.6 kg, respectively. In-flight precision of the fixed-wing Vertical profiles were presented in the daytime over a coastal area up to 600 m a.g.l. and during night-time at a grassland site in southern Germany. The night-time vertical profiles were compared to available observations from two heights at a 9 m tower. Additionally, a same-gas reference cell was described providing additional evidence of stability during the flights.
The ability of the sensor to operate on small UAS opens up new opportunities for studies of trace gases in the atmospheric boundary layer, due to their low cost, ability to operate down to ground-level, highly flexible flight patterns, and reduced regulatory scrutiny. The combination of small UAS with a high-precision methane sensor will provide a more-detailed knowledge of vertical and horizontal heterogeneity of methane, as well as applications in localizing and quantifying emission sources. Further work will improve the long-term stability of the sensor to reduce the need for ground-based reference instruments and ruggedize the sensor for repeated field use. Fig. 8 Methane measurements taken from the tower and from the UAS during the same period as in Fig. 7. The open circles are data measured by the Picarro instrument and connected with lines, with black and blue corresponding to the 1 and 9 m inlets, respectively. The 3 m inlet is omitted for clarity and showed similar trends. Cross marks use the same colors and correspond to measurements from the UAS averaged at 1, 9, and 49 m