In this paper, we present the operational forecasting system developed to assist in the response to the 2002 Prestige oil spill in Cantabria. The objective of the system developed was to forecast the wave climate, tidal and wind currents, and oil spill trajectories to provide a technical assessment to decision makers for a response to the oil spill. The two main components of the system were data collection and processing and integration with numerical models for forecasting. The information from overflights received daily became essential in achieving a correct initial position of the oil slicks. Meteorological and oceanographic data were also received daily by means of an emergency protocol established between Puertos del Estado (Spain), the Naval Research Laboratory (USA), and the University of Cantabria (Spain). These data were used to run the trajectory model, the wave propagation model, and the shallow depth-integrated flow model. The information generated by the numerical simulations was presented to the decision makers every day in the form of maps that were easy and quick to interpretation as a tool to help in the response planning.In addition, to develop a defensive or protection strategy for sensitive areas like estuaries and marshes, a hydrodynamic study was carried out by the University of Cantabria in all the estuaries of the region. The result of this study consisted of a boom deployment plan for each.
The analysis of three different approximations of the turbulent diffusion terms, widely used to simulate shallow water flows, is carried out both analytically and experimentally. Based on the eddy viscosity concept, the terms are solved for steady, uniform, turbulent flow in a simplified geometry, which may represent a tidal estuary or a compound channel. It is shown that, although the three approximations are identical in constant depth, they behave differently if strong depth gradients exist and, consequently, the transverse velocity profile obtained varies depending on the turbulence term used. It is also shown that the relative depth (flood plain depth-to-main channel depth ratio) has an important influence on the lateral momentum transfer and, consequently, depending on the approximation adopted, different dimensionless eddy viscosity coefficient λ values must be used to best fit the experimental data. A tentative relationship between the relative depth and the dimensionless eddy viscosity coefficient is presented for each expression. Comparison of analytical results with experimental data shows that not all the widely used expressions for the turbulence terms can adequately represent the velocity profile if strong depth gradients exist.
Abstract Salinity is a key variable used to explain the distribution of species within estuaries and the functioning of estuarine ecosystems. Despite the relevance of average conditions and extreme events, the existing schemes of zonation are only based on average values; the extreme values obtained through high‐resolution long‐term data have not been still integrated in an appropriate methodology capable of recognizing salinity types. Therefore, the background variability of salinity has been frequently underestimated in ecological studies. The primary goal of this research is the identification of ecologically significant salinity zones capable of encompassing the entire estuarine regime. A two‐step methodological approach was developed: (1) the reconstruction of long‐term salinity series using the Deflt3D hydrodynamic model, the analog method and a subset of short‐term representative states of river flow and tidal level and (2) the identification of different zones using eight descriptors of the salinity regime (i.e., the median, the range and the intensity, duration and frequency of extreme events) and the application of a combination of two clustering techniques of unsupervised learning (Self‐Organizing Maps (SOM) and K‐Means). Thus, five ecologically significant salinity types that are representative of estuarine variability were identified based on a salinity series reconstructed using a validated method. Differences in the mean values of salinity among typologies allow explaining patterns in the general descriptors of benthic macroinvertebrates assemblages (i.e., richness and diversity). If extreme salinity conditions are also considered, typologies increase their ecological significance and they are able to recognize differences in species composition.
The experience acquired in the Prestige crisis management has demonstrated the importance of forecasting oil slick trajectories to plan an effective oil spill response. To have a reliable prediction system, we need to perform a detailed calibration and validation of the oil spill transport model. In this work, the Lagrangian transport model, PICHI, developed by the University of Cantabria during the Prestige accident, is calibrated by means of an automatic calibration methodology. The shuffled complex evolution method, developed by the University of Arizona (SCE-UA), is applied to estimate the optimal coefficients of the model. The calibration of the model has been carried out using 13 buoys deployed in the Bay of Biscay during the Prestige accident as well as coetaneous meteorological and oceanographic data. Moreover, reanalysis data collected in the Spanish ESEOO project framework has also been used. Results suggest that buoys outside the continental slope were mainly driven by wind, whereas ocean currents played an important role in the motion of the buoys located over the continental slope and shelf. According to these findings, the final calibration of the coefficients is performed considering different buoy data. The methodology applied to this broad buoy database, has allowed us to calibrate the model, taking into account the relative importance of the forcings in buoy movement as well as the dynamics associated with each area.