In Situ Monitoring of Free-Phase Gas Accumulation and Release in Peatlands Using Ground Penetrating Radar (GPR)

We tested a set of surface common mid-point (CMP) ground penetrating radar (GPR) surveys combined with elevation rods ( to monitor surface deformation) and gas flux measurements to investigate in-situ biogenic gas dynamics and ebullition events in a northern peatland ( raised bog). The main findings are: ( 1) changes in the two-way travel time from the surface to prominent reflectors allow estimation of average gas contents and evolution of free-phase gas (FPG); ( 2) peat surface deformation and gas flux measurements are strongly consistent with GPR estimated changes in FPG content over time; ( 3) rapid decreases in atmospheric pressure are associated with increased gas flux; and ( 4) single ebullition events can induce releases of methane much larger ( up to 192 g/m(2)) than fluxes reported by others. These results indicate that GPR is a useful tool for assessing the spatial distribution, temporal variation, and volume of biogenic gas deposits in peatlands.

Stratigraphic Continuity in 400 Mhz Short-Pulse Radar Profiles of Firn in West Antarctica

We track dated firn horizons within 400 MHz short-pulse radar profiles to find the continuous extent over which they can be used as historical benchmarks to study past accumulation rates in West Antarctica. The 30-40 cm pulse resolution compares with the accumulation rates of most areas. We tracked a particular set that varied from 30 to 90 m in depth over a distance of 600 km. The main limitations to continuity are fading at depth, pinching associated with accumulation rate differences within hills and valleys, and artificial fading caused by stacking along dips. The latter two may be overcome through multi-kilometer distances by matching the relative amplitude and spacing of several close horizons, along with their pulse forms and phases. Modeling of reflections from thin layers suggests that the - 37 to - 50 dB range of reflectivity and the pulse waveforms we observed are caused by the numerous thin ice layers observed in core stratigraphy. Constructive interference between reflections from these close, high-density layers can explain the maintenance of reflective strength throughout the depth of the firn despite the effects of compaction. The continuity suggests that these layers formed throughout West Antarctica and possibly into East Antarctica as well.

Spatial Variability in Biogenic Gas Accumulations in Peat Soils Is Revealed By Ground Penetrating Radar (GPR)

We performed surface and borehole ground penetrating radar (GPR) tests, together with moisture probe measurements and direct gas sampling to detect areas of biogenic gas accumulation in a northern peatland. The main findings are: (1) shadow zones (signal scattering) observed in surface GPR correlate with areas of elevated CH4 and CO2 concentration; (2) high velocities in zero offset profiles and lower water content inferred from moisture probes correlate with surface GPR shadow zones; (3) zero offset profiles depict depth variable gas accumulation from 0-10% by volume; (4) strong reflectors may represent confining layers restricting upward gas migration. Our results have implications for defining the spatial distribution, volume and movement of biogenic gas in peatlands at multiple scales.

Phase Structure of Radar Stratigraphic Horizons Within Antarctic Firn

We have recorded reflection profiles of firn through large areas of West Antarctica and part of the East Antarctic plateau using 400 MHz short-pulse radar. The locations show accumulation rates that vary from well above to well below the vertical radar resolution. Most reflection horizons have extensive lateral continuity, and are composed of distinctive wavelets with a consistent phase polarity sequence within their successive half-cycles. We modeled these waveforms, and conclude that they arise from thin, double layers of ice over hoar, which is consistent with the standard model of firn stratification. In addition, we conclude that ice/hoar layers are extensive throughout West Antarctica and also present (although more sparsely) beneath the Antarctic Plateau.

Stratigraphic Variation Within Polar Firn Caused by Differential Accumulation and Ice Flow: Interpretation of a 400 Mhz Short-Pulse Radar Profile from West Antarctica

We investigate causes of the stratigraphic variation revealed in a 177 km, 400 MHz short-pulse radar profile of firn from West Antarctica. The profile covers 56 m depth, and its direction was close to those of the ice flow and mean wind. The average, near-surface accumulation rates calculated from the time delays of one radar horizon consistently show minima on leeward slopes and maxima on windward slopes, confirming an earlier study based on stake observations. The stratigraphic variation includes up to 30 m depth variation in individual horizons over tens of km, fold limbs that become progressively steeper with depth, and fold-hinge loci that change direction or propagate down-ice with depth over distances far less than predicted by the ice speeds. We use an accumulation rate model to show how local rate anomalies and the effect of ice speed upon a periodic variation in accumulation rate cause these phenomena, and we reproduce two key features seen in the stratigraphic variations. We conclude that the model provides an explanation of changes in spatial stratigraphy and local measures of accumulation history given the constraints of surface topography, ice and wind velocities, and a general accumulation rate for an area.

