Real-time remote detection and measurement for airborne imaging spectroscopy: a case study with methane

Localized anthropogenic sources of atmospheric CH4 are highly uncertain and temporally variable. Airborne remote measurement is an effective method to detect and quantify these emissions. In a campaign context, the science yield can be dramatically increased by real-time retrievals that allow operators to coordinate multiple measurements of the most active areas. This can improve science outcomes for both singleand multiple-platform missions. We describe a case study of the NASA/ESA CO2 and MEthane eXperiment (COMEX) campaign in California during June and August/September 2014. COMEX was a multi-platform campaign to measure CH4 plumes released from anthropogenic sources including oil and gas infrastructure. We discuss principles for real-time spectral signature detection and measurement, and report performance on the NASA Next Generation Airborne Visible Infrared Spectrometer (AVIRISNG). AVIRIS-NG successfully detected CH4 plumes in realtime at Gb s 1 data rates, characterizing fugitive releases in concert with other in situ and remote instruments. The teams used these real-time CH4 detections to coordinate measurements across multiple platforms, including airborne in situ, airborne non-imaging remote sensing, and groundbased in situ instruments. To our knowledge this is the first reported use of real-time trace-gas signature detection in an airborne science campaign, and presages many future applications. Post-analysis demonstrates matched filter methods providing noise-equivalent (1 ) detection sensitivity for 1.0% CH4 column enhancements equal to 141 ppmm.


Introduction
Airborne imaging spectrometers have been deployed for a wide range of scientific, regulatory, and disaster response objectives. Traditionally these campaigns wait for favorable environmental conditions and then fly pre-arranged survey patterns (typically "mowing the lawn"), recording data for post-flight radiometric calibration and geolocation. Significant time can pass before data are analyzed fully, and results often arrive too late for mid-course corrections during the campaign. However, improvements in computing power, communication, and telemetry are changing this situation. Tactical remote measurement generates in-flight calibrated data products to inform a real-time adaptive survey strategy. This can be coordinated to direct other platforms in multiplatform campaigns. We use the term "tactical" to emphasize environmental awareness and real-time decision making, with no military connotation. Its applications include the following.
i. Detection of transient or rare targets -Many airborne missions hunt isolated or nonstationary phenomena. Examples include trace-gas emissions Gerilowski et al., 2015), algal blooms (Karaska et al., 2004), invasive species (Ustin et al., 2002), isolated microhabitats (Thompson et al., 2013b), and hurricane intensity (Braun et al., 2013). Aircraft use radar to hunt extreme weather, and lidar to find cirrus, thunderstorms or biomass burning (Rolph, 2003). In each case, tactical remote measurement can identify desired iii. Data quality assurance -Tactical remote measurement adds flexibility and confidence to flight management decisions. Currently, mid-campaign flight planning often occurs without knowing the quality of data already collected. This risks wasting resources if, for example, the mission continues under marginal environmental conditions. On the other hand, conservative planning can miss opportunities. Tactical science products can inform flight plans and mid-day scrub decisions to avoid spending flight hours on low-value or redundant data. For example, it may reveal interference such as cirrus clouds (Gao et al., 1993a), sun glint (Kay et al., 2009), and unacceptable aerosol scattering (Bojinski et al., 2002). This also allows instrument subsystem failures to be recognized and addressed immediately.
iv. Robotic exploration -Real-time analysis can improve autonomous operations when communication opportunities are rare and bandwidth is limited, such as in space exploration. Remote spacecraft that are out of touch with ground control can autonomously detect highvalue spectral signatures that guide prioritized downlink or trigger additional measurements (Thompson et al., 2013a). Operators can generate compact map products onboard the spacecraft and downlink them to supplement raw spectra, expanding spatial coverage at a low bandwidth cost. Onboard cloud screening is one example of data volume reduction; it can improve yields by a factor of 2 or more for Earth orbiting instruments .
v. Hypothesis formation and testing -Real-time data analysis and visualization in a mapping environment, like Google Earth (Google Earth, 2015), is common in surface and airborne in situ applications. Many systems allow the scientist to visualize spatial relationships between measured parameters, forming hypotheses on the fly for immediate testing. Adaptive surveying can address new hypotheses during the campaign, while instruments are deployed and environmental conditions are favorable. Telemetering live data allows remote investigators to observe and participate in operational decisions .
These techniques require high-performance data telemetry and communication. As the technologies proliferate, unanticipated applications are likely to appear -just as instant results from the digital CCD transformed chemical photography in dramatic and unforeseen ways.
This study demonstrates tactical remote measurement with imaging spectroscopy during a multi-aircraft, multi-platform campaign, CO 2 and MEthane eXperiment (COMEX). The COMEX campaign was funded by NASA and ESA to explore synergies between NASA's proposed HyspIRI (Hyperspectral Infrared Imager) mission and ESA's CarbonSat Earth Explorer 8 candidate mission. Greenhouse gas emissions were measured from a range of important anthropogenic sources. Investigators surveyed landfills, husbandry, and fossil-fuel production sites in southern California during summer and fall, 2014. A multi-scale experimental design combined airborne and surface measurements to characterize CH 4 sources on scales of meters to tens of kilometers. Ground-validated airborne imaging spectroscopy identified sources and their heterogeneity. This was followed by downwind surface surveys together with airborne sounding and in situ observations transecting plumes at different upwind and downwind distances. Surface mobile survey teams carried sensors to specific locations of interest. Finally, repeated surface in situ surveys studied longer term temporal variability and larger spatial context. COMEX exploited tactical remote measurements from multiple platforms. We focus on one participating instrument, the Airborne Visible Infrared Spectrometer -Next Generation (AVIRIS-NG) (Hamlin et al., 2011;Green et al., 1998), which mapped CH 4 enhancements in real time. A simple detection method based on a band ratio (BR) was sufficient to detect several sources and enhance the COMEX campaign. These initial results motivated the development of a more sophisticated matched filter detection approach, described in this paper, which was developed after COMEX and has been adopted by subsequent CH 4 monitoring campaigns. Although prior studies quantified CH 4 anomalies using Visible Shortwave Infrared (VSWIR) imaging spectrometers Bradley et al., 2011), we believe the COMEX campaign to be the first real-time tactical deployment for remote trace-gas imaging.
Section 2 describes the real-time algorithms, system architecture, and implementation decisions. Section 3 reviews the campaign results including an AVIRIS-NG sensitivity analysis and discussion of lessons learned. We conclude with a discussion of future directions for tactical remote sensing.

