A Model Study of the Seasonal Circulation in the Gulf of Maine

The Princeton Ocean Model is used to study the circulation in the Gulf of Maine and its seasonal transition in response to wind, surface heat flux, river discharge, and the M-2 tide. The model has an orthogonal-curvature linear grid in the horizontal with variable spacing from 3 km nearshore to 7 km offshore and 19 levels in the vertical. It is initialized and forced at the open boundary with model results from the East Coast Forecast System. The first experiment is forced by monthly climatological wind and heat flux from the Comprehensive Ocean Atmosphere Data Set; discharges from the Saint John, Penobscot, Kennebec, and Merrimack Rivers are added in the second experiment; the semidiurnal lunar tide (M-2) is included as part of the open boundary forcing in the third experiment. It is found that the surface heat flux plays an important role in regulating the annual cycle of the circulation in the Gulf of Maine. The spinup of the cyclonic circulation between April and June is likely caused by the differential heating between the interior gulf and the exterior shelf/slope region. From June to December the cyclonic circulation continues to strengthen, but gradually shrinks in size. When winter cooling erodes the stratification, the cyclonic circulation penetrates deeper into the water column. The circulation quickly spins down from December to February as most of the energy is consumed by bottom friction. While inclusion of river discharge changes details of the circulation pattern, the annual evolution of the circulation is largely unaffected. On the other hand, inclusion of the tide results in not only the anticyclonic circulation on Georges Bank but also modifications to the seasonal circulation.

Size-Fractionated Nitrogen Uptake Measurements in the Equatorial Pacific and Confirmation of the Low Si-High-Nitrate Low-Chlorophyll Condition

The equatorial Pacific Ocean is the largest natural source of CO(2) to the atmosphere, and it significantly impacts the global carbon cycle. Much of the large flux of upwelled CO(2) to the atmosphere is due to incomplete use of the available nitrate (NO(3)) and low net productivity. This high-nutrient low-chlorophyll (HNLC) condition of the equatorial upwelling zone (EUZ) has been interpreted from modeling efforts to be due to low levels of silicate ( Si( OH) 4) that limit the new production of diatoms. These ideas were incorporated into an ecosystem model, CoSINE. This model predicted production by the larger phytoplankton and the picoplankton and effects on air-sea CO(2) fluxes in the Pacific Ocean. However, there were no size-fractionated rates available for verification. Here we report the first size-fractionated new and regenerated production rates (obtained with (15)N - NO(3) and (15)N - NH(4) incubations) for the EUZ with the objective of validating the conceptual basis and functioning of the CoSINE model. Specifically, the larger phytoplankton ( with cell diameters > 5 mu m) had greater rates of new production and higher f-ratios (i.e., the proportion of NO(3) to the sum of NO(3) and NH(4) uptake) than the picoplankton that had high rates of NH(4) uptake and low f-ratios. The way that the larger primary producers are regulated in the EUZ is discussed using a continuous chemostat approach. This combines control of Si(OH)(4) production by supply rate (bottom-up) and control of growth rate ( or dilution) by grazing ( top-down control).

Modeling the Effects of Size on Patch Dynamics of an Inert Tracer

Mesoscale iron enrichment experiments have revealed that additional iron affects the phytoplankton productivity and carbon cycle. However, the role of initial size of fertilized patch in determining the patch evolution is poorly quantified due to the limited observational capability and complex of physical processes. Using a three-dimensional ocean circulation model, we simulated different sizes of inert tracer patches that were only regulated by physical circulation and diffusion. Model results showed that during the first few days since release of inert tracer, the calculated dilution rate was found to be a linear function with time, which was sensitive to the initial patch size with steeper slope for smaller size patch. After the initial phase of rapid decay, the relationship between dilution rate and time became an exponential function, which was also size dependent. Therefore, larger initial size patches can usually last longer and ultimately affect biogeochemical processes much stronger than smaller patches.