Thermokarst lakes are important emitters of methane, a potent greenhouse
gas. However, accurate estimation of methane flux from thermokarst lakes is
difficult due to their remoteness and observational challenges associated
with the heterogeneous nature of ebullition. We used high-resolution (9–11 cm)
snow-free aerial images of an interior Alaskan thermokarst lake acquired
2 and 4 days following freeze-up in 2011 and 2012, respectively, to detect
and characterize methane ebullition seeps and to estimate whole-lake
ebullition. Bubbles impeded by the lake ice sheet form distinct white
patches as a function of bubbling when lake ice grows downward and around
them, trapping the gas in the ice. Our aerial imagery thus captured a
snapshot of bubbles trapped in lake ice during the ebullition events that
occurred before the image acquisition. Image analysis showed that low-flux
A- and B-type seeps are associated with low brightness patches and are
statistically distinct from high-flux C-type and hotspot seeps associated
with high brightness patches. Mean whole-lake ebullition based on optical
image analysis in combination with bubble-trap flux measurements was
estimated to be 174
Soils in the northern permafrost region contain 1300–1370 Pg organic
carbon with an uncertainty range of 930–1690 Pg (Hugelius et al., 2014). A
large amount of soil carbon in the Yedoma permafrost region (
Thermokarst lakes are a prominent landscape feature in the high northern latitudes (Smith et al., 2007; Grosse et al., 2013). They are formed in closed depressions following the thawing of ice-rich permafrost or melting of massive ice. Once initiated, the presence of a water body on permafrost serves as a positive feedback to permafrost degradation. Depending on the amount of excess ice content in permafrost, this positive feedback accelerates the growth of thermokarst lakes in both lateral and vertical directions (Jorgenson and Shur, 2007; Plug and West, 2009; Kokelj and Jorgenson, 2013). Over many years, taliks (thaw bulbs) of perennially thawed soil develop beneath thermokarst lakes (Hinzman et al., 2005; West and Plug, 2008; Rowland et al., 2011), creating conditions favorable for year-round methane production through anaerobic decomposition of organic matter by microbes (Zimov et al., 1997; Walter et al., 2006, 2008; Kessler et al., 2012). During lateral expansion, thermal erosion along the lake margin also releases both Holocene and Pleistocene organic matter from adjacent soils into anaerobic lake bottoms, further enhancing methanogenesis (Zimov et al., 1997; Walter Anthony et al., 2014).
Ebullition (bubbling) is considered the dominant pathway of methane release
from lakes to the atmosphere (Keller and Stallard, 1994; Bastviken et al.,
2011). Methane produced in dense lake sediments and thaw bulbs emerges
primarily through intrasedimentary bubble tubes as point-source seeps on the
lake bed (Walter Anthony et al., 2010). In the high northern latitude
region, where lake surfaces freeze throughout the winter, most bubbles
emerging from the lake bed ascend through the water column and get trapped
by ice as gas pockets (Walter et al., 2008; Greene et al., 2014). Ongoing
ice growth can separate ice-trapped bubbles from an individual seep by thin
films of ice, resulting in vertically oriented bubble columns in the ice.
Walter et al. (2006) took advantage of this phenomenon to reveal locations
and relative strength of point sources of methane seep ebullition across
lake ice. They identified four major types of methane ebullition seeps based
on ice-bubble cluster morphology and they measured daily mean ebullition
rates (mean
Photos show four distinct patterns of point-source ebullition
seeps in early winter lake ice:
Accounting for methane ebullition from northern thermokarst lakes can significantly improve estimates of lake contributions to regional and global atmospheric carbon budgets (Walter et al., 2007; Bastviken et al., 2011). However, due to challenges associated with the logistics of fieldwork in remote locations as well as spatial and temporal heterogeneity of ebullition, accurate estimation of methane flux from thermokarst lakes is difficult (Casper et al., 2000; Bastviken, 2004; Wik et al., 2011). Most studies have been carried out using field measurements to understand the spatial and temporal variability of methane ebullition. However, the insufficiency of field data is a recurring issue since it is difficult to sample the entire lake area, particularly when lakes have remote locations. This may lead to an unrealistic characterization of variability of ebullition bubbles and a less accurate estimation of methane flux at a regional scale. Recently, Walter Anthony and Anthony (2013) combined point-process modeling with field-measured data to understand the drivers of ebullition spatial variability in thermokarst lakes and provided ways to reduce uncertainty in regional-scale lake ebullition estimates based on limited field data; nonetheless spatially limited field sampling remains a hindrance to whole-lake ebullition quantification.
Remote sensing methods combined with field observations can help overcome
some of the limitations that exist in a sole field-survey method. One of the
major advantages of remote sensing tools is that they may provide the
possibility to map the entire population of methane ebullition bubbles on a
lake. Moreover, remote sensing can overcome the logistical difficulties that
exist in accessing methane-bubbling lakes in the remote regions of the
Arctic and Subarctic. Walter et al. (2008) and Engram et al. (2012)
demonstrated the potential application of synthetic aperture radar (SAR) satellite imagery to estimate
whole-lake ebullition from spatially limited field measurements of
ebullition along survey transects. These studies showed correlation of radar
backscatter values with the percent of the surface area of lake ice covered with
bubbles and field-measured methane ebullition rates based on bubble-trap
measurements from lakes. Additionally, Walter Anthony et al. (2012) used
aerial surveys to identify, photograph, and map large (
Since open holes induced by bubbling in snow-covered lake ice are visually distinct, and since lower-flux ebullition bubble clusters trapped inside ice appear as bright white features that have a strong contrast against dark, bubble-, and snow-free congelation ice (Fig. 1), there is the potential and need to detect and quantify methane bubbles with optical remote sensing. In this study, we explored high-resolution optical remote sensing images to characterize methane ebullition seeps on Goldstream Lake, an interior Alaska thermokarst lake, and to study their spatiotemporal dynamics.
Goldstream Lake (64.91
Ebullition accounts for a total of 83 % of methane emission from Goldstream Lake (Greene et al., 2014). The concentration of methane in Goldstream Lake's bubbles is 82–89 % (Greene et al., 2014). Ice formation on the lake usually occurs between the end of September and mid-October, reaches maximum thickness by mid-March, and ice break-up occurs around the end of April or early May. Vertically oriented layers of methane ebullition bubbles (Fig. 1), representing point-source seeps, are widespread in the lake ice, particularly along the eastern margin (Walter Anthony and Anthony, 2013). Many hotspot seeps are also found near the eastern eroding shore and are seen as open holes in lake ice during early winter and spring.
We used three sets of data in our study: (1) high-resolution snow-free early winter lake ice images from fall 2011 and 2012, (2) high-resolution snow-covered early winter lake images from fall 2012, and (3) field-based surveys of ebullition bubble seep types and their fluxes.
We first processed the aerial images of snow-free early winter lake ice acquired in the fall of 2011 and 2012 and then mapped ebullition bubbles trapped in the ice. We refer to the bubble features seen in our snow-free images as “bubble patches” henceforth since the image resolution was not sufficient to fully resolve smaller individual bubbles (Fig. 1). Then we characterized imaged bubble patches based on field-collected ebullition bubble seep data collected approximately 1–2 weeks after image acquisition when ice was safe to walk on in the fall of 2011 and 2012 and again in spring of the following year. We hypothesized that the brightness of bubble patches correlates with the strength of methane flux associated with four classes of ebullition bubble seeps (A, B, C, and hotspot) identified by Walter Anthony et al. (2010). We estimated from aerial photos the bubble-patch density for each seep class and the mean whole-lake seep ebullition. Finally, we examined the spatial patterns of seep locations in the lake with respect to eroding thermokarst shores, and analyzed interannual variability of seep occurrences by comparing imagery from different years.
Due to similar spectral characteristics of open water and dark lake ice, open-hole hotspots are difficult to map on snow-free lake ice images; but they are easily identifiable on snow-covered lake ice images. Therefore, we also collected high-resolution snow-covered aerial images in the fall of 2012 to map open-hole hotspots on the lake. We compared the locations of hotspots in aerial images with maps of hotspot locations determined by field measurements over multiple years to assess hotspot regularity.
We scheduled low-altitude, high-resolution aerial image acquisitions to map
and characterize methane ebullition bubble patches (A, B, C, and ice-covered
hotspots) during a narrow time window in the early winter, when first ice
had formed but was still snow-free. Images were acquired in nadir with a
Navion L17a plane using a Nikon D300 camera system mounted in a belly port on
14 October 2011 and 13 October 2012, 2 and 4 days following freeze-up,
respectively. Flight altitude for the acquisitions was
We collected images of the snow-covered lake in fall 14 October 2012 using an unmanned aerial vehicle (UAV) mounted with an Aptina MT9P031 board camera to map open-hole hotspot seep locations. The images were acquired from a flying height of approximately 230 m a.s.l., corresponding to an image scale 1 : 30 000 and GSD of 6 cm. All the images consisted of three visible bands: red, green, and blue (RGB).
We collected the point location data of (1) the lake perimeter and permanently installed reference markers as ground control points (GCPs) to perform rectification of aerial images and (2) methane ebullition seeps on Goldstream Lake in the fall and spring of 2011 and 2012 using a survey-grade LEICA VIVA™ real-time kinematic Differential Global Positioning System (DGPS) with centimeter accuracy.
We surveyed the lake perimeter and measured several permanently installed
reference markers. We conducted detailed ebullition ice-bubble surveys in
October 2011 2 weeks after image acquisition when lake ice was safe to
walk on. The surveys were performed within two large polygons that are
identified in Fig. 1f: one about 7 m from the eastern
thermokarst shore and a second near the center of the lake. The surveyed
polygons in the east and center of the lake covered
In October 2012, we performed bubble surveys 6 days after image acquisition
in three other polygons (total area
While field-based estimations of A-, B-, and C-type seeps were limited to survey plots covering about 13 % of the lake area, the locations of hotspot seeps were mapped across the whole lake using detailed DGPS surveys of open holes in October and April 2011 and 2012. Hotspots were detected visually at these times of the year as open water holes in lake ice.
In April 2013, we extracted several blocks of the full lake ice column at seep locations to investigate the temporal ebullition patterns that developed throughout the winter season.
We conducted the following image pre-processing. (1) We performed mosaicking of multiple images of Goldstream Lake to construct a complete image of the lake. This was achieved by using Agisoft PhotoScan Professional Software™ Version 0.9.0. (2) We then performed geometric image rectification with 22 DGPS-collected GCPs using a second-order polynomial transformation with bilinear resampling. The GCPs were distributed mostly around the lake perimeter. Some of them were identifiable reference points on the lake such as cattail vegetation and a LiCor methane analyzer installed on a platform on the lake. (3) For image enhancement we applied a feature linear transformation on all three visible spectral bands of the lake images using unstandardized principal component analysis (PCA). Both geometric and spectral image transformations were performed in ENVI™ image processing software, Version 4.8. PCA spectral transformation produced three independent principal component (PC) bands. The first band (PC 1 band) consisted of the variables that explained the most variance (> 98 %) in the data set attributing to bubble patches (Fig. 2, Supplement Sect. S1).
Panel
We applied a classification technique based on object-based image analysis (OBIA) to semi-automatically identify and map methane ebullition bubble patches in the PCA-transformed images using eCognition Developer™ 8 (Lindgren et al., 2015). Our object-based classification method comprised two steps: (1) image segmentation, i.e., aggregation of homogenous image pixels based on their spatial and spectral homogeneity into meaningful clusters known as image objects, and (2) classification of image objects (Navulur, 2007; Blaschke and Strobl, 2001). Varying ice conditions on the lake such as (a) clear, dark congelation ice, (b) milky white snow–ice, and (c) ice with shadows from neighboring trees added challenges to identifying ebullition bubble patches. We were able to resolve these challenges by integrating semantic information associated with image objects in classification (Lindgren et al., 2015). For this, the scene is first decomposed into meaningful regions that represent different areas of lakes such as vegetation and shadow, dark, and white ice. These regions are then organized in a conceptual image object hierarchy, creating a semantic network between different sized image objects – large-scale objects in the upper level called “super-objects”, and small-scale objects in the lower level called “sub-objects” (Supplement Fig. S1; Lindgren et al., 2015). For example, the lake area is a super-object composed of sub-objects associated with various lake ice characteristics (e.g., shadow, dark black ice), whereas areas of specific lake ice characteristics are super-objects of our final target feature, ebullition bubble patches. At each level, image segmentation and classification are performed to delineate and label target regions. For example in the first level, segmentation is performed on the whole lake image to identify lake shore and lake. In the second level, only the lake region is segmented and image objects derived from the lake are classified into different lake ice characteristics. The process continued as it proceeded towards lower and finer classification levels until bubble patches were identified in the lake ice in the final stage. This approach of detecting image objects from coarser to finer scale has been described as an effective way to classify images in OBIA (Blaschke et al., 2008).
More information on this hierarchical approach of bubble-patch identification can be found in the Supplement (Supplement Sect. S2) and in Lindgren et al. (2015).
In images acquired after the first snow fall, hotspots appeared dark against the snow-covered lake. Hence, they were mapped using a simple contrast and split segmentation technique in eCognition Developer™ (Definiens, 2007). This approach involves choosing a threshold value on the RGB image bands for the algorithm to maximize the contrast between hotspots and snow-covered lake pixels that separates the image content into dark objects (consisting of pixels below the threshold, i.e., hotspots) and bright objects (consisting of pixels above the threshold, i.e., snow-covered lake ice).
We extracted PC 1 gray values of individual ebullition bubble patches mapped in images from the year 2011 and 2012. PC 1 values of the image ranged between 0 and 255. Bubble patches are visually bright in true-color composite (RGB composite) images but appeared darker (i.e., low PC 1 values) than the surrounding lake ice in the PC 1 band (Supplement Fig. S2). Therefore, we inverted PC 1 values of bubble patches to make it visually intuitive, i.e., bubbles that appeared bright in natural color composite also appeared bright in the PC 1 band (Fig. 2). Henceforth, we refer to this brightness obtained in the inverted PC 1 as “PC 1 brightness”.
We assessed the relationship of ebullition bubble-patch PC 1 brightness values with four distinct types of ebullition seeps that we identified during our field surveys. We performed an analysis of variance (ANOVA) to test the null hypothesis that the mean PC 1 values (and thus true bubble brightnesses via its inverse relationship with the PC 1) of four types of seeps are not significantly different. We applied a post-hoc Tukey's honest significant difference (HSD) test, in case the null hypothesis was rejected, to identify significantly distinct seeps. Results of this analysis are shown in Sect. 4.1.
We applied a supervised classification using a maximum-likelihood classifier (MLC) on the three original visible bands and the extracted PC 1 band to classify mapped bubble patches into four distinct seep classes. The MLC calculates a Bayesian probability function from the input training classes and then assigns each pixel in the image to the class of highest membership probability (Mather, 2009). We collected 98 random samples, 35 for training and 63 for validation, on the 2011 image and similarly 181 random samples, 50 for training and 131 for validation, on the 2012 image. The samples were located at seep locations identified using field-collected DGPS data points.
The MLC approach categorized bubble patches solely based on the pixel spectral characteristics, i.e., only using the brightness values of the training samples. Since the size of bubble patches is also an additional important indicator of seep class and methane flux (Walter Anthony et al., 2010), in a subsequent step, we further refined our classification results by integrating size as an additional feature to more accurately assign bubble patches with a seep type (Supplement Sect. S3, Table S1). Finally, we estimated the seep density and mean whole-lake ebullition rate by assigning the mean long-term flux values for seep types provided by Walter Anthony and Anthony (2013) to our classified bubble patches.
We studied the spatial distribution of ebullition bubble patches as a function of distance from the eroding eastern thermokarst shore. For this, we divided the lake area into multiple 5 m wide zones starting from the eastern eroding margin as mapped in a 1949 aerial image (Fig. 1). Lake zones were created on both sides of the 1949 lake margin to cover the present-day lake area. We calculated the percent of lake ice area covered with ebullition bubble patches for each zone and then analyzed its relationship to the distance from the eastern shore lines of the lake observed in 1949, 1978, and 2012.
We evaluated the multi-temporal (year 2011 and 2012) variability of ebullition bubble patches and assessed their regularities in space and time. We utilized a marked point-process model to analyze spatial seep patterns in our multi-year bubble-patch data set derived from the images. Point-process modeling was performed on a set of bubble-patch centroids with their respective location and year information, which served as a marked point data set for the model, to derive and test the spatial characteristics of bubble-patch distribution against a null hypothesis based on complete spatial randomness. The null hypothesis suggests that the bubble patches are the results of a spatially random process over the study area and thus the difference of spatial patterns between years is random, i.e., the locations of bubble patches are independent when comparing between years (Bivand et al., 2008). For this, we generated a multi-type nearest-neighbor distance function derived from the locations of the bubbles mapped in the images using Gcross from the spatstat statistical package in R (Baddeley and Turner, 2005). Gcross first determines clustering parameters for the data set in the first year. These clustering parameters are then used to model the expected number of the second-year point given a certain distance from the first-year points if the second-year point placement is random relative to the first-year point placement. Based on the deviation between observed empirical value and expected theoretical value estimated by the model, we determined the stability of seep locations between 2011 and 2012. Similarly, we performed the multi-type nearest-neighbor distance function analysis using Gcross on the field data set of hotspot locations collected in year 2011 and 2012 to check regularity of hotspots.
We also considered that the centroid of a bubble patch, representing an ebullition bubble-patch point location, could move from one year to another due to changes in the shape and size of a bubble patch or changes in bubble-tube configuration in the sediment. We compared the overlap area between ebullition patches mapped in 2011 and 2012 images. If some area of a 2011 bubble patch appeared within the area of a 2012 bubble patch or vice versa, then we considered the bubble patch to be stable in location (i.e., reappearing). We assumed that the overlapping bubble patches originated from the same point-source seep. We checked location stability among four classes of overlapping patches that were defined by setting thresholds on area overlap: “all overlapping bubble patches”, “more than 25 % area overlap”, “more than 50 % area overlap”, and “more than 75 % area overlap”.
We used a map of open-hole hotspot seeps derived from UAV images to compare the frequency of hotspots with hotspot occurrences observed during multiple years of fieldwork by Greene et al. (2014).
We found that PC 1 brightness values of bubble patches correlated with the
strength of field-measured methane flux of ebullition seeps (A, B, C and
hotspot seeps). The lowest mean PC 1 brightness belonged to A-type seeps
followed by B-type seeps in both 2011 and 2012 (Fig. 3). The highest mean
PC 1 brightness was demonstrated by hotspots in 2011 but C-type seeps had
slightly higher mean brightness in 2012 (Fig. 3). Our ANOVA test rejected
the null hypothesis that the mean PC 1 values of the different seep types
are the same suggesting significant distinctions between mean PC 1
brightness values of different seep classes. Further post-hoc analysis using
Tukey's HSD test demonstrated that C- and A-type, hotspot and A-type, and
hotspot and B-type seeps are significantly distinct based on their mean PC 1
with
An absolute discrimination of individual seep type based on brightness was not supported by the post-hoc tests due to overlapping brightness ranges between different seep types (Fig. 3). This is likely because ebullition is episodic with varying bubbling rates over time and because individual low-flux methane seeps were not resolved given the spatial resolution of the image. A possible explanation for low PC 1 brightness of some hotspots is that fresh thin nighttime ice temporarily covered some hotspots on the image acquisition day, allowing the formation of few small white gas bubbles, while much of the remaining gas escaped through cracks in the thin ice, resulting in low true brightness for these high-flux seeps. We have observed this phenomenon on several occasions during our field visits in early winter and spring, particularly on days when temperatures stayed low and hotspots were covered with a few millimeters of ice with small bubbles beneath (Fig. 4); these hotspots usually open up when atmospheric temperature rises again during the day. Conversely, hotspots that remained open could not be identified in our snow-free lake ice imagery due to spectral similarities between open water and clear black ice (Fig. 1). ANOVA analysis was only performed on ice-covered hotspots.
Box plots are shown of PC 1 brightness values for bubble patches with different classes of seeps in 2011
Panels
We found that a large number of A-type seeps clustered together were not mapped as individual bubble patches but rather as a single large bubble patch. A-type seeps and high-flux seeps that were close together were also mapped in a single feature associated with a brighter bubble patch. Therefore, some A-type seeps showed high PC 1 brightness values. Similar to A-type seeps, occasionally individual B-type seeps were also not distinct. In a time series analysis of bubbling frequency by A- and B-type seeps, Walter Anthony et al. (2010) showed that bubbling from these shallow-sourced seeps is highly seasonal. Bubbling rates are high in summer when surface sediments are warmer, and low in winter when sediments cool down. Bubble traps left in place over these seep types year round revealed that low-flux seeps can have periods of no bubbling for up to several months. Ice blocks harvested by us in spring over seeps marked as A-type seeps in October confirm this pattern (Supplement Fig. S3). It is very likely that A- and B-type seep conduits were present in the sediments, but not actively bubbling during the 2- and 4-day periods after ice formation captured by the 2011 and 2012 imagery. Thus they did not appear under the given spatial resolution of the image and its specific acquisition time. Also, bubble traps placed over C-type seeps year round revealed that these seeps can also undergo long periods (weeks to months) of no bubbling, but when they bubble, the bubbling rates are usually very high (Walter Anthony et al., 2010). This intermittent flux behavior probably contributed to some discrepancies in the relationship between bubble-patch brightness derived from images that captured a snapshot of ebullition activity and methane flux values of seeps estimated from long-term field observations (Table 1).
Seep density and estimated mean whole-lake ebullition flux of Goldstream Lake, Fairbanks, Alaska, was derived from 2011 and 2012 optical aerial image data set and from ground surveys. The ground survey estimates are from previously published studies by Walter Anthony and Anthony (2013) and Greene et al. (2014) based on ground surveys conducted over multiple years (2007–2011) at Goldstream Lake.
In other parts of Goldstream Lake, especially along the eastern shore, we found large patches of many individual ebullition bubbles (each typically 3 to 10 mm diameter) that formed large diffuse patches rather than clustering as tightly packed bubbles the way A, B, C, and hotspot seeps do. In our optical images, these ebullition bubbles appeared as irregular patches of fuzzy, white-colored bright ice with some bright regular bubble spots (Fig. 4). Therefore, the brightness values corresponding to the surrounding diffuse patches were assigned to other seeps, particularly to low-flux seeps that were within the patch and had not expressed completely when the images were acquired. Until recently, these millimeter-scale ebullition bubbles were only recorded in transect-survey data as “tiny-type” seep but never assigned a mean daily flux value or included in whole-lake ebullition estimates due to a lack of associated flux field measurements. Recent flux measurements made continuously year round with submerged bubble traps on the tiny-type seep class in Goldstream Lake and other lakes suggest that fluxes from these seeps may also be important (Walter Anthony et al., unpublished). Analysis of bubbles collected with bubble traps placed over tiny-type seeps revealed that these bubbles were 60–80 % methane in volume (Walter Anthony et al., unpublished). When we extracted an ice block in spring 2013, we observed that tiny-type ebullition had been frequent throughout winter, resulting in long, vertically oriented stacks of tiny ebullition bubbles trapped in ice (Fig. 4).
The overall MLC classification accuracy for differentiating seep types was
Generally higher densities of A-type seeps (and also slightly in B- and C-type seeps) in ground surveys (Walter Anthony and Anthony, 2013; Greene et al., 2014) compared to aerial images (Table 1) can be explained by the time in which observations were made and image resolution. Results reported in Walter Anthony and Anthony (2013) and Greene et al. (2014) are based on ground surveys conducted over multiple years (2007–2011) at Goldstream Lake usually 1 to 2 weeks following freeze-up when ice was safe to walk on (Walter Anthony et al., 2010). Since our aerial surveys were conducted only 2–4 days after ice formation, and the frequency of bubbling events from A-type seeps is often weeks to months in winter, it is not surprising that the field surveys several weeks after ice formation capture an order of magnitude more A-type seep bubbles. Additionally, it is very likely that some active A-type seeps that occurred in very small patches were not distinct under the given resolution of the aerial images. Relatively more frequent bubbling in B- and C-type seeps allows for similar seep density values between ground surveys and aerial images; however, as expected, the 2012 seep densities are closer to the ground-ice survey values due to (a) more time since freeze-up and (b) a much higher barometric pressure drop preceding the aerial image acquisition in October 2012 compared to October 2011. It is well established that ebullition dynamics are related to changes in barometric pressure (Mattson and Likens, 1990; Fechner-Levy and Hemond, 1996; Scandella et al., 2011).
The comparison of hotspot densities in optical images vs. ground surveys in
Table 1 also shows the expected pattern. The ground-survey data of hotspots
reflect multiple years of surveys of whole-lake hotspots when ice is thick
enough to safely walk on. When ice is very thin a few days after freeze-up,
more open holes are present on the lake and classified as hotspot seeps in
aerial images. A week or more later, many holes freeze over and will be
classified as C-type seeps in ground surveys. This could have also led to a
high classification error for C-type seeps. The total density of C-type and
hotspot seeps combined remains consistent (
Our image-based analysis shows the whole-lake flux to be 174
High methane production in response to thermokarst activity on the
Goldstream Lake is evident from the distribution pattern of ebullition bubble
patches at the eroding margins in different years. We found a strong inverse
relationship (
Panel
We observed fewer ebullition bubble patches in the center of the lake, which we interpret as a sign that labile Pleistocene-aged organic carbon in the talik under this area has been largely depleted, and unlike at the edge along the active erosion margin, there is no significant additional accumulation of ancient labile carbon in the lake center (Brosius et al., 2012). Radiocarbon dating of bubble patches found in the lake center showed that these seeps originate from Holocene-aged and more recent organic matter that is found in the upper lake sediments (Brosius et al., 2012). Generally, methane bubbling was the lowest along the 1949 eastern lake margin and the highest along the 2012 eastern lake margin (Fig. 5), indicating that depletion of labile carbon progressed since these areas were included in the lake, and the active thermo-erosion margin migrated eastward. This shows that optical remote sensing is a powerful tool to understand the spatial variability of methane ebullition on thermokarst lakes.
A comparison of bubble patches visible in thin lake ice 2 days after
freeze-up in October 2011 (left-side images) and 4 days after freeze-up in
October 2012 (right side images) is shown. Image pairs in
We observed four possible characteristics of bubble patch dynamics in our
images (Fig. 6). (i) Bubble patches may move horizontally; (ii) bubble
patches do not maintain the same morphology between years (e.g., single
bubble patches re-appear in a cluster of multiple patches the next year or
vice-versa); (iii) bubble patches appear in an image in one year and not
another; and (iv) bubble patches maintain the location and shape but patch
size is different between the years. It is important to note that these
observations are made during the two very short windows of time 2–4 days
after freeze-up. Our analysis does not take into account the changes in
long-term bubble patch morphology. Hence, it is important to highlight that
the characteristics of bubble patches are driven by the dynamics of bubble
formation and transport, hydrostatic pressure, and ice growth. Other changes
in the characteristics of bubble patches could be because of evolution of
point sources or changes in point-source conduits (bubble tubes) in the
sediment (Walter et al., 2008; Scandella et al., 2011). Atmospheric
pressure dynamics can also strongly impact bubbling over short timescales,
resulting in different ice-bubble patterns one year from the next if
insufficient time passes to allow all seeps to be expressed in the lake ice
cover. Field measurements have shown that ebullition is inversely related to
changes in hydrostatic pressure (Mattson and Likens, 1990; Varadharajan,
2009; Casper et al., 2000; Glaser et al., 2004; Tokida et al., 2005;
Scandella et al., 2011). A significant air pressure drop during the week
preceding image acquisition in October 2012 may have allowed methane that
previously accumulated in the sediment during high-pressure days to rise up
into the water column, manifesting itself as larger numbers of bubbles (Fig. 7),
and larger and brighter bubble patches in the lake ice (Fig. 6).
Conversely, air pressure change in October 2011 was not large enough to
enhance ebullition before the image was acquired. As a result, bubble patch
density was 55 % higher in 2012 (0.185 m
We rejected the null hypothesis of complete spatial randomness to show that the difference in spatial patterns of bubble patches and hotspots between years is not random, i.e., that the locations of seeps in the years 2011 and 2012 are not independent. The Gcross distribution function showed that a statistically significant number of second-year bubble-patch center points are less than 2 m away from the first-year center points and that there are far fewer than expected that are 3 m or more apart (Fig. 8a). For the hotspots, a statistically significant number of seeps moved less than a meter (Fig. 8b). Since our image rectification accounted for a geolocation error of less than 20 cm and DGPS geolocation error is even smaller and considered negligible, we conclude that the seep locations are consistent between years 2011 and 2012.
Cumulative distribution function of distances (
Based on our DGPS data, the number of hotspots was relatively stable among the various surveys with about 105 hotspots for the whole lake as the average of various measurements during different years and spring and fall field seasons (Greene et al., 2014). UAV-based aerial images taken 5 days after ice formation when snow covered the lake also demonstrated close agreement with the hotspot seep numbers and locations. We were able to identify 78 dark open water holes in the white, snow-covered UAV lake image acquired in early winter of 2012. Among these 78 locations there was a total of about 95–100 active open-hole hotspot seeps since some large, irregularly shaped holes consisted of multiple, coalesced holes produced by hotspot seeps of close proximity (Supplement Fig. S4).
When we compared the location of bubble patches in 2011 and 2012, we found that 47.2 % of a total of 1195 ebullition bubble patches mapped in 2011 reappeared in 2012, which is 35.7 % of a total of 1860 ebullition bubble patches mapped in 2012. We found that 37.5, 30, and 17.7 % of bubble patches mapped in 2011 reappeared in 2012 with an overlap area of “more than 25 %”, “more than 50 %”, and “more than 75 % area”, respectively. We expect that if more time passed between the time of freeze-up and aerial image acquisition date we would see an even higher percentage of seep location re-occurrences because more seeps would be actively expressed.
We also observed a relationship between bubble-patch brightness and location stability of bubble patches. Very bright patches in 2012 seemed to appear at locations where bubble patches were already observed in 2011. This could indicate locations of high-flux seeps where methane was able to rise through the sediment even under relatively high hydrostatic pressure conditions that we observed in October 2011. Based on our bubble patch classification results (Table 1), we also noticed that seep density of high-flux C- and hotspot-type seeps is less variable during our study period compared to low-flux A- and B-type seeps. However, long-term remote sensing and ground-based observations are required to further test our hypothesis of seep regularity that high-flux seeps are temporally more stable in their location than low-flux seeps. Additionally, long-term data may also help to account for the difference in pressure and look at possible changes in seep type over the years.
The regularity of bubble patches observed despite the differences in atmospheric pressure conditions following the lake freeze-up events in 2011 and 2012 as well as the location stability of hotspots indicates that a large number of point-source seeps in thermokarst lakes are stable over at least annual timescales. Walter Anthony et al. (2010) also found seeps to maintain stable locations in Goldstream Lake when submerged bubble traps were placed over individual seeps to monitor their ebullition dynamics for periods of up to 700 days. In Siberia one hotspot seep location was marked and found stable for at least 8 years (Walter Anthony et al., 2010).
We found numerous significant benefits of using aerial images for characterizing ebullition seeps on lake ice. Aerial images of early winter lake ice without snow cover allowed us to map and characterize bubble patches on the entire lake surface as well as assess their spatial distribution more accurately. Snow-covered lake ice images allowed us to map open-hole hotspots. We were able to differentiate high-methane-emitting seeps from low-methane-emitting seeps on the lake based on PC 1 brightness values of bubble patches. Image-derived estimates of seep densities by class agreed with those of field-based survey methods, except for some overestimation of hotspots and underestimation of A-type seeps. We were able to differentiate lake areas with high seep densities versus low seep densities; having this ability is especially useful for quantifying methane ebullition on larger lakes that are harder to survey extensively by foot.
Our results also imply a potential to apply high-resolution optical images at a regional scale to quantify relative methane flux from many lakes, which at a minimum should allow for the classification of high-ebullition versus low-ebullition lakes and their distribution in a region. It is important to note that while image analysis is useful to comprehensive mapping of lake ice bubbles, for the estimation of whole-lake methane emissions, this technique should be coupled with bubble-trap field measurements and laboratory measurements of methane concentration in bubbles.
But because ebullition is a temporally dynamic phenomenon, our ability to accurately identify the distinct seep type of bubble patches on a snapshot of ebullition activity during only 2 and 4 days since lake ice formation is limited. The morphology and distribution of bubbles can undergo significant changes in response to freeze/thaw cycles during winter (Jeffries et al., 2005). Furthermore, ebullition is highly controlled by the balance between atmospheric pressure and sediment strength, making it an episodic phenomenon (Varadharajan, 2009; Scandella et al., 2011). Ebullition is triggered following the fall of hydrostatic pressure or after a sufficient volume of gas is produced in the sediment that allows bubble tubes or gas conduits in lake sediments to open or dilate (Scandella et al., 2011). Bubbles previously trapped in lake sediment then break out through these open bubble tubes and rise up in the water column. Moreover, microbial activity of methane-producing bacteria is temperature-dependent. As a result, seep ebullition slows down when the lake surface sediments cool down in winter and it increases as lake sediment warms up in summer (Walter Anthony and Anthony, 2010). Therefore, discrepancies arise in estimates of the number of seeps and seep morphology derived from observations made at different times of the ice cover season (Wik et al., 2011). Ideally, optical image acquisition would occur at least several weeks following freeze-up of lakes to allow more time for seep expression in lake ice. Unfortunately, snow-free conditions several weeks after freeze-up is rare in many regions of the Arctic; and early snow cover inhibits the mapping of bubble patches with optical data.
SAR data have an advantage over optical remote sensing data in detecting methane bubbles trapped in lake ice under snow cover conditions (Walter et al., 2008; Engram et al., 2012). Engram et al. (2012) showed that particularly L-band SAR data acquired in the fall have the potential to estimate whole-lake methane ebullition since a longer-wavelength L-band is able to detect bubbles under other conditions such as the presence of snow, a thin layer of white ice, and aquatic vegetation. However, the moderate spatial resolution of current L-band SAR systems can be a limiting factor to estimate methane emission from small lakes and to capture delicate spatial patterns of ebullition seeps on lakes. SAR lake images further tend to have false backscatter signals from the lake shore (Walter et al., 2008; Engram et al., 2012), therefore limiting its usability in proximity to shores (about 1 pixel around lake shores is excluded in SAR analyses) where we showed an important component of ebullition may take place on eroding thermokarst margins. Thus, high-resolution optical images can supplement SAR-based studies by revealing the location of methane ebullition seeps and their types on the lake more precisely. Our study shows that optical high-resolution remote sensing methods have the potential, given the caveats raised above, to improve understanding of spatial and temporal variability of ebullition and therefore the dynamics of microbial processing of organic matter within an individual lake.
It is important to understand the dynamics of methane ebullition from thermokarst lakes to estimate the amount of carbon release from thawing permafrost and to evaluate its feedback to the global carbon cycle. Our study focusing on Goldstream Lake, Interior Alaska, shows that high-resolution optical remote sensing is a promising tool to map the distribution of point-source methane ebullition seeps across an entire thermokarst lake surface, a task that is difficult to achieve through field-based surveys alone. This method helps to reveal the location and relative sizes of high- and low-flux seepage zones within lakes. We also demonstrated that a large proportion of ebullition seeps in the study lake were location-stable over at least two winter seasons in the 2011–2012 observation period. Such observations may be used to indirectly characterize permafrost carbon mobilization in a lake since lake portions with greater numbers of high-flux seeps likely indicate either the presence of rapidly thawing organic-rich permafrost deposits or eroding lake margins. Our approach is also applicable to other regions and will help to characterize methane ebullition emissions from seasonally ice-covered lakes, including thermokarst and non-thermokarst lakes in tundra and boreal zones. It will help to differentiate lakes in a region based on methane emission by estimating ebullition seep density, and their relative methane flux. This differentiation could potentially be used to identify presence or absence of organic-rich permafrost deposits such as yedoma in the area. For example, yedoma-type thermokarst lakes, such as Goldstream Lake, where large amounts of labile carbon is readily available for microbes to decompose, emit more methane than non-yedoma-type thermokarst lakes. This can be a useful supplement to surveying soil carbon pools and yedoma distribution at a regional scale. Multi-temporal spatial information derived from remotely sensed optical data allows the identification of variables that control methane ebullition dynamics and spatial patterns. However, the timing of optical image acquisitions is a critical and a potentially limiting factor, with respect to both atmospheric pressure changes and snow/no-snow conditions during early lake freeze-up. Therefore, high-resolution remotely sensed optical images in combination with SAR and field data could be a very valuable tool to improve the estimation of methane emission from lakes at the regional scale.
G. Grosse and K. M. Walter Anthony conceived this study. P. R. Lindgren developed the method, performed data analysis, and wrote the manuscript with significant input from all co-authors. P. R. Lindgren, G. Grosse, and K. M. Walter Anthony were responsible for the field work.
We thank A. Bondurant for harvesting ice blocks; P. Anthony for statistical coding; and A. Strohm, M. Engram, and J. Lenz for assistance with other fieldwork. We thank J. Cherry for aerial image acquisitions and B. Crevensten and G. Walker for UAV image acquisitions. This research was funded by NASA Carbon Cycle Science, grant no. NNX11AH20G. Additional support came from NSF no. 1107892. Edited by: I. Laurion