Modelling the Mediterranean Sea ecosystem at high spatial resolution to inform the ecosystem-based management in the region.

Cumulative pressures are rapidly expanding in the Mediterranean Sea with consequences for marine biodiversity and marine resources, and the services they provide. Policy makers urge for a marine ecosystem assessment of the region in space and time. This study evaluates how the whole Mediterranean food web may have responded to historical changes in the climate, environment and fisheries, through the use of an ecosystem modelling over a long time span (decades) at high spatial resolution (8 × 8 km), to inform regional and sub-regional management. Results indicate coastal and shelf areas to be the sites with highest marine biodiversity and marine resources biomass, which decrease towards the south-eastern regions. High levels of total catches and discards are predicted to be concentrated in the Western sub-basin and the Adriatic Sea. Mean spatial-temporal changes of total and commercial biomass show increases in offshore waters of the region, while biodiversity indicators show marginal changes. Total catches and discards increase greatly in offshore waters of the Western and Eastern sub-basins. Spatial patterns and temporal mean changes of marine biodiversity, community biomasses and trophic indices, assessed in this study, aim at identifying areas and food web components that show signs of deterioration with the overall goal of assisting policy makers in designing and implementing spatial management actions for the region.

Resilience dynamics and productivity-driven shifts in the marine communities of the Western Mediterranean Sea.

Ecological resilience has become a conceptual cornerstone bridging ecological processes to conservation needs. Global change is increasingly associated with local changes in environmental conditions that can cause abrupt ecosystem reorganizations attending to system-specific resilience fluctuations with time (i.e. resilience dynamics). Here we assess resilience dynamics associated with climate-driven ecosystems transitions, expressed as changes in the relevant contribution of species with different life-history strategies, in two benthopelagic systems. We analysed data from 1994 to 2019 coming from a scientific bottom trawl survey in two environmentally contrasting ecosystems in the Western Mediterranean Sea-Northern Spain and Alboran Sea. Benthopelagic species were categorized according to their life-history strategies (opportunistic, periodic and equilibrium), ecosystem functions and habitats. We implemented an Integrated Resilience Assessment (IRA) to elucidate the response mechanism of the studied ecosystems to several candidate environmental stressors and quantify the ecosystems' resilience. We demonstrate that both ecosystems responded discontinuously to changes in chlorophyll-a concentration more than any other stressor. The response in Northern Spain indicated a more overarching regime shift than in the Alboran Sea. Opportunistic fish were unfavoured in both ecosystems in the recent periods, while invertebrate species of short life cycle were generally favoured, particularly benthic species in the Alboran Sea. The study illustrates that the resilience dynamics of the two ecosystems were mostly associated with fluctuating productivity, but subtle and long-term effects from sea warming and fishing reduction were also discernible. Such dynamics are typical of systems with wide environmental gradient such as the Northern Spain, as well as systems with highly hydrodynamic and of biogeographical complexity such as the Alboran Sea. We stress that management should become more adaptive by utilizing the knowledge on the systems' productivity thresholds and underlying shifts to help anticipate both short-term/less predictable events and long-term/expected effects of climate change.

Non-Redfieldian dynamics driven by phytoplankton phosphate frugality explain nutrient and chlorophyll patterns in model simulations for the Mediterranean Sea.

The relative abundance of nitrate (N) over phosphate (P) measured as a molar ratio (N:P) is typically considered to determine the macronutrient limiting marine primary production. In low-complexity biogeochemical models, a simple threshold value is usually applied based on the canonical Redfield ratio (N:P = 16). However, the N:P ratio is not constant in many oceanic areas, especially marginal, semi-enclosed seas, such as the Mediterranean basin. In this work, a flexible definition of the N:P ratio based on the capacity of phytoplankton to modulate phosphate uptake according to its availability in seawater, the so-called Line of Frugality, is incorporated into the biogeochemical model MedERGOM. This modification allows the acquisition of a more realistic representation of the stoichiometry of nutrients in the Mediterranean basin and allows to better reproduce the observed phytoplankton biomass in productive areas such as the Gulf of Gabes and the Adriatic Sea. This approach is, thus, especially suitable for coastal areas in which basin-scale biogeochemical models fail to reproduce patterns observed by remote sensing or in situ measurements. Our results show that implementation of the stoichiometric flexibility of phytoplankton in a low-complexity biogeochemical model enhances the reproducibility of ecosystem dynamics without increasing the computational demand, representing a simple approximation easily implemented in models aiming to describe regions with a Non-Redfieldian stoichiometry.

