The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part I - System overview and formulation

The Regional Ocean Modeling System (ROMS) is one of the few community ocean general circulation models for which a 4-dimensional variational data assimilation (4D-Var) capability has been developed. The ROMS 4D-Var capability is unique in that three variants of 4D-Var are supported: a primal formulation of incremental strong constraint 4D-Var (I4D-Var), a dual formulation based on a physical-space statistical analysis system (4D-PSAS), and a dual formulation representer-based variant of 4D-Var (R4D-Var). In each case, ROMS is used in conjunction with available observations to identify a best estimate of the ocean circulation based on a set of a priori hypotheses about errors in the initial conditions, boundary conditions, surface forcing, and errors in the model in the case of 4D-PSAS and R4D-Var. In the primal formulation of I4D-Var the search for the best circulation estimate is performed in the full space of the model control vector, while for the dual formulations of 4D-PSAS and R4D-Var only the sub-space of linear functions of the model state vector spanned by the observations (i.e. the dual space) is searched. In oceanographic applications, the number of observations is typically much less than the dimension of the model control vector, so there are clear advantages to limiting the search to the space spanned by the observations. In the case of 4D-PSAS and R4D-Var, the strong constraint assumption (i.e. that the model is error free) can be relaxed leading to the so-called weak constraint formulation. This paper describes the three aforementioned variants of 4D-Var as they are implemented in ROMS. Critical components that are common to each approach are conjugate gradient descent, preconditioning, and error covariance models, which are also described. Finally, several powerful 4D-Var diagnostic tools are discussed, namely computation of posterior errors, eigenvector analysis of the posterior error covariance, observation impact, and observation sensitivity.     

Drivers of Change in Shallow Coastal Photic Systems: An Introduction to a Special Issue

Coastal ecosystems are characterized by relatively deep, plankton-based estuaries and much shallower systems where light reaches the bottom. These latter systems, including lagoons, bar-built estuaries, the fringing regions of deeper systems, and other systems of only a few meters deep, are characterized by a variety of benthic primary producers that augment and, in many cases, dominate the production supplied by phytoplankton. These “shallow coastal photic systems” are subject to a wide variety of both natural and anthropogenic drivers and possess numerous natural “filters” that modulate their response to these drivers; in many cases, the responses are much different from those in deeper estuaries. Natural drivers include meteorological forcing, freshwater inflow, episodic events such as storms, wet/dry periods, and background loading of optically active constituents. Anthropogenic drivers include accelerated inputs of nutrients and sediments, chemical contaminants, physical alteration and hydrodynamic manipulation, climate change, the presence of intensive aquaculture, fishery harvests, and introduction of exotic species. The response of these systems is modulated by a number of factors, notably bathymetry, physical flushing, fetch, sediment type, background light attenuation, and the presence of benthic autotrophs, suspension feeding bivalves, and fringing tidal wetlands. Finally, responses to stressors in these systems, particularly anthropogenic nutrient enrichment, consist of blooms of phytoplankton, macroalgae, and epiphytic algae, including harmful algal blooms, subsequent declines in submerged aquatic vegetation and loss of critical habitat, development of hypoxia/anoxia particularly on short time scales (i.e., “diel-cycling”), fish kills, and loss of secondary production. This special issue of Estuaries and Coasts serves to integrate current understanding of the structure and function of shallow coastal photic systems, illustrate the many drivers that cause change in these systems, and synthesize their varied responses.     

Verification of mid-ocean ballast water exchange using naturally occurring coastal tracers

We examined methods for verifying whether or not ships have performed mid-ocean ballast water exchange (BWE) on four commercial vessels operating in the Pacific and Atlantic Oceans. During BWE, a ship replaces the coastal water in its ballast tanks with water drawn from the open ocean, which is considered to harbor fewer organisms capable of establishing in coastal environments. We measured concentrations of several naturally occurring chemical tracers (salinity, six trace elements, colored dissolved organic matter fluorescence and radium isotopes) along ocean transects and in ballast tanks subjected to varying degrees of BWE (0-99%). Many coastal tracers showed significant concentration changes due to BWE, and our ability to detect differences between exchanged and unexchanged ballast tanks was greatest under multivariate analysis. An expanded dataset, which includes additional geographic regions, is now needed to test the generality of our results.     

