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.

Surface "Waves" on Byrd Glacier, Antarctica

Byrd Glacier has one of the largest ice catchment areas in Antarctica, delivers more ice to the Ross Ice Shelf than any other ice stream, and is the fastest of these ice streams. A force balance, combined with a mass balance, demonstrates that stream flow in Byrd Glacier is transitional from sheet flow in East Antarctica to shelf flow in the Ross Ice Shelf. The longitudinal pulling stress, calculated along an ice flowband from the force balance, is linked to variations of ice thickness, to the ratio of the basal water pressure to the ice overburden pressure where Byrd Glacier is grounded, and is reduced by an ice-shelf buttressing stress where Byrd Glacier is floating. Longitudinal tension peaks at pressure-ratio maxima in grounded ice and close to minima in the ratio of the pulling stress to the buttressing stress in floating ice. The longitudinal spacing of these tension peaks is rather uniform and, for grounded ice, the peaks occur at maxima in surface slope that have no clear relation to the bed slope. This implies that the maxima in surface slope constitute a "wave train" that is related to regular variations in ice-bed coupling, not primarily to bed topography. It is unclear whether these surface "waves" are "standing waves" or are migrating either upslope or downslope, possibly causing the grounding line to either retreat or advance. Deciding which is the case will require obtaining bed topography in the map plane, a new map of surface topography, and more sophisticated modeling that includes ice flow linked to subglacial hydrology in the map plane.

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.

Effects of Surface Forcing on Interannual Variability of the Fall Phytoplankton Bloom in the Gulf of Maine Revealed Using a Process-Oriented Model

SeaWiFS (Sea-viewing Wide Field-of-view Sensor) chlorophyll data revealed strong interannual variability in fall phytoplankton dynamics in the Gulf of Maine, with 3 general features in any one year: (1) rapid chlorophyll increases in response to storm events in fall; (2) gradual chlorophyll increases in response to seasonal wind-and cooling-induced mixing that gradually deepens the mixed layer; and (3) the absence of any observable fall bloom. We applied a mixed-layer box model and a 1-dimensional physical-biological numerical model to examine the influence of physical forcing (surface wind, heat flux, and freshening) on the mixed-layer dynamics and its impact on the entrainment of deep-water nutrients and thus on the appearance of fall bloom. The model results suggest that during early fall, the surface mixed-layer depth is controlled by both wind-and cooling-induced mixing. Strong interannual variability in mixed-layer depth has a direct impact on short-and long-term vertical nutrient fluxes and thus the fall bloom. Phytoplankton concentrations over time are sensitive to initial pre-bloom profiles of nutrients. The strength of the initial stratification can affect the modeled phytoplankton concentration, while the timing of intermittent freshening events is related to the significant interannual variability of fall blooms.