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.

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.

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.

Topographic Control of Regional Accumulation Rate Variability at South Pole and Implications for Ice-Core Interpretation

Snow-accumulation rates are known to be sensitive to local changes in ice-sheet surface slope because of the effect of katabatic winds. These topographic effects can be preserved in ice cores that are collected at non-ice-divide locations. The trajectory of an ice-core site at South Pole is reconstructed using measurements of ice-sheet motion to show that snow was probably deposited at places of different surface slope during the past 1000 years. Recent accumulation rates, derived from shallow firn cores, vary along this trajectory according to surface topography, so that on a relatively steep flank mean annual accumulation is similar to 18% smaller than on a nearby topographic depression. These modern accumulation rates are used to reinterpret the cause of accumulation rate variability with time in the long ice-core record as an ice-dynamics effect and not a climate-change signal. The results highlight the importance of conducting ancillary ice-dynamics measurements as part of ice-coring programs so that topographic effects can be deconvolved from potential climate signals.

Climate Variability in West Antarctica Derived from Annual Accumulation-Rate Records from ITASE Firn/Ice Cores

Thirteen annually resolved accumulation-rate records covering the last similar to 200 years from the Pine Island-Thwaites and Ross drainage systems and the South Pole are used to examine climate variability over West Antarctica. Accumulation is controlled spatially by the topography of the ice sheet, and temporally by changes in moisture transport and cyclonic activity. A comparison of mean accumulation since 1970 at each site to the long-term mean indicates an increase in accumulation for sites located in the western sector of the Pine Island-Thwaites drainage system. Accumulation is negatively associated with the Southern Oscillation Index (Sol) for sites near the ice divide, and periods of sustained negative Sol (1940-42, 1991-95) correspond to above-mean accumulation at most sites. Correlations of the accumulation-rate records with sea-level pressure (SLP) and the SOI suggest that accumulation near the ice divide and in the Ross drainage system may be associated with the midlatitudes. The post-1970 increase in accumulation coupled with strong SLP-accumulation-rate correlations near the coast suggests recent intensification of cyclonic activity in the Pine Island-Thwaites drainage system.

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.

Snow Accumulation Rate on Qomolangma (Mount Everest), Himalaya: Synchroneity With Sites Across the Tibetan Plateau on 50-100 Year Timescales

Annual-layer thickness data, spanning AD 1534-2001, from an ice core from East Rongbuk Coll on Qomolangma (Mount Everest, Himalaya) yield an age-depth profile that deviates systematically from a constant accumulation-rate analytical model. The profile clearly shows that the mean accumulation rate has changed every 50-100 years. A numerical model was developed to determine the magnitude of these multi-decadal-scale rates. The model was used to obtain a time series of annual accumulation. The mean annual accumulation rate decreased from similar to 0.8 m ice equivalent in the 1500s to similar to 0.3 m in the mid-1800s. From similar to 1880 to similar to 1970 the rate increased. However, it has decreased since similar to 1970. Comparison with six other records from the Himalaya and the Tibetan Plateau shows that the changes in accumulation in East Rongbuk Col are broadly consistent with a regional pattern over much of the Plateau. This suggests that there may be an overarching mechanism controlling precipitation and mass balance over this area. However, a record from Dasuopu, only 125 km northwest of Qomolangma and 700 m higher than East Rongbuk Col, shows a maximum in accumulation during the 1800s, a time during which the East Rongbuk Col and Tibetan Plateau ice-core and tree-ring records show a minimum. This asynchroneity may be due to altitudinal or seasonal differences in monsoon versus westerly moisture sources or complex mountain meteorology.

Post-Depositional Movement of Methanesulphonic Acid at Law Dome, Antarctica, and the Influence of Accumulation Rate

A series of ice cores from sites with different snow-accumulation rates across Law Dome, East Antarctica, was investigated for methanesulphonic acid (MSA) movement. The precipitation at these sites (up to 35 km apart) is influenced by the same air masses, the principal difference being the accumulation rate. At the low-accumulation-rate W20k site (0.17 in ice equivalent), MSA was completely relocated from the summer to winter layer. Moderate movement was observed at the intermediate-accumulation-rate site (0.7 in ice equivalent), Dome Summit South (DSS), while there was no evidence of movement at the high-accumulation-rate DE08 site (1.4 in ice equivalent). The main DSS record of MSA covered the epoch AD 1727-2000 and was used to investigate temporal post-depositional changes. Co-deposition of MSA and sea-salt ions was observed of the surface layers, outside of the main summer MSA peak, which complicates interpretation of these peaks as evidence of movement in deeper layers. A seasonal study of the 273 year DSS record revealed MSA migration predominantly from summer into autumn (in the up-core direction), but this migration was suppressed during the Tambora (1815) and unknown (1809) volcanic eruption period, and enhanced during an epoch (1770-1800) with high summer nitrate levels. A complex interaction between the gradients in nss-sulphate, nitrate and sea salts (which are influenced by accumulation rate) is believed to control the rate and extent of movement of MSA.

