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

Physical forcing and biological response within the California Current System (CCS) are highly variable over a wide range of scales. Satellite remote sensing offers the only feasible means of quantifying this variability over the full extent of the CCS. Using six years (1997-2003) of daily SST and chlorophyll imagery, we map the spatial dependence of dominant temporal variability at resolutions sufficient to identify recurrent mesoscale circulation and local pattern associated with coastal topography. Here we describe mean seasonal cycles and interannual variation; intraseasonal variability is left to a companion paper ( K. R. Legaard and A. C. Thomas, manuscript in preparation, 2006). Coastal upwelling dictates seasonality along north-central California, where weak cycles of SST fluctuate between spring minima and late summer maxima and chlorophyll peaks in early summer. Off northern California, chlorophyll maxima are bounded offshore by the seasonally recurrent upwelling jet. Seasonal cycles differ across higher latitudes and in the midlatitude Southern California Bight, where upwelling winds are less vigorous and/or persistent. Seasonality along south-central Baja is strongly affected by processes other than upwelling, despite year-round upwelling-favorable winds. Interannual variation is generally dominated by El Nino and La Nina conditions. Interannual SST variance is greatest along south-central Baja, although interannual variability constitutes a greater fraction of total variance inshore along southern Oregon and much of California. Patterns of interannual chlorophyll variance are consistent with dominant forcing through the widespread depression and elevation of the nutricline during El Nino and La Nina, respectively. Interannual variability constitutes a greater fraction of total chlorophyll variance offshore.

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.

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.

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.

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.

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.

Oxygen Isotopic and Soluble Ionic Composition of a Shallow Firn Core, Inilchek Glacier, Central Tien Shan

