A Model Study of the Circulation in the Pearl River Estuary (PRE) and Its Adjacent Coastal Waters: 1. Simulations and Comparison with Observations

The Princeton Ocean Model is used to study the circulation features in the Pearl River Estuary and their responses to tide, river discharge, wind, and heat flux in the winter dry and summer wet seasons. The model has an orthogonal curvilinear grid in the horizontal plane with variable spacing from 0.5 km in the estuary to 1 km on the shelf and 15 sigma levels in the vertical direction. The initial conditions and the subtidal open boundary forcing are obtained from an associated larger-scale model of the northern South China Sea. Buoyancy forcing uses the climatological monthly heat fluxes and river discharges, and both the climatological monthly wind and the realistic wind are used in the sensitivity experiments. The tidal forcing is represented by sinusoidal functions with the observed amplitudes and phases. In this paper, the simulated tide is first examined. The simulated seasonal distributions of the salinity, as well as the temporal variations of the salinity and velocity over a tidal cycle are described and then compared with the in situ survey data from July 1999 and January 2000. The model successfully reproduces the main hydrodynamic processes, such as the stratification, mixing, frontal dynamics, summer upwelling, two-layer gravitational circulation, etc., and the distributions of hydrodynamic parameters in the Pearl River Estuary and coastal waters for both the winter and the summer season.

A Model Study of the Circulation in the Pearl River Estuary (PRE) and Its Adjacent Coastal Waters: 2. Sensitivity Experiments

The Princeton Ocean Model is used to study the circulation in the Pear River Estuary (PRE) and the adjacent coastal waters in the winter and summer seasons. Wong et al. [2003] compares the simulation results with the in situ measurements collected during the Pearl River Estuary Pollution Project (PREPP). In this paper, sensitivity experiments are carried out to examine the plume and the associated frontal dynamics in response to seasonal discharges and monsoon winds. During the winter, convergence between the seaward spreading plume water and the saline coastal water sets up a salinity front that aligns from the northeast to the southwest inside the PRE. During the summer the plume water fills the PRE at the surface and spreads eastward in the coastal waters in response to the prevailing southwesterly monsoon. The overall alignment of the plume is from the northwest to the southeast. The subsurface front is similar to that in the winter and summer except that the summer front is closer to the mouth and the winter front closer to the head of the estuary. Inside the PRE, bottom flows are always toward the head of the estuary, attributed to the density gradient associated with the plume front. In contrast, bottom flows in the shelf change from offshore in winter to onshore in summer, reflecting respectively the wintertime downwelling and summertime upwelling. Wind also plays an essential role in controlling the plume at the surface. An easterly wind drives the plume westward regardless winter or summer. The eastward spreading of the plume during the summer can be attributed to the southerly component of the wind. On the other hand, the surface area of the plume is positively proportional to the amount of discharge.

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.