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.

Derived Bedrock Elevations, Strain Rates and Stresses from Measured Surface Elevations and Velocities - Jakobshavns-Isbrae, Greenland

Jakobshavns Isbrae (69 degrees 10'N, 49 degrees 5'W) drains about 6.5% of the Greenland ice sheet and is the fastest ice stream known. The Jakobshavns Isbrae basin of about 10 000 km(2) was mapped photogrammetrically from four sets of aerial photography, two taken in July 1985 and two in July 1986. Positions and elevations of several hundred natural features on the ice surface were determined for each epoch by photogrammetric block-aerial triangulation, and surface velocity vectors were computed from the positions. The two flights in 1985 yielded the best results and provided most common points (716) for velocity determinations and are therefore used in the modeling studies. The data from these irregularly spaced points were used to calculate ice elevations and velocity vectors at uniformly spaced grid paints 3 km apart by interpolation. The field of surface strain rates was then calculated from these gridded data and used to compute the field of surface deviatoric stresses, using the flow law of ice, for rectilinear coordinates, X, Y pointing eastward and northward. and curvilinear coordinates, L, T pointing longitudinally and transversely to the changing ice-flow direction. Ice-surface elevations and slopes were then used to calculate ice thicknesses and the fraction of the ice velocity due to basal sliding. Our calculated ice thicknesses are in fair agreement with an ice-thickness map based on seismic sounding and supplied to us by K. Echelmeyer. Ice thicknesses were subtracted from measured ice-surface elevations to map bed topography. Our calculation shows that basal sliding is significant only in the 10-15 km before Jakobshavns Isbrae becomes afloat in Jakobshavns IsfJord.

Observations of the Sensitivity of Beam Attenuation to Particle Size in a Coastal Bottom Boundary Layer

The goal of this study was to test the hypothesis that the aggregated state of natural marine particles constrains the sensitivity of optical beam attenuation to particle size. An instrumented bottom tripod was deployed at the 12-m node of the Martha's Vineyard Coastal Observatory to monitor particle size distributions, particle size-versus-settling-velocity relationships, and the beam attenuation coefficient (c(p)) in the bottom boundary layer in September 2007. An automated in situ filtration system on the tripod collected 24 direct estimates of suspended particulate mass (SPM) during each of five deployments. On a sampling interval of 5 min, data from a Sequoia Scientific LISST 100x Type B were merged with data from a digital floc camera to generate suspended particle volume size distributions spanning diameters from approximately 2 mu m to 4 cm. Diameter-dependent densities were calculated from size-versus-settling-velocity data, allowing conversion of the volume size distributions to mass distributions, which were used to estimate SPM every 5 min. Estimated SPM and measured c(p) from the LISST 100x were linearly correlated throughout the experiment, despite wide variations in particle size. The slope of the line, which is the ratio of c(p) to SPM, was 0.22 g m(-2). Individual estimates of c(p):SPM were between 0.2 and 0.4 g m(-2) for volumetric median particle diameters ranging from 10 to 150 mu m. The wide range of values in c(p):SPM in the literature likely results from three factors capable of producing factor-of-two variability in the ratio: particle size, particle composition, and the finite acceptance angle of commercial beam-transmissometers.