Effects of Climate Crisis on Marine Ecosystem: Mass Mortality Event of the Mediterranean Mussels (Mytilus galloprovincialis) in the Middle Adriatic

1. Introduction

The climate crisis is a global threat to species, biodiversity, and ecosystems that affects individual organisms and the way they interact with other organisms and their habitats, altering ecosystem structure and functions [1]. In this regard, the Intergovernmental Panel on Climate Change claims that the rapid global warming over the past few decades has had consequences for weather, climate, ecosystems, human society, and the economy [2]. The impacts of the climate crisis are widespread but not uniform over time and space, as the responses of species and habitats vary according to their relative vulnerability, degree of exposure, sensitivity, and capacity to adapt [3-6]. The climate crisis affects not only the average temperature but also the duration, magnitude, and frequency of extreme events such as droughts, forest fires, and heat waves [7]. This synergy implies strong pressures on ecosystems, reducing their resilience [8].

On a global scale, recent decades have seen several high-impact marine heatwaves (MHWs) (anomalously warm water events) [9], that can affect local environmental dynamics. Significant events were observed in the northern Mediterranean Sea in 2003 [10,11], in Australia (west coast) in 2011 [12], in the northwestern Atlantic Ocean in 2012 [13], in the north-eastern Pacific Ocean between 2013 and 2015 [14,15], off South East Australia in 2015/16 [16] and across Northern Australia in 2016 [17]. MHWs can generate significant impacts on ecological communities, such as the loss of kelp forests [18], coral bleaching [19], reduction of surface chlorophyll levels (Bond et al., 2015), mass mortality of marine invertebrates [20,21] and changes in distribution and structure of species and communities respectively [20,18,22].

Based on their characteristics, the Mediterranean regions, including the Adriatic ones, are considered hot spots of climate crisis [23]. In particular, the Adriatic ecosystems suffer the combined effect of climate change in the Northern Hemisphere and local climate variations [24]. The Adriatic Sea is a basin into which most of the rivers of northern and central Italy flow and which communicates with the Mediterranean Sea only through the Otranto channel (just 70 km wide); therefore, it could be considered an important example of a marine system strongly and dynamically influenced by both anthropic pressures led by rivers, and climate conditions [25-30]. A number of MHW events have been reported in the Adriatic Sea [31]. However, there is no report regarding the latest MHW events, nor even their potential association with mass mortality events (MMEs) of marine organisms.

The aim of this work was to analyze the potential causes and the ecological consequences of an MME of the Mediterranean mussel (Mytilus galloprovincialis) observed along the Piceno Adriatic coast (Marche Region, central Italy) at the end of the summer of 2022. Specifically, we integrated a climatological analysis with observational data in order to explore the potential relationship between MHW and mussel MME. It is worth mentioning that in the Mediterranean Sea, M. galloprovincialis is ecological and economic importance [32-34]. As far as we know, this report is the first observation of a mussel MME within this coastal area, in which there are two marine Sites of Community Importance (SCI) and three mussel farms.

2. Materials and Methods

2.1. Study area

The study area extends for about 16 km along the coast of the middle Adriatic, locally known as the Piceno coast (Figure 1). The climate has warm and moderately stormy summer, fresh winter with rare snowfall and mild and humid equinoctial season; the precipitation evidenced an absolute maximum in autumn and a slight secondary maximum in the month of May [35]. Summers are not usually too dry. The rainfall regime is bimodal - sublittoral adriatic and the annual rainfall and temperatures average (1997-2022), are 593 mm and 15.7°C respectively [36]. The wave regime is mostly characterized by storms from NNE to ESE. The tides are very low, with average amplitudes around 0.4-0.5 m with a maximum of around 0.75 m.

The investigated area is characterized by shallow sandy bottoms that slowly degrade towards the open sea (medium deep of about -4.5 mt at 300 mt from the coast). Continuity solutions of this type of bottom are represented by artificial reefs (linear barriers placed to defend the coast against marine erosion) and by a natural reef of about 7 hectares, which is included in the boundaries of the SCI IT534022 "Costa del Piceno - San Nicola a mare" [37]. These hard substrates are here usually and strongly colonized by mussel beds (i.e., Mytilus galloprovincialis), which perform various ecosystem functions such as the increase of the substrate three-dimensional complexity (ecosystem engineering function), the control of seston abundance in the water column by filtration (control particle abundance function), the transfer of organic matter and nutrients from the water column to the seabed through the production of biodeposits such as feces and pseudofeces (benthic-pelagic coupling and nutrient cycling function), the removal of waste from the water (bioremediation function) or being a food resource for marine organisms (feeding function) [34,38-40]. The observation sites insist on the following reef environments: Site 1, Grottammare, in the Marine Site of Community Importance (4761502N, 408104E) (natural reef); Site 2, San Benedetto del Tronto harbor (4756933N, 409848E) (artificial reef); Sites 3, San Benedetto del Tronto beach (4754238N, 409930E) (artificial reef); Site 4, Martinsicuro beach (4746627N, 412413E) (artificial reef) (Figure 1).

