Investigating Peatland Stratigraphy and Hydrogeology Using Integrated Electrical Geophysics

Hydrology has been suggested as the mechanism controlling vegetation and related surficial pore-water chemistry in large peatlands. Peatland hydrology influences the carbon dynamics within these large carbon reservoirs and will influence their response to global warming. A geophysical survey was completed in Caribou Bog, a large peatland in Maine, to evaluate peatland stratigraphy and hydrology. Geophysical measurements were integrated with direct measurements of peat stratigraphy from probing, fluid chemistry, and vegetation patterns in the peatland. Consistent with previous field studies, ground-penetrating radar (GPR) was an excellent method for delineating peatland stratigraphy. Prominent reflectors from the peat-lake sediment and lake sediment-mineral soil contacts were precisely recorded up to 8 m deep. Two-dimensional resistivity and induced polarization imaging were used to investigate stratigraphy beneath the mineral soil, beyond the range of GPR. We observe that the peat is chargeable, and that IP imaging is an alternative method for defining peat thickness. The chargeability of peat is attributed to the high surface-charge density on partially decomposed organic matter. The electrical conductivity imaging resolved glaciomarine sediment thickness (a confining layer) and its variability across the basin. Comparison of the bulk conductivity images with peatland vegetation revealed a correlation between confining layer thickness and dominant vegetation type, suggesting that stratigraphy exerts a control on hydrogeology and vegetation distribution within this peatland. Terrain conductivity measured with a Geonics EM31 meter correlated with confining glaciomarine sediment thickness and was an effective method for estimating variability in glaciomarine sediment thickness over approximately 18 km(2). Our understanding of the hydrogeology, stratigraphy, and controls on vegetation growth in this peatland was much enhanced from the geophysical study.

Spatial Patterns in Mass Balance of the Siple Coast and Amundsen Sea Sectors, West Antarctica

Local rates of change in ice-sheet thickness were calculated at IS sites in West Antarctica using the submergence velocity technique. This method entails a comparison of the vertical velocity of the ice sheet, measured using repeat global positioning system surveys of markers, and local long-term rates of snow accumulation obtained using firn-core stratigraphy. Any significant difference between these two quantities represents a thickness change with time. Measurements were conducted at sites located similar to 100-200 km apart along US ITASE traverse routes, and at several isolated locations. All but one of the sites are distributed in the Siple Coast and the Amundsen Sea basin along contours of constant elevation, along flowlines, across ice divides and close to regions of enhanced flow. Calculated rates of thickness change are different from site to site. Most of the large rates of change in ice thickness (similar to 10 cm a(-1) or larger) are observed in or close to regions of rapid flow, and are probably related to ice-dynamics effects. Near-steady-state conditions are calculated mostly at sites in the slow-moving ice-sheet interior and near the main West Antarctic ice divide. These results are consistent with regional estimates of ice-sheet change derived from remote-sensing measurements at similar locations in West Antarctica.

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.

Dating the Siple Dome (Antarctica) Ice Core By Manual and Computer Interpretation of Annual Layering

The Holocene portion of the Siple Dome (Antarctica) ice core was dated by interpreting the electrical, visual and chemical properties of the core. The data were interpreted manually and with a computer algorithm. The algorithm interpretation was adjusted to be consistent with atmospheric methane stratigraphic ties to the GISP2 (Greenland Ice Sheet Project 2) ice core, (BE)-B-10 stratigraphic ties to the dendrochronology C-14 record and the dated volcanic stratigraphy. The algorithm interpretation is more consistent and better quantified than the tedious and subjective manual interpretation.

Mass Balance and Accumulation Rate Across Siple Dome, West Antarctica

Snow-accumulation rates and rates of ice-thickness change (mass balance) are studied at several sites on Siple Dome, West Antarctica. Accumulation rates are derived from analyses of gross beta radioactivity in shallow firn cores located along a 60 km transect spanning both flanks and the crest of the dome. There is a north-south gradient in snow-accumulation rate across the dome that is consistent with earlier radar mapping of internal stratigraphy. Orographic processes probably control this distribution. Mass balance is inferred from the difference between global positioning system (GPS)-derived vertical velocities and snow-accumulation rates for sites close to the firn-core locations. Results indicate that there is virtually no net thickness change at four of the five sites. The exception is at the northernmost site where a small amount of thinning is detected, that appears to be inconsistent with other studies. A possible cause of this anomalous thinning is recent retreat of the grounding line of Ice Stream D.

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.

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.

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.