Quantity and distribution of methane entrapped in sediments of calcareous, Alpine glacier forefields

Aside from many well-known sources, the greenhouse gas methane (CH4) was recently discovered entrapped in the sediments of Swiss Alpine glacier forefields derived from calcareous bedrock. A first study performed in one glacial catchment indicated that CH4 was ubiquitous in sediments and rocks and was largely of thermogenic origin. Here, we present the results of a follow-up study that aimed at (1) determining the occurrence and origin of sedimententrapped CH4 in other calcareous glacier forefields across Switzerland and (2) providing an inventory of this sedimententrapped CH4, i.e., determining the contents and total mass of CH4 present, and its spatial distribution within and between five different Swiss glacier forefields situated on calcareous formations of the Helvetic nappes in the Central Alps. Sediment and bedrock samples were collected at high spatial resolution from the forefields of Im Griess, Griessfirn, Griessen, Wildstrubel, and Tsanfleuron glaciers, representing different geographic and geologic regions of the Helvetic nappes. We performed geochemical analyses on gas extracted from sediments and rocks, including the determination of CH4 contents, stable carbon-isotope analyses (δCCH4), and the determination of gas-wetness ratios (ratio of CH4 to ethane and propane contents). To estimate the total mass of CH4 entrapped in glacier-forefield sediments, the total volume of sediment was determined based on the measured forefield area and either literature values of mean sediment thickness or direct depth measurements using electrical resistivity tomography. Methane was found in all sediments (0.08– 73.81 μgCH4 g−1 dry weight) and most rocks (0.06– 108.58 μgCH4 g−1) collected from the five glacier forefields, confirming that entrapped CH4 is ubiquitous in these calcareous formations. Geochemical analyses further confirmed a thermogenic origin of the entrapped CH4 (average δCCH4 of sediment of −28.23 (± 3.42) ‰; average gas-wetness ratio of 75.2 (± 48.4)). Whereas sediment-entrapped CH4 contents varied moderately within individual forefields, we noted a large, significant difference in the CH4 content and total CH4 mass (range of 200–3881 t CH4) between glacier forefields at the regional scale. The lithology and tectonic setting within the Helvetic nappes appeared to be dominant factors determining rock and sediment CH4 contents. Overall, a substantial quantity of CH4 was found to be entrapped in Swiss calcareous glacier forefields. Its potential release and subsequent fate in this environment is the subject of ongoing studies.

Abstract. Aside from many well-known sources, the greenhouse gas methane (CH 4 ) was recently discovered entrapped in the sediments of Swiss Alpine glacier forefields derived from calcareous bedrock. A first study performed in one glacial catchment indicated that CH 4 was ubiquitous in sediments and rocks and was largely of thermogenic origin. Here, we present the results of a follow-up study that aimed at (1) determining the occurrence and origin of sedimententrapped CH 4 in other calcareous glacier forefields across Switzerland and (2) providing an inventory of this sedimententrapped CH 4 , i.e., determining the contents and total mass of CH 4 present, and its spatial distribution within and between five different Swiss glacier forefields situated on calcareous formations of the Helvetic nappes in the Central Alps.
Sediment and bedrock samples were collected at high spatial resolution from the forefields of Im Griess, Griessfirn, Griessen, Wildstrubel, and Tsanfleuron glaciers, representing different geographic and geologic regions of the Helvetic nappes. We performed geochemical analyses on gas extracted from sediments and rocks, including the determination of CH 4 contents, stable carbon-isotope analyses (δ 13 C CH4 ), and the determination of gas-wetness ratios (ratio of CH 4 to ethane and propane contents). To estimate the total mass of CH 4 entrapped in glacier-forefield sediments, the total volume of sediment was determined based on the measured forefield area and either literature values of mean sediment thickness or direct depth measurements using electrical resistivity tomography.
Methane was found in all sediments (0.08-73.81 µg CH 4 g −1 dry weight) and most rocks (0.06-108.58 µg CH 4 g −1 ) collected from the five glacier forefields, confirming that entrapped CH 4 is ubiquitous in these cal-careous formations. Geochemical analyses further confirmed a thermogenic origin of the entrapped CH 4 (average δ 13 C CH4 of sediment of −28.23 (± 3.42) ‰; average gas-wetness ratio of 75.2 (± 48.4)). Whereas sediment-entrapped CH 4 contents varied moderately within individual forefields, we noted a large, significant difference in the CH 4 content and total CH 4 mass (range of 200-3881 t CH 4 ) between glacier forefields at the regional scale. The lithology and tectonic setting within the Helvetic nappes appeared to be dominant factors determining rock and sediment CH 4 contents. Overall, a substantial quantity of CH 4 was found to be entrapped in Swiss calcareous glacier forefields. Its potential release and subsequent fate in this environment is the subject of ongoing studies.
Accelerated melting of glaciers and ice sheets as a result of global warming (Haeberli et al., 2007;Paul et al., 2004;UNEP and WGMS, 2008) has prompted intense research activities in glacial environments, including investigations on their role in the turnover of greenhouse gases. Several studies have identified subglacial environments as habitats for methanogens and, consequently, as a potentially important CH 4 source (Wadham et al., , 2013Christner et al., 2012;Souchez et al., 1995;Stibal et al., 2012). Methane emissions from these locations have recently been confirmed and quantified in field measurements (Burns et al., 2018;Christiansen and Jørgensen, 2018;Lamarche-Gagnon et al., 2019). Conversely, other studies have provided evidence that aerated glacier-forefield sediments can act as a sink for atmospheric CH 4 (Chiri et al., 2015;Nauer et al., 2012;Bárcena et al., 2010;Hofmann et al., 2013). This function is mediated by a group of aerobic methane-oxidizing bacteria (MOB), which catalyze CH 4 oxidation at near-atmospheric concentrations (Curry, 2009;Zhuang et al., 2013;Dunfield, 2007). The strength of this CH 4 sink appears to vary between different glacier-forefield landforms and increases with sediment age (Bárcena et al., 2010;Chiri et al., 2017;Hofmann et al., 2013). In this context, sediment age refers to the number of years sediment has been exposed to the atmosphere following glacier retreat. Note that both terms, sediment age and landform, serve as proxies for all edaphic variations present in these sediments at different locations within the glacier forefield. We will adopt this convention and use the terms sediment age and landform in this fashion throughout the paper.
An alternative potential CH 4 source was recently detected in sediments of Swiss glacier forefields, in particular in those derived from calcareous bedrock (Nauer et al., 2012). Subsequent laboratory experiments revealed that CH 4 was released from glacier-forefield sediments upon mechanical impact and during acidification (Nauer et al., 2014). In a recent study focusing on one particular Swiss glacial catchment (Wildstrubel catchment, Canton of Valais), we established that entrapped CH 4 was present in nearly all sediment and bedrock samples collected throughout this area, but CH 4 contents exhibited substantial variation between sampling locations (Zhu et al., 2018). We also provided robust evidence based on stable-isotope and other geochemical data that CH 4 entrapped in sediment and bedrock samples was predominantly of thermogenic origin and that microbial CH 4 production was likely of minor importance at this site. However, as the focus of the previous study was on the occurrence and origin of entrapped CH 4 in different regions of the catchment, the number of samples collected was insufficient to rigorously assess the spatial distribution and total quantity (here defined in terms of content, i.e., concentration, and total mass) of entrapped CH 4 within the forefield sediments (Zhu et al., 2018). Nevertheless, to better characterize this potential CH 4 source, it is important to assess its spatial distribution and total quantity, particularly in glacier-forefield sediments, as we expect the potential for CH 4 release from these sediments to far exceed that from large bedrock surfaces due to the much higher specific surface area of the former (André et al., 2009;Michel and Courard, 2014). Moreover, as calcareous glacier-forefield sediments throughout the Swiss Alps are of similar origin (Weissert and Stössel, 2015), sediment-entrapped CH 4 may be a feature common to most if not all Swiss glacier forefields derived from calcareous bedrock. While this hypothesis remains to be tested, its confirmation would greatly increase the magnitude of this potential CH 4 source. Therefore, the overall goal of this study was to extend the work of Zhu et al. (2018) to other calcareous glacier forefields located in different regions of the Swiss Alps as well as to assess the distribution of entrapped CH 4 contents within a glacier forefield and compare the total mass of entrapped CH 4 among all sampled glacier forefields. Specific objectives included (1) testing the occurrence and origin of sediment-entrapped CH 4 in four additional calcareous glacier forefields; (2) assessing the spatial distribution of the sediment-entrapped CH 4 contents in detail within one glacier forefield, testing for dependencies on sediment depth, sediment age, and glacier-forefield landforms; based on results obtained from the previous step, (3) efficiently sampling sediments of the other glacier forefields to quantify the contents and total mass of sediment-entrapped CH 4 ; and, finally, (4) upscaling these results and deriving a first estimate of the total mass of the sediment-entrapped CH 4 contained in all Swiss glacier forefields situated on calcareous bedrock. the five glacier forefields from which sediment and bedrock samples were collected: Central Switzerland, including Im Griess (IMG), Griessfirn (GRF), and Griessen (GRI); Canton of Valais, including Wildstrubel (WIL) and Tsanfleuron (TSA). All forefields are located within the Helvetic nappes (green shaded area), which consist largely of Mesozoic limestones, shales, and marls (map modified from Weissert and Stössel, 2015).