Moreover, a mussel farm faces the town of San Benedetto del Tronto, about three nautical miles from the coast; two others are located just south and north of the study area.

2.2. Observations

For all observation sites, the pre and post-summer 2022 status was detected by authors’ direct observations; for sites 1 and 3 the changes are also documented by photos. Moreover, about Site 1, there is a detailed description of the pre-summer 2022 condition [32,34] since it is included in the SCI-IT534022. Direct interviews with mussels’ farmers were conducted to quantify the loss of productivity.

2.3. Climatic data

In order to analyze the possible causes of the massive death of Mytilus galloprovincialis observed within the study area, we carried out an analysis of the available climatic data. In detail, on the basis of the ecological characteristics of this species, the climatic signal was characterized by focusing on air and sea temperature (time trend and heat waves) and precipitation/flow rivers.

Validated and analyzed meteo-marine data related to the measuring station of San Benedetto del Tronto (4756389N, 409228E) with a variable time step between 10 and 60 minutes according to the 2008 World Meteorological Organization guide [41] have been recovered from the dedicated website of marine weather observations of ISPRA [42]. This signal is homogeneous and almost continuous (except for very short recording gaps) and the data were analyzed from 1 January 2011. Evidently, this dataset is not a sufficiently extended historical series to supply statistically meaningful indications on the climatic course of the main marine weather features (12 complete years) but it falls anyway in a temporal phase in which the climatic crisis in existence is evidencing its greater effects in terms of thermodynamic variations of the Adriatic basin. Finally, we paid particular attention to both the climatic conditions of the last summers (2022-2019) and those of the summers of 2017 and 2015 in which the other two heat waves occurred.

3. Results

3.1. Climatic data analysis

The historical series relating to thermal parameters for water and air, for a total of more than 1.2 mln of numerical data (period 2011-2022), was analyzed in detail. A rather significant thermal increase of the sea surface temperature (SST), quantifiable at about 0.37°C and therefore in line with the indications deriving from the 17th Report on the climate in Italy of the “Higher Institute for Environmental Protection and Research” [43], highlights these data. The air temperature recorded in San Benedetto del Tronto during the same period evidenced an increase of about 1.08°C. In order to analyze summer temperature variations of the sea, Figure 2 shows the hourly details of July and August for the last few years.

In the 2022 year, the observed threshold of the SST of 29°C was continuously exceeded from 24th of July to the 10th of August, therefore for 17 days and the maximum temperature detected was 30.6°C on the 6th of August. However, comparing this intense wave of sea heat with the signals of recent years, it is evident that thermal values of the SST higher than 29°C are recurrent in the study area during the summer season. The highest absolute temperature values of the limited series were recorded in the summer of 2017 and precisely on the 5th of August, with 31.2°C; during this year, the SST above 29°C extended from day 1 to 14 of August, therefore less extensive than in 2022. In the summer of 2015, the maximum temperature of 31°C was detected on July 24th, and the period with SST>29°C extended from July 19th to August 15th, for 25 days (in this period the SST drops below 29°C for two days).

The status of the sea waters of the study area depends not only on the conditions at the mesoscale (basin of Tronto and Aso rivers) but especially on what happens in the upper Adriatic basin, where the main Italian rivers (in terms of water flow: Po, Adige, and Brenta) flow their freshwater; indeed, a dominant sea current from north to south, runs along the Italian Adriatic coast, transporting southward the runoff.

Analyzing the data of the above-cited ISPRA Report [43], it is possible to observe that the flow rates of the aforementioned rivers and especially the rainfall that occurred in the catchment areas of the Po and the Adige were extremely low even in the first half of the year, due to a long phase of climatic drought, extended temporally from December 2021 and lasted until July 2022; moreover, this phase was also characterized by a strong nivometric deficit on the Alpine and Apennine mountains. The negative anomalies related to the first 6 months of the 2021-22 meteorological years were about 55% in the regions of Piemonte and Valle d'Aosta and about 50% in Lombardy and Tridentine Venice (in northern Italy).