Tactical imaging spectroscopy
The visible/shortwave infrared (VSWIR) imaging spectrometers serve diverse applications including mineralogical mapping (Kruse, 2012), characterization of coastal and terrestrial ecosystems (Ustin et al., 2004), and atmospheric studies (Popp et al., 2012). Imaging spectrometers are valuable for tactical operations because they can map and localize targets over wide areas, providing reconnaissance for other instruments along with spatial and spectral context. Realtime airborne imaging spectroscopy has been deployed in a few previous instances. For example, Thompson et al. (2014) demonstrated real-time cloud screening for future space missions. Bue et al. (2015) calculated reflectance products using the model-based ATREM atmospheric correction. Eismann et al. (2009) demonstrated the ARCHER system which provided real-time processing for search and rescue applications. They performed matched filter spectral signature detection and change detection using the chronochrome method (Schaum and Stocker, 1998). They also demonstrated spectral anomaly detection using the Reed-Xiaoli (RX) detector (Stein et al., 2002), an anomaly score based on the Mahalanobis distance (Chang and Chiang, 2002). These methods detected artificial objects in wilderness scenes, such as parts of aircrafts and vehicles near crash sites. Another real-time airborne investigation used non-imaging spectroscopy for detecting dangerous volcanic plumes (Vogel et al., 2011). On the ground, real-time analysis has imaged these plumes' SO 2 and SiF 4 absorption (Lübcke et al., 2013;Stremme et al., 2012), enabling subsequent analysis to infer emission rates Krueger et al. (2013).
This section describes the real-time system used by AVIRIS-NG during the COMEX campaign. AVIRIS-NG measures reflected sunlight in the 0.38-2.5 µm range with 0.005 µm spectral resolution. Its 1 mrad instantaneous field of view (iFOV) provides sub-meter ground sampling distance (GSD). The real-time system characterizes CH 4 plumes by analyzing absorption features from 2.1-2.4 µm (Clark et al., 2009). Its design must balance the competing needs of speed and algorithm sensitivity, and several guiding requirements drive our decisions. First, the system must provide a sensitivity floor sufficient to detect the phenomena of interest reliably. In other words, it must have a signal to noise ratio (SNR) sufficient to find sources under relevant wind, illumination, and substrate conditions. Only then can planners safely act on a null detection result. For similar reasons, it must minimize false positives. Prior studies of CH 4 with the "Classic" AVIRIS instrument by Thorpe et al. (2014) detected local enhancements of 1 ppm within a kilometer-thick atmospheric model layer. Later stud-ies by Thorpe et al. (2015) using AVIRIS-NG found similar enhancements in the distal regions of plumes associated with CH 4 fluxes of 14.2 m 3 h 1 (500 standard cubic feet per hour, scfh) under moderate (5 m s 1 ) winds. Resolving plumes of this magnitude under similar conditions should be possible with a sensitivity of 1000 ppm m. Better performance would further reduce ambiguity and improve the detail of diffuse plumes.
A second requirement is high spatial resolution. The phenomena should subtend multiple pixels with sufficient resolution for the operator to identify typical morphologies. The diagnostic shape of atmospheric plumes can be corroborated with ancillary wind information , while in the case of oil slicks thickness asymmetry and shape are useful cues (Leifer et al., 2012). For plumes, resolution can enhance detection sensitivity due to non-uniformity: many plumes are initially buoyant, rising abruptly in a column for tens of meters before dispersing downwind. High spatial resolution avoids diluting this feature, which may be only a few meters in diameter. Fine spatial resolution also helps exclude false positives caused by artificial features with obvious geometric shapes. For these reasons, we desired that the system would process AVIRIS-NG at native resolution without subsampling.
A third requirement is speed. Speed follows a "threshold" utility function: the system must operate at the instrument data rate, but additional performance provides no extra benefit. Real-time operation avoids a confusing temporal association puzzle where a detection appears at a location passed seconds or minutes ago. In addition, keeping pace with data collection simplifies operations by permitting the system to operate whenever the instrument collects data. We find it possible to implement many of the most common detection methods from literature within this requirement, though speed considerations determine the physical quantities that we retrieve. Specifically, we focus on measuring plume absorption along the optical path. This is sufficient to indicate the relative strengths of different sources. We do not estimate the vertical structure or flux; such products generally require complex iterative retrievals involving ancillary data such as wind speed, and are less critical for real-time decisions.
The real-time system first executes the standard AVIRIS-NG ground data pipeline to create calibrated radiance measurements. It then matches these spectra to the known gas absorption signature of CH 4 . The procedure successfully operates at instrument data rates of approximately 500 Mb s 1 and allows operators to detect CH 4 emissions in real time.
The following sections detail specific design choices for the software architecture and detection algorithms.

