Decadal Variability in the Northeast Pacific in a Physical-Ecosystem Model: Role of Mixed Layer Depth and Trophic Interactions

A basin-wide interdecadal change in both the physical state and the ecology of the North Pacific occurred near the end of 1976. Here we use a physical-ecosystem model to examine whether changes in the physical environment associated with the 1976-1977 transition influenced the lower trophic levels of the food web and if so by what means. The physical component is an ocean general circulation model, while the biological component contains 10 compartments: two phytoplankton, two zooplankton, two detritus pools, nitrate, ammonium, silicate, and carbon dioxide. The model is forced with observed atmospheric fields during 1960-1999. During spring, there is a similar to 40% reduction in plankton biomass in all four plankton groups during 1977-1988 relative to 1970-1976 in the central Gulf of Alaska (GOA). The epoch difference in plankton appears to be controlled by the mixed layer depth. Enhanced Ekman pumping after 1976 caused the halocline to shoal, and thus the mixed layer depth, which extends to the top of the halocline in late winter, did not penetrate as deep in the central GOA. As a result, more phytoplankton remained in the euphotic zone, and phytoplankton biomass began to increase earlier in the year after the 1976 transition. Zooplankton biomass also increased, but then grazing pressure led to a strong decrease in phytoplankton by April followed by a drop in zooplankton by May: Essentially, the mean seasonal cycle of plankton biomass was shifted earlier in the year. As the seasonal cycle progressed, the difference in plankton concentrations between epochs reversed sign again, leading to slightly greater zooplankton biomass during summer in the later epoch.

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.

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.

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.

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.

Influence of Surface Oceanographic Variability on Abundance of the Western Winter-Spring Cohort of Neon Flying Squid Ommastrephes Bartramii in the Nw Pacific Ocean

Abundance of the Ommastrephes bartramii winter-spring cohort fluctuated greatly from 1995 to 2004. To understand how abundance was influenced by sea surface conditions, we examined the variations in the proportion of thermal habitats with favourable sea surface temperature (SST). The SST data of both the spawning and feeding grounds were used to calculate the monthly proportion of favourable-SST areas (PFSSTA). Catch per fishing day per fishing boat (catch per unit effort, CPUE) of the Chinese mainland squid-jigging fleet was used as squid abundance index. The relationships between CPUE and monthly PFSSTA at spawning and feeding grounds were analyzed, and the relationship between CPUE and selected PFSSTA was quantified with a multiple linear regression model. Results showed that February PFSSTA at the spawning ground and August to November PFSSTA at the feeding ground could account for about 60% of the variability in O. bartramii abundance between 1995 and 2004, that February was the most important period influencing squid recruitment during the spawning season, and that feeding ground PFSSTA during the fishing season would influence CPUE by causing squid to aggregate. Our forecast model was found to perform well when we compared the model-predicted CPUEs and the average CPUEs observed during August to November in 2005 and 2006 from the Chinese squid-jigging fishery.

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.