Thermal Convection in Ice Sheets: We Look But Do Not See

Thermal convection in the Antarctic and Greenland ice sheets has been dismissed on the grounds that radio-echo stratigraphy is undisturbed for long distances. However, the undisturbed stratigraphy lies, for the most part, above the density inversion in polar ice sheets and therefore does not disprove convection. An echo-free zone is widespread below the density inversion, yet nobody has cited this as a strong indication that convection is indeed present at d�pth. A generalized Rayleigh criterion for thermal convection in e1astic-viscoplastic polycrystalline solids heated from below is developed and applied to ice-sheet convection. An infinite Rayleigh number at the onset of primary creep decreases with time and becomes constant when secondary creep dominates, suggesting that any thermal buoyancy stress can initiate convection but convection cannot be sustained below a buoyancy stress of about 3 kPa. An analysis of the temperature profile down the Byrd Station core hole suggests that about 1000 m of ice below the density inversion will sustain convection. Creep along the Byrd Station strain network, radar sounding in East Antarctica, and seismic sounding in West Antarctica are examined for evidence of convective creep superimposed on advective creep. It is concluded that the evidence for convection is there, if we look for it with the intention offinding it.

Variability in Accumulation Rates from GPR Profiling on the West Antarctic Plateau

lsochronal layers in firn detected with ground-penetrating radar (GPR) and dated using results from ice-core analyses are used to calculate accumulation rates along a 100 km across-flow profile in West Antarctica. Accumulation rates are shown to be highly variable over short distances. Elevation measurements from global positioning system surveys show that accumulation rates derived from shallow horizons correlate well with surface undulations, which implies that wind redistribution of snow is the leading cause of this variability. Temporal changes in accumulation rate over 25-185 year intervals are smoothed to along-track length scales comparable to surface undulations in order to identify trends in accumulation that are likely related to changes in climate. Results show that accumulation rates along this profile have decreased in recent decades, which is consistent with core-derived time series of annual accumulation rates measured at the two ends of the radar profile. These results suggest that temporal variability observed in accumulation-rate records from ice cores and GPR profiles can be obscured by spatial influences, although it is possible to resolve temporal signals if the effects of local topography and ice flow are quantified and removed.

Modeling Englacial Radar Attenuation at Siple Dome, West Antarctica, Using Ice Chemistry and Temperature Data

The radar reflectivity of an ice-sheet bed is a primary measurement for discriminating between thawed and frozen beds. Uncertainty in englacial radar attenuation and its spatial variation introduces corresponding uncertainty in estimates of basal reflectivity. Radar attenuation is proportional to ice conductivity, which depends on the concentrations of acid and sea-salt chloride and the temperature of the ice. We synthesize published conductivity measurements to specify an ice-conductivity model and find that some of the dielectric properties of ice at radar frequencies are not yet well constrained. Using depth profiles of ice-core chemistry and borehole temperature and an average of the experimental values for the dielectric properties, we calculate an attenuation rate profile for Siple Dome, West Antarctica. The depth-averaged modeled attenuation rate at Siple Dome (20.0 +/- 5.7 dB km(-1)) is somewhat lower than the value derived from radar profiles (25.3 +/- 1.1 dB km(-1)). Pending more experimental data on the dielectric properties of ice, we can match the modeled and radar-derived attenuation rates by an adjustment to the value for the pure ice conductivity that is within the range of reported values. Alternatively, using the pure ice dielectric properties derived from the most extensive single data set, the modeled depth-averaged attenuation rate is 24.0 +/- 2.2 dB km(-1). This work shows how to calculate englacial radar attenuation using ice chemistry and temperature data and establishes a basis for mapping spatial variations in radar attenuation across an ice sheet.

Seasonal Geophysical Monitoring of Biogenic Gases in a Northern Peatland: Implications for Temporal and Spatial Variability in Free Phase Gas Production Rates