The Po River reached the absolute minimum of the last 200 years during 2022, with a flow of 104 mc/sec on the 24th of July (hydrometric station of Ferrara-Pontelagoscuro) (Figure 3); other exceptional dry periods for the Po River were in June 2022 with a range of 303 mc/sec, in June 2006 with 320 mc/sec, in June 2005 with 444 mc/sec and in June 2003 with 521 mc/sec [44].

In addition, during the summer of 2017, there was a slightly lower negative anomaly of winter and spring precipitation than in 2022 but the hydrological data of the flows of the Po in spring and summer were characterized by values much higher than those of 2022. Regarding the summer 2015, we observe that the first half of the year is characterized by rainfall around the average of the last thirty years, with moderate snowfall in the sectors of the central and eastern Alps and alpine river flows were significantly higher than in 2022 [44].

3.2. Field observation of mussel MME

In the present study, we recorded a massive absence of Mediterranean mussel beds on rocky substrates that were almost completely covered before the summer 2022 (Figure 4). This phenomenon was observed after the heat wave hit the central Adriatic coast in July-August 2022, when water temperatures reached 30.6°C and remained above 29°C for 17 days (24th July - 10th August) (Figure 2). The drop in mussel production recorded in the farms located in front of the study area three nautical miles from the coast, was around 30%.

4. Discussion

4.1. Climatic data

The air thermal increase of about 1.08°C recorded in San Benedetto del Tronto harbor station, during the 2011-2022 timespan, is confirmed by the data recorded in the nearby harbors of Ancona and Ortona, respectively 1.27°C and 1.22 °C. At meso-scale level we can observe a strong acceleration of climate forcing. The increase of 0.37°C of the SST, in the same period, is comparable to the data acquired from satellite observations within the Copernicus Marine Environment Monitoring Service (CMEMS) project - Mediterranean SST [45] and with those acquired directly in the Croatian surface waters of the Central Adriatic Sea [23]. These authors highlighted two fundamental points: i) contrary to what happened in the previous 20 years, the period after 2005 would be characterized by a constant increasing thermometric trend (about 0.04ºC/year) without cooling or breaking phases, ii) there is a change in the width of the seasonal cycle, characterized by a sharp increase in the uptrend of SST in summer (about 0.056°C/year), and a slowdown in winter (~0.029°C/ year). The increase in SST on the eastern coast of the Adriatic Sea, with a total increase exceeding 1°C between 1979 to 2015, was also shown by [46] and [47]. The seasonal variability changes could affect the SST trend during solstice seasons - with summers warmer than the winters - and could amplify the effects of the sea heat waves in terms of frequency, magnitude, and duration, with negative impacts on marine ecosystems as a consequence.

The 2022 sea heat wave resulted from the high air temperatures recorded at the end of June along the middle and upper Adriatic Italian coast. The month of July has been indeed the hottest of the last 50 years [48]. These sea heat waves confirm a trend of increasing temperatures in the Mediterranean basin in the last twenty years [31].

The reduced river runoff is another crucial aspect that characterized the Central Adriatic waters during 2022. The Po River is the largest Italian river and the main contributor to the nutrient load in the central northern Adriatic waters [26]. In summer, the typical flow of the Po River falls below the range of 550 and 700 mc/sec. Interestingly, the slowest flow (i.e. 104 mc/sec) was recorded exactly in July 2022 [44]. In 2022, the fluvial regime of the Po River (Figure 3) and, simultaneously, of the Adige, Piave and Brenta rivers, which are normally of the nivo-glacial type, showed a typical Mediterranean regime, with a slight maximum in winter and a slight recovery in the flow rate in May. This regime was due to the rapid melting of the very scarce residual snowpack present on the Alpine chain, and a strong summarizing minimum in summer.

4.2. Mussels MME

The disappearance of the mussels bed observed in the study area is not attributable to phenomena of predation or illegal fishing (Figure 5), as it is widespread throughout the investigated area and is totalitarian in the observation sites.

Our findings align with previous studies demonstrating that HMWs correlate with MMEs in the Mediterranean Sea [21,49]. However, there is still limited knowledge of this association in the Adriatic Sea area. In the Northwestern Mediterranean Sea, mortality incidence between 20% and 100% was registered according to MHW duration ranging from 25 to 50 days [50]. This finding suggests the existence of high variability in the level of mortality impact for similar MHWs and demonstrates the need to consider further factors in analyzing the relationship between heat exposure conditions and MMEs. Among others, these factors include the ecological memory of the recurrence of MHWs over the same geographic areas [51], and potential species or populations-related differences [52,53]. For instance, benthic organisms colonizing intertidal zones, such as M. galloprovincialis, are exposed to a variety of harsh marine and terrestrial conditions that could affect their survival. Indeed [54] demonstrated that the Mediterranean mussel exposed to a low food load and high temperatures (30ºC) had a survival of 33%, while individuals exposed to the same thermal conditions but with a high food load had a survival of 57%.

