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.

A High-Resolution Record of Atmospheric Dust Composition and Variability Since Ad 1650 from a Mount Everest Ice Core

A Mount Everest ice core analyzed at high resolution for major and trace elements (Sr, Cs, Ba, La, Ce, Pr, Nd, Sm, Eu, Tb, Dy, Ho, Er, Tm, Yb, Lu, Bi, U, Tl, Al, S, Ca, Ti, V, Cr, Mn, Fe, Co) and spanning the period A. D. 1650- 2002 is used to investigate the sources of and variations in atmospheric dust through time. The chemical composition of dust varies seasonally, and peak dust concentrations occur during the winter-spring months. Significant correlations between the Everest dust record and dust observations at stations suggest that the Everest record is representative of regional variations in atmospheric dust loading. Back-trajectory analysis in addition to a significant correlation of Everest dust concentrations and the Total Ozone Mapping Spectrometer (TOMS) aerosol index indicates that the dominant winter sources of dust are the Arabian Peninsula, Thar Desert, and northern Sahara. Factors that contribute to dust generation at the surface include soil moisture and temperature, and the long-range transport of dust aerosols appears to be sensitive to the strength of 500-mb zonal winds. There are periods of high dust concentration throughout the 350-yr Mount Everest dust record; however, there is an increase in these periods since the early 1800s. The record was examined for recent increases in dust emissions associated with anthropogenic activities, but no recent dust variations can be conclusively attributed to anthropogenic inputs of dust.

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.

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.

Spatial and Temporal Melt Variability at Helheim Glacier, East Greenland, and Its Effect on Ice Dynamics

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.

Satellite-Derived Variability in Chlorophyll, Wind Stress, Sea Surface Height, and Temperature in the Northern California Current System

Satellite-derived data provide the temporal means and seasonal and nonseasonal variability of four physical and biological parameters off Oregon and Washington ( 41 degrees - 48.5 degrees N). Eight years of data ( 1998 - 2005) are available for surface chlorophyll concentrations, sea surface temperature ( SST), and sea surface height, while six years of data ( 2000 - 2005) are available for surface wind stress. Strong cross-shelf and alongshore variability is apparent in the temporal mean and seasonal climatology of all four variables. Two latitudinal regions are identified and separated at 44 degrees - 46 degrees N, where the coastal ocean experiences a change in the direction of the mean alongshore wind stress, is influenced by topographic features, and has differing exposure to the Columbia River Plume. All these factors may play a part in defining the distinct regimes in the northern and southern regions. Nonseasonal signals account for similar to 60 - 75% of the dynamical variables. An empirical orthogonal function analysis shows stronger intra-annual variability for alongshore wind, coastal SST, and surface chlorophyll, with stronger interannual variability for surface height. Interannual variability can be caused by distant forcing from equatorial and basin-scale changes in circulation, or by more localized changes in regional winds, all of which can be found in the time series. Correlations are mostly as expected for upwelling systems on intra-annual timescales. Correlations of the interannual timescales are complicated by residual quasi-annual signals created by changes in the timing and strength of the seasonal cycles. Examination of the interannual time series, however, provides a convincing picture of the covariability of chlorophyll, surface temperature, and surface height, with some evidence of regional wind forcing.

Spatial Patterns of Intraseasonal Variability of Chlorophyll and Sea Surface Temperature in the California Current

