Weak and strong constraint data assimilation in the inverse Regional Ocean Modeling System (ROMS): Development and application for a baroclinic coastal upwelling system

We describe the development and preliminary application of the inverse Regional Ocean Modeling System (ROMS), a four dimensional variational (4DVAR) data assimilation system for high-resolution basin-wide and coastal oceanic flows. Inverse ROMS makes use of the recently developed perturbation tangent linear (TL), representer tangent linear (RP) and adjoint (AD) models to implement an indirect representer-based generalized inverse modeling system. This modeling framework is modular. The TL, RP and AD models are used as stand-alone sub-models within the Inverse Ocean Modeling (IOM) system described in [Chua, B.S., Bennett, A.F., 2001. An inverse ocean modeling system. Ocean Modell. 3, 137–165.]. The system allows the assimilation of a wide range of observation types and uses an iterative algorithm to solve nonlinear assimilation problems. The assimilation is performed either under the perfect model assumption (strong constraint) or by also allowing for errors in the model dynamics (weak constraints). For the weak constraint case the TL and RP models are modified to include additional forcing terms on the right hand side of the model equations. These terms are needed to account for errors in the model dynamics.Inverse ROMS is tested in a realistic 3D baroclinic upwelling system with complex bottom topography, characterized by strong mesoscale eddy variability. We assimilate synthetic data for upper ocean (0–450 m) temperatures and currents over a period of 10 days using both a high resolution and a spatially and temporally aliased sampling array. During the assimilation period the flow field undergoes substantial changes from the initial state. This allows the inverse solution to extract the dynamically active information from the synthetic observations and improve the trajectory of the model state beyond the assimilation window. Both the strong and weak constraint assimilation experiments show forecast skill greater than persistence and climatology during the 10–20 days after the last observation is assimilated.Further investigation in the functional form of the model error covariance and in the use of the representer tangent linear model may lead to improvement in the forecast skill.     

From strike-slip faulting to oblique subduction: A survey of the Alpine Fault-Puysegur Trench transition,New Zealand,results of cruise Geodynz-sud leg 2

The Geodynz-sud cruise on board the R/V l'Atalante collected bathymetric, side-scan sonar and seismic reflection data along the obliquely convergent boundary between the Australian and Pacific plates southwest of the South Island, New Zealand. The survey area extended from 44°05 S to 49°40 S, covering the transition zone between the offshore extension of the Alpine Fault and the Puysegur Trench and Puysegur Ridge. Based on variations in the nature and structure of the crust on either side of the margin, the plate boundary zone can be divided into three domains with distinctive structural and sedimentary characteristics. The northern domain involves subduction of probably thinned continental crust of the southern Challenger Plateau beneath the continental crust of Fiordland. It is characterized by thick sediments on the downgoing slab and a steep continental slope disrupted by fault scarps and canyons. The middle domain marks the transition between subduction of likely continental and oceanic crust defined by a series of en echelon ridges on the downgoing slab. This domain is characterized by a large collapse terrace on the continental slope which appears to be due to the collision of the en echelon ridges with the plate margin. The southern domain involves subduction of oceanic crust beneath continental and oceanic crust. This domain is characterized by exposed fabric of seafloor spreading on the downgoing slab, a steep inner trench wall and linear ridges and valleys on the Puysegur ridge crest. The data collected on this cruise provide insights into the nature and history of both plates, and factors influencing the distribution of strike-slip and compressive strain and the evolution of subduction processes along a highly oblique convergent margin.     

An acidification model (MAGIC) with organic acids evaluated using whole-catchment manipulations in Norway

The results from two whole-catchment manipulation experiments in Norway are used to evaluate MAGIC (Model of Acidification of Groundwater In Catchments), a model of ecosystem biogeochemical response to acid deposition. MAGIC is an aggregated catchment-scale model of acidification that has been widely used in assessment activities in Europe and North America. The experiments involved artificial decrease (Reversing Acidification In Norway, RAIN) and increases (Humid Lake Acidification Experiment, HUMEX) in acidic deposition. Runoff from both sites is influenced by moderate levels of organic acids. At each site the model was calibrated to the control catchment and then applied to the manipulated catchments with only minor adjustments. The major responses in runoff chemistry to the manipulations were closely simulated by the model. Differences between simulated and observed volume-weighted annual average concentrations of all major ions were less than ±6 μequiv 1⁻¹ for the entire 4–8 year period of prediction. Trends in response to the manipulations were correctly simulated. Most of the residual error resulted from an inability of the model to reproduce the year-to-year variability (noise) around the trends. Although such model evaluations cannot ‘prove’ the correctness of the model structure, good fits to experimental data increase confidence in model applications for assessment and management purposes. Evaluations of this sort can also identify aspects of the model that need further development. For MAGIC, these are primarily a need for improvement in the calibration of aluminium solubility and a better process basis for nitrogen dynamics.     

Soufrière of guadeloupe 1976–1977 eruption — mass and energy transfer and volcanic health hazards

The Soufrière volcano in Guadeloupe island delivered a phreatic eruption that commenced on July 8th, 1976 and lasted until March 1st, 1977. This eruption was similar to the 1797, 1798, 1809 and 1956 outbreaks. Phreatic activity ejected blocks derived from the fissure walls and fine pyroclasts produced by hydrothermal alteration of the old dome’s rocks. Field observations and measurements allowed the present authors to calculate the mass and energy transfer of steam and ashes: 10⁷ tons of water (very low considering that on the mountain summit the annual precipitation is 10 tons m)²,10⁶ m³ of ashes. The most important energy transfers was thermal: about 5 × 10²⁰ ergs for each phreatic eruption. The total kinetic energy output was 2 × 10¹⁹ ergs for a total thermal energy output of 64 × 10²⁰ ergs. The gases and fine pyroclasts did pollute the atmosphere, waters and soils and consequently affected the population living on the slopes of the volcano.     

Use of On-Site High Performance Liquid Chromatography to Evaluate the Magnitude and Extent of Organic Contaminants in Aquifers

Scientists engaged in the appraisal of ground water contaminated by organic substances are faced with the problems of difficult sample collection and timely chemical analysis. High-performance liquid chromatography was evaluated for on-site determination of specific organic contaminants in ground water samples and was used at three study sites. Organic solutes were determined directly in water samples, with little or no preparation, and usually in less than an hour after collection. This information improved sampling efficiency and was useful in screening for subsequent laboratory analysis. On two occasions, on-site analysis revealed that samples were undergoing rapid change, with major solutes being upgraded and alteration products being formed. In addition to sample stability, this technique proved valuable for monitoring other sampling factors such as compositional changes with respect to pumping, filtration, and cross contamination.     

Sublimation of the Martian CO2 Seasonal South Polar Cap

The polar condensation/sublimation of CO₂, that involve about one fourth of the atmosphere mass, is the major Martian climatic cycle. Early observations in visible and thermal infrared have shown that the sublimation of the Seasonal South Polar Cap (SSPC) is not symmetric around the geographic South Pole.Here we use observations by OMEGA/Mars Express in the near-infrared to detect unambiguously the presence of CO₂ at the surface, and to estimate albedo. Second, we estimate the sublimation of CO₂ released in the atmosphere and show that there is a two-step process. From Ls=180° to 220°, the sublimation is nearly symmetric with a slight advantage for the cryptic region. After Ls=220° the anti-cryptic region sublimation is stronger. Those two phases are not balanced such that there is 22% ± 9 more mass the anti-cryptic region, arguing for more snow precipitation. We compare those results with the MOLA height measurements. Finally we discuss implications for the Martian atmosphere about general circulation and gas tracers, e.g. Ar.