John C. Kindle

Atmospheric forcing in the Arabian Sea during 1994–1995: observations and comparisons with climatology and models

Accurate, year-long time series of winds, incoming shortwave and longwave radiation, air and sea temperatures, relative humidity, barometric pressure, and precipitation were collected from a surface mooring deployed off the coast of Oman along the climatological axis of the Findlater Jet from October 1994 to October 1995. Wind stress, heat flux, and freshwater flux were computed using bulk formulae. The Northeast Monsoon was characterized by steady but moderate winds, clear skies, relatively dry air, and two months, December and January, in which the ocean, on average, lost 45 W m-2 to the atmosphere. The Southwest Monsoon had strong winds, cloudy skies, and moist air. Because of reduced latent and longwave heat loss, it was accompanied by sustained oceanic heat gain, with the strongest monthly mean warming, 147 W m-2, in August. Large differences are found between the observations and older climatologies. Recent climatologies agree better with the observations. The means of the Southampton Oceanography Center climatology for 1980–1995 are close to the buoy monthly means. Monthly means from that climatology show that 1994–1995 was in general a typical year, with surface meteorology and air–sea fluxes within one standard deviation of the long term means. Concurrent data from the NCEP, ECMWF, and FNMOC show that the models provide realistic surface winds. FNMOC winds show that the timing and character of the onset of the Southwest Monsoon in 1995 differed from 1994 and 1996 when variability within one month is resolved. The models fail to replicate other observed surface meteorology and to produce realistic heat fluxes. Annual and monsoonal mean net heat fluxes from the models differed from those of the buoy by 50 to 80 W m-2. Because of these differences, some care is warranted in selecting and using air-sea flux fields in studies of the Arabian Sea.

Euphotic zone depth: Its derivation and implication to ocean‐color remote sensing

[1] Euphotic zone depth, z1%, reflects the depth where photosynthetic available radiation (PAR) is 1% of its surface value. The value of z1% is a measure of water clarity, which is an important parameter regarding ecosystems. Based on the Case-1 water assumption, z1% can be estimated empirically from the remotely derived concentration of chlorophyll-a ([Chl]), commonly retrieved by employing band ratios of remote sensing reflectance (Rrs). Recently, a model based on water’s inherent optical properties (IOPs) has been developed to describe the vertical attenuation of visible solar radiation. Since IOPs can be nearanalytically calculated from Rrs, so too can z1%. In this study, for measurements made over three different regions and at different seasons (z1% were in a range of 4.3–82.0 m with [Chl] ranging from 0.07 to 49.4 mg/m 3 ), z1% calculated from Rrs was compared with z1% from in situ measured PAR profiles. It is found that the z1% values calculated via Rrs-derived IOPs are, on average, within � 14% of the measured values, and similar results were obtained for depths of 10% and 50% of surface PAR. In comparison, however, the error was � 33% when z1% is calculated via Rrs-derived [Chl]. Further, the importance of deriving euphotic zone depth from satellite ocean-color remote sensing is discussed.

Evidence for eddy formation in the eastern Arabian Sea during the northeast monsoon

The seasonal formation of a large (500–800 km diameter) anticyclonic eddy in the upper 300–400 m of the eastern Arabian Sea during the northeast monsoon period (December-April) is indicated from hydrographic and satellite altimetry sea level observations, as well as from numerical model experiments. The center of the eddy circulation is approximately 10°N, 70°E, just to the west of the north-south Laccadive Island chain. In this paper the eddy is called the Laccadive High (LH). In some ways it is a mirrorlike counterpart to the Great Whirl, which develops during the southwest monsoon off the Somali coast (western Arabian Sea). The LH occurs at the same latitude but on the opposite side of the basin during the reversed monsoon. It is different from the Great Whirl, however, in its formation process, its intensity, and its decay. The hydrographic data obtained from surveys all during a single season give sufficiently close station spacing to allow reasonable contouring of the geopotential surfaces and of the properties within and around the LH region with minimum time aliasing. The Geosat altimeter record extends over 4 years, during which the seasonal variability of the LH indicates a dynamic relief of approximately 15–20 cm, which is in good agreement with the hydrographic observations. The altimetry time series also suggests a westward translation of the LH by January with a subsequent dissipation in midbasin. The model used is a wind-forced three-layer primitive equation model which depicts a LH in agreement with the timing, position, and amplitude of both the hydrographic and altimetric measurements. The numerical simulation includes a passive tracer located in the western Bay of Bengal; the western advection of the tracer around the south coasts of Sri Lanka and India in December and January is consistent with the appearance of low-salinity water observed to extend into the Arabian Sea during this period. The modeling studies suggest that both local and remote forcing are important in formation of the LH.

Influences of diurnal and intraseasonal forcing on mixed‐layer and biological variability in the central Arabian Sea