Six years of daily satellite data are used to quantify and map intraseasonal variability of chlorophyll and sea surface temperature (SST) in the California Current. We define intraseasonal variability as temporal variation remaining after removal of interannual variability and stationary seasonal cycles. Semivariograms are used to quantify the temporal structure of residual time series. Empirical orthogonal function (EOF) analyses of semivariograms calculated across the region isolate dominant scales and corresponding spatial patterns of intraseasonal variability. The mode 1 EOFs for both chlorophyll and SST semivariograms indicate a dominant timescale of similar to 60 days. Spatial amplitudes and patterns of intraseasonal variance derived from mode 1 suggest dominant forcing of intraseasonal variability through distortion of large scale chlorophyll and SST gradients by mesoscale circulation. Intraseasonal SST variance is greatest off southern Baja and along southern Oregon and northern California. Chlorophyll variance is greatest over the shelf and slope, with elevated values closely confined to the Baja shelf and extending farthest from shore off California and the Pacific Northwest. Intraseasonal contributions to total SST variability are strongest near upwelling centers off southern Oregon and northern California, where seasonal contributions are weak. Intraseasonal variability accounts for the majority of total chlorophyll variance in most inshore areas save for southern Baja, where seasonal cycles dominate. Contributions of higher EOF modes to semivariogram structure indicate the degree to which intraseasonal variability is shifted to shorter timescales in certain areas. Comparisons of satellite-derived SST semivariograms to those calculated from co-located and concurrent buoy SST time series show similar features.

Phytoplankton Scales of Variability in the California Current System: 1. Interannual and Cross-Shelf Variability

In the California Current System, strong mesoscale variability associated with eddies and meanders of the coastal jet play an important role in the biological productivity of the area. To assess the dominant timescales of variability, a wavelet analysis is applied to almost nine years (October 1997 to July 2006) of 1-km-resolution, 5-day-averaged, Sea-viewing Wide Field-of-view Sensor (SeaWiFS) chlorophyll a (chl a) concentration data. The dominant periods of chlorophyll variance, and how these change in time, are quantified as a function of distance offshore. The maximum variance in chlorophyll occurs with a period of similar to 100-200 days. A seasonal cycle in the timing of peak variance is revealed, with maxima in spring/summer close to shore (20 km) and in autumn/winter 200 km offshore. Interannual variability in the magnitude of chlorophyll variance shows maxima in 1999, 2001, 2002, and 2005. There is a very strong out-of-phase correspondence between the time series of chlorophyll variance and the Pacific Decadal Oscillation (PDO) index. We hypothesize that positive PDO conditions, which reflect weak winds and poor upwelling conditions, result in reduced mesoscale variability in the coastal region, and a subsequent decrease in chlorophyll variance. Although the chlorophyll variance responds to basin-scale forcing, chlorophyll biomass does not necessarily correspond to the phase of the PDO, suggesting that it is influenced more by local-scale processes. The mesoscale variability in the system may be as important as the chl a biomass in determining the potential productivity of higher trophic levels.

Interannual Variability in Timing of Bloom Initiation in the California Current System

In the California Current System the spring transition from poleward to equatorward alongshore wind stress heralds the beginning of upwelling-favorable conditions. The phytoplankton response to this transition is investigated using 8 years ( 1998-2005) of daily, 4-km resolution, Sea-viewing Wide Field of view Sensor ( SeaWiFS) chlorophyll a concentration data. Cluster analysis of the chlorophyll a time series at each location is used to separate the inshore upwelling region from offshore and oligotrophic areas. An objective method for estimating the timing of bloom initiation is used to construct a map of the mean bloom start date. Interannual variability in bloom timing and magnitude is investigated in four regions: 45 degrees N - 50 degrees N, 40 degrees N - 45 degrees N, 35 degrees N - 40 degrees N and 20 degrees N - 35 degrees N. Daily satellite derived wind data ( QuikSCAT) allow the timing of the first episode of persistently upwelling favorable winds to be estimated. Bloom initiation generally coincides with the onset of upwelling winds ( +/- 15 days). South of similar to 35 degrees N, where winds are southward year-round, the timing of increased chlorophyll concentration corresponds closely to timing of the seasonal increase in upwelling intensity. A 1-D model and satellite derived photosynthetically available radiation data are used to estimate time series of depth- averaged irradiance. In the far north of the region (> 46 degrees N) light is shown to limit phytoplankton growth in early spring. In 2005 the spring bloom in the northern regions (> 35 degrees N) had a "false start''. A sharp increase in chl a in February quickly receded, and a sustained increase in biomass was delayed until July. We hypothesize that this resulted in a mismatch in timing of food availability to higher trophic levels.