Instrument and software architecture
The AVIRIS-NG instrument acquires 598 cross-track spectra at 100 Hz. Frames are captured with a custom field programmable gate array (FPGA) frame grabber over a dedi- cated Camera Link interface at 500 Mb s 1 data rate. Data are synchronized with an onboard inertial measurement unit (IMU)/GPS system , and finally stored in a solid-state RAID array. The AVIRIS-NG console records un-orthorectified raw data and displays it for the operator. A backup computer records a second copy in parallel, and the detection system runs on this machine. Real-time analysis requires that detection algorithms keep pace with the data recording rate, while leaving enough CPU cycles for the backup data recorder.
Our solution exploits parallelism with multi-core CPUs ( Fig. 1). A watchdog process waits for a new image to appear on disk. As the instrument writes to this file, an executive real-time process begins reading from the other end and buffering 1000 lines at a time. The real-time analysis applies radiometric calibrations to each 1000-line block and partitions the resulting data into spatially independent regions for multi-thread detection. When all threads have finished, the results are reassembled and recorded to storage, where they are immediately available on the operator display. The system processes 10 s intervals of data in well under 10 s, achieving the real-time speed requirement.

Onboard radiance processing
The detection pipeline first transforms the frame (a single cross-track slice of data) to a calibrated radiance product (Green, 1998;Hamlin et al., 2011). We calculate radiance in W cm 2 nm 1 sr 1 at each cross-track spatial location c and  Figure 2. A comparison of spectral shapes between the CH 4 transmission spectrum, resampled to AVIRIS-NG wavelengths, and the target signature t used for detection. The vertical axis plots two different quantities as noted in the legend. Both signatures were calculated from a 20-layer atmosphere based on HITRAN 2012 absorption cross sections  wavelength using: where R(c, ) are the raw digital numbers from the instrument. R d (c, ) is the electronic dark current estimated from a closed-shutter segment at the beginning of each flight line. R p (c) are electronic "pedestal shift" effects, in which a spatially compact signal depresses the signal at other spatial locations. The onboard system estimates the pedestal shift of each spectrum based on the residual dark current in nonilluminated edges of the detector. r( ) and f (c, ) are the radiometric calibration coefficients and flat field corrections, respectively. Both are calculated from laboratory calibration sequences using a known spatially uniform illuminant under fixed imaging geometry. Appendix A discusses wavelength calibration. Figure 2 shows a typical CH 4 transmission signature, calculated using a model atmosphere with absorption coefficients of Brown et al. (2013). The detection algorithm calculates a scalar score to estimate any local enhancement of this background. We evaluated several detection algorithms based on their sensitivity and speed. At one extreme, an iterative nonlinear or "optimal estimation" solution such as iterative maximum a posteriori differential optical absorption spectroscopy (IMAP-DOAS) (Frankenberg et al., 2005;Thorpe et al., 2014) is more quantitative, but somewhat slow for real-time operation. At the other extreme, an absorption band depth score uses simple arithmetic, but its low SNR can detect only the strongest signatures. This paper focuses on a third approach: a novel matched filter variant with a good balance of sensitivity, stability, and speed, and which also permits a quantitative interpretation.