Deep winter convection and phytoplankton dynamics in the NW Mediterranean Sea under present climate and future (horizon 2030) scenarios.

Deep water convection (DC) in winter is one of the major processes driving open-ocean primary productivity in the Northwestern Mediterranean Sea. DC is highly variable in time, depending on the specific conditions (stratification, circulation and ocean-atmosphere interactions) of each specific winter. This variability also drives the interannual oscillations of open-ocean primary productivity in this important region for many commercially-important fish species. We use a coupled model system to 1) understand to what extent DC impacts phytoplankton seasonality in the present-day and 2) to explore potential changes in future scenarios (~2030). Our model represents quite accurately the present-day characteristics of DC and its importance for open-ocean phytoplankton blooms. However, for the future scenarios the importance of deep nutrients in fertilizing the euphotic layer of the NW Mediterranean decreases. The model simulates changes in surface density and on the levels of kinetic energy that make mesoscale activity associated with horizontal currents to become a more important fertilization mechanism, inducing subsequently phenological changes in seasonal plankton cycles. Because of our focus on the open-sea, an exact quantification of the impact of those changes on the overall biological production of the NW Mediterranean cannot be made at the moment.

Hydrological and biogeochemical response of the Mediterranean Sea to freshwater flow changes for the end of the 21st century.

We evaluate the changes on the hydrological (temperature and salinity) and biogeochemical (phytoplankton biomass) characteristics of the Mediterranean Sea induced by freshwater flow modifications under two different scenarios for the end of the 21st century. An ensemble of four regional climate model realizations using different global circulation models at the boundary and different emission scenarios are used to force a single ocean model for the Mediterranean Sea. Freshwater flow is modified according to the simulated changes in the precipitation rates for the different rivers' catchment regions. To isolate the effect resulting from a change in freshwater flow, model results are evaluated against a 'baseline' simulation realized assuming a constant inflow equivalent to climatologic values. Our model results indicate that sea surface salinity could be significantly altered by freshwater flow modification in specific regions and that the affected area and the sign of the anomaly are highly dependent on the used climate model and emission scenario. Sea surface temperature and phytoplankton biomass, on the contrary, show no coherent spatial pattern but a rather widespread scattered response. We found in open-water regions a significant negative relationship between sea surface temperature anomalies and phytoplankton biomass anomalies. This indicates that freshwater flow modification could alter the vertical stability of the water column throughout the Mediterranean Sea, by changing the strength of vertical mixing and consequently upper water fertilization. In coastal regions, however, the correlation between sea temperature anomalies and phytoplankton biomass is positive, indicating a larger importance of the physiological control of growth rates by temperature.

On the causal structure between CO2 and global temperature.

We use a newly developed technique that is based on the information flow concept to investigate the causal structure between the global radiative forcing and the annual global mean surface temperature anomalies (GMTA) since 1850. Our study unambiguously shows one-way causality between the total Greenhouse Gases and GMTA. Specifically, it is confirmed that the former, especially CO2, are the main causal drivers of the recent warming. A significant but smaller information flow comes from aerosol direct and indirect forcing, and on short time periods, volcanic forcings. In contrast the causality contribution from natural forcings (solar irradiance and volcanic forcing) to the long term trend is not significant. The spatial explicit analysis reveals that the anthropogenic forcing fingerprint is significantly regionally varying in both hemispheres. On paleoclimate time scales, however, the cause-effect direction is reversed: temperature changes cause subsequent CO2/CH4 changes.

Application of the singular spectrum analysis technique to study the recent hiatus on the global surface temperature record.

