A global ocean temperature and altimeter data assimilation system for studies of climate variability

 An ocean data assimilation (ODA) system which can assimilate both temperature and altimeter observations has been applied to the global ocean and tested between January 1993–October 1996. A statistical method has been used to convert sea surface height (SSH) anomalies observations from TOPEX/POSEIDON into synthetic temperature profiles. The innovative aspect of this method is the introduction of time dependency in the correlations used to transform the altimeter observations into temperature corrections. The assimilation system is based on a univariate variational optimal interpolation scheme applied to assimilate both in situ and synthetic temperature profiles. In addition, a longer global analysis for the upper-ocean temperature starting from January 1979 and ending November 1997, has been produced to examine the skill of sea temperature assimilation with a rather simple and practical method. The temperature analysis shows encouraging improvement over a corresponding ocean simulation when compared to independent (not assimilated) temperature data both at seasonal and interannual time scales. However, the univariate data assimilation of hydrographic data does not result in an improvement of the velocity field. In fact the assimilation of sparse in situ data can introduce unrealistic spatial variability in the temperature field which affects the velocity field in a negative way. This deficiency is partially overcome when we also assimilate altimeter observations since the coverage is complete and uniform for this data. In particular, our study shows that temperature corrections due to the altimeter signal have a positive impact on the current system in the tropical Pacific. Received: 28 May 2000 / Accepted: 6 November 2000     

Space- and ground-based measurements of sulphur dioxide emissions from Turrialba Volcano (Costa Rica)

Remotely sensed measurements of sulphur dioxide (SO₂) emitted by Turrialba Volcano (Costa Rica) are reported for the period September 2009–January 2011. These measurements were obtained using images from Advanced Spaceborne Thermal Emission and Reflexion radiometer, Ozone Monitoring Instrument and a ground-based UV camera. These three very different instruments provide flux measurements in good agreement with each other, which demonstrate that they can be integrated for monitoring SO₂ fluxes. Fluxes from Turrialba increased fourfold in January 2010, following a phreatic explosion that formed a degassing vent in the W crater of Turrialba. Since then, the SO₂ flux has remained high (30–50?kg/s) but seems to be showing a slowly decreasing trend. We interpret this evolution as the start of open vent degassing from a recently intruded magma body. The opening of the degassing vent decreased the confining pressure of the magma body and allowed the gases to bypass the hydrothermal system.     

Landslide HotSpot Mapping by means of Persistent Scatterer Interferometry

Landslide detection and mapping represent fundamental requirements for every hazard and risk evaluation and consequent improvement of the management strategies for such natural hazards. Optical and radar remote sensing can be used to observe landslide-induced ground deformation, ranging from regional to local scales. This work presents a methodology called Landslide HotSpot Mapping; this approach integrates cartographic, thematic and optical data with Persistent Scatterer Interferometry for the identification of extremely slow to very slow moving landslides, and for the evaluation of their state of activity and intensity. This methodology scans wide areas to detect hotspots, which are narrow unstable zones characterized by higher landslide hazard. To these hotspots, priority has to be given when planning field surveys and in situ validation campaigns, so that field work time and effort can be optimized and significantly reduced. The approach is tested in Central Calabria, over a 4,470?km² area located in southern Italy. ENVISAT ascending images acquired between 2003 and 2009 and processed with the Persistent Scatterer Pairs (PSP) technique are used to analyse deformation patterns. Combining conventional photo-interpretation with the analysis of PSP data, 64 new landslides are identified and the spatial (boundaries) and temporal (activity) information of 980 pre-mapped phenomena (23.6% of updated inventory) are updated. 1,012 active (continuous or reactivated) landslides are identified and 4 hotspot areas selected: San Fili, Rende, Lago, Catanzaro. Urgent field checks have to be organized for these hotspots to validate the satellite-based observations and to design appropriate mitigation measures to reduce impacts on the elements at risk.     

On the assessment of Argo float trajectory assimilation in the Mediterranean Forecasting System