A three-dimensional, physical-biological model of the Indian Ocean is used to study the influences of diurnal and intraseasonal forcing on mixed-layer and biological variability in the central Arabian Sea, where a mooring was deployed and maintained from October 1994 to October 1995 by the Woods Hole Oceanographic Institution Upper Ocean Processes group. The physical model consists of four active layers overlying an inert deep ocean, namely, a surface mixed layer of thickness h1, diurnal thermocline layer, seasonal thermocline, and main thermocline. The biological model consists of a set of advective-diffusive equations in each layer that determine nitrogen concentrations in four compartments: nutrients, phytoplankton, zooplankton, and detritus. Both monthly climatological and “daily” fields are used to force solutions, the latter being a blend of daily-averaged fields measured at the mooring site and other products that include intraseasonal forcing. Diurnal forcing is included by allowing the incoming solar radiation to have a daily cycle. In solutions forced by climatological fields, h1 thickens steadily throughout both monsoons. When h1 detrains at their ends, short-lived, intense blooms develop (the models spring and fall blooms) owing to the increase in depth-averaged light intensity sensed by the phytoplankton in layer 1. In solutions forced by daily fields, h1 thins in a series of events associated with monsoon break periods. As a result, the spring and fall blooms are split into a series of detrainment blooms, broadening them considerably. Diurnal forcing alters the mixed-layer and biological responses, among other things, by lengthening the time that h1 is thick during the northeast monsoon, by strengthening the spring and fall blooms and delaying them by 3 weeks, and by intensifying phytoplankton levels during intermonsoon periods. Solutions are compared with the mixed-layer thickness, phytoplankton biomass, and phytoplankton production fields estimated from mooring observations. The solution driven by daily fields with diurnal forcing reproduces the observed fields most faithfully.

Monsoons boost biological productivity in Arabian Sea

Monsoons over the Arabian Sea—the oceanic basin that separates the Arabian peninsula from the Indian subcontinent—follow seasonal cycles, reversing directions twice a year, in summer and winter. Rather than spreading across the expanse of the sea, the southwest (summer) monsoon is often concentrated into a jet over the central Arabian Sea. Evidence suggests that variations in wind stress force substantial upwelling in the ocean to the west of the jet, and weaker upwelling or even downwelling to the east. This upwelling provides nutrients to the euphotic zone and enhances biological productivity.

High resolution modeling and data assimilation in the Monterey Bay area

Abstract A high resolution, data assimilating ocean model of the Monterey Bay area (ICON model) is under development within the framework of the project “An Innovative Coastal-Ocean Observing Network” (ICON) sponsored by the National Oceanographic Partnership Program. The main objective of the ICON model development is demonstration of the capability of a high resolution model to track the major mesoscale ocean features in the Monterey Bay area when constrained by the measurements and nested within a regional larger-scale model. This paper focuses on the development of the major ICON model components, including grid generation and open boundary conditions, coupling with a larger scale, Pacific West Coast (PWC) model, atmospheric forcing etc. Impact of these components on the Models predictive skills in reproducing major hydrographic conditions in the Monterey Bay area are analyzed. Comparisons between observations and the ICON model predictions with and without coupling to the PWC model, show that coupling with the regional model improves significantly both the correlation between the ICON model and observed ADCP currents, and the ICON models skill in predicting the location and intensity of observed upwelling events. Analysis of the ICON model mixed layer depth predictions show that the ICON model tends to develop a thicker than observed mixed layer during the summer time, and while assimilation of sea surface temperature data is enough for development of observed thin mixed layer in the regional larger-scale model, the fine-resolution ICON model needs variable heat fluxes as surface boundary conditions for the accurate prediction of the vertical thermal structure. The paper targets researchers involved in high-resolution numerical modeling of coastal areas in which the dynamics are determined by the complex geometry of a coastline, variable bathymetry and by the influence of complex water masses from a complicated hydrographic system (such as the California Current system).

The ocean response to operational westerly wind bursts during the 1991 - 1992 El Nino

Numerical simulations of the remotely forced ocean response to westerly wind bursts prior to and during the 1991–1992 El Nino are examined; the models are forced by wind stress from the U.S. Navys atmospheric global operational analysis/forecast system. The study focuses on (1) the relative response of the first and second internal modes to a single episode of westerly bursts; (2) the role of westerly bursts in producing the eastern Pacific sea level variations from October 1990 to February 1992; and (3) the relative importance of the remotely forced sea level responses generated by central and western Pacific wind anomalies. The simulations use the Naval Research Laboratory global multilayer formulation; the suite of experiments includes hydrodynamic simulations that use both one- and three-active-layer reduced gravity models as well as an experiment that also includes thermodynamic effects. The models are executed from January 1, 1990, to March 1, 1992, a period that includes 10 significant westerly wind bursts or burst clusters and the 1991–1992 El Nino event. The numerical experiments reveal an ability to accurately simulate the eastern Pacific sea level variations during this period. In response to a single burst the three-layer hydrodynamic simulation reveals that the second internal mode Kelvin wave yields a sea level change at the eastern boundary that is approximately one third that of the first mode and a surface velocity signature that is equivalent to the first mode. During the onset of the El Nino event the inclusion of higher modes also produces a more realistic representation of the observed eastern boundary sea level signal. Furthermore, by comparing the model response to particular wind bursts with the observed sea level at Baltra, Galapagos, a value of 2.5–2.6 m/s is suggested as the most appropriate mean speed for the first internal mode Kelvin wave during the onset phase. The models reveal the following scenario for the onset of the 1991–1992 El Nino. The eastern boundary sea level exhibited three, distinct pulselike events superimposed on a general rise; the observed sea level pulses occurred in October–November, December–January, and February. The numerical simulations indicate that the pulses in December–January and February were the result of three powerful, westerly wind bursts in the western Pacific preceded by westerly anomalies in the central Pacific. The Kelvin wave pulse generated in the western Pacific by the wind burst in late December to early January was reinforced by a strong central Pacific westerly wind anomaly in mid-January. These westerly wind events were associated with the development of western Pacific tropical cyclones near the equator in one or both hemispheres.