Regulation of Phytoplankton Carbon to Chlorophyll Ratio By Light, Nutrients and Temperature in the Equatorial Pacific Ocean: A Basin-Scale Model

The complex effects of light, nutrients and temperature lead to a variable carbon to chlorophyll (C:Chl) ratio in phytoplankton cells. Using field data collected in the Equatorial Pacific, we derived a new dynamic model with a non-steady C:Chl ratio as a function of irradiance, nitrate, iron, and temperature. The dynamic model is implemented into a basin-scale ocean circulation-biogeochemistry model and tested in the Equatorial Pacific Ocean. The model reproduces well the general features of phytoplankton dynamics in this region. For instance, the simulated deep chlorophyll maximum (DCM) is much deeper in the western warm pool (similar to 100 m) than in the Eastern Equatorial Pacific (similar to 50 m). The model also shows the ability to reproduce chlorophyll, including not only the zonal, meridional and vertical variations, but also the interannual variability. This modeling study demonstrates that combination of nitrate and iron regulates the spatial and temporal variations in the phytoplankton C:Chl ratio in the Equatorial Pacific. Sensitivity simulations suggest that nitrate is mainly responsible for the high C:Chl ratio in the western warm pool while iron is responsible for the frontal features in the C:Chl ratio between the warm pool and the upwelling region. In addition, iron plays a dominant role in regulating the spatial and temporal variations of the C:Chl ratio in the Central and Eastern Equatorial Pacific. While temperature has a relatively small effect on the C:Chl ratio, light is primarily responsible for the vertical decrease of phytoplankton C:Chl ratio in the euphotic zone.

Modeling Responses of Diatom Productivity and Biogenic Silica Export to Iron Enrichment in the Equatorial Pacific Ocean

Using a three-dimensional physical-biogeochemical model, we have investigated the modeled responses of diatom productivity and biogenic silica export to iron enrichment in the equatorial Pacific, and compared the model simulation with in situ (IronEx II) iron fertilization results. In the eastern equatorial Pacific, an area of 540,000 km(2) was enhanced with iron by changing the photosynthetic efficiency and silicate and nitrogen uptake kinetics of phytoplankton in the model for a period of 20 days. The vertically integrated Chl a and primary production increased by about threefold 5 days after the start of the experiment, similar to that observed in the IronEx II experiment. Diatoms contribute to the initial increase of the total phytoplankton biomass, but decrease sharply after 10 days because of mesozooplankton grazing. The modeled surface nutrients (silicate and nitrate) and TCO(2) anomaly fields, obtained from the difference between the "iron addition'' and "ambient'' (without iron) concentrations, also agreed well with the IronEx II observations. The enriched patch is tracked with an inert tracer similar to the SF6 used in the IronEx II. The modeled depth-time distribution of sinking biogenic silica (BSi) indicates that it would take more than 30 days after iron injection to detect any significant BSi export out of the euphotic zone. Sensitivity studies were performed to establish the importance of fertilized patch size, duration of fertilization, and the role of mesozooplankton grazing. A larger size of the iron patch tends to produce a broader extent and longer-lasting phytoplankton blooms. Longer duration prolongs phytoplankton growth, but higher zooplankton grazing pressure prevents significant phytoplankton biomass accumulation. With the same treatment of iron fertilization in the model, lowering mesozooplankton grazing rate generates much stronger diatom bloom, but it is terminated by Si(OH)(4) limitation after the initial rapid increase. Increasing mesozooplankton grazing rate, the diatom increase due to iron addition stays at minimum level, but small phytoplankton tend to increase. The numerical model experiments demonstrate the value of ecosystem modeling for evaluating the detailed interaction between biogeochemical cycle and iron fertilization in the equatorial Pacific.

Influence of Equatorial Diatom Processes on Si Deposition and Atmospheric CO(2) Cycles at Glacial/Interglacial Timescales