The Mediterranean Forecasting System (MFS) has been operational for a decade, and is continuously providing forecasts and analyses for the region. These forecasts comprise local- and basin-scale information of the environmental state of the sea and can be useful for tracking oil spills and supporting search-and-rescue missions. Data assimilation is a widely used method to improve the forecast skill of operational models and, in this study, the three-dimensional variational (OceanVar) scheme has been extended to include Argo float trajectories, with the objective of constraining and ameliorating the numerical output primarily in terms of the intermediate velocity fields at 350 m depth. When adding new datasets, it is furthermore crucial to ensure that the extended OceanVar scheme does not decrease the performance of the assimilation of other observations, e.g., sea-level anomalies, temperature, and salinity. Numerical experiments were undertaken for a 3-year period (2005–2007), and it was concluded that the Argo float trajectory assimilation improves the quality of the forecasted trajectories with ~15%, thus, increasing the realism of the model. Furthermore, the MFS proved to maintain the forecast quality of the sea-surface height and mass fields after the extended assimilation scheme had been introduced. A comparison between the modeled velocity fields and independent surface drifter observations suggested that assimilating trajectories at intermediate depth could yield improved forecasts of the upper ocean currents.     

Baroclinic wind adjustment processes in the Mediterranean Sea

The wind-driven general circulation of the Mediterranean Sea is studied using a primitive equation model. The model uses a 0.25° horizontal resolution and eight or 16 levels in the vertical. The model uses the Mediterranean basin geometry, and the Strait of Gibraltar is closed. The vertical density structure is initialized with annual average data, and the temperature and salinity values are fixed at the surface to simulate perpetual annual mean conditions. The wind forcing consists of monthly mean climatological stresses.The results show that the general circulation of the Mediterranean Sea has a multiple time-scale character (seasonal excursions and steady state amplitudes are comparable) and it is composed by sub-basin scale gyres corresponding to the scale of the wind stress curl centers. The steady state circulation (annual mean average) is determined by a Sverdrup balacne modified by viscous effects.The unsteady vertically integrated transport circulation consists of sub-basin scale gyres similar to the steady state transport components, which amplify seasonally and the partial or total reversal of the currents in many subportions of the basin. The gyres can be stationary in position or propagating. This seasonal ocean response is partly constituted by Rossby modes due to the wind stress curl annual harmonic. The baroclinic circulation shows the seasonal shift of the North African Current from a position along the African coasts during winter to the center of the Balearic and Ionian basin during summer.     

Impact of data assimilation of glider observations in the Ionian Sea (Eastern Mediterranean)

Glider observations of temperature, salinity and vertically averaged velocity in the Ionian Sea (Eastern Mediterranean Sea), made in the period October 2004–December 2004, were assimilated into an operational forecasting model together with other in situ and satellite observations. The study area has a high spatial and temporal variability of near surface dynamics, characterized by the entrance of the Atlantic Ionian Stream (AIS) into the Northern Ionian Sea. The impact of glider observations on the estimation of the circulation is studied, and it is found that their assimilation locally improves the prediction of temperature, salinity, velocity and surface elevation fields. However, only the assimilation of temperature and salinity together with the vertically averaged velocity improves the forecast of all observed parameters. It is also found that glider observations rapidly impact the analyses even remotely, and the remote impacts on the analyses remain several months after the presence of the glider. The study emphasizes the importance of assimilating as much as possible all available information from gliders, especially in dynamically complex areas.     

Linking 1D coastal ocean modelling to environmental management: an ensemble approach

The use of a one-dimensional interdisciplinary numerical model of the coastal ocean as a tool contributing to the formulation of ecosystem-based management (EBM) is explored. The focus is on the definition of an experimental design based on ensemble simulations, integrating variability linked to scenarios (characterised by changes in the system forcing) and to the concurrent variation of selected, and poorly constrained, model parameters. The modelling system used was previously specifically designed for the use in “data-rich” areas, so that horizontal dynamics can be resolved by a diagnostic approach and external inputs can be parameterised by nudging schemes properly calibrated. Ensembles determined by changes in the simulated environmental (physical and biogeochemical) dynamics, under joint forcing and parameterisation variations, highlight the uncertainties associated to the application of specific scenarios that are relevant to EBM, providing an assessment of the reliability of the predicted changes. The work has been carried out by implementing the coupled modelling system BFM-POM1D in an area of Gulf of Trieste (northern Adriatic Sea), considered homogeneous from the point of view of hydrological properties, and forcing it by changing climatic (warming) and anthropogenic (reduction of the land-based nutrient input) pressure. Model parameters affected by considerable uncertainties (due to the lack of relevant observations) were varied jointly with the scenarios of change. The resulting large set of ensemble simulations provided a general estimation of the model uncertainties related to the joint variation of pressures and model parameters. The information of the model result variability aimed at conveying efficiently and comprehensibly the information on the uncertainties/reliability of the model results to non-technical EBM planners and stakeholders, in order to have the model-based information effectively contributing to EBM.     