Global surface temperature has been increasing since the beginning of the 20th century but with a highly variable warming rate, and the alternation of rapid warming periods with 'hiatus' decades is a constant throughout the series. The superimposition of a secular warming trend with natural multidecadal variability is the most accepted explanation for such a pattern. Since the start of the 21st century, the surface global mean temperature has not risen at the same rate as the top-of-atmosphere radiative energy input or greenhouse gas emissions, provoking scientific and social interest in determining the causes of this apparent discrepancy. Multidecadal natural variability is the most commonly proposed cause for the present hiatus period. Here, we analyze the HadCRUT4 surface temperature database with spectral techniques to separate a multidecadal oscillation (MDV) from a secular trend (ST). Both signals combined account for nearly 88% of the total variability of the temperature series showing the main acceleration/deceleration periods already described elsewhere. Three stalling periods with very little warming could be found within the series, from 1878 to 1907, from 1945 to 1969 and from 2001 to the end of the series, all of them coincided with a cooling phase of the MDV. Henceforth, MDV seems to be the main cause of the different hiatus periods shown by the global surface temperature records. However, and contrary to the two previous events, during the current hiatus period, the ST shows a strong fluctuation on the warming rate, with a large acceleration (0.0085°C year(-1) to 0.017°C year(-1)) during 1992-2001 and a sharp deceleration (0.017°C year(-1) to 0.003°C year(-1)) from 2002 onwards. This is the first time in the observational record that the ST shows such variability, so determining the causes and consequences of this change of behavior needs to be addressed by the scientific community.

Biogeochemical control of marine productivity in the Mediterranean Sea during the last 50 years.

The temporal dynamics of biogeochemical variables derived from a coupled 3-D model of the Mediterranean Sea are evaluated for the last 50 years (1960-2010) against independent data on fisheries catch per unit effort (CPUE) for the same time period. Concordant patterns are found in the time series of all of the biological variables (from the model and from fisheries statistics), with low values at the beginning of the series, a later increase, with maximum levels reached at the end of the 1990s, and a posterior stabilization. Spectral analysis of the annual biological time series reveals coincident low-frequency signals in all of them. The first, more energetic signal peaks around the year 2000, while the second, less energetic signal peaks near 1982. Almost identical low-frequency signals are found in the nutrient loads of the rivers and in the integrated nutrient levels in the surface marine ecosystem. Nitrate concentration shows a maximum level in 1998, with a later stabilization to present-day values, coincident with the first low-frequency signal found in the biological series. Phosphate shows maximum concentrations around 1982 and a posterior sharp decline, in concordance with the second low-frequency signal observed in the biological series. That result seems to indicate that the control of marine productivity (plankton to fish) in the Mediterranean is principally mediated through bottom-up processes that could be traced back to the characteristics of riverine discharges. The high sensitivity of CPUE time series to environmental conditions might be another indicator of the overexploitation of this marine ecosystem. KEY POINTS: Biogeochemical evolution of the Mediterranean over the past 50 yearsRiver nutrient loads drive primary and secondary productionsStrong link between low trophic levels and fisheries.

Understanding the causes of recent warming of mediterranean waters. How much could be attributed to climate change?

During the past two decades, Mediterranean waters have been warming at a rather high rate resulting in scientific and social concern. This warming trend is observed in satellite data, field data and model simulations, and affects both surface and deep waters throughout the Mediterranean basin. However, the warming rate is regionally different and seems to change with time, which has led to the question of what causes underlie the observed trends. Here, we analyze available satellite information on sea surface temperature (SST) from the last 25 years using spectral techniques and find that more than half of the warming tendency during this period is due to a non-linear, wave-like tendency. Using a state of the art hydrodynamic model, we perform a hindcast simulation and obtain the simulated SST evolution of the Mediterranean basin for the last 52 years. These SST results show a clear sinusoidal tendency that follows the Atlantic Multidecadal Oscillation (AMO) during the simulation period. Our results reveal that 58% of recent warming in Mediterranean waters could be attributed to this AMO-like oscillation, being anthropogenic-induced climate change only responsible for 42% of total trend. The observed acceleration of water warming during the 1990s therefore appears to be caused by a superimposition of anthropogenic-induced warming with the positive phase of the AMO, while the recent slowdown of this tendency is likely due to a shift in the AMO phase. It has been proposed that this change in the AMO phase will mask the effect of global warming in the forthcoming decades, and our results indicate that the same could also be applicable to the Mediterranean Sea. Henceforth, natural multidecadal temperature oscillations should be taken into account to avoid underestimation of the anthropogenic-induced warming of the Mediterranean basin in the future.