Oxygen isotopic and soluble ionic measurements made on snow-pit (2 in depth) and firn-core (12.4 m depth samples recovered from the accumulation zone 5100 m) of Inilchek glacier 43degrees N, 79degrees E) provide information on recent (1992-98) climatic and environmental conditions in the central Tien Shan region of central Asia. The combined 14.4 m snow-pit/firn-core profile lies within the firn zone, arid contains only one observed melt feature (10 m temperature = - 12 degreesC), Although some post-depositional attenuation of the sub-seasonal delta(18)O record is possible, annual cycles are apparent throughout the isotope profile. We therefore use the preserved delta(18)O record to establish a depth/age scale for the core. Mean delta(18)O values for the entire core and for summer periods are consistent with delta(18)O/temperature observations, and suggest the delta(18)O record provides a means to reconstruct past changes in summer surface temperature at the site. Major-ion (Na(+), K(+), Mg(2+), Ca(2+), NH(4)(+), Cl(-), NO(3)(-), SO(4)(2-)) data from the core demonstrate the dominant influence of dust deposition on the soluble chemistry at the site, arid indicate significant interannual variability in atmospheric-dust loading during the 1900s. Anthropogenic impacts oil NH(4)(+) concentrations are observed at the site, and suggest a summer increase in atmospheric NH(4)(+) that may be related to regional agricultural (nitrogen-rich fertilizer use activities.

An Aquaculture-Based Method for Calibrated Bivalve Isotope Paleothermometry

To quantify species- specific relationships between bivalve carbonate isotope geochemistry ( delta O-18(c)) and water conditions ( temperature and salinity, related to water isotopic composition [delta O-18(w)]), an aquaculture-based methodology was developed and applied to Mytilus edulis ( blue mussel). The four- by- three factorial design consisted of four circulating temperature baths ( 7, 11, 15, and 19 degrees C) and three salinity ranges ( 23, 28, and 32 parts per thousand ( ppt); monitored for delta O-18(w) weekly). In mid- July of 2003, 4800 juvenile mussels were collected in Salt Bay, Damariscotta, Maine, and were placed in each configuration. The size distribution of harvested mussels, based on 105 specimens, ranged from 10.9 mm to 29.5 mm with a mean size of 19.8 mm. The mussels were grown in controlled conditions for up to 8.5 months, and a paleotemperature relationship based on juvenile M. edulis from Maine was developed from animals harvested at months 4, 5, and 8.5. This relationship [ T degrees C = 16.19 (+/- 0.14) - 4.69 (+/- 0.21) {delta O-18(c) VPBD - delta O-18(w) VSMOW} + 0.17 (+/- 0.13) {delta O-18(c) VPBD - delta O-18(w) VSMOW}(2); r(2) = 0.99; N = 105; P < 0.0001] is nearly identical to the Kim and O'Neil ( 1997) abiogenic calcite equation over the entire temperature range ( 7 - 19 degrees C), and it closely resembles the commonly used paleotemperature equations of Epstein et al. ( 1953) and Horibe and Oba ( 1972). Further, the comparison of the M. edulis paleotemperature equation with the Kim and O'Neil ( 1997) equilibrium- based equation indicates that M. edulis specimens used in this study precipitated their shell in isotopic equilibrium with ambient water within the experimental uncertainties of both studies. The aquaculture- based methodology described here allows similar species- specific isotope paleothermometer calibrations to be performed with other bivalve species and thus provides improved quantitative paleoenvironmental reconstructions.

The Effects of Joint ENSO-Antarctic Oscillation Forcing on the McMurdo Dry Valleys, Antarctica

Stable oxygen analyses and snow accumulation rates from snow pits sampled in the McMurdo Dry Valleys have been used to reconstruct variations in summer temperature and moisture availability over the last four decades. The temperature data show a common interannual variability, with strong regional warmings occurring especially in 1984/85, 1995/96 and 1990/91 and profound coolings during 1977/78, 1983/84, 1988/89, 1993/94, and 1996/97. Annual snow accumulation shows a larger variance between sites, but the early 1970s, 1984, 1997, and to a lesser degree 1990/91 are characterized overall by wetter conditions, while the early and late 1980s show low snow accumulation values. Comparison of the reconstructed and measured summer temperatures with the Southern Oscillation Index (SOI) and the Antarctic Oscillation (AAO) yield statistically significant correlations, which improve when phaserelationships are considered. A distinct change in the phase relationship of the correlation is observed, with the SOI-AAO leading over the temperature records by one year before, and lagging by one year after 1988. These results suggest that over the last two decades summer temperatures are influenced by opposing El Niho Southern Oscillation and AAO forcings and support previous studies that identified a change in the Tropical-Antarctic teleconnection between the 1980s and 1990s.

Model-Based Estimates of Calanus Finmarchicus Abundance in the Gulf of Maine

Ocean observing systems and satellites routinely collect a wealth of information on physical conditions in the ocean. With few exceptions, such as chlorophyll concentrations, information on biological properties is harder to measure autonomously. Here, we present a system to produce estimates of the distribution and abundance of the copepod Calanus finmarchicus in the Gulf of Maine. Our system uses satellite-based measurements of sea surface temperature and chlorophyll concentration to determine the developmental and reproductive rates of C. finmarchicus. The rate information then drives a population dynamics model of C. finmarchicus that is embedded in a 2-dimensional circulation field. The first generation of this system produces realistic information on interannual variability in C. finmarchicus distribution and abundance during the winter and spring. The model can also be used to identify key drivers of interannual variability in C. finmarchicus. Experiments with the model suggest that changes in initial conditions are overwhelmed by variability in growth rates after approximately 50 d. Temperature has the largest effect on growth rate. Elevated chlorophyll during the late winter can lead to increased C. finmarchicus abundance during the spring, but the effect of variations in chlorophyll concentrations is secondary to the other inputs. Our system could be used to provide real-time estimates or even forecasts of C. finmarchicus distribution. These estimates could then be used to support management of copepod predators such as herring and right whales.

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.

Dispersal Modeling of Fish Early Life Stages: Sensitivity with Application to Atlantic Cod in the Western Gulf of Maine

As an initial step in establishing mechanistic relationships between environmental variability and recruitment in Atlantic cod Gadhus morhua along the coast of the western Gulf of Maine, we assessed transport success of larvae from major spawning grounds to nursery areas with particle tracking using the unstructured grid model FVCOM (finite volume coastal ocean model). In coastal areas, dispersal of early planktonic life stages of fish and invertebrate species is highly dependent on the regional dynamics and its variability, which has to be captured by our models. With state-of-the-art forcing for the year 1995, we evaluate the sensitivity of particle dispersal to the timing and location of spawning, the spatial and temporal resolution of the model, and the vertical mixing scheme. A 3 d frequency for the release of particles is necessary to capture the effect of the circulation variability into an averaged dispersal pattern of the spawning season. The analysis of sensitivity to model setup showed that a higher resolution mesh, tidal forcing, and current variability do not change the general pattern of connectivity, but do tend to increase within-site retention. Our results indicate strong downstream connectivity among spawning grounds and higher chances for successful transport from spawning areas closer to the coast. The model run for January egg release indicates 1 to 19 % within-spawning ground retention of initial particles, which may be sufficient to sustain local populations. A systematic sensitivity analysis still needs to be conducted to determine the minimum mesh and forcing resolution that adequately resolves the complex dynamics of the western Gulf of Maine. Other sources of variability, i.e. large-scale upstream forcing and the biological environment, also need to be considered in future studies of the interannual variability in transport and survival of the early life stages of cod.