A high resolution operational oil spill model at Santander Bay (Spain): Implementation and validation

A High Resolution Operational Oceanography System that provides decision makers with short-term (within 48 hours) oil spill trajectory forecasting at local scale, has been developed in the Bay of Santander (Spain). The system is based on process models applied in a set of nested grids. Hydrodynamics in the study area are calculated with the COAWST modelling system which uses daily boundary conditions and meteorological forcing obtained from the European network MYOCEAN (http://www.myocean.eu/) and from the Spanish met office, AEMET, respectively. Daily COAWST’s outputs and meteorological forecast are ready to be used by the oil spill transport and fate model, TESEO. A web service that manages the operational system and allows the user to run hypothetical as well as real oil spill trajectories has been implemented. Data from two hydrodynamics field campaigns and from experimental tests carried out with two types of oils have been used to validate the hydrodynamic and the oil spill models.This work has been co-funded by project SPRES in the framework of the EU-Atlantic Area Programme and by project PLVMA (TRA2011-28900-Spanish Ministry of Economy and Competitiveness). Authors B.P. and M.C. would like to thank the Spanish Ministry for its support within the FPI Program. The authors would like to thank AEMET, AZTI-Tecnalia, Cedre, MyOcean and Puertos del Estado for the data and collaboration provided in this project

Historic 2005 toxic bloom of Alexandrium fundyense in the western Gulf of Maine : 2. Coupled biophysical numerical modeling

Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 113 (2008): C07040, doi:10.1029/2007JC004602.A coupled physical/biological modeling system was used to hindcast a massive Alexandrium fundyense bloom that occurred in the western Gulf of Maine in 2005 and to investigate the relative importance of factors governing the bloom's initiation and development. The coupled system consists of a state-of-the-art, free-surface primitive equation Regional Ocean Modeling System (ROMS) tailored for the Gulf of Maine (GOM) using a multinested configuration, and a population dynamics model for A. fundyense. The system was forced by realistic momentum and buoyancy fluxes, tides, river runoff, observed A. fundyense benthic cyst abundance, and climatological nutrient fields. Extensive comparisons were made between simulated (both physical and biological) fields and in situ observations, revealing that the hindcast model is capable of reproducing the temporal evolution and spatial distribution of the 2005 bloom. Sensitivity experiments were then performed to distinguish the roles of three major factors hypothesized to contribute to the bloom: (1) the high abundance of cysts in western GOM sediments; (2) strong ‘northeaster' storms with prevailing downwelling-favorable winds; and (3) a large amount of fresh water input due to abundant rainfall and heavy snowmelt. Model results suggest the following. (1) The high abundance of cysts in western GOM was the primary factor of the 2005 bloom. (2) Wind-forcing was an important regulator, as episodic bursts of northeast winds caused onshore advection of offshore populations. These downwelling favorable winds accelerated the alongshore flow, resulting in transport of high cell concentrations into Massachusetts Bay. A large regional bloom would still have happened, however, even with normal or typical winds for that period. (3) Anomalously high river runoff in 2005 resulted in stronger buoyant plumes/currents, which facilitated the transport of cell population to the western GOM. While affecting nearshore cell abundance in Massachusetts Bay, the buoyant plumes were confined near to the coast, and had limited impact on the gulf-wide bloom distribution.Research support was provided through the Woods Hole Center for Oceans and Human Health, National Science Foundation (NSF) grant OCE-0430723 and National Institute of Environmental Health Science (NIEHS) grant 1-P50-ES012742-01, ECOHAB program through NSF grant OCE-9808173 and NOAA grant NA96OP0099, and GOMTOX program through NOAA grant NA06NOS4780245

Oceanic three-dimensional Lagrangian Coherent Structures: A study of a mesoscale eddy in the Benguela ocean region

We study three dimensional oceanic Lagrangian Coherent Structures (LCSs) in the Benguela region, as obtained from an output of the ROMS model. To do that we first compute Finite-Size Lyapunov exponent (FSLE) fields in the region volume, characterizing mesoscale stirring and mixing. Average FSLE values show a general decreasing trend with depth, but there is a local maximum at about 100 m depth. LCSs are extracted as ridges of the calculated FSLE fields. They present a "curtain-like" geometry in which the strongest attracting and repelling structures appear as quasivertical surfaces. LCSs around a particular cyclonic eddy, pinched off from the upwelling front are also calculated. The LCSs are confirmed to provide pathways and barriers to transport in and out of the eddy

Development of a Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) Modeling System