Onboard signature detection
Our first algorithm uses a continuum interpolated band ratio (CIBR), defined as the depth of an absorption feature relative to a local linear continuum (Green et al., 1989;Bruegge et al., 1990). It is written as follows: where center , left , and right are wavelengths in the middle and either side of the absorption feature. The weighting coefficients w sum to unity, and make the denominator the linearly interpolated continuum at the location of the absorption center center . We find the 2.37 µm feature provides the best overall contrast. The CIBR method is simple to implement and fast to execute. For the sources studied during COMEX, its sensitivity is sufficient to detect strong local CH 4 enhancements. The second detection strategy is a classical matched filter (Manolakis et al., 2014), a variant of which was used previously for CH 4 detection by Thorpe et al. (2013). The matched filter tests each spectrum against a target signature t while accounting for the background covariance. Here t is a vector with one element per wavelength. If the background spectra are distributed as a multivariate Gaussian N with mean vector µ and covariance matrix 6, the matched filter is equivalent to a hypothesis test between the case H 0 where the target is absent and H 1 where it is present.
Here, t is the target signature. The matched filter estimates the scalar value ↵, the fraction of the target (potentially larger than unity) which perturbs the background. Larger values of ↵ signify a stronger match. The matched filter is written.
The hat symbols indicate that the background mean vector and covariance matrix are estimated using samples from the scene. One typically draws samples from a rectangular region near the target. However, most push-broom sensors have a slightly nonuniform behavior at different cross-track positions, which violates the Gaussian background assumption. The cross-track push-broom elements are separate detectors, so it often is better to model their noise distributions independently. Thus, we apply an independent matched filter to each column of the (non-orthorectified) image, calculating a new µ and 6 for each cross-track element. This columnwise matched filter dramatically reduces the number of samples available for estimating each 6. We compensate by estimating a stable, low-rank approximation of the inverse sample covariance as in Manolakis et al. (2009). The covariance matrix 6 decomposes as product of p column eigenvectors q and p eigenvalues : The top d eigenvalues approximate the inverse. With the identity matrix I and trace operator tr, we havê We typically estimate 30 eigenvalues for vertical blocks consisting of 1000-2000 samples per column.