Satellite-Measured Chlorophyll and Temperature Variability Off Northern Chile During the 1996-1998 La Niña and El Niño

Time series of satellite measurements are used to describe patterns of surface temperature and chlorophyll associated with the 1996 cold La Nina phase and the 1997-1998 warm El Nino phase of the El Nino - Southern Oscillation cycle in the upwelling region off northern Chile. Surface temperature data are available through the entire study period. Sea-viewing Wide Field-of-view Sensor (SeaWiFS) data first became available in September 1997 during a relaxation in El Nino conditions identified by in situ hydrographic data. Over the time period of coincident satellite data, chlorophyll patterns closely track surface temperature patterns. Increases both in nearshore chlorophyll concentration and in cross-shelf extension of elevated concentrations are associated with decreased coastal temperatures during both the relaxation in El Nino conditions in September-November 1997 and the recovery from EI Nino conditions after March 1998. Between these two periods during austral summer (December 1997 to March 1998) and maximum El Nino temperature anomalies, temperature patterns normally associated with upwelling were absent and chlorophyll concentrations were minimal. Cross-shelf chlorophyll distributions appear to be modulated by surface temperature frontal zones and are positively correlated with a satellite-derived upwelling index. Frontal zone patterns and the upwelling index in 1996 imply an austral summer nearshore chlorophyll maximum, consistent with SeaWiFS data from I 1998-1999, after the El Nino. SeaWiFS retrievals in the data set used here are higher than in situ measurements by a factor of 2-4; however, consistency in the offset suggests relative patterns are valid.

Twentieth Century Delta δC Variability in Surface Water Dissolved Inorganic Carbon Recorded by Coralline Algae in the Northern North Pacific Ocean and the Bering Sea

The oxygen isotopic composition and Mg/Ca ratios in the skeletons of long-lived coralline algae record ambient seawater temperature over time. Similarly, the carbon isotopic composition in the skeletons record delta(13)C values of ambient seawater dissolved inorganic carbon. Here, we measured delta(13)C in the coralline alga Clathromorphum nereostratum to test the feasibility of reconstructing the intrusion of anthropogenic CO(2) into the northern North Pacific Ocean and Bering Sea. The delta(13)C was measured in the high Mgcalcite skeleton of three C. nereostratum specimens from two islands 500 km apart in the Aleutian archipelago. In the records spanning 1887 to 2003, the average decadal rate of decline in delta(13)C values increased from 0.03% yr(-1) in the 1960s to 0.095% yr(-1) in the 1990s, which was higher than expected due to solely the delta(13)C-Suess effect. Deeper water in this region exhibits higher concentrations of CO(2) and low delta(13)C values. Transport of deeper water into surface water (i.e., upwelling) increases when the Aleutian Low is intensified. We hypothesized that the acceleration of the delta(13)C decline may result from increased upwelling from the 1960s to 1990s, which in turn was driven by increased intensity of the Aleutian Low. Detrended delta(13)C records also varied on 4-7 year and bidecadal timescales supporting an atmospheric teleconnection of tropical climate patterns to the northern North Pacific Ocean and Bering Sea manifested as changes in upwelling.

Seasonal Variability in the Bacteriolytic Capacity of the Deposit Feeder Arenicola Marina: Environmental Correlates