This paper is not subject to U.S. copyright. The definitive version was published in Ocean Modelling 35 (2010): 230-244, doi:10.1016/j.ocemod.2010.07.010.Understanding the processes responsible for coastal change is important for managing our coastal resources, both natural and economic. The current scientific understanding of coastal sediment transport and geology suggests that examining coastal processes at regional scales can lead to significant insight into how the coastal zone evolves. To better identify the significant processes affecting our coastlines and how those processes create coastal change we developed a Coupled Ocean–Atmosphere–Wave–Sediment Transport (COAWST) Modeling System, which is comprised of the Model Coupling Toolkit to exchange data fields between the ocean model ROMS, the atmosphere model WRF, the wave model SWAN, and the sediment capabilities of the Community Sediment Transport Model. This formulation builds upon previous developments by coupling the atmospheric model to the ocean and wave models, providing one-way grid refinement in the ocean model, one-way grid refinement in the wave model, and coupling on refined levels. Herein we describe the modeling components and the data fields exchanged. The modeling system is used to identify model sensitivity by exchanging prognostic variable fields between different model components during an application to simulate Hurricane Isabel during September 2003. Results identify that hurricane intensity is extremely sensitive to sea surface temperature. Intensity is reduced when coupled to the ocean model although the coupling provides a more realistic simulation of the sea surface temperature. Coupling of the ocean to the atmosphere also results in decreased boundary layer stress and coupling of the waves to the atmosphere results in increased bottom stress. Wave results are sensitive to both ocean and atmospheric coupling due to wave–current interactions with the ocean and wave growth from the atmosphere wind stress. Sediment resuspension at regional scale during the hurricane is controlled by shelf width and wave propagation during hurricane approach

Persistent ocean monitoring with underwater gliders: Adapting sampling resolution

Ocean processes are dynamic and complex and occur on multiple spatial and temporal scales. To obtain a synoptic view of such processes, ocean scientists collect data over long time periods. Historically, measurements were continually provided by fixed sensors, e.g., moorings, or gathered from ships. Recently, an increase in the utilization of autonomous underwater vehicles has enabled a more dynamic data acquisition approach. However, we still do not utilize the full capabilities of these vehicles. Here we present algorithms that produce persistent monitoring missions for underwater vehicles by balancing path following accuracy and sampling resolution for a given region of interest, which addresses a pressing need among ocean scientists to efficiently and effectively collect high-value data. More specifically, this paper proposes a path planning algorithm and a speed control algorithm for underwater gliders, which together give informative trajectories for the glider to persistently monitor a patch of ocean. We optimize a cost function that blends two competing factors: maximize the information value along the path while minimizing deviation from the planned path due to ocean currents. Speed is controlled along the planned path by adjusting the pitch angle of the underwater glider, so that higher resolution samples are collected in areas of higher information value. The resulting paths are closed circuits that can be repeatedly traversed to collect long-term ocean data in dynamic environments. The algorithms were tested during sea trials on an underwater glider operating off the coast of southern California, as well as in Monterey Bay, California. The experimental results show improvements in both data resolution and path reliability compared to previously executed sampling paths used in the respective regions.United States. National Oceanic and Atmospheric Administration. Monitoring and Event Response for Harmful Algal Blooms (NA05NOS4781228)National Science Foundation (U.S.). Center for Embedded Networked Sensing (CCR-0120778)National Science Foundation (U.S.). (Grant number CNS-0520305)National Science Foundation (U.S.). (Grant number CNS-0540420)United States. Office of Naval Research. Multidisciplinary University Research Initiative (N00014-09-1-1031)United States. Office of Naval Research. Multidisciplinary University Research Initiative (N00014-08-1-0693)United States. Office of Naval Research. Service-Oriented Architectur

Caribbean Current variability and the influence of the Amazon and Orinoco freshwater plumes