Field sites and field work stages
Field work was conducted in five different glacier forefields: the Im Griess (IMG), Griessfirn (GRF), and Griessen (GRI) glaciers located in Central Switzerland in the Cantons of Uri (IMG and GRF) and Obwalden (GRI), and Tsanfleuron (TSA) and the previously investigated Wildstrubel (WIL; Zhu et al., 2018) glaciers located in the Canton of Valais (Figs. 1 and S1 in the Supplement). These forefields were selected for two main reasons. Foremost, their sediments are mainly derived from the calcareous bedrock of the Helvetic nappes (green shaded area in Fig. 1), which consist of a series of nappes (sheets of thrusted rocks) largely composed of Mesozoic limestones, shales, and marls of Jurassic to Eocene age (Pfiffner, 2014;Weissert and Stössel, 2015). They were originally deposited on the shallow northern shelf of the ancient Alpine Tethys Ocean (Weissert and Mohr, 1996) and were subsequently deformed, folded, and stacked on top of each other during Alpine orogeny (Herwegh and Pfiffner, 2005). Whereas individual nappes within the Helvetic nappe system therefore share a similar origin, the lithology and tectonic settings between individual nappes can be quite diverse (Weissert and Stössel, 2015). This was suggested to be a dominant factor determining rock CH 4 contents in the WIL catchment (Zhu et al., 2018). Consequently, we chose to investigate glacier forefields distant from one another (e.g., the distance from TSA to IMG is ∼ 136 km), for which sediments are derived from different nappes, but also glacier forefields in close proximity to each other (e.g., the distance from IMG to GRF is ∼ 3.8 km, and the distance from TSA to WIL is ∼ 24 km), for which sediments are derived, at least in part, from the same nappe. A second important reason for the selection was that all five glacier forefields are relatively easy to access, facilitating sample collection and transport to the laboratory.
We conducted our field work in two stages. During stage I in summer 2016, we performed a detailed investigation on the spatial distribution of sediment-entrapped CH 4 within a designated sampling zone at the GRF glacier forefield, using high spatial resolution sampling to determine variations in entrapped CH 4 contents in relation to sediment depth, sediment age, and glacier-forefield landforms. The GRF forefield was chosen for this purpose mainly because it features well-defined sediment-age classes and well-developed, clearly distinguishable landforms within a previously discretized and characterized sampling zone (Chiri et al., 2015(Chiri et al., , 2017. We also conducted measurements of sediment thickness (distance between the ground surface and the underlying bedrock) to estimate sediment volumes and, thus, the total mass of entrapped CH 4 present in these sediments. Results of the GRF field work were then used to adapt our sampling strategies for field work stage II, which was performed in summer 2017 to quantify the contents and total mass of sediment-entrapped CH 4 in the IMG, GRI, TSA, and WIL glacier forefields. During both field work stages, selected sediment and rock samples were used to identify the origin of the entrapped CH 4 based on CH 4 stable carbonisotope analyses and analyses of entrapped gas composition (see below).

Field work stage I (GRF glacier forefield)
Sampling and measurements during stage I in the GRF forefield was conducted in three steps: first, we tested the effect of sediment depth; second, we tested the effects of sediment age and glacier-forefield landforms on entrapped CH 4 contents; and, finally, we estimated the total mass of the sediment-entrapped CH 4 based on the measured CH 4 content, sediment thickness, and sediment-covered area.
To study the effect of sediment depth on entrapped CH 4 contents, we implemented a randomized design, selecting 14 random locations within our sampling zone (not shown). We collected a total of 52 sediment samples (each ∼ 500 g) by excavation from depths ranging from 20 to 70 cm b.s. (centimeters below the ground surface). All sediment samples were stored in clean plastic bags, transferred to the laboratory, and kept in the dark at 4 • C before further treatment. Following the extraction of entrapped gas and the subsequent quantification of CH 4 contents in sediment samples (see below), the effect of sediment depth on entrapped CH 4 contents was studied using a one-way ANOVA. The result of this analysis was then used to adapt the sampling scheme for the following step.
To study the effect of sediment age and glacier-forefield landforms on entrapped CH 4 contents, we implemented a randomized block-sampling design. We first divided the GRF sampling zone into nine blocks (a combination of three sediment-age classes and three landforms; Fig. 2a), using the spatial discretization of Chiri et al. (2017). The three sediment-age classes were A (0-20 years), B (20-50 years), and C (50-100 years). The three forefield landforms at GRF were floodplain, terrace, and sandhills. A floodplain refers to the frequently flooded area in the immediate vicinity of the glacial stream, which commonly consists of sediments of fine particle size (mostly silt) and a lack of vegetation. A terrace refers to an elevated, previously flooded area, i.e., a former floodplain, usually featuring some vegetation coverage. Finally, sandhills consist of unoriented, hummocky glacial-debris deposits, typically featuring poorly sorted, well-aerated sediments of sandy loam to sandy clayloam texture. We collected a total of 78 sediment samples (each ∼ 500 g) by excavation from a depth of 20 cm b.s., with 8-12 samples collected at random locations from within each block (Fig. 2a). The sampling depth of 20 cm b.s. was chosen based on our results from the previous step. Following laboratory analyses (see below), the impact of sediment age and landforms on entrapped CH 4 contents was studied using a two-way ANOVA.
In addition to sediments, we also collected a total of 17 bedrock samples from outcrops and large boulders within the GRF glacier forefield. These samples were used to determine the CH 4 content of the parent material (Zhu et al., 2018). All bedrock samples were stored in plastic bags, transferred to the laboratory, and stored in the dark at 4 • C before further treatment.
Estimation of the total mass of CH 4 entrapped in glacierforefield sediments also requires information on sediment thickness. For the GRF sampling zone we employed the elec-trical resistivity tomography (ERT) method (e.g., Kneisel, 2006;Reynolds, 1997;Scapozza et al., 2011). Five twodimensional, vertical ERT profiles (ERT1-ERT5) were measured during two field campaigns, covering the three sediment-age classes and the three landforms (Fig. 2b). Two profiles were measured parallel to the glacier stream (ERT2 and ERT5), and three profiles were measured perpendicular to the glacier stream (ERT1, 3, and 4). For each profile, 48 stainless-steel electrodes (30 cm long, 1.2 cm diameter) were hammered into the sediment to a depth of ∼ 15 cm and connected to two 24-core copper cables, which were linked to the ERT instrument (Syscal Pro; IRIS Instruments, Orléans, France) at the profile's midpoint. To improve electrical coupling of the electrodes with the skeleton-rich glacierforefield sediments, water-soaked sponges were positioned at the sediment surface surrounding each electrode. Profile ERT1 was measured with an electrode interspacing of 2.5 m (total profile length 120 m), whereas the other four profiles were measured with a 5 m distance between the electrodes (240 m profile length). Using a so-called Wenner-Schlumberger configuration (Loke, 2001), an electrical current was sent to the subsurface using a pair of electrodes. The voltage difference measured across the other pairs of electrodes was used to calculate the electrical resistivity of the subsurface. To infer the location of the sediment-bedrock interface, inversion of apparent resistivities was performed using the two-dimensional program RES2DINV (Loke and Barker, 1996). The average sediment thickness and its uncertainty within the GRF forefield was then analyzed in R. Electrical resistivities > 2000 m were considered indicative of solid bedrock, whereas resistivities < 2000 m were considered indicative of unconsolidated sediment (Kneisel, 2006;Reynolds, 1997). The portions of the ERT profiles for which the sediment-bedrock interface could not be detected were omitted from further analyses.

Field work stage II (IMG, GRI, WIL, and TSA glacier forefields)
During stage II, we collected a total of 111 sediment samples at a depth of 20 cm b.s. from the glacier forefields: 25 samples from IMG, 25 samples from GRI, 33 samples from WIL, and 28 samples from TSA (sampling locations shown in Fig. 3). Based on results obtained during field work stage I and given that glacier-forefield landforms were much less prominent at IMG, GRI, WIL, and TSA, we divided each of the four forefields into six blocks and collected four to eight sediment samples (each ∼ 500 g) from each block at random locations. We also collected 55 bedrock samples from outcrops and boulders at the glacier forefields: 13 samples from IMG, 14 samples from GRI, 12 samples from WIL, and 16 samples from TSA (sampling locations shown in Fig. 3).

Extraction of entrapped gas
We extracted entrapped gas from sediments and rocks using the acidification method described in Nauer et al. (2014) and Zhu et al. (2018). Before acid treatment, sediments were sieved with a clean 20 mm mesh sieve. Particles > 20 mm were excluded from subsequent analyses. For each sample, ∼ 3-5 g of sediment was weighed and transferred into a 117 mL serum bottle, sealed with a butyl rubber stopper, and crimped with an aluminum cap. The vial's headspace was then flushed with N 2 gas. Thereafter, 5 mL of deionized water was added to the vial followed by ∼ 50 mL of 6 N HCl to dissolve carbonate minerals. The headspace of each vial was connected to one or multiple 1 L gas bags (Tesseraux Spezialverpackungen GmbH, Bürstadt, Germany). Sediment samples released large amounts of gas immediately after the acid was added. When bubbling stopped, an additional 2 mL of 6 N HCl was added to each vial to confirm that the carbonate minerals were fully dissolved. Full dissolution of all carbonate minerals took ∼ 4 h. After gas extraction, ∼ 200 mL of gas was removed from the gas bags with syringes and stored in glass vials for further analysis. The total volume of gas remaining in gas bags was measured with a mass-flow meter (Bronkhorst, Reinach, Switzerland). Rocks were first hammered or sawed into ∼ 1 cm diameter fragments and then dissolved in the same way as sediments. Initial tests, in which we compared hammering with sawing to obtain rock fragments from both the surface and the core of larger rocks, showed insignificant effects on the fragments' entrapped CH 4 contents and other geochemical parameters.
As the duration of the respective mechanical treatment varied greatly between the collected fragments, we consider this to be evidence that neither hammering nor sawing had an adverse effect on measured geochemical parameters.

Quantification of methane, ethane, and propane
Concentrations of CH 4 were measured with a gas chromatograph equipped with a flame ionization detector (GC-FID; Trace GC Ultra, Thermo Fisher Scientific, Rodano, Italy) and a Porapak N100-120 column. The column-oven temperature was 30 • C, and the runtime was 36 s. Nitrogen carriergas flow was set to 26 mL min −1 . The FID was operated at 150 • C in high-sensitivity mode. Concentrations of ethane (C 2 H 6 ) and propane (C 3 H 8 ) were quantified in selected gas samples using the same GC-FID system but with an oven temperature of 40 • C for 2 min; this temperature was then increased to 140 • C at a rate of 25 • C min −1 before being held at a constant oven temperature of 140 • C for another 9 min. Gas contents were calculated as the mass of CH 4 , C 2 H 6 , and C 3 H 8 released during acidification, normalized to the dry weight of the sample. The dry weight of sediments was determined by the oven-drying of subsamples at 60 • C for 72 h. The computed entrapped gas contents C CH4 , C C2H6 , and C C3H8 were subsequently used to calculate the gas-wetness ratio as follows (Jackson et al., 2013): The gas-wetness ratio is a commonly used indicator of CH 4 origin: a value of > 1000 is considered evidence of microbial CH 4 , whereas a value of 1000 is considered indicative of thermogenic CH 4 (Rowe and Muehlenbachs, 1999).

Stable carbon-isotope analysis of entrapped methane
A total of 31 sediment and bedrock samples from the five glacier forefields were selected for stable carbon-isotope analysis of entrapped CH 4 (δ 13 C CH4 ). To determine δ 13 C CH4 we used a modified acidification protocol for gas extraction, which consisted of flushing the vials' headspace with He instead of N 2 to remove ambient air. Gas released during the acidification treatment was passed through two 1 M NaOH solutions to remove the majority of CO 2 , an Ascarite trap to remove final traces of CO 2 , a Drierite trap to remove H 2 O vapor, and a 1 M ZnCl 2 trap to remove potential H 2 S (all chemicals from Sigma Aldrich, Buchs, Switzerland). The purified gas samples were subsequently analyzed by GC-IRMS (gas chromatography isotope ratio mass spectrometry; Isoprime, Elementar UK Ltd., Stockport, UK).
2.3 Estimation of total mass of CH 4 entrapped in glacier-forefield sediments 2.3.1 Estimation of entrapped CH 4 mass for the GRF sampling zone The mass of CH 4 (m CH4 ) entrapped in a specific volume of porous sediment may be calculated using the following equation: where C CH4 is the sediment-entrapped CH 4 content (mass of CH 4 per mass of sediment); ρ sed is the sediment-particle density; A sed and T sed are the sediment-covered area and sediment thickness in the glacier forefield, respectively; and θ t,sed is the total interparticle sediment porosity, hereafter referred to as sediment porosity. To determine m CH4 for the GRF sampling zone, we applied Eq. (1) separately to each landform but also used averaged values for entrapped CH 4 contents (from laboratory analyses), sediment thickness (from ERT field measurements), and sediment-covered area estimated from aerial maps (https://map.geo.admin.ch, last access: 12 December 2019). In Eq. (1), the term in brackets represents the sediment's solid volume. To compute the sediment's solid volume, we assumed a mean θ t,sed = 0.42 ± 0.02, as determined for this site by Nauer et al. (2012).
To convert solid volume to sediment mass, a mean value of ρ sed = 2.71 ± 0.15 g cm −3 was used, as derived by Daly (1935) from measurements of a variety of calcite rock samples. The total uncertainty in the estimated mean CH 4 mass m CH4 , expressed as the standard error (SE) of the mean (σ m CH4 ), was computed using the following equation: where σ x represents the SE associated with any parameter's mean value x. The individual contribution of any parameter x (frac x , in %) to the total uncertainty in m CH4 was then computed using We note that SE values (reported as x ± σ x ) are used as a measure of uncertainty of any parameter's mean value x throughout this paper, whereas standard deviations (SD, reported as x (± σ x )) are used as a measure of general parameter variability. 4 mass for the five  glacier forefields (IMG, GRF, GRI, WIL, and  TSA) To compute the total mass and associated uncertainty of sediment-entrapped CH 4 for all five glacier forefields, we employed Eqs. (1) and (2), although with partially modified parameters. For C CH4 we used mean values of sedimententrapped CH 4 contents determined for each glacier forefield. In addition, we determined mean values A sed from estimates of the maximum and minimum extents of the sediment-covered area within each glacier forefield. We used the areas exposed as a result of glacier retreat since the last glacial maximum (Little Ice Age, ∼ 1850) as a maximum. These areas were estimated from the difference in glacial extent as taken from the most current (2018) and historic (∼ 1850) topographic maps (Swisstopo; https://map. geo.admin.ch, last access: 12 December 2019; Fig. S2). Minimum areas were directly estimated from the 2018 aerial maps. Moreover, data on sediment thickness were unavailable for the IMG, GRI, WIL, and TSA glacier forefields as well as for the GRF forefield outside of the designated sampling zone. Therefore, we used the average value of T sed = 10.0 ± 3.0 m obtained from our ERT measurements in the GRF sampling zone (see below) as an average T sed for all five glacier forefields. We note that our average T sed value agrees well with previous measurements performed in another Swiss glacier forefield, in which T sed ∼ 8 m was obtained by borehole drilling (Kneisel and Kääb, 2007). Finally, we used values of θ t,sed for the GRF, GRI, and WIL forefields as determined for these sites by Nauer et al. (2012). As estimates of sediment porosity were unavailable for IMG and TSA, we used θ t,sed = 0.44 ± 0.05 for these forefields, averaged from data reported for the five calcareous glacier forefields (Nauer et al., 2012).

Estimation of entrapped CH 4 mass for sediments in all Swiss glacier forefields derived from calcareous bedrock
We again used Eqs. (1) and (2) to upscale results and to compute a first estimate of the total mass of sediment-entrapped CH 4 contained in all Swiss glacier forefields derived from calcareous bedrock. In this case, we used the mean C CH4 of the five glacier forefields. Calcareous glacier-forefield surface area in Switzerland (A sed in Eq. 1) was estimated from available data on the decrease in glaciated area in the Swiss Alps between the Little Ice Age (∼ 1850; Zemp et al., 2008) and the year 2010 (Fischer et al., 2014) as well as from an estimate of the fraction of calcareous bedrock area to the total area of the Swiss Alps taken from the Tectonic Map of Switzerland 1 : 500 000 (Federal Office of Topography, swisstopo). Mean values for ρ sed , T sed , and θ t,sed were used as described above.

Geochemistry of gas entrapped in sediment and bedrock samples
From the 271 sediment samples from the five glacier forefields, we analyzed 256 samples for entrapped CH 4 contents. All analyzed sediments contained detectable amounts of CH 4 ranging from 0.08 to 73.81 µg CH 4 g −1 dry weight (d.w.; Fig. 3), with an average of 14.9 (± 17.0) µg CH 4 g −1 d.w. Gas released from 225 samples was analyzed for C 2 H 6 and C 3 H 8 contents, of which 215 contained detectable amounts of C 2 H 6 ranging from 0.002 to 1.67 µg C 2 H 6 g −1 d.w., with an average of 0.25 (± 0.32) µg C 2 H 6 g −1 d.w. In addition, 146 out of 225 samples contained detectable amounts of C 3 H 8 ranging from 0.001 to 0.82 µg C 3 H 8 g −1 d.w., with an average of 0.11 (± 0.15) µg C 3 H 8 g −1 d.w. (not shown). The average gas-wetness ratio for all sediment samples was 75.2 (± 48.4), and the average δ 13 C CH4 was −28.23 (± 3.42) ‰. Plotting δ 13 C CH4 values vs. gas-wetness ratios in a so-called Bernard diagram ( Fig. 4; Bernard et al., 1978) indicated a thermogenic origin for sediment-entrapped CH 4 , derived from ancient terrestrial or marine organic matter (kerogen types III and II; Fig. 4). Although CH 4 extracted from sediments collected in the IMG glacier forefield showed a higher variability in the gas-wetness ratios than CH 4 extracted from the sediments of other glacier forefields, it still fell into the same origin type in the Bernard diagram.
All 72 bedrock samples were analyzed for CH 4 content, and 64 contained detectable amounts of CH 4 ranging from 0.06 to 108.58 µg CH 4 g −1 , with an average of 11.4 (± 20.0) µg CH 4 g −1 (Fig. 3). The average δ 13 C CH4 value of −29.21 (± 2.77) ‰ was similar to that of sedimententrapped CH 4 . Likewise, the average gas-wetness ratio of gas extracted from rocks was 78.45 (± 121.84), which is similar in value but with higher variability than gas-wetness ratios for sediment-entrapped CH 4 (Fig. 4). Together, these data suggest a common, thermogenic origin of entrapped CH 4 in sediments and rocks, with little apparent alteration from physical and/or chemical weathering. Moreover, our data suggest that entrapped CH 4 is of similar origin in all five glacier forefields.

Spatial distribution of sediment-entrapped CH 4 contents in the GRF sampling zone
Methane contents in 52 samples collected at depths from 20 to 70 cm b.s. ranged from 1.19 to 11.24 µg CH 4 g −1 d.w., with one exceptionally high value at 40 cm b.s. (Fig. 5). Based on these data, there was no clear correlation between sediment depth and entrapped CH 4 contents (one-way ANOVA, Figure 4. Adapted Bernard diagram (Bernard et al., 1978) showing the gas-wetness ratio (C CH4 /(C C2H6 + C C3H8 )) vs. δ 13 C CH4 for gas released from selected sediment and rock samples collected from the Im Griess (IMG), Griessfirn (GRF), Griessen (GRI), Wildstrubel (WIL), and Tsanfleuron (TSA) glacier forefields. Gray shaded areas indicate different CH 4 origins (microbial vs. thermogenic). p = 0.9). Thus, we subsequently proceeded to collect sediments from 20 cm b.s. only, and we assumed that these samples were representative in terms of entrapped CH 4 content for the entire sediment thickness. The effects of the sediment age and landform proxies on the entrapped CH 4 contents were tested using sediments collected from 20 cm b.s. at 99 locations (Fig. 2a). The CH 4 contents in these samples ranged from 0.59 to 34.82 µg CH 4 g −1 d.w. (Fig. 3b), with an average of 5.30 (± 4.86) µg CH 4 g −1 d.w. Two-way ANOVA analysis indicated that landform had a significant effect on sedimententrapped CH 4 contents (p = 0.03), whereas the effects of sediment age (p = 0.19) and the combined effects of sediment age and landform on entrapped CH 4 contents (p = 0.37) were insignificant. Post hoc comparisons using the Tukey HSD test indicated that mean values for the sedimententrapped CH 4 content (Table 1) were significantly different between floodplain and sandhill landforms (p = 0.03) and were weakly different between floodplain and terrace landforms (p = 0.10). The difference between terrace and sandhill landforms with respect to mean sediment-entrapped CH 4 content was insignificant (p = 0.88).

Mass of sediment-entrapped CH 4 in the GRF sampling zone
To estimate the mass of sediment-entrapped CH 4 stored within the GRF sampling zone, we used Eq. (1) with mean values of entrapped CH 4 contents, sediment thickness, and sediment-covered area determined for each of the three landforms (Table 1). While C CH4 varied by a factor < 1.4 be- Floodplain 6.37 ± 0.55 11.8 ± 3.0 2.07 × 10 4 ± 2.0 × 10 2 3.84 × 10 5 ± 1.0 × 10 5 2.4 ± 0.7 Terrace 4.72 ± 0.97 12.5 ± 4.0 2.06 × 10 4 ± 2.0 × 10 2 4.04 × 10 5 ± 1.3 × 10 5 1.9 ± 0.7 Sandhill 5.04 ± 0.78 6.4 ± 3.2 1.05 × 10 5 ± 1.0 × 10 3 1.06 × 10 6 ± 5.4 × 10 5 5.4 ± 2.8 Combined 5.38 ± 0.49 10.2 ± 3.0 1.47 × 10 5 ± 1.4 × 10 3 2.36 × 10 6 ± 7.1 × 10 5 9.7 ± 3.0 b 12.7 ± 4.0 c a Standard error of the mean (SE). b Calculated by summing the estimated mass of entrapped CH 4 from each landform. c Calculated using average values for entrapped CH 4 contents, sediment thickness, and sediment-covered area. Figure 5. Sediment-entrapped CH 4 contents as a function of sediment depth for samples collected in three sediment-age classes in the Griessfirn (GRF) sampling zone. The box and whiskers plot (top panel) shows the range of the entrapped CH 4 contents displayed below, with the box representing the 25th, 50th (median), and 75th percentiles; the whiskers indicating the 5th and 95th percentiles; and the outliers marked as dots. tween landforms, sediment thickness was highly variable along the five measured ERT profiles (range 1.0-31.5 m; Figs. 6, S3), and T sed varied by a factor of ∼ 2 between landforms (Table 1). Sediment-covered area also showed substantial variation between the different landforms. Within the GRF sampling zone, the sandhill landform comprised the largest sediment-covered area with A sed ≈ 10 5 m 2 , which was about 5 times larger than the values for floodplain and terrace landforms. Consequently, the largest sediment mass was contained in the sandhill landform (a factor 2-3 larger Figure 6. Vertical, two-dimensional electrical resistivity tomography (ERT) cross sections of the (a) ERT1 and (b) ERT2 profiles collected in the sampling zone of Griessfirn (GRF) glacier forefield. Solid black lines indicate the approximate location of the interface between unconsolidated sediment and the bedrock underneath. Lines were omitted at locations where the sediment-rock interface was too deep to be detected. than floodplain and terrace; Table 1). All three landforms combined comprised a surface area of ∼ 1.5 × 10 5 m 2 and contained an estimated mass of ∼ 2.4 × 10 6 t sediment. Summing the masses of sediment-entrapped CH 4 for each landform yielded a total m CH4 of 9.7 ± 3.0 t CH 4 . When calculated using average values for entrapped CH 4 contents and sediment thickness, as well as combined sediment-covered area, the estimated m CH4 within the GRF sampling zone was 12.7 ± 4.0 t CH 4 (last row in Table 1). Uncertainties in individual m CH4 values up to ∼ 50% mostly arose from uncertainties in T sed (dominated by the large variability in sediment thickness across the GRF sampling zone) and to a smaller degree from uncertainties in C CH4 .

Contents and total mass of sediment-entrapped CH 4 in five glacier forefields
Methane contents varied substantially between different glacier forefields (Table 2), with the distance between the forefields playing an apparently important role. Specifically, the IMG, GRF, and GRI glacier forefields are located in the northeast of the Helvetic nappes relatively close to each other (Fig. 1) and featured similar, low sediment-entrapped CH 4 contents. Likewise, the WIL and TSA glacier forefields are located close to each other in the southwest of the Helvetic nappes and featured similar but high sediment-entrapped CH 4 contents. Indeed, our ANOVA results indicated that differences in the sediment-entrapped CH 4 contents were insignificant between the IMG, GRF, and GRI glacier forefields (p = 0.36) and between the WIL and TSA glacier forefields (p = 0.18). Conversely, differences in entrapped CH 4 contents between the two groups of glacier forefields were highly significant (p < 0.0001). The total mass of CH 4 entrapped in sediments of the five glacier forefields was calculated using estimated values for sediment thickness (10.0 ± 3.0 m; the thickness measured in the GRF sampling zone, see above) and sedimentparticle density (2.71 ± 0.15 g cm −3 ; Daly, 1935) that were assumed to be identical for all five forefields, as well as specific data for each glacier forefield on entrapped CH 4 contents, sediment-covered area, and sediment porosity (Table 2). Whereas θ t,sed values only varied slightly between the five forefields, A sed varied by up to a factor of ∼ 3 (IMG vs. WIL), and C CH4 varied by up to a factor of ∼ 7 (GRF vs. WIL). This led to substantial differences in the estimated total mass of sediment-entrapped CH 4 between the five forefields, which ranged from 200 ± 74 t CH 4 for the GRF glacier forefield to 3881 ± 1367 t CH 4 for the WIL forefield (Fig. 7a). Estimates of sediment-entrapped CH 4 mass for the WIL and TSA glacier forefields were significantly larger than for IMG, GRF, and GRI. For all five forefields, sediment thickness and sediment-covered area contributed most to uncertainties in the quantification (Fig. 7b). Conversely, entrapped CH 4 contents, sediment porosity, and sediment-particle density together contributed ≤ 16% to the calculated uncertainties.

Mass of sediment-entrapped CH 4 in all Swiss glacier forefields on calcareous bedrock
The first estimate of the total mass of sediment-entrapped CH 4 in all calcareous glacier forefields in Switzerland was based on published data on glacier retreat in the Swiss Alps, an estimation of the fraction of calcareous glacierforefield surface area, mean values for sediment thickness, sediment-particle density and sediment porosity, as well as a mean value for sediment-entrapped CH 4 content obtained from the five investigated glacier forefields (18.5 ± 4.4 µg CH 4 g −1 d.w.; Table 3). Between the end of the Little Ice Age (∼ 1850) and 2010, the glaciated area within the Swiss Alps has decreased by ∼ 676 km 2 to less than 60% of its original value (data sources given in Table 3). When multiplied by the fraction of the calcareous bedrock area in the Swiss Alps (54.6 ± 1.7%), this yielded an exposed calcareous glacier-forefield area of ∼ 369 km 2 . The total sediment mass contained within this exposed calcareous glacier-forefield area was then computed as 5.62 × 10 9 ± 1.46 × 10 9 t. From these numbers, we estimated the total mass of sediment-entrapped CH 4 in all Swiss glacier forefields derived from calcareous bedrock to be 1.04 × 10 5 ± 3.7 × 10 4 t CH 4 .

The widespread occurrence of sediment-entrapped, thermogenic CH 4 in calcareous glacier forefields
We detected substantial quantities of sediment-entrapped CH 4 in all sampled glacier forefields. Entrapped CH 4 was ubiquitously encountered at different sediment depths as well as in different forefield landforms and sediment-age classes. We also detected entrapped CH 4 in most bedrock samples obtained from these glacial catchments. Furthermore, our data indicated that both sediment-and rock-entrapped CH 4 are of thermogenic origin. Thus, the results presented here extend our previous studies (Nauer et al., 2012;Zhu et al., 2018) by providing a more detailed survey on entrapped CH 4 contained in glacier-forefield sediments across the Helvetic nappes, and they support our hypothesis on its widespread occurrence and thermogenic origin in calcareous, Swiss Alpine glacier forefields. On the other hand, we cannot entirely reject the possibility of the presence of microbial CH 4 sources in certain parts of glacier forefields, particularly in water-logged sediments. Methanogenic potential in isolated hot spots of water-logged sediments has previously been confirmed for the WIL glacier forefield, but it was considered to be of minor importance under field conditions (Zhu et al., 2018). In the present study, no attempt was made to specifically identify potential methanogenic hot spots in the sediments of the other four glacier forefields. Methane is commonly found in organic-rich sedimentary rocks such as shales, marls, and limestones as a product of the thermal maturation of buried organic matter (Etiope, 2017;Horsfield and Rullkötter, 1994). Previous studies on fluid inclusions in quartz and calcite minerals collected from Alpine fissures and veins within the Helvetic nappes have revealed the existence of four fluid zones, including a large thermogenic CH 4 zone (Gautschi et al., 1990;Mazurek et al., 1998;Mullis et al., 1994;Tarantola et al., 2007). The five glacier forefields that we sampled in this study were all located within or near the border of this thermogenic CH 4 zone (see Fig. 1 in Tarantola et al., 2007). Thus, our results agree with previous findings on the occurrence of thermo- Table 2. Mean values and uncertainties of sediment-entrapped CH 4 content, sediment-covered area, and total interparticle sediment porosity for the Im Griess (IMG), Griessfirn (GRF), Griessen (GRI), Wildstrubel (WIL), and Tsanfleuron (TSA) glacier forefields located within the Helvetic nappes of Switzerland. Also listed are individual nappes and major geological formations, from which glacier-forefield sediments are derived.

Glacier
Entrapped CH 4 Sediment-covered Sediment forefield content area (km 2 ) porosity a (µg CH 4 g −1 d.w.)  Contribution of sediment-entrapped CH 4 contents, sediment thickness and area, sedimentparticle density, and sediment porosity to total uncertainties in the estimation of the mass of sediment-entrapped CH 4 (error bars given in panel a) for each of the five glacier forefields.
genic CH 4 in this region, including the occurrence of thermogenic CH 4 detected in gas seeps near Giswil, Central Switzerland, which lies on Penninic Flysch underlain by the Helvetic nappes (Etiope et al., 2010). On the other hand, our results also show that CH 4 entrapment within the Helvetic nappes is not restricted to fluid inclusions in fissure minerals and that substantial quantities of CH 4 are entrapped within the matrix of the sedimentary bedrock and sediment particles themselves, presumably within inter-and intragranular micro-and macropores (Hashim and Kaczmarek, 2019;Moshier, 1989;Léonide et al., 2014;Abrams, 2017). Our geochemical data further indicate a common origin for CH 4 entrapped in bedrock and glacier-forefield sediments, derived from ancient terrestrial and marine organic matter (kerogen types III and II, respectively; Fig. 4). This provides further evidence that CH 4 entrapped in the forefield sediments of the Helvetic nappes has its origin in the calcareous parent bedrock. Moreover, terrestrial and marine organic matter as the ultimate source of sediment-and rockentrapped CH 4 agrees with the origin of the Helvetic nappes: their sediments and organic matter were originally deposited under highly variable climatic conditions on the shallow northern shelf of the ancient Alpine Tethys Ocean (Weissert and Mohr, 1996;Weissert and Stössel, 2015).  (Chesworth et al., 2008;Fu and Harbor, 2011), the observed variability in the sediment-entrapped CH 4 contents reflects the variability in the entrapped CH 4 contents of the various geological formations and the associated mineralogy present in each catchment (Fig. 3f). Entrapped CH 4 contents in sedimentary bedrock is typically affected by three main factors: the quantity and quality of organic matter buried during sediment deposition, the thermal history during sediment diagenesis and subsequent organic matter catagenesis, and the resulting permeability of the calcareous bedrock, which affects potential gas migration (e.g., Horsfield and Rullkötter, 1994;Mani et al., 2017). Whereas geological formations contained within the same nappe are expected to possess a similar thermal history, the quantity and quality of organic matter buried may vary substantially between individual formations depending on prevailing conditions during the period of sediment deposition (Weissert and Mohr, 1996;Weissert et al., 1985). Thus, variability in rock-and sediment-entrapped CH 4 contents is to be expected for glacial catchments featuring geological formations from different time periods, as was observed for all of the glacier forefields sampled in this study (Table 2). Our study in the GRF forefield sampling zone indicated that sediment-entrapped CH 4 content varied little with sediment depth (Fig. 5) and sediment age. However, we cannot exclude the possibility that such variations could be somewhat larger outside of the sampled depth interval, e.g., in top-layer sediments at depths < 5 cm as a result of enhanced chemical, physical, or biological weathering (Bernasconi et al., 2011;van der Meij et al., 2016;Lazzaro et al., 2009). We refrained from collecting top-layer sediments because they were generally much coarser in all five glacier forefields and, thus, did not appear representative of bulk sediments present at greater depth. We assume that sediment fines are continuously removed from the top layer as a result of physical (wind and water) erosion.
On the other hand, we consider the lack of significant variation with sediment age as an indication that CH 4 in glacier-forefield sediments is relatively stable in its entrapped state. This hypothesis is supported by the results of our geochemical analyses for all five glacier forefields, which mostly indicated high similarity between sediment-and rockentrapped CH 4 in terms of the range of measured CH 4 contents (Fig. 3f), as well as gas-wetness ratios and δ 13 C CH4 values (Fig. 4). Thus, although sediments have undergone erosion from the parent bedrock and subsequent weathering, changes in entrapped CH 4 geochemical characteristics appeared negligible. This indicates that a potential release of entrapped CH 4 from sediment particles by molecular diffusion or from the oxidation of CH 4 in its entrapped state within sediment particles may be small in magnitude, as these processes would be expected to cause a noticeable change in CH 4 geochemical characteristics (Schloemer and Krooss, 2004;Whiticar, 1999;Zhang and Krooss, 2001). Therefore, our findings suggest that CH 4 entrapped in bedrock and sediment matrices resides largely in inaccessible, occluded rather than connected pore spaces. However, a potential release of entrapped CH 4 from occluded pore spaces may yet occur via sediment erosion processes, in particular by means of physical and/or chemical weathering of calcareous minerals (Emmanuel and Levenson, 2014;Ryb et al., 2014;Trudgill and Viles, 1998). As these processes act on rock surfaces, they are of great importance to sediments with large specific surface areas, with the specific surface area being inversely related to particle size (Michel and Courard, 2014). Although we are aware that similar erosion processes will act upon large bedrock surfaces, e.g., rock walls and other outcrops within glacial catchments, we have refrained from considering CH 4 release from these locations (thus far) due to the much smaller specific surface areas involved. Unfortunately, the release of entrapped CH 4 as a result of sediment erosion may not be detectable in the sedi-ment's entrapped CH 4 contents, as both CH 4 and sediment mass is lost as a result of erosion. Hence, our present data set yields no information on the relevance of erosion processes for CH 4 release.
In contrast to sediment depth and sediment age, we detected a small but significant difference in the mean sediment-entrapped CH 4 content between landforms within the GRF sampling zone. Specifically, mean entrapped CH 4 content in floodplain sediments was significantly higher than in terrace and sandhill sediments (Table 1). We can only speculate about the possible reasons for this observation. One reason could be that floodplain sediments, which are intermittently removed and deposited by the glacial stream during and after flooding events, originate from locations outside of our sampling zone, i.e., from different parent bedrock (Fig. S2). Therefore, it is possible that we missed sampling parent bedrock types (e.g., from steep rock walls) with different (in this case higher) entrapped CH 4 contents in this or any of the other glacial catchments. An improved mineralogical investigation of sediments in the various landforms would aid in clarifying this issue.
Finally, our data revealed large regional differences in mean sediment-entrapped CH 4 contents between glacier forefields (Table 2). This may be explained by the fact that sediments in glacier forefields located in close proximity to one another are, at least in part, derived from the same individual nappes and geological formations contained therein. For example, both the WIL and TSA glacier forefields harbor sediments derived from the Wildhorn Nappe, featuring several identical geological formations. Hence, this result supports our previous hypothesis that differences in lithology, mineralogy, and tectonic settings between individual nappes play an important role in determining bedrock-and, thus, sediment-entrapped CH 4 contents (Zhu et al., 2018). Regional differences in the entrapped CH 4 contents paired with differences in the sediment-covered area led to significant variation in the estimates for the total mass of CH 4 stored in sediments of the five glacier forefields (Fig. 7a). Uncertainties associated with these estimates were reasonably small and arose largely from uncertainties in sediment thickness and sediment-covered area (Fig. 7b). To further reduce these uncertainties, measurements of these parameters across entire glacier forefields would be of help using techniques such as geophysical methods for sediment thickness (such as the ERT method used in the GRF sampling zone) and field mapping of sediment-covered area in combination with GIS-based methods utilizing digital elevation models (e.g., Geilhausen et al., 2012;Smith and Clark, 2005;Zemp et al., 2005). Unfortunately, field measurements in the rugged alpine environment are typically time-consuming, expensive, and challenging to perform.

A substantial quantity of sediment-entrapped CH 4 with a currently unknown fate
Our first, rough estimate of the total quantity of CH 4 entrapped in sediments of all calcareous Swiss glacier forefields combined yielded a substantial mass of 1.04 × 10 5 ± 3.7 × 10 4 t CH 4 , contained within a solid volume of ∼ 2.1 km 3 glacier-forefield sediments. At first glance, this number appears large when compared with an estimate of the annual CH 4 emissions to the atmosphere (5.7 × 10 3 t CH 4 ) from natural and seminatural sources in Switzerland, including emissions from lakes, reservoirs, wetlands, and wild animals (Hiller et al., 2014). However, while the latter data represent annual CH 4 fluxes, the fate of sediment-entrapped CH 4 still remains elusive (see below). On the other hand, our number is in good agreement with a previous estimate of the CH 4 content for Valanginian marl, a geological formation within the Helvetic nappes, containing calcite fracture fill (∼ 0.7 × 10 5 -2.1 × 10 5 t CH 4 km −3 bedrock; Gautschi et al., 1990).
Our estimate for the total sediment-entrapped CH 4 mass is subject to substantial uncertainty. The two largest contributors to the calculated uncertainty are the sediment-entrapped CH 4 content and sediment thickness. In addition, there is considerable uncertainty in the exposed calcareous glacierforefield area, as the latter was only roughly estimated based on glacier retreat and the fraction of calcareous bedrock area in the Swiss Alps. As discussed above for individual glacier forefields, field measurements and GIS-based methods may help to reduce uncertainties related to sediment thickness and exposed area. An important way to reduce uncertainty related to the entrapped CH 4 contents would be to generate a database of CH 4 contents for different geological formations present within the Helvetic nappes, as lithology and tectonic settings appear to control CH 4 contents. Determination of the areal extent of different geological formations would likely help to reduce uncertainties in the sediment-entrapped CH 4 mass.
Whether or not sediment-entrapped CH 4 plays a role as an emission source to the atmosphere will largely depend upon its rate of release from sediment particles and its potential consumption by MOB in aerated sediments. While we produced some evidence that CH 4 is stable in its entrapped state (see discussion above), further investigations will be required to specifically elucidate possible mechanisms and fluxes of CH 4 release in forefield sediments, in particular during periods of enhanced physical and/or chemical weathering, e.g., during rainstorms or snow melt (Winnick et al., 2017). Notably, atmospheric CH 4 oxidation has previously been detected in several glacier forefields including our GRF site (Bárcena et al., 2011;Chiri et al., 2015;Hofmann et al., 2013). These studies indicated that MOB activity in forefield sediments establishes quickly (within the first 10 years after glacier retreat) and that fluxes of CH 4 uptake from the atmosphere increase to values comparable with mature soils 3626 B. Zhu et al.: Methane in glacier-forefield sediments within a few decades (Chiri et al., 2015). Nonetheless, intermittent CH 4 emissions to the atmosphere were also observed in GRF floodplain sediments (Chiri et al., 2017). Hence, we hypothesize that traces of CH 4 released from sediment particles may be consumed by MOB, at least under favorable environmental conditions, and serve as an additional source of energy and carbon to this group of microorganisms. This hypothesis, of course, awaits experimental confirmation.

Summary and conclusions
Our results provide new evidence for the widespread occurrence of sediment-entrapped, thermogenic CH 4 in Swiss calcareous glacier forefields. As entrapped CH 4 with highly similar geochemical characteristics was also detected in most bedrock samples collected from nearby geological formations, we conclude that CH 4 entrapped in forefield sediments of the Helvetic nappes has its origin in the calcareous parent bedrock. Hence, the spatial variability in the sedimententrapped CH 4 contents within glacier forefields largely reflects the variability in the entrapped CH 4 contents of the surrounding bedrock types.
Within glacier forefields, sediment-entrapped CH 4 contents and other geochemical characteristics showed little systematic variation with sediment age and, thus, time of exposure to the atmosphere following glacier retreat. In combination with the noted similarity in geochemical characteristics, we took this finding to be evidence of the fact that CH 4 in glacier-forefield sediments is relatively stable in its entrapped state, presumably because it resides in occluded pore spaces within bedrock and sediment matrices. This further indicates that CH 4 entrapment within the Helvetic nappes is not restricted to fluid inclusions in fissure minerals and that substantial quantities of CH 4 are entrapped within the matrix of the sedimentary bedrock and sediment particles themselves. Moreover, our results revealed large regional differences in the mean sediment-entrapped CH 4 contents between glacier forefields, supporting our previous hypothesis that differences in lithology and tectonic settings between individual nappes play an important role in determining bedrockand, thus, sediment-entrapped CH 4 contents.
Data availability. The data used in this paper will be made available on the ETH Zurich Research Collection repository after the article is published.
Author contributions. BZ, MR, MK, and MHS helped with sample collection and/or geochemical measurements and also substantially contributed to the interpretation of data. DB helped with the electrical resistivity tomography measurements and the subsequent electrical resistivity tomography data analyses. BZ and MHS wrote the paper. MHS designed the study and acquired the funding for the project. All authors commented on the paper and approved the final version of the article.