A set of high resolution surface ground penetrating radar (GPR) surveys, combined with elevation rod ( to monitor surface deformation) and gas flux measurements, were used to investigate in situ biogenic gas dynamics within a northern peatland (Caribou Bog, Maine). Gas production rates were directly estimated from the time series of GPR measurements. Spatial variability in gas production was also investigated by comparing two sites with different geological and ecological attributes, showing differences and/or similarities depending on season. One site characterized by thick highly humified peat deposits (5-6 m), wooded heath vegetation and open pools showed large ebullition events during the summer season, with estimated emissions (based on an assumed range of CH(4) concentration) between 100 and 172 g CH(4) m(-2) during a single event. The other site characterized by thinner less humified peat deposits (2-3 m) and shrub vegetation showed much smaller ebullition events during the same season (between 13 and 23 g CH(4) m(-2)). A consistent period of free-phase gas (FPG) accumulation during the fall and winter, enhanced by the frozen surficial peat acting as a confining layer, was followed by a decrease in FPG after the snow/ice melt that released estimated fluxes between 100 and 200 g CH(4) m(-2) from both sites. Estimated FPG production rates during periods of biogenic gas accumulation ranged between 0.22 and 2.00 g CH(4) m(3) d(-1) and reflected strong seasonal and spatial variability associated with differences in temperature, peat soil properties, and/or depositional attributes (e. g., stratigraphy). Periods of decreased atmospheric pressure coincided with short-period increases in biogenic gas flux, including a very rapid decrease in FPG content associated with an ebullition event that released an estimated 39 and 67 g CH(4) m(-2) in less than 3.5 hours. These results provide insights into the spatial and seasonal variability in production and emission of biogenic gases from northern peatlands.

Geophysical and Hydrological Evaluation of Two Bog Complexes in a Northern Peatland: Implications for the Distribution of Biogenic Gases at the Basin Scale

Ground penetrating radar (GPR) was used to determine peat basin geometry and the spatial distribution of free-phase biogenic gasses in two separate units of a northern peatland (Central and Southern Unit of Caribou Bog, Maine). The Central Unit is characterized by a deep basin structure (15 m maximum depth) and a raised (eccentric) bog topographic profile (up to 2 m topographic variation). Here numerous regions of electromagnetic (EM) wave scattering are considered diagnostic of the presence of extensive free-phase biogenic gas. In contrast, the Southern Unit is shallower (8 m maximum depth) and has a slightly convex upwards bog profile (less than 1 m topographic variation), and areas of EM wave scattering are notably absent. The biogenic gas zones interpreted from GPR in the Central Unit are associated with: (1) wooded heath vegetation at the surface, (2) open pools at the surface, (3) high water table elevations near the center of the basin, and (4) a region of overpressure (at approximately 5 m depth) immediately below the zone of free-phase gas accumulation. The latter suggests (1) a transient pressure head associated with low hydraulic conductivity resulting from the biogenic gasses themselves or confining layers in the peat that restrict both gas release and groundwater flow and/or (2) overpressure in the peat column as a result of the gas buildup itself. In contrast, the Southern Unit, where zones of EM scattering are absent, is characterized by: (1) predominantly shrub vegetation, (2) a lack of open pools, (3) only minor variations (less than 1 m) in water table elevation throughout the entire unit; and (4) generally upward groundwater flow throughout the basin. The results illustrate the nonuniformity of free-phase biogenic gas distribution at the peat basin scale and provide insights into the processes and controls associated with CH4 and CO2 accumulation in peatlands.

Controls on the Basal Water Pressure in Subglacial Channels Near the Margin of the Greenland Ice Sheet

Assuming a channelized drainage system in steady state, we investigate the influence of enhanced surface melting on the water pressure in subglacial channels, compared to that of changes in conduit geometry, ice rheology and catchment variations. The analysis is carried out for a specific part of the western Greenland ice-sheet margin between 66 degrees N and 66 degrees 30' N using new high-resolution digital elevation models of the subglacial topography and the ice-sheet surface, based on an airborne ice-penetrating radar survey in 2003 and satellite repeat-track interferometric synthetic aperture radar analysis of European Remote-sensing Satellite 1 and 2 (ERS-1/-2) imagery, respectively. The water pressure is calculated up-glacier along a likely subglacial channel at distances of 1, 5 and 9 km from the outlet at the ice margin, using a modified version of Rothlisberger's equation. Our results show that for the margin of the western Greenland ice sheet, the water pressure in subglacial channels is not sensitive to realistic variations in catchment size and mean surface water input compared to small changes in conduit geometry and ice rheology.