Author Posting. © Elsevier, 2007. This is the author's version of the work. It is posted here by permission of Elsevier for personal use, not for redistribution. The definitive version was published in Deep Sea Research Part I: Oceanographic Research Papers 54 (2007): 1451-1473, doi:10.1016/j.dsr.2007.04.021.The variability of the Caribbean Current is studied in terms of the influence on its dynamics of the freshwater inflow from the Orinoco and Amazon rivers. Sea-surface salinity maps of the eastern Caribbean and SeaWiFS color images show that a freshwater plume from the Orinoco and Amazon Rivers extends seasonally northwestward across the Caribbean basin, from August to November, 3 to 4 months after the peak of the seasonal rains in northeastern South America. The plume is sustained by two main inflows from the North Brazil Current and its current rings. The southern inflow enters the Caribbean Sea south of Grenada Island and becomes the main branch of the Caribbean Current in the southern Caribbean. The northern inflow (14°N) passes northward around the Grenadine Islands and St. Vincent. As North Brazil Current rings stall and decay east of the Lesser Antilles, between 14°N and 18°N, they release freshwater into the northern part of the eastern Caribbean Sea merging with inflow from the North Equatorial Current. Velocity vectors derived from surface drifters in the eastern Caribbean indicate three westward flowing jets: (1) the southern and fastest at 11°N; (2) the center and second fastest at 14°N; (3) the northern and slowest at 17°N. The center jet (14°N) flows faster between the months of August and December and is located near the southern part of the freshwater plume. Using the MICOM North Atlantic simulation, it is shown that the Caribbean Current is seasonally intensified near 14°N, partly by the inflow of river plumes. Three to four times more anticyclonic eddies are formed during August-December, which agrees with a pronounced rise in the number of anticyclonic looper days in the drifter data then. A climatology-forced regional simulation embedding only the northern (14°N) Caribbean Current (without the influence of the vorticity of the NBC rings), using the ROMS model, shows that the low salinity plume coincides with a negative potential vorticity anomaly that intensifies the center jet located at the salinity front. The jet forms cyclones south of the plume, which are moved northwestward as the anticyclonic circulation intensifies in the eastern Caribbean Sea, north of 14°N. Friction on the shelves of the Greater Antilles also generates cyclones, which propagate westward and eastward from 67°W.The study was supported by National Science Foundation grants OCE 03-271808 and OCE 01-36477

Modulation of Wind Work by Oceanic Current Interaction with the Atmosphere

In this study uncoupled and coupled ocean-atmosphere simulations are carried out for the California Upwelling System to assess the dynamic ocean-atmosphere interactions, viz.,the ocean surface current feedback to the atmosphere. We show the current feedback by modulating the energy transfer from the atmosphere to the ocean, controls the oceanic Eddy Kinetic Energy (EKE). For the first time, we demonstrate the current feedback has an effect on the surface stress and an counteracting effect on the wind itself. The current feedback acts as an oceanic eddy killer, reducing by half the surface EKE, and by 27% the depth-integrated EKE. On one hand, it reduces the coastal generation of eddies by weakening the surface stress and hence the near-shore supply of positive wind work (i.e., the work done by the wind on the ocean). On the other hand, by inducing a surface stress curl opposite to the current vorticity, it deflects energy from the geostrophic current into the atmosphere and dampens eddies. The wind response counteracts the surface stress response. It partly re-energizes the ocean in the coastal region and decreases the offshore return of energy to the atmosphere. Eddy statistics confirm the current feedback dampens the eddies and reduces their lifetime, improving the realism of the simulation. Finally, we propose an additional energy element in the Lorenz diagram of energy conversion, viz., the current-induced transfer of energy from the ocean to the atmosphere at the eddy scale

Towards an integrated observation and modeling system in the New York Bight using variational methods. Part I : 4DVAR data assimilation