Net autotrophy in a fluvial lake: the relative role of phytoplankton and floating-leaved macrophytes

This study combined water- and sediment flux measurements with mass balances of dissolved gas and inorganic matter to determine the importance of pelagic and benthic processes for whole-system metabolism in a eutrophic fluvial lake. Mass balances of dissolved O₂, inorganic carbon (DIC), nitrogen (DIN), phosphorous (SRP), particulate N (PN) and P (PP) and Chl a were calculated at a nearly monthly frequency by means of repeated sampling at the lake inlet and outlet. Simultaneously, benthic fluxes of gas and nutrients, including denitrification rates, and the biomass of the dominant pleustophyte (Trapa natans) were measured, and fluxes of O₂ and CO₂ across the water–atmosphere interface were estimated from diel changes in outlet concentrations. On an annual scale, Middle Lake exhibited CO₂ supersaturation, averaging 313% (range 86–562%), but was autotrophic with a net O₂ production (6.35 ± 2.05 mol m⁻² y⁻¹), DIC consumption (−31.18 ± 18.77 mol m⁻² y⁻¹) and net export of Chl a downstream (8.38 ± 0.95 mol C m⁻² y⁻¹). Phytoplankton was the main driver of Middle Lake metabolism, with a net primary production estimated at 33.24 mol O₂ m⁻² y⁻¹, corresponding to a sequestration of 4.18 and 0.26 mol m⁻² y⁻¹ of N and P, respectively. At peak biomass, T. natans covered about 18% of Middle Lake’s surface and fixed 2.46, 0.17 and 0.02 mol m⁻² of C, N and P, respectively. Surficial sediments were a sink for O₂ (−14.47 ± 0.65 mol O₂ m⁻² y⁻¹) and a source of DIC and NH₄ ⁺ (18.84 ± 2.80 mol DIC m⁻² y⁻¹ and 0.83 ± 0.16 mol NH₄ ⁺ m⁻² y⁻¹), and dissipated nitrate via denitrification (1.44 ± 0.11 mol NO₃ ⁻ m⁻² y⁻¹). Overall, nutrient uptake by primary producers and regeneration from sediments were a minor fraction of external loads. This work suggests that the creation of fluvial lakes can produce net autotrophic systems, with elevated rates of phytoplanktonic primary production, largely sustained by allochtonous nutrient inputs. These hypereutrophic aquatic bodies are net C sinks, although they simultaneously release CO₂ to the atmosphere.     

IT-OSRA: applying ensemble simulations to estimate the oil spill risk associated to operational and accidental oil spills

Oil Spill Risk Assessments (OSRAs) are widely employed to support decision making regarding oil spill risks. This article adapts the ISO-compliant OSRA framework developed by Sepp Neves et al. (J Environ Manag 159:158–168, 2015) to estimate risks in a complex scenario where uncertainties related to the meteo-oceanographic conditions, where and how a spill could happen exist and the risk computation methodology is not yet well established (ensemble oil spill modeling). The improved method was applied to the Algarve coast, Portugal. Over 50,000 simulations were performed in 2 ensemble experiments to estimate the risks due to operational and accidental spill scenarios associated with maritime traffic. The level of risk was found to be important for both types of scenarios, with significant seasonal variations due to the the currents and waves variability. Higher frequency variability in the meteo-oceanographic variables were also found to contribute to the level of risk. The ensemble results show that the distribution of oil concentrations found on the coast is not Gaussian, opening up new fields of research on how to deal with oil spill risks and related uncertainties.