Target signatures
The signature t should match the spectrum of the target feature. A reasonable approach is to use the transmission shape itself (the red curve in Fig. 2). However, this is inaccurate when absorption is strong; further attenuation becomes nonlinear as absorption lines saturate. The matched filter assumes a linear perturbation, so the Jacobian of the radiance spectrum is an appropriate signature. We calculate it by modeling local CH 4 enhancement as a uniform cell. The airborne instrument measures absorption along a path transecting the CH 4 cloud. For thin, uniform plumes, the unknown quantities of absorption length and concentration are interchangeable, so we consider the combined quantity, the mixing ratio length expressed in ppm m (Sandsten et al., 2000). Our derivation is similar to that of Theiler et al. (2005). Following Eq. (4) the matched filter estimates ↵, which is a multiplicative scaling of the target signal that perturbs the mean background radiance µ. This background includes absorption by ambient CH 4 . Under hypothesis H 1 , a local enhancement acts as a concentration-dependent absorption coefficient ( ) and absorption path length`. For clarity, we write this relation with a functional form µ( ) to represent a single wavelength of the mean vector µ.
For x near zero, the first-order Taylor expansion exp(x) ⇡ 1 + x permits Combining all wavelengths using vector notation, and folding unknowns into ↵, yields a near-surface plume of unit concentration and unit length, multiplied by the background mean radiance. The resulting matched filter estimates ↵, the scaling of the unit concentration path length. The target signature represents the perturbation, in radiance units, of the background radiance by an additional unit mixing ratio length of CH 4 absorption, which acts as a thin Beer-Lambert attenuation of the (already attenuated) background µ. Evaluating the partial derivative of Eq. (9) at`= 0 gives We can estimate the enhancement of CH 4 using the linear scaling of a target signature that perturbs the mean radiance; that signature is defined as the negative absorption coefficient scaled by the (wavelength-dependent) radiance. Figure 2 compares the shape of the Jacobian target signature to the typical transmission signature of ambient CH 4 . The linearization works for thin plumes even when the background is saturated, because deviations are small and can be modeled linearly to permit a fast yet accurate quantitative retrieval . It ignores scattering effects, which is a reasonable compromise at low flight altitudes; spectral features caused by actual CH 4 enhancements by far exceed typical retrieval biases that could be induced by atmospheric scattering . In addition, Rayleigh and aerosol scattering is much lower in the 2.3 µm region than in the UV and visible spectral range. The linearized approach is complementary to other more complete retrieval algorithms, such as the IMAP-DOAS approach (Frankenberg et al., 2005;Thorpe et al., 2014).

Operator display
After the detection stage maps plume intensities, an interface displays this information to the operator in a more tactically useful format. Specifically, the display overlays detected plumes on RGB images for visual interpretation and localization (Fig. 3). It supports variable detection thresholds so the operators can set the cutoff concentration according to their tolerance for false positives. This is important because source strength varies and the detection sensitivity changes with solar zenith angle. Also, it is important that the system preserve the overlay images in memory until explicitly reset, so that the operator has time to consider ambiguous detections. Figure 3 was produced by playing back a June flight line using a recent iteration of the detection software. A simple, intuitive interface minimizes unnecessary controls. A vertical slider scrolls the flight line to review previously collected data. During recording, the system appends data to the end of this image. A horizontal slider adjusts the detection threshold, allowing the operator to change the overlay sensitivity based on their tolerance for error. Bright red pixels signify a strong signal well above the threshold, while dark red pixels signify an ambiguous signal.