Author Posting. © The Author(s), 2010. This is the author's version of the work. It is posted here by permission of Elsevier B.V. for personal use, not for redistribution. The definitive version was published in Ocean Modelling 35 (2010): 119-133, doi:10.1016/j.ocemod.2010.08.003.Four-dimensional Variational data assimilation (4DVAR) in the Regional Ocean Modeling System (ROMS) is used to produce a best-estimate analysis of ocean circulation in the New York Bight during spring 2006 by assimilating observations collected by a variety of instruments during an intensive field program. An incremental approach is applied in an overlapped cycling system with 3-day data assimilation window to adjust model initial conditions. The model-observation mismatch for all observed variables is reduced substantially. Comparisons between model forecast and independent observations show improved forecast skill for about 15 days for temperature and salinity, and 2 to 3 days for velocity. Tests assimilating only certain subsets of the data indicate that assimilating satellite sea surface temperature improves the forecast of surface and subsurface temperature but worsens the salinity forecast. Assimilating in situ temperature and salinity from gliders improves the salinity forecast but has little effect on temperature. Assimilating HF-radar surface current data improves the velocity forecast by 1-2 days yet worsens the forecast of subsurface temperature. During some time periods the convergence for velocity is poor as a result of the data assimilation system being unable to reduce errors in the applied winds because surface forcing is not among the control variables. This study demonstrates the capability of 4DVAR data assimilation system to reduce model-observation mismatch and improve forecasts in the coastal ocean, and highlights the value of accurate meteorological forcing.This work was funded by National Science Foundation grant OCE-0238957

Variational analysis of drifter positions and model outputs for the reconstruction of surface currents in the central Adriatic during fall 2002

Author Posting. © American Geophysical Union, 2008. This article is posted here by permission of American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research 113 (2008): C04004, doi:10.1029/2007JC004148.In this paper we present an application of a variational method for the reconstruction of the velocity field in a coastal flow in the central Adriatic Sea, using in situ data from surface drifters and outputs from the ROMS circulation model. The variational approach, previously developed and tested for mesoscale open ocean flows, has been improved and adapted to account for inhomogeneities on boundary current dynamics over complex bathymetry and coastline and for weak Lagrangian persistence in coastal flows. The velocity reconstruction is performed using nine drifter trajectories over 45 d, and a hierarchy of indirect tests is introduced to evaluate the results as the real ocean state is not known. For internal consistency and impact of the analysis, three diagnostics characterizing the particle prediction and transport, in terms of residence times in various zones and export rates from the boundary current toward the interior, show that the reconstruction is quite effective. A qualitative comparison with sea color data from the MODIS satellite images shows that the reconstruction significantly improves the description of the boundary current with respect to the ROMS model first guess, capturing its main features and its exchanges with the interior when sampled by the drifters.Four of the authors are supported by the Office of Naval Research, V.T. and A.G. under grants N00014-05-1-0094 and N00014-05-1-0095, P.M.P. under grant N00014-03-1-0291, and S.C. under grant N00014-05-1-0730. CNR-ISMAR activity was partially supported by P.O.R. ‘‘CAINO’’ (Regione Puglia), VECTOR (Italian MIUR) project, and ECOOP (EU project)

Interannual and seasonal variabilities in air-sea CO2 fluxes along the US eastern continental shelf and their sensitivity to increasing air temperatures and variable winds

Uncertainty in continental shelf air-sea CO2 fluxes motivated us to investigate the impact of interannual and seasonal variabilities in atmospheric forcing on the capacity of three shelf regions along the U.S. eastern continental shelf to act as a sink or source of atmospheric CO2. Our study uses a coupled biogeochemical-circulation model to simulate scenarios of present-day and future-perturbed mesoscale forcing variability. Overall, the U.S. eastern continental shelf acts as a sink for atmospheric CO2. There is a clear gradient in air-sea CO2 flux along the shelf region, with estimates ranging from -0.6MtCyr(-1) in the South Atlantic Bight (SAB) to -1.0MtCyr(-1) in the Mid-Atlantic Bight (MAB) and -2.5MtCyr(-1) in the Gulf of Maine (GOM). These fluxes are associated with considerable interannual variability, with the largest interannual signal exhibited in the Gulf of Maine. Seasonal variability in the fluxes is also evident, with autumn and winter being the strongest CO2 sink periods and summer months exhibiting some outgassing. In our future-perturbed scenario spatial differences tend to cancel each other out when the fluxes are integrated over the MAB and GOM, resulting in only minor differences between future-perturbed and present-day air-sea CO2 fluxes. This is not the case in the SAB where the position of the along-shelf gradient shifts northward and the SAB becomes a source of CO2 to the atmosphere (0.7MtCyr(-1)) primarily in response to surface warming. Our results highlight the importance of temperature in regulating air-sea CO2 flux variability