H2O–CO2 solubility in mafic alkaline magma: applications to volatile sources and degassing behavior at Erebus volcano, Antarctica

We present new equilibrium mixed-volatile (H₂O–CO₂) solubility data for a phonotephrite from Erebus volcano, Antarctica. H₂O–CO₂-saturated experiments were conducted at 400–700 MPa, 1,190 °C, and ~NNO + 1 in non-end-loaded piston cylinders. Equilibrium H₂O–CO₂ fluid compositions were determined using low-temperature vacuum manometry, and the volatile and major element compositions of the glassy run products were determined by Fourier transform infrared spectroscopy and electron microprobe. Results show that the phonotephrite used in this study will dissolve ~0.8 wt% CO₂ at 700 MPa and a fluid composition of $ X_{{{\text{H}}_{ 2} {\text{O}}}} $ ~0.4, in agreement with previous experimental studies on mafic alkaline rocks. Furthermore, the dissolution of CO₂ at moderate to high $ X_{{{\text{H}}_{ 2} {\text{O}}}}^{\text{fluid}} $ in our experiments exceeds that predicted using lower-pressure experiments on similar melts from the literature, suggesting a departure from Henrian behavior of volatiles in the melt at pressures above 400 MPa. With these data, we place new constraints on the modeling of Erebus melt inclusion and gas emission data and thus the interpretation of its magma plumbing system and the contributions of primitive magmas to passive and explosive degassing from the Erebus phonolite lava lake.     

Clinopyroxene–liquid thermometers and barometers specific to alkaline differentiated magmas

We present new thermometers and barometers based on clinopyroxene–liquid equilibria specific to alkaline differentiated magmas. The new models were calibrated through the regression analyses of experimental datasets obtained by merging phase equilibria experiments from the literature with new experiments performed by using trachytic and phonolitic starting compositions. The regression strategy was twofold: (1) we have tested previous thermometric and barometric equations and recalibrated these models using the new datasets; (2) we have calibrated a new thermometer and a new barometer including only regression parameters that closely describe the compositional variability of the datasets. The new models yield more precise estimates than previous thermometers and barometers when used to predict temperatures and pressures of alkaline differentiated magmas. We have tested the reliability of the new equations by using clinopyroxene–liquid pairs from trachytes and phonolites erupted during major explosive eruptions at the Phlegrean Fields and Mt. Vesuvius (central Italy). The test yielded crystallization conditions comparable to those determined by means of melt and fluid inclusion analyses and phase equilibria studies; this validates the use of the proposed models for precise estimates of crystallization temperatures and pressures in differentiated alkaline magmas. Because these magmas feed some of the most voluminous, explosive, and threatening volcanic eruptions in the world, a better understanding of the environmental conditions of their reservoirs is mandatory and this is now possible with the new models provided here.     

The Morphological Evolution of Galaxy Satellites

We study the evolution of galaxy satellites with high resolutionN-body simulations. Satellites are modelled as replicas of typical low and high surface brightness galaxies (LSBs and HSBs). Encounters on high eccentricity orbits (as typical in hierarchical models of galaxy formation) strip LSBs of most of their stars and tend to decrease their surface brightness. In contrast, bar instability in HSBs leads to substantial loss of angular momentum of the stellar component and to an increase in central surface brightness. In both cases the remnant resembles a spheroidal galaxy with an exponential surface brightness profile. A simple modelling of colour evolution and interaction-driven star formation gives M/L ratios for the remnants that are roughly consistent with observations. These results suggest an evolutionary scenario for the dwarf galaxies in our Local Group, faint dSphs being the descendants of LSBs and brighter dSphs/dEs being the final state of HSB satellites. This revised version was published online in July 2006 with corrections to the Cover Date.