Results from the COMEX campaign
COMEX field data collection included the Kern River, Kern Front, and Poso Creek Oil Fields, located to the north and northwest of Bakersfield, CA in the San Joaquin Valley (Fig. 4). Along with AVIRIS-NG, the COMEX campaign deployed a second aircraft: the CIRPAS Twin Otter, which carried the Methane Airborne MAPper (MAMAP) Krings et al., 2011), a non-imaging spectrometer, and an in situ Picarro CH 4 sensor sponsored by NASA, Ames. The campaign also deployed the AutoMObile greenhouse Gas Surveyor (AMOG) car-mounted system for in situ CH 4 and wind measurement . These platforms used several real-time displays and communications links. AMOG used a map overlay (Google Earth, 2015) displaying CH 4 measurements along with the wind direction. MAMAP was modified for COMEX to deliver real-time retrieved CH 4 information using a WFM-DOAS algorithm described in Krings et al. (2011Krings et al. ( , 2013. These data were displayed on the MAMAP instrument scientist aboard the aircraft, and overlaid on a map (Google Earth, 2015) for tactical decision making. MAMAP also transmitted its real-time CH 4 measurements together with telemetry and data from other CIRPAS sensors to the CIRPAS data acquisition and assimilation system, where they were downlinked by satellite to the command center. All aircraft were tracked during the mission using the Airborne Science Mission Tools Suite (Duley et al., 2011), the ground segment to the NASA Airborne Science Sensor Network (Van Gilst and . This system provided a web-based service for real-time communications between aircraft operators and the science team (Fig. 5). It also integrated real-time aircraft position and state information through a common map display. We focus here on 3 days when all sensors and platforms were active in the field: 13 June, 2 September, and 4 September 2014. During the June investigation only MAMAP was analyzed in real time; AVIRIS-NG CH 4 detection took place offline. We later installed the band ratio algorithm and used it onboard AVIRIS-NG in September, where it operated successfully. Finally, the full columnwise matched filter for CH 4 was developed and installed after the campaign, and used in post-analysis.

Operational implementation of tactical remote measurement
AVIRIS-NG flew along six neighboring flight lines. Tactical remote measurement was implemented during the 2 September data acquisition, and detected many "hot spots" of high CH 4 concentrations. Operators noted strong plumes on flight lines 2 and 3, on the west side of the study area, and relatively weak activity on the east side. AVIRIS-NG operators transmitted plume coordinates to the ground team by text message via the command center. The MAMAP and AVIRIS-NG aircraft also shared their observations using direct radio communications. Flight line 2, recorded in flight as having the largest number of hits, was revealed by post-analysis to contain patches of high activity (Fig. 6). The final data acquisitions focused on this flight line, which comprised three of the final nine images. Focusing on the west side significantly increased the total plumes imaged for the day, and also provided improved data on temporal variability. Operations on 4 September made further use of real-time AVIRIS-NG and MAMAP data. As before, the flight lines contained many active plumes. However, the ground team's initial data collection at the Kern River oil field did not show significant CH 4 . The ground team coordinator suggested shifting data collection west to the Kern Front Oil Field. AVIRIS-NG confirmed the presence of CH 4 plumes in this area, and the experiment coordinator rerouted the surface teams. The surface teams relocated and identified an exceptionally strong CH 4 plume coming from the direction of a drilling rig, approximately 1 km from the road. Subsequent analysis of AVIRIS-NG data indicated a plume in this  area that originated at a small structure about 100 m from the road. The real-time MAMAP data also observed a large scale plume originating from that area (Gerilowski et al., 2014).
Thanks to the tactical remote sensing, AVIRIS-NG continued covering this area with repeat overflights at 10-20 min intervals that revealed the temporal evolution of many plume features. Figure 7 shows a sequence of revisits to one of the stronger targets. Commercial satellite imagery from Google Earth reveals the source, a pump jack (Google Earth, 2015).
The AVIRIS-NG data shows strong CH 4 absorption near the source, which would be expected for a vertical rise by a buoyant column. Turbulence causes the plume structure to become discontinuous as it disperses downwind. The matched filter resolves the plume at concentrations as low as 500 ppm m, approximately 3 standard deviations above the background noise. Figure 8 shows the agreement between the different instruments on 4 September. Colored pixels indicate CH 4 mixing ratio lengths from AVIRIS-NG. The monochrome dots show MAMAP retrievals: black signifies the background up Table 1. Selected flight days in the COMEX campaign. Often AVIRIS-NG overflew the same plume multiple times. Here the "Plumes" column records the total number of instances that a plume appears in the data, rather than the number of physical plumes.

Date
Real With these caveats in mind, we directly compared the two measurements. We calculated a spatially weighted average of the AVIRIS-NG estimated enhancement above background, matching it to the MAMAP response while simultaneously searching over the sounding's position uncertainty radius of 50 m. The resulting estimates were as high as 103.5 ppm m (white dot) and 101.3 ppm m (grey dot), of similar magnitude to the MAMAP retrievals. Despite differences in observing conditions, both data sets evidence similar-scale enhancements at this site. In other cases, small plumes visible in AVIRIS-NG were sometimes invisible in MAMAP data. Table 1 summarizes the total number of plume instances appearing in each day's data, as revealed in post-analysis of AVIRIS-NG data. We record only the unambiguous detections with a long axis greater than five pixels, and exclude compact detections from becalmed image segments -their bright, concentrated plumes show a strong signal but lack the morphological cues needed for unambiguous attribution.

CH 4 detection sensitivity
We use the 13 June flight lines as a control case to evaluate sensitivity, because the tactical remote sensing system was not operational; i.e., data collection was "blind." Figure 9 shows a typical map of plume thicknesses in ppm m. Retrieved values approach 2500 ppm m. A defined plume is evident, along with turbulent structures that disperse 100-200 m downwind of the source. The insert shows a highresolution visible image, which reveals the source to be a pump jack. Figure 10 shows the same scene analyzed with alternative algorithms: band ratio, matched filter, and columnwise matched filter detection strategies respectively, with intensities scaled to the maximum on-plume pixel. The band ratio barely reveals the largest plumes with many contiguous "hot" pixels. A classical matched filter improves performance and the columnwise version is cleaner still. One gleans a final SNR benefit using the Jacobian rather than the transmission as a target signature.
We evaluated detection sensitivity by exhaustively labeling all CH 4 plumes in the 13 June flight lines by manual inspection. There were 29 obvious plumes, some of which were repeat overflights of the same physical location. For each plume, we identified 3-10 on-plume pixels having the highest estimated concentration. We then calculated SNR using a large rectangular region of pixels 100 m upwind as the background. Figure 11 shows the relationship between plume strength (in units of ppm m) and SNR. The dark lines of best fit are constrained to intersect the origin. Table 2 reports these slopes ↵ and the reciprocal, the noise equivalent mixing ratio length (NEMRL), defined here as 1 standard deviation above the background. The first three rows show a CIBR method, a classical matched filter based on the transmission spectrum, but applied columnwise, and a more traditional matched filter with rectangular support but using the Jacobian spectrum. The fourth row shows the columnwise Jacobian matched filter. This combination achieves a NEMRL of 140 ppm m the best overall performance of any algorithm.
We have defined sensitivities in terms of ppm m at standard temperature and pressure (STP), which is a reasonable approximation for the Kern fields. However, one could also express sensitivity as a fractional enhancement of the entire atmospheric column. To calculate this value, we average the seasonal and geographic CH 4 profiles of Matricardi (2003), which incorporate tropospheric models from Müller and Brasseur (1995) and Clerbaux et al. (1998). The resulting column has an STP-equivalent mixing ratio length of 14 600 ppm m for which the noise-equivalent enhancement is 1 %. A larger total-column enhancement of 3.4 % would be equivalent to 500 ppm m, producing a strong > 3 detection.
The SNRs reported here understate the effective system sensitivity when paired with a trained operator. SNR estimates assume pixel independence, but in practice plumes show multiple contiguous pixels, none of which need to be above the standard 3-5 detection limit to distinguish the structure as a whole. Consequently plumes with strengths near the noise floor were visually apparent. The SNR disregards contaminating clutter, such as painted structures with similar absorption features, which could cause localized false positives. However, such features were rare; they occurred at most a few times per flight line. Moreover, they were generally easy to ignore by eye because their sizes and morphologies were very distinct from plumes. Consequently false positives were not a significant problem during actual use.   Figure 12. Typical ratio of radiances in the plume and out of plume. The ratio is continuum-removed, with an offset and vertical scaling to highlight the similar shapes.
We used several methods to verify that detections were actually caused by CH 4 , and not false positives due to interfering surface features or gases. First, a visual assessment verified that the source location of each plume lay near artificial structures, as expected, and that its main structure was uncorrelated with visible surface features or changes in albedo. In other words, the plumes were contiguous phenomena that crossed rather than followed the boundaries of surface features. Second, for a subset of plumes, we verified that the ratio of in-plume and background radiances showed a clear CH 4 signature. Figure 12 illustrates this process for the plume in Fig. 9. The figure shows the ratio of the average radiance spectra, further divided by the linear continuum stretching across the two endpoints of the spectral interval. The plot compares this ratio to the modeled transmittance of CH 4 . The model fits three parameters consisting of a linear continuum and the absorption given by the mixing ratio length. The retrieved path length of 2438 ppm m is a close match to the linear solution from the matched filter, and the resulting spectrum is a good fit to the empirical ratio. This provides additional confidence that the detection was due to the actual presence of CH 4 and not (for example) a false positive artifact. Motivated by the use of steam injection for enhanced oil recovery at Bakersfield, we checked the ratios for H 2 O vapor absorption and failed to find any excess concentration of this potential interferant. Finally, corroborative Atmos. Meas. Tech., 8, 4383-4397, 2015 www.atmos-meas-tech.net/8/4383/2015/ measurements by other COMEX in situ and remote sensing instruments confirmed the presence of CH 4 , as illustrated in Fig. 8.

Discussion
These experiments underscore the value of merging diverse instruments with partial overlap in measurement capabilities. This overlap permits the instruments to cross-check each other, and allows fast mapping platforms to provide reconnaissance to in situ teams. For example, COMEX used a multi-tier strategy of nested measurement scales. At the top level, MAMAP provided high accuracy retrievals over very wide areas, and resolved CH 4 plumes on the scale of hundreds of meters. Airborne in situ measurements provided validation of remotely observed large scale plumes. AVIRIS-NG provided unambiguous images of specific plume sources, locating them to within a few meters. The ground team augmented this remote view with in situ point samples. The arrangement proved effective, but relied on a relatively slowmoving ground team to close the loop with unambiguous in situ measurements. Tactical remote measurement with efficient airborne reconnaissance makes full use of these in situ instruments.
There are several logical next steps for both hardware and software development. For hardware, a future instrument designed specifically for the purposes of CH 4 would offer far better detection performance. For example, an imaging spectrometer that focused just on the 2.1-2.4 µm range would permit much finer spectral resolution for highly sensitive and accurate retrievals. A slower, low-flying platform would provide higher spatial resolution and improved SNR. On software, COMEX revealed a need for a real-time geolocalization capability to display accurate plume locations. This capability has since been implemented using a geometric camera model, geometric ray-tracing, and an onboard digital elevation model (DEM). One could further mitigate clutterrelated noise with more sophisticated background modeling (Funk et al., 2001;Thorpe et al., 2013) or explicit outlier rejection (DiPietro et al., 2012;Boardman and Kruse, 2011). Finally, the system could be made more general by pairing it with the real-time reflectance processing (Bue et al., 2015) to recognize materials at the surface in addition to gaseous absorbers. This would expand missions to applications including search and rescue, tracking of spatially variable phytoplankton or algal blooms for water color studies, fire response, and oil spill response.
During COMEX, communicating tactical science data across platforms improved science outcomes and produced a robust data set that could be used to validate interpretations. Advanced communications infrastructure can further streamline data transfer between ground and flight. The NASA Airborne Science Data and Telemetry system  functions as a flight data recorder, onboard payload network server, and low bandwidth telemetry system. A Payload Telemetry Link Module augments this system with higher level data products and high bandwidth satellite telemetry. Similar systems can be used to more seamlessly share data across participants. As the number of realtime data products increases, it becomes important to avoid information overload. Displaying data as optional map layers provides user control, while chat capabilities allow different team members to flag important observations.

Conclusions
We have demonstrated tactical signature detection onboard AVIRIS-NG in service of a coordinated science campaign, COMEX. Informing the operators when CH 4 plumes were detected improved the science yield of this investigation, both by confirming data quality in real time and by enabling more flexible asset deployment. Post-analysis demonstrated that the strongest detected signatures were related to real sources. The latest iteration of the system provides 3 detection sensitivity of 500 ppm m or weaker.

Appendix A: Wavelength calibration
Accurate wavelength calibration is critical for detecting narrow spectral absorption features. Our wavelength calibration uses laser sources to characterize the center wavelengths and and full width half maximum (FWHM) of the detector array. This initial wavelength calibration is derived from six sources (at 0.4067, 0.532, 0.632, 1.064, 1.551, and 2.064 µm). We refine this result with flight data, optimizing a single uniform shift to match atmospheric absorption features in the top of atmosphere (TOA) reflectance spectrum ⇢(c, ). Following Gao et al. (1993b), the TOA reflectance represents the radiance measurements at wavelengths , after normalizing for extra-terrestrial solar irradiance F and solar zenith ✓ : ⇢(c, ) = ⇡L m (c, ) F ( ) cos(✓) .
We model this spectrum as a locally linear continuum attenuated by gaseous absorption of the 0.74 µm oxygen band and the 0.96 and 1.14 µm water vapor bands. The attenuation is governed by a Beer-Lambert relation based on the gas absorption coefficient ( ) obtained from a 20 layer model atmosphere: where ⌧ represents free parameters we optimize during the fit. Specifically, h(⌧ 1 ) is convolution with a Gaussian spectral response function with a FWHM given by ⌧ 1 . ⌧ 2 is the continuum level at 100 %, ⌧ 3 represents the absorption path length, ⌧ 4 the wavelength shift and ⌧ 5 a linear slope. We fit these free parameters using a Nelder-Mead simplex algorithm. Figure A1 (Right) shows an example from a bright, spectrally smooth playa. Here the model matches the measured spectrum with residual error under 1 % -within the limits of the spectrometer's radiometric accuracy. The empirical calibration procedure is an independent check of crosstrack spectral uniformity. Figure A2 shows the wavelength calibration shift for different cross-track elements, after averaging 500 downtrack samples in each column. The average shift is less than 0.1 nm, or 2 % of the full width at half maximum.