Physical processes affecting availability of dissolved silicate for diatom production in the Arabian Sea

A passive tracer to represent dissolved silicate concentrations, with biologically realistic uptake kinetics, is successfully incorporated into a three-dimensional, eddy-resolving, ocean circulation model of the Indian Ocean. Hypotheses are tested to evaluate physical processes which potentially affect the availability of silicate for diatom production in the Arabian Sea. An alternative mechanism is offered to the idea that open ocean upwelling is primarily responsible for the high, vertical nutrient flux and consequent large-scale phytoplankton bloom in the northwestern Arabian Sea during the southwest monsoon. Model results show that dissolved silicate in surface waters available for uptake by diatoms is primarily influenced by the intensity of nearshore upwelling from southwest monsoonal wind forcing and by the offshore advective transport of surface waters. The upwelling, which in the model occurs within 200±50 km of the coast, appears to be a result of a combination of coastal upwelling, Ekman pumping, and divergence of the coastal flow as it turns offshore. Localized intensifications of silicate concentrations appear to be hydrodynamically driven and geographically correlated to coastal topographic features. The absence of diatoms in sediments of the eastern Arabian Basin is consistent with modeled distributional patterns of dissolved silicate resulting from limited westward advection of upwelled coastal waters from the western continental margin of India and rapid uptake of available silicate by diatoms. Concentrations of modeled silicate become sufficiently low to become unavailable for diatom production in the eastern Arabian Sea, a region between 61°E and 70°E at 8°N on the south, with the east and west boundaries converging on the north at ∼67°E, 20°N.

Modeling of upwelling/relaxation events with the Navy Coastal Ocean Model

Prediction System (COAMPS TM )(COAMPS is a registered trademark of the Naval Research Laboratory). Issues investigated in this study are: NCOM-based model simulations of upwelling and relaxation events, coupling to COAMPS, use of sigma versus hybrid (sigma-z) vertical grids, and coupling with a larger-scale model on the open boundaries. The NCOM simulations were able to reproduce the observed sequence of the upwelling and relaxation events, which can be attributed, in part, to the good agreement between the observed and COAMPS winds. Comparisons with the mooring observations show that COAMPS overestimates shortwave radiation values, which makes the NCOM modeled SSTs too warm in comparison with observations. The NCOM runs forced with different resolution atmospheric forcing (3 versus 9 km) do not show significant differences in the predicted SSTs and mixed-layer depths at the mooring locations. At the same time, during the extended upwelling event, the model runs forced with 3 and 9 km resolution COAMPS fields show differences in the surface circulation patterns, which are the most distinct in the southern portion of the model domain. The model run with 9-km forcing develops a northward flow along the coast, which is not present in the run with 3-km forcing and in observations (for example, HF radar-derived radials). Comparison of the wind patterns of the 3- and 9-km products shows a weakening of the 9-km wind stress along the southern coast of the NCOM model domain, which is responsible for the development of the artificial northward flow in the NCOM run with 9-km forcing.

Bay of Bengal currents during the northeast monsoon

Velocity and property observations were made during February and March 1995 as part of the World Ocean Circulation Experiment (WOCE) Hydrographic Program (WHP) expedition in the Indian Ocean. The observed circulation in the upper 300 m of the ocean during the northeast monsoon is compared to the output of a high-resolution, 3-layer, nonlinear model forced by European Center for Medium Range Weather Forecasting (ECMWF) winds. The data identify several new features in the Bay of Bengal: a shelf-break coastal current in the northeast corner, an eddy and subsurface jet near the South Preparis Channel, and an eastward countercurrent extending from 80°E to the eastern boundary. The countercurrent transports high-salinity surface water from the northwest Indian Ocean into the bay, and separates the historically observed North Equatorial Current into two currents with separate sources and water properties. The data also detail the spatial structure of the South Equatorial Current and Countercurrent, and show the Equatorial Undercurrent in the eastern half of the Indian Ocean. The model compares well enough with the data to suggest that such realistic models may provide a useful temporal context for the WHP snapshot.