Previous experimental studies focusing on measuring the death levels of M. galloprovincialis in response to increasing temperature showed that individuals exposed to 30ºC reached a 40% mortality after 10 days of exposure, peaking at 100% on the 12th day [55]. Indeed, [56,57] identified a temperature of 24–25°C as an upper limit for optimal physiological processes, indicating that Mediterranean mussels inhabiting natural populations already live close to their thermal acclimation limits. Physiological parameters that may contribute to increased animal stress include the heart rate. [58] observed that increasing temperatures promote the rise of heartbeats until a certain critical temperature, after which the heart rate drops. Therefore, the effects of the climate crisis are challenging their physiological limits, making these bivalves particularly vulnerable to extreme thermal variations, especially considering that rocky intertidal coasts are one of the most thermally variable and stressful habitats [59,60].

Another aspect to consider for both natural and farmed mussel populations is the status of the byssus. Byssal threads are composed of extracellular proteinaceous (collagenous) fibers and are synthesized along the mussel's foot starting from a mussel foot protein [61]. The byssus' role includes preventing dislodgement of the bivalve in harsh hydrodynamic conditions, resisting mobile competitors, and competing for space [62]. Environmental conditions like the warming of waters and ocean acidification may alter byssal production processes causing a reduction in strength, stiffness, and extension [62-64]. [65] further demonstrated that in Mytilus coruscus the byssus may be highly susceptible to the adverse effects of ocean warming because the expression of mussel foot protein genes was affected by elevated temperature, leading to a risk of dislodgement at high sea temperature (up to 27°C). Similarly, [66] found that in M. galloprovincialis, the attachment strength was negatively correlated with sea surface temperature.

One further consequence of the increased ocean warming is decreased marine productivity [67]. [68] have suggested that depending on the timing of phytoplankton blooms, marine bivalves may face periods of limited food availability and elevated temperatures, especially in the high intertidal zone. These authors indeed recorded a decrease in the enzyme kinetics with a controlled increase of water temperature in the Mediterranean mussels in anoxia or normoxia conditions. Interestingly, the authors found that food availability can buffer the detrimental effects of harsh environmental conditions.

Our observations also show that during the summer of 2022, the study area experienced a dry period (i.e. strong nivometric deficit in winter, low rainfall in springer, and an absolute minimum of rivers’ flow in springer and summer), which resulted in a low concentration of nutrients released into the sea by rivers. In fact, local mussel farmers estimated a reduction of 30% in their productivity in this period. We thus hypothesize that the mussel MMEs observed in our study are due to a vicious cycle of heat stressors (i.e. MHW) which increased the metabolic rate, and to the low nutrient availability, which did not meet the metabolic demand. Indeed, heatwaves with temperatures reaching 31ºC also occurred in the summer seasons of 2015 and 2017. However, the Po River runoffs during the spring and summer of these years were characterized by values much higher than those of 2022 [44], and as consequence, neither mussel MMEs nor mussel bed reduction were observed in the same study area, during our monthly surveys.

5. Conclusions

Based on our field observation in the study area, climatic data analysis, and also considering previous findings on M. galloprovincialis’ tolerance under certain stress conditions, we suggest that the MME observed in summer 2022 was likely due to prolonged anomalously warm water combined with reduced nutrient supply. These findings clearly indicate that the Adriatic Sea is undergoing an evident acceleration of climate change impact, whose socio-ecological implications need to be addressed.

Probably, M. galloprovincialis from our study area lives in a defined range of summer temperatures and food availability; the knowledge of the climatic conditions that led to the observed massive die-off can represent an important contribution both to the monitoring of ecosystems and the management of mussel farms. The consequences of this stressor combination (i.e., MHW and low nutrient load) on mussel farms are extremely different to what is possible to observe in natural settlements because the good farming practices adopted by farmers (e.g., anticipating the sale of the product before the hottest periods and using fine knit socks (even in cotton), can reduce the damage in certain periods allowing production to resume. On the contrary, climate crisis-driven alterations of natural habitats can be dramatic and require a longer time to restore, with the risk that once altered, natural habitat can reassemble itself with components and structures different from the initial condition.

Considering that these extreme climatological events will probably intensify [69], our preliminary analysis stresses the need for adopting specific mitigation and restoration measures in order to improve the management of marine natural ecosystems and fisheries.