Although deposit-feeding macrofauna consume and digest sedimentary bacteria, it is unclear whether feeding rates and digestion efficiencies are high enough to significantly impact the composition and abundance of bacteria in marine sediments. It is likely that both feeding rates and efficiency of digestion vary markedly through space and time. We used a turbidimetric assay to compare the rate of bacteriolysis by digestive fluids collected seasonally from the deposit-feeding polychaete Arenicola marina. Under standardized, experimental conditions, bacteriolytic rates represent concentrations of lytic agents. This concentration was found to vary significantly throughout the year (p = 0.001), showing greater than a 2x range. Lytic agent concentration was positively correlated with bioavailable amino acid concentrations in the surface sediment (r = 0.85, p = 0.03) but showed no apparent relationship to other proxies for food resources (e.g, chl a), sediment temperature, or gut throughput time. In vitro, temperature has been shown to have a strong positive influence on bacteriolytic rate. Temperature has no influence, however, on the in situ concentration of lytic agent in gut fluids, thus it appears that compensation for this temperature dependence is unimportant. These findings, combined with previous kinetics studies with A. marina gut fluids, predict that the quantitative influence of deposit feeding on the microbial ecology of sediments will exhibit clear seasonal variation.

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.

North Pacific Carbon Cycle Response to Climate Variability on Seasonal to Decadal Timescales

Climate variability drives significant changes in the physical state of the North Pacific, and there may be important impacts of this variability on the upper ocean carbon balance across the basin. We address this issue by considering the response of seven biogeochemical ocean models to climate variability in the North Pacific. The models' upper ocean pCO(2) and air-sea CO(2) flux respond similarly to climate variability on seasonal to decadal timescales. Modeled seasonal cycles of pCO(2) and its temperature- and non-temperature-driven components at three contrasting oceanographic sites capture the basic features found in observations (Takahashi et al., 2002, 2006; Keeling et al., 2004; Brix et al., 2004). However, particularly in the Western Subarctic Gyre, the models have difficulty representing the temporal structure of the total pCO(2) seasonal cycle because it results from the difference of these two large and opposing components. In all but one model, the air-sea CO(2) flux interannual variability (1 sigma) in the North Pacific is smaller ( ranges across models from 0.03 to 0.11 PgC/yr) than in the Tropical Pacific ( ranges across models from 0.08 to 0.19 PgC/yr), and the time series of the first or second EOF of the air-sea CO(2) flux has a significant correlation with the Pacific Decadal Oscillation (PDO). Though air-sea CO(2) flux anomalies are correlated with the PDO, their magnitudes are small ( up to +/- 0.025 PgC/yr ( 1 sigma)). Flux anomalies are damped because anomalies in the key drivers of pCO(2) ( temperature, dissolved inorganic carbon (DIC), and alkalinity) are all of similar magnitude and have strongly opposing effects that damp total pCO(2) anomalies.

Distributions and Variability of Particulate Organic Matter in a Coastal Upwelling System

In this study we examined the spatial and temporal variability of particulate organic material (POM) off Oregon during the upwelling season. High-resolution vertical profiling of beam attenuation was conducted along two cross-shelf transects. One transect was located in a region where the shelf is relatively uniform and narrow (off Cascade Head (CH)); the second transect was located in a region where the shelf is shallow and wide (off Cape Perpetua (CP)). In addition, water samples were collected for direct analysis of chlorophyll, particulate organic carbon (POC), and particulate organic nitrogen (PON). Beam attenuation was highly correlated with POC and PON. Striking differences in distribution patterns and characteristics of POM were observed between CH and CP. Off CH, elevated concentrations of chlorophyll and POC were restricted to the inner shelf and were highly variable in time. The magnitude of the observed short-term temporal variability was of the same order as that of the seasonal variability reported in previous studies. Elevated concentrations of nondegraded chlorophyll and POM were observed near the bottom. Downwelling and rapid sinking are two mechanisms by which phytoplankton cells can be delivered to the bottom before being degraded. POM may be then transported across the shelf via the benthic nepheloid layer. Along the CP transect, concentrations of POM were generally higher than they were along the CH transect and extended farther across the shelf. Characteristics of surface POM, namely, C: N ratios and carbon: chlorophyll ratios, differed between the two sites. These differences can be attributed to differences in shelf circulation.