The causes of the glacial cycle remain unknown, although the primary driver is changes in atmospheric CO(2), likely controlled by the biological pump and biogeochemical cycles. The two most important regions of the ocean for exchange of CO(2) with the atmosphere are the equatorial Pacific and the Southern Ocean ( SO), the former a net source and the latter a net sink under present conditions. The equatorial Pacific has been shown to be a Si(OH)(4)-limited ecosystem, a consequence of the low source Si(OH)(4) concentrations in upwelled water that has its origin in the SO. This teleconnection for nutrients between the two regions suggests an oscillatory relationship that may influence or control glacial cycles. Opal mass accumulation rate (MAR) data and delta(15)N measurements in equatorial cores are interpreted with predictions from a one- dimensional Si(OH)(4)-limited ecosystem model (CoSINE) for the equatorial Pacific. The results suggest that equatorial Pacific surface CO(2) processes are in opposite phase to that of the global atmosphere, providing a negative feedback to the glacial cycle. This negative feedback is implemented through the effect of the SO on the equatorial Si(OH)(4) supply. An alternative hypothesis, that the whole ocean becomes Si(OH)(4) poor during cooling periods, is suggested by low opal MAR in cores from both equatorial and Antarctic regions, perhaps as a result of low river input. terminations in this scenario would result from blooms of coccolithophorids triggered by low Si(OH)(4) concentrations.

Physical and Biological Controls on the Latitudinal Asymmetry of Surface Nutrients and pCO(2) in the Central and Eastern Equatorial Pacific

Surface nutrients and dissolved inorganic carbon (DIC) in the central (CEP) and eastern equatorial Pacific (EEP) show much higher concentrations to the south than to the north of the equator. In this study, the physical and biological controls on this asymmetry are investigated using a coupled physical-biogeochemical model. Two numerical experiments are conducted to examine the effects of asymmetrical photosynthetic efficiency (a) due to asymmetrical iron supply about the equator. The experiment with asymmetrical photosynthesis produces improved results as compared with historical observations. A nitrate budget analysis suggests that in the EEP the divergence of upwelling waters controls the surface nitrate asymmetry with additional contribution from the South Equatorial Current (SEC) carrying nutrient-rich Peru upwelling water. The changes of a affect the surface nitrate distribution but not the overall asymmetry. The SEC further carries excess nitrate to the west and thus extends the asymmetry in the east to the CEP. In the CEP, however, stronger northward than southward transport tends to reduce the nitrate asymmetry, while the asymmetrical photosynthesis would help to maintain it. Similar processes also control the distributions of surface silicate and DIC in the equatorial Pacific, which is also affected by the air-sea CO(2) exchange. The asymmetrical photosynthesis influences the distribution of surface DIC, pCO(2), and the air-sea CO(2) flux, by redistributing about 20% CO(2) flux from the north to the south of the equator. Owing to the adjustment of air-sea CO(2) flux, however, the net surface DIC change is smaller than the direct change associated with primary production.

The Role of Forest Soils in a Northern New England Effluent Management System

The Carrabassett Valley Sanitary District in Carrabassett Valley, Maine has utilized both a forest spray irrigation system and a Snowfluent™ system for the treatment of their wastewater effluent. This study was designed to evaluate potential changes in soil properties after approximately 20 years of treatment in the forested spray irrigation site and three years of treatment in the field Snowfluent™ site. In addition, grass yield and composition were evaluated on the field study sites. After treatment with effluent or Snowfluent™, soils showed an increase in soil exchangeable Ca, Mg, Na, and K, base saturation, and pH. While most constituents were higher in treated soils, available P was lower in treated soils compared to the controls. This difference was attributed to higher rates of P mineralization from soil organic matter due to an irrigation effect of the treatment, depleting available P pools despite the P addition with the treatment. Most of the differences due to treatment were greatest at the surface and diminished with depth. Depth patterns in soil properties mostly reflected the decreasing influence of organic matter and its decomposition products with depth as evidenced by significantly higher total C in the surface compared to lower horizons. There were decreasing concentrations of total N, and exchangeable or extractable Ca, Mg, Na, K, Mn, Zn, and P with depth. In addition, there was decreasing BS with depth, driven primarily by declining exchangeable Ca and Mg. Imgation with Snowfluent™ altered the chemical composition of the grass on the site. All element concentrations were significantly higher in the grass foliage except for Ca. The differences were attributed to the additional nutrients and moisture derived from the Snowfluent™. The use of forest spray imgation and Snowfluent™ as a wastewater treatment strategy appears to work well. The soil and vegetation were able to retain most of the applied nutrients, and do not appear to be moving toward saturation. Vegetation management may be a key tool for managing nutrient accumulation on the grass sites as the system ages.