Heatwave impacts on intertidal seagrass reflectance: from laboratory experiment to satellite mapping of seagrass heat shock index

Authors
Affiliations

Simon Oiry10

Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France

Bede Ffinian Rowe Davies¹

Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France

Philippe Rosa¹

Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France

Augustin Debly¹

Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France

Maria Laura Zoffoli11

Consiglio Nazionale delle Ricerche, Istituto di Scienze Marine (CNR-ISMAR), 00133 Rome, Italy

Anne-Laure Barillé12

Bio-littoral, Immeuble Le Nevada, 2 Rue du Château de l’Eraudière, 44300 Nantes, France

Nicolas Harin³

Bio-littoral, Immeuble Le Nevada, 2 Rue du Château de l’Eraudière, 44300 Nantes, France

Marta Román¹⁴

Centro de Investigación Mariña, Departamento de Ecoloxía e Bioloxía Animal, Facultade de Ciencias do Mar, Universidade de Vigo, Campus Universitario As Lagoas-Marcosende, 36310 Vigo, Spain

Pierre Gernez¹

Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France

Laurent Barillé¹

Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France

Published

February 26, 2025

Other Formats
Abstract
Seagrasses play a vital role in coastal ecosystems, providing habitat, stabilizing sediments, and contributing to carbon sequestration. However, climate change has increased the frequency and intensity of heatwaves, posing a significant threat to seagrass health. This study investigates the effects of marine and atmospheric heatwaves on the spectral reflectance of the intertidal seagrass Zostera noltei. Laboratory experiments were conducted under controlled heatwave conditions, where hyperspectral reflectance measurements were taken to assess the impacts over time. Heatwaves caused a substantial decline in seagrass reflectance, particularly in the green and near-infrared regions, corresponding to the browning of green leaves. Key vegetation indices, including the Normalized Difference Vegetation Index (NDVI) and Green Leaf Index (GLI), showed pronounced reductions under heatwave stress, with NDVI values decreasing by up to 34% and GLI by 57%. A novel metric, the Seagrass Heat Shock Index (SHSI), was developed to quantify the transition of seagrass leaves from green to brown, demonstrating a strong ability to capture the effects of heatwave exposure on seagrass coloration. Multispectral satellite observations corroborated the laboratory results, revealing widespread browning of seagrass leaves during marine and atmospheric heatwave events in South Brittany, France. Notably, darkened seagrass patches were observed in intertidal areas exposed to temperatures exceeding 32°C for over 13.5 hours per day. These findings highlight the potential of spectral reflectance as a tool for detecting early signs of heatwave-induced stress in seagrasses, offering a valuable method for remote sensing-based habitat assessment. The present study underscores the potential of remote sensing to capture rapid environmental changes in intertidal zones, enabling for continuous monitoring of seagrass meadows under the current and future climate regimes.
Keywords

Remote Sensing, Pigment Composition, Seagrass, Coastal Ecosystems, Heatwaves

1 Abstract Long version (For LPS Submission)

Seagrasses meadows are important coastal ecosystems, serving as crucial habitats for marine biodiversity, stabilizing sediments to mitigate erosion, and acting as significant carbon sinks in global climate regulation. However, these vital ecosystems face mounting threats from climate change, with the intensification and increased frequency of marine and atmospheric heatwaves posing profound risks to their health and functionality. This study investigates the impact of these extreme thermal events on the intertidal seagrass, employing a combination of laboratory-controlled experiments and satellite-based remote sensing to capture changes in spectral reflectance and to assess its ecological implications.

During laboratory experiments, seagrass of the species Zostera noltei were exposed to controlled simulated heatwave conditions during three days to assess the physiological and structural impacts of extreme thermal stress on their spectral reflectance. The experimental design involved placing seagrass samples in intertidal chambers that simulated natural tidal cycles, allowing to precisely regulate temperature conditions during both high and low tides. One chamber served as a control, maintaining typical seasonal temperatures, while the other was used to simulate heatwave conditions by progressively increasing air and water temperatures to mimic an actual heatwave event observed in the field. Hyperspectral reflectance measurements were recorded at regular intervals during low-tide to monitor changes in the seagrass over time. After two days of heatwave exposure the reflectance decreased significantly, particularly in the green (around 560 nm) and near-infrared (NIR) regions of the spectrum. This decline in reflectance was closely linked to visible leaf browning, suggesting alterations in pigment content, in the internal structure of seagrass leaves and their overall physiological performance. The green reflectance decline indicated a reduction in the plant’s photosynthetic capacity, while changes in NIR reflectance were related to the heat-indiced damage to internal arrangement of cells and air spaces. Vegetation indices like the Normalized Difference Vegetation Index (NDVI) and the Green Leaf Index (GLI), which are indicative of the overall health and structural integrity of vegetation, showed marked decreases under heatwave conditions, with NDVI dropping by up to 34% and GLI by 57%. This reduction highlights the adverse effects on the seagrass’s ability to maintain its normal physiological processes under heat stress. To quantify these changes more effectively, a novel Seagrass Heat Shock Index (SHSI) was developed. The SHSI, applicable on emerged seagrass, was particularly effective in detecting the transition from green leaves to darkened, stressed leaves, providing a sensitive and reliable tool for assessing thermal stress effects on seagrass. By focusing on specific changes in reflectance across certain spectral bands, the SHSI allowed for a clear differentiation between unimpacted and impacted vegetation, capturing the onset of heat-induced stress with high accuracy. This sensitivity makes the SHSI valuable for early intervention, enabling managers and researchers to identify vulnerable seagrass meadows before substantial damage occurs, thereby facilitating more timely conservation measures.

Complementing the experimental data, Sentinel-2 observations provided clear evidence of the effects of a documented heatwave event in South Brittany, France, on natural seagrass meadows in September 2021. These intertidal zones were exposed to extreme air temperatures up to 32°C for more than 13.5 hours per day, leading to significant leaf darkening, which affected up to 24% of the meadow’s area. The satellite-derived SHSI indicated a strong spatial correlation between prolonged heat exposure and areas experiencing spectral darkening, which was consistent with the laboratory experiment. Spatial analysis showed that darkening was especially pronounced in the higher intertidal regions, where seagrasses were exposed to air for longer durations during low tide, underlining the interaction between tidal exposure and thermal stress. Using a Sentinel-2 image acquired one month after the heatwave, SHSI showed partial recovery of darkened patches. Despite this recovery, the seagrasses that had experienced the most severe darkening did not fully return to their original state, indicating that while Z. noltei meadows have some capacity for resilience, prolonged or repeated thermal stress can have lasting impacts, particularly in the more exposed intertidal regions. This highlights the need for focused conservation efforts to support their recovery and enhance their resilience in the face of increasing climate-driven thermal extremes.

The combination of laboratory-controlled experiments and satellite-based remote sensing provided a comprehensive understanding of the impacts of heatwaves on seagrass meadows at multiple scales, from individual leaf-level responses to entire meadow-wide effects, offering valuable insights for the management and conservation of these vulnerable habitats. The study underscores the critical role of spectral reflectance as an early warning indicator of heatwave-induced stress, laying the foundation for remote sensing-based monitoring and conservation efforts. By employing innovative indices like the Seagrass Heat Shock Index (SHSI), we captured ecologically significant changes in the physiological performance of Z. noltei, advancing the precision of habitat assessments in intertidal zones. As climate scenarios predict more frequent and intense heatwaves, the need for continuous monitoring of intertidal seagrass meadows becomes increasingly urgent. This research demonstrates the efficacy of remote sensing in capturing rapid environmental changes, providing a framework for mitigating the impacts of climate-driven stressors and calling for adaptive conservation strategies. These strategies should integrate advancements in remote sensing technologies with targeted field-based interventions to preserve the resilience of intertidal seagrass meadows, thereby addressing the escalating challenges posed by climate change and ensuring the continued health of these indispensable coastal habitats.

2 Introduction

Seagrasses play a crucial role in coastal ecosystems by providing habitats and feeding grounds for various marine species, supporting marine biodiversity, and contributing to primary production and carbon sequestration (Sousa et al., 2019; Unsworth et al., 2022). Seagrasses are essential for several ecological functions, such as sediment stabilization (Infantes et al., 2022) or eutrophication mitigation by consuming nutrients (Gladstone-Gallagher et al., 2018). This justifies their use as indicators of environmental changes due to their sensitivity to water quality variations (Zoffoli et al., 2021). The interactions between seagrass meadows and their associated herbivores further enhance the delivery of ecosystem services, including coastal protection, fisheries support and provision of habitat and resources for birds (Gardner and Finlayson, 2018; Jankowska et al., 2019; Unsworth and Butterworth, 2021; Zoffoli et al., 2023). Understanding and preserving seagrass is vital for maintaining the biodiversity and productivity of coastal regions (Ramesh and Mohanraju, 2020; Scott et al., 2018).

Despite their crucial role in marine ecosystems, seagrasses face numerous threats that compromise their health and functionality. Intertidal seagrasses are subjected to a combination of aquatic and aerial conditions linked with tidal cycles, and face disturbance from terrestrial and aquatic stressors. Coastal development and human activities are primary threats, reducing the available habitat for seagrasses and increasing water turbidity, limiting light penetration, and photosynthesis (Waycott et al., 2009). Seagrasses are also threatened by runoff from agricultural fields and urban areas leading to nutrient enrichment. Eutrophication promotes the growth of seaweed in coastal waters, causing macroalgal blooms that compete with seagrasses for light and nutrients (Brun et al., 2003; Oiry et al., 2024; Thomsen et al., 2023). Pollution from industrial and agricultural sources introduces harmful chemicals and heavy metals into coastal waters, posing toxic risks to seagrass health (Bastos et al., 2023; Green et al., 2021; Zahoor and Mushtaq, 2023). Among manifold anthropogenic stressors, heatwaves (HWs), exacerbated by climate change, pose a severe threat to seagrasses, with catastrophic dieback events observed woldwide (Carlson et al., 2018; Marbà and Duarte, 2010; Moore and Jarvis, 2008; Strydom et al., 2020; Thomson et al., 2015).

Marine Heatwaves (MHWs) are defined by Hobday et al. (2016) as prolonged discrete anomalously warm water events, while Atmospheric Heatwaves (AHW) are defined by Perkins and Alexander (2013) as periods of at least three consecutive days with temperatures exceeding the 90th percentile of a time series covering at least 30 years. Subtidal seagrass meadows are exposed to MHWs, whereas at the interface between land and ocean, intertidal seagrasses are exposed to both MHWs and AHWs. HWs profoundly impact seagrass physiology, with effects varying between species and geographic location. Widespread seagrass species such as Zostera marina exhibits high susceptibility to elevated sea surface temperatures during winter and spring, leading to advanced flowering, high mortality rates, and reduced biomass (Sawall et al., 2021). Similarly, Cymodocea nodosa shows increased photosynthetic activity during HWs but suffers negative effects on photosynthetic performance and leaf biomass during recovery (Deguette et al., 2022). Additionally, different populations of Zostera marina along the European thermal gradient exhibit varied photophysiological responses during the recovery phase of HWs, indicating differential adaptation capabilities among populations (Winters et al., 2011). High-latitude populations exhibited prolonged declines in photophysiological performance even after temperatures returned to control levels, whereas the low-latitude Adriatic population showed full recovery (Winters et al., 2011). These events intensify other stressors, such as overgrazing and seed burial, compromising recruitment (Guerrero-Meseguer et al., 2020). Although extensive research exists on marine heatwaves’ effects on subtidal seagrasses (Arias-Ortiz et al., 2018; Deguette et al., 2022; Strydom et al., 2020), less attention has been given to intertidal habitats and even less to the effect of atmospheric extreme events on intertidal seagrass. Nonetheless, recent research showed that the low tide exposure of Zostera noltei to a simulated four-day atmospheric HW caused significant decreases in its photosynthetic efficiency, resulting in leaf necrosis and decay (Román et al., 2023).

The increased occurrence of extreme climate events calls for the implementation of monitoring strategies able to provide detailed and spatially explicit assessments of HWs effects on seagrass meadows. In such context, remote sensing, whose ability to map seagrass distribution over a variety of spatio-temporal scales has been demonstrated (Davies et al., 2024a, 2024b; Oiry et al., 2024; Román et al., 2021), proved useful to study the changes in seagrass coverage caused by extreme HW event (Strydom et al., 2020). The pigment composition of plants, such as chlorophylls, carotenoids, and anthocyanins, significantly influences their spectral signature in the visible range due to their specific light absorption properties (Davies et al., 2023; Douay et al., 2022; Olmedo-Masat et al., 2020; Ustin and Jacquemoud, 2020). During the senescence phase of seagrass’ life-cycle, the degradation of chlorophyll and the unmasking of accessory pigments result in noticeable changes in leaf coloration and reflectance, including increased reflectance in the red and green wavelengths and shifts in the red-edge position (Boyer et al., 1988; Mariën et al., 2019; Peñuelas et al., 2004). Leaf browning, often observed after stress events, produces reflectance changes similar to those caused by senescence, enabling the detection of vegetation stress through remote sensing (Boyer et al., 1988; Peñuelas et al., 2004). Spectral indices such as the Brown Pigment Index (BPI) and the Photochemical Reflectance Index (PRI) have been developed to assess changes in terrestrial plant physiological status, including oxidative and drought stress (Garbulsky et al., 2011; Skendzic, 2023). While these effects are well-documented in terrestrial plants, the spectral reflectance changes associated with senescence and stress events such as MHWs or AHWs remain poorly studied on intertidal seagrasses.

This study aims to experimentally test the hypothesis that HWs alter the reflectance of the intertidal seagrass Zostera noltei. Controlled experiments in intertidal chambers were conducted to evaluate the direct impact of heat stress on seagrass reflectance. The findings will then be applied to satellite remote sensing images, providing critical insights into the spatial extent and temporal dynamics of HW effects on seagrass meadows. By linking experimental results with large-scale observations of seagrass leaves’ browning, the study aims to underscore the potential of remote sensing to enhance our understanding of seagrass responses to extreme thermal events across diverse settings and timescales.

3 Material & Methods

3.1 Laboratory Experiment

3.1.1 Sampling and acclimation of seagrasses

Seagrass samples were taken in summer 2024, at low tide, from a Zostera noltei (dwarf eelgrass) meadow located in Bourgneuf Bay, France (46°57’32.0”N, 2°10’37.0”W). A metal coring devices was used to sample seagrass from an area of 30x15x5 cm (length x width x depth, respectively), maintaining the sediment structure and avoiding damage to seagrass rhizomes and leaves (Figure 1 A). These coring device enabled the collection of sediment samples at a consistent depth, minimizing variability between samples. A total of six samples were collected. Samples including seagrass, sediment, meiofauna, and macrofauna, were placed in plastic trays. Keeping the entire biota allowed for natural interactions between components and reduced stress on the seagrass. Seawater was added to each tray to avoid hydric stress caused by insufficient moisture during transportation (1h drive from the laboratory). Simultaneously, seawater was sampled 4km away from the seagrass sampling site and transported to the lab, where it was filtered using a 0.22 µm nitrocellulose filter to remove suspended particulate matter. The filtered seawater was used in the acclimation tank and the intertidal chambers. The seagrasses were acclimated for one week with a water temperature of 17°C, matching the in situ temperature during sampling, and a photosynthetically active radiation (PAR) of 150 µmol.s-1.m-2.

Figure 1: Illustration of the experiment. A: Seagrass field sampling using a coring device; B: Intertidal chamber used during the experiment; C: Seagrass sample inside a chamber during the experiment at high tide; D: Treatment sample at the start of the experiment; E: Treatment sample at the end of the experiment, 3 days after the start of the HW event.

3.1.2 Experimental design

A tidal cycle (i.e. regularly alternating 6h of low-tide and 6h of high-tide) was simulated in the laboratory using an intertidal chamber system from ElectricBlue® (Figure 1 B ; Electric Blue (2023)). The experimental setup allowed only two tidal states: high tide or low tide, with no intermediate stages. The transition between these states took about 15 minutes to complete after initiation. During the phase of high tide, a volume of 30 L of filtered seawater was pumped and circulated through the chamber (Figure 1 B, C). During low tide, the seagrass sample was emerged. The acclimated seagrasses were split into two subsets and placed in two independent chambers used in parallel, with one chamber used for control and the other for the experimental treatment. The intertidal chambers were equipped with LED lights that emitted a low mount of red and infrared radiation. To achieve a Photosynthetically Active Radiation (PAR) intensity of up to 400 μmol·m⁻²·s⁻¹, a filament bulb was added inside the chambers. During the diurnal phase of the experiment, the PAR was kept constant in both intertidal chambers. To follow the circadian cycle, these lights (both LED and filament bulb) were turned on and off each day, at the time of sunrise and sunset, respectively.

Air temperature and water temperature were controlled inside the experiment chambers in order to reproduce the range of variability observed in the field (Figure 2). Field temperature was measured using in situ sensors (T7.3 EnvLoggers from ElectricBlue®) deployed at the sampling site in August 2024. The loggers were positioned along a transect from the upper to the lower intertidal zone, attached to pre-existing wooden poles at the sediment surface. In complement, the daily temperature maxima recorded in situ were compared with measurements from the nearest Météo France weather station (Annexe A1, Section 7.1). The control chamber was kept at temperatures representing typical summer conditions, with water temperatures at 18°C and air temperatures from 19°C to 23°C, following natural daily temperature fluctuations (Figure 2). In the treatment chamber, the air temperature was adjusted to mimic an AHW that affected the seagrass meadow in Quiberon, South Brittany, France (47°35’40.0”N, 3°07’30.0”W), from September 2 to September 6, 2021. Air temperature in the experimental chamber was set to vary from 23°C (at night) to 35°C (daytime) during the first day of the experiment, and increase by 1°C daily during three consecutive days. Water temperature in the experimental chamber was also adjusted to mimic MHW conditions, starting at the seasonal baseline (18°C) and rising incrementally by 0.5°C daily to simulate the increasing temperatures during the event. This aimed to reproduce the thermal stress experienced by the seagrass meadow during a MHW (Figure 2). The experiment, with both treatment and control chambers, was repeated three times to obtain replicates (hereafter referred to as “Run”).

Figure 2: Temperature variation in the control (left) and treatment (right) intertidal chambers, during the HW experiment. The red line indicates air temperature, and the blue line water temperature. Due to the tidal cycle of immersion / emersion, the seagrasses experienced the temperatures represented by solid lines.

3.1.3 Optical measurements

3.1.3.1 Hyperspectral reflectance measurements

Throughout the experiment, the hyperspectral reflectance, \(R(\lambda)\), of both the control and treatment seagrasses was measured using an ASD HandHeld 2 equipped with a fiber optic extension placed inside the chamber. The measurement set up made it possible to automatically acquire \(R(\lambda)\) without opening the chamber. An average of five \(R(\lambda)\) spectra, each with an integration time of 544 ms, was taken every minute during daytime (Malvern Panalytical, 2023). Every 10 minutes, the fiber optic was switched from one intertidal chamber to the other, in order to measure \(R(\lambda)\) in both the treatment and control. Light conditions were controlled inside of the chambers and the reflectance calibration was performed each morning at the very first moment of low tide using a Spectralon white reference with 99% Lambertian reflectivity.

3.1.3.2 Spectrum post-processing

A Savitzky-Golay smoothing function with a 5 nm moving window was applied to each spectrum using the “hsdar” package in R (Lehnert et al., 2017). The second derivative at 665 nm (δδ665), showing the highest variability between the control and the treatment, was tested as an indicator of the spectral changes following HWs.

The effect the HW on \(R(\lambda)\) was also quantified using two radiometric indices:

  • The Normalized Difference Vegetation Index (NDVI, Rouse et al. (1974)), a proxy of chlorophyll-a concentration (Equation 1)

\[ NDVI = \frac{R(840)-R(668)}{R(840)+R(668)} \tag{1}\]

where \(R(840)\) and \(R(668)\) are the reflectance at 840 and 668 nm respectively.

  • The Green Leaf Index (GLI, Louhaichi et al. (2001)), a quantification of the seagrass leaves greenness (Equation 2)

\[ GLI = \frac{[R(550)-R(668)]+[R(550)-R(450)]}{(2 \times R(550) )+ R(668) + R(450) } \tag{2}\]

where \(R(550)\) and \(R(450)\) are the reflectance in the green (at 550 nm) and in the blue (at 450 nm) spectral bands, respectively.

Based on the observed spectral changes in seagrasses exposed to HWs, we developed a new radiometric index to better detect the radiometric change caused by the HW. The browning of the leaves was characterized by substantial radiometric changes in both the green and red-edge spectral regions. The Seagrass Heat Shock Index (SHSI) was introduced as the reflectance line height at 740 nm, compared to the 560 - 842 nm baseline (Figure 3). Namely the SHSI subtracts the reflectance observed at 740 nm to the linearly interpolated reflectance between 560 and 842 nm, so that the index is positive in the case of brown, HW-impacted seagrass leaves, and negative in the case of green, non-impacted leaves:

\[ \text{SHSI} = I_{SHSI} - R(740) \tag{3}\]

where :

\[ I_{SHSI} = R(560) + \tau [R(842) - R(560)] \] and :

\[ \tau = \frac{740 - 560}{842 - 560} \]

where \(R(560)\), \(R(740)\), and \(R(842)\) represent the reflectance at 560, 740, and 842 nm, respectively, and \(\tau\) is equal to a contant of 0.64. These wavelengths were selected to align with the spectral resolution of satellites missions such as Sentinel-2, for broader remote sensing application.

Figure 3: Computation of the reflectance Seagrass Heat Shock Index (SHSI) for Impacted (A) and Unimpacted (B) seagrass leaves. The dashed line represents the reflectance interpolation between 560 and 842 nm. The red vertical at 740 nm represents the SHSI line height.

3.2 Observation of a seagrass bed impacted by a HWs

Field measurements were taken the 10th of September 2021 after an atmospheric and marine HW in order to assess the impact of heat stress on seagrass. The study site was a seagrass meadow near Quiberon (France : 46°57’32.0”N, 2°10’37.0”W, Figure 4). Brown seagrass leaves were observed over large patches of the meadow alongside areas covered by green seagrass (Figure 5). A total of 96 Quadrat Points (QPs) were collected as georeferenced quadrat images across the meadow. These images allowed for visual assessment of vegetation type, density, and coloration. The quadrats were then divided into two categories: green seagrasses (henceforth: unimpacted QPs) and brown seagrasses (henceforth: impacted QPs), based on a visual estimation of the leaf coloration (Figure 4).

Figure 4: Location of field observations in a seagrass meadow impacted by a HW that occurred on the 10th of September 2021 in Quiberon, South Brittany, France. The red line indicates the intertidal zone (Zone between high tide and low tide, exposed during low tide), the dark green area indicates the extent of the seagrass meadow and the olive area indicate saltmarshes. Green points indicate the location of quadrat pictures over unimpacted seagrasses (i.e. showing a green colour on the field), and orange points indicate the location of quadrats taken over impacted seagrasses (i.e. showing a brown color on the field).
Figure 5: Illustration of the two colorations of seagrass leaves observed in situ on the 10th of September 2021 after a heatwave in Quiberon, South Brittany (France). A: Picture of a zone with both green and brown seagrass; B: Seagrass quadrat with green leaves; C: Seagrass quadrat with brown leaves; D: Picture of a zone where all leaves turned brown.

3.2.1 Temperature data and HW detection

3.2.1.1 Air temperature

Hourly air temperature data from 1952 to 2024 (more than 395,000 observation) from a nearby weather station (Lorient-Lann Bihoue, 47°45’46”N 3°26’11”W) was retrieved from Meteo France (https://portail-api.meteofrance.fr).

3.2.1.2 Water temperature

Sea Surface Temperature (SST) data from 1982 -2022 over the Quiberon coastal area was downloaded from the Copernicus Marine Data Store (Copernicus Marine Environment Monitoring Service, CMEMS (2024)). An area of 2700 km² was extracted and analyzed. This area was large enough to minimize missing values caused by cloud cover and small enough to limit the influence of offshore SST stability.

3.2.1.3 Heatwave detection and characterization

MHW and AHW detection was performed using the HeatwaveR package in R (Schlegel and Smit, 2018). This package utilizes the methodology proposed by Hobday et al. (2016) to detect HW events. The annual climatology (i.e. the average temperature of each day of the year since the start of the time serie) of both air and water temperature was computed. HWs were defined as events when the temperature exceeded the 90th percentile of the climatology during three consecutive days. Furthermore, the severity of each event was assessed using the methodology proposed by Hobday et al. (2018).

3.2.2 Satellite observations

Three 2021 Sentinel-2 images of the study site were selected to assess the effect of thecombined AHW and MHW (“HW event”, henceforth) on the seagrass meadow: the first image was taken 5 days before the HW (1st of September 2021), the second image during the HW (6th of September 2021) and the third image one month later (8th of October 2021). Level-2 data were downloaded from the Copernicus open access hub (European Space Agency, 2024a) provided by the European Space Agency (ESA). Level-2 images consist of orthorectified surface reflectance corrected from the effect of the atmosphere using ESA’s standard correction (i.e., Sen2cor, European Space Agency (2024b)).

The SHSI (Equation 3) was computed and mapped for each image. For the pixel containing a field QP (Figure 4), the satellite-derived reflectance was extracted, and compared before and after the HW event.

3.2.3 Emersion time of the seagrass meadow

The spatial distribution of seagrass emersion time during low tide was estimated using bathymetric and water level data. High resolution bathymetry data (Litto3D® product) for the Quiberon intertidal meadow were sourced from the “Service Hydrographique et Océanographique de la Marine” (SHOM, 2021), while one-minute interval water level data were downloaded from Intergovernmental Oceanographic Commission data portal (IOC, 2024), using measurements from the nearest tide gauge at Le Crouesty. A 2.85 m vertical correction was applied to the Litto3D data to align its zero reference with that of the water level data (RAM, SHOM (2022))

Once aligned, the corrected elevation was compared to water height for each pixel and each time step during the HW event. The emersion time was then calculated as the daily total time each pixel remained exposed along the duration of the AHW.

3.3 Statistics

General Linear Models (GLMs) were used to assess relative differences over time in response variables (Spectral Indices) with different treatments (Impacted vs Unimpacted). To analyze the effect of HW on the reflectance indices observed during the lab experiment, the relative change was modeled as a function of Days (1-3: Discrete) with Runs (1-3: Factor) and Timestep within Run (1-6: Factor) as cross random factors. Satellite-derived SHSI were modeled as a function of Date (1-3: Discrete) and Treatment (Impacted vs Unimpacted: Categorical). A General Additive Model (GAM) was used to assess the relationship between relative SHSI change with emersion time. SHSI was modeled as a function of emersion time with a basis spline. All model parameters were estimated within a Bayesian framework using the “brms” and “RStan” packages in R to leverage the stan language (Bürkner, 2021; Carpenter et al., 2017; R Core Team, 2023; Stan Development Team et al., 2020). The response variables were modeled assuming a Gaussian distribution, with weakly informative priors (Student-T(3,0,2.5)). Model parameters were estimated using Markov Chain Monte Carlo (MCMC) sampling, with 4 chains of 5000 iterations and a warm-up of 500.

4 Results

4.1 Laboratory Experiment

4.1.1 Heatwave effect on seagrass reflectance

During the laboratory HW experiment, the seagrass reflectance was drastically impacted by the increase in both air and water temperature (Figure 6). The Control \(R(\lambda)\), displayed a spectral shape typical of seagrass, with a green peak around 560 nm, a valley associated with chlorophyll-a absorption around 665 nm, and a high near-infrared (NIR) plateau beyond 705 nm. This remained stable over time, with only minimal changes in magnitude and spectral features. In contrast, the Treatment \(R(\lambda)\) showed severe changes throughout the experiment. During day 1, the Treatment \(R(\lambda)\) was generally similar to the Control \(R(\lambda)\), despite a slightly less marked peak around 560 nm and slightly lower NIR values from 750 – 900 nm. A drastic decrease was then observed during days 2 and 3 across all wavelengths, particularly in the green – yellow spectral region (from 500 – 650 nm) and in the NIR (from 750 - 900 nm). During day 3, the collapse in \(R(\lambda)\) appeared to stabilize in the NIR, whereas it slightly continued in the green spectral region. At the end of the HW experiment, the \(R(\lambda)\) valley around 665 nm was also less pronounced, suggesting a decrease of chlorophyll-a concentration.

Figure 6: Standardized hyperspectral reflectance of Z. noltei leaves during the HW experiment, showing the Control (Left) and Treatment (Right) measurements. The color indicates the progression along the experiment from the beginning (Day 1: Green), middle (Day 2: Yellow) and end (Day 3: Brown). A min-max standardization was applied to each individual spectrum.

4.1.2 Heatwave effect on radiometric indices

All radiometric indices, \(R''_{665 \, \text{nm}}\), NDVI, GLI and SHSI changed after the experimental heatwave (Figure 7 ; Section 7.2).

At the start of the experiment (day 1), the difference between Treatment and Control was not significant in \(R''_{665 \, \text{nm}}\), NDVI and GLI (Figure 7 A, B & C). During days 2 and 3 the radiometric indices all decreased significantly, with an overall decline of 68%, 31% and 54 for \(R''_{665 \, \text{nm}}\), NDVI and GLI.

Unlike the other metrics, the SHSI of the Treatment was on average 55 % higher than that of the Control in day 1 (Figure 7 D). By day 2, the SHSI exhibited a rapid increase of approximately 241 %, eventually reaching an overall rise of 420 % by day 3.

With a maximum deviation of 420 %, SHSI emerges as the most sensitive index for detecting seagrass browning. Consequently, only this index was considered for the next steps of the study.

Figure 7: Comparison of spectral metrics for detecting reflectance changes of seagrass leaves after a HW. A: Relative difference between the Treatment and the Control over time for A) the second derivative at 665 nm B) the NDVI C) the GLI and D) the SHSI. Points indicate raw data, the line represents a GLM estimates, while the shaded area is the model’s 89 % confidence interval. The dashed lines represent no difference between the Control and the Treatment.

Looking at raw SHSI values revealed clear distinctions between the Control and Treatment groups (Figure 8). On day 1, the SHSI of the Control and Treatment groups were comparable, with median values of -0.11 and -0.08, respectively. By the end of the experiment, seagrasses in the Treatment group exhibited a median SHSI of 0.15, consistent with their visibly darkened appearance. In contrast, the Control group retained a green appearance throughout the experiment, with a median SHSI of -0.07. A negative SHSI was considered indicative of non-impacted seagrasses, while a positive SHSI was used as a marker for impacted seagrasses.

Figure 8: Median of the SHSI across experimental runs, on each day of the experiment. Error structure represent the 89 % confidence interval (Quantiles 0.055 and 0.945 for lower and upper ribbon, respectively). The green line shows values of the Control group while the orange line indicates values of the Treatment group.

4.2 HW of September 2021 in Quiberon, South Brittany

4.2.1 Spectral changes

Sentinel-2 images acquired the 1st and 6th of September 2021 were analysed to assess the short-term impact of a HW on seagrass leaves in South Brittany, France (Figure 9 A and C). Both an AHW and a MHW occurred during the 2021 summer, from the 4th to the 7th and from the 3rd to the 8th of September, respectively (Figure 9 B). Within just a few days, the temperature experienced a sharp increase, from 22.2 to 30.8°C in air, and from 17.7 to 19.3 °C in water temperature. During this period, the 90th percentile of the air temperatures was 25.3 °C and 18.8 °C for the water temperature. The air temperature anomaly of 9.9 °C classified the AHW as a strong event, whereas the 1.7 °C anomaly in water temperature classified the MHW as a moderate event.

Two days after the start of the AHW, the Sentinel-2 image from the 6th of September revealed patches of brown seagrass in the true-color composition (Figure 9 C). Such brown patches were absent from the image from the 1st of September, taken before the HW began (Figure 9 A). Before the HW, all QPs appeared green on the Sentinel-2 image, with similar reflectance spectra, typical of green seagrass leaves (Figure 9 A and D). Their reflectance spectra showed a peak at 560 nm (in the green part of the spectra), low values at 665 nm and a high infrared plateau (> 705 nm). However, on the 6th of September, QPs classified as impacted during the field campaign, showed significant differences in their reflectance spectral shape compared to unimpacted QPs (Figure 9 C and E). The reflectance spectra of brown seagrass were characterized by the loss of the reflectance peak at 560 nm and a decrease in the infrared plateau, which was replaced by a steadily increasing slope up to 940 nm. The darkening of large seagrass patches could also be observed in the true color composition (Figure 9 C)

Figure 9: Intertidal seagrass meadow in South Brittany (France) observed before and during a heatwave (HW). A: RGB color composition of the Sentinel-2 image of the 1st of September 2021 before the HW; C: RGB color composition of the Sentinel-2 image of the 6th of September 2021 on the second day of a strong AHW. The circles correspond to QPs collected on the 10th of September 2021, with unimpacted seagrass in green and impacted seagrass in orange; B: Detection of HW events based on both Air Temperature and Sea Surface Temperature (SST). The solid line represents the daily average temperature, while the dashed line indicates the 90th percentile of the climatology. Coloured areas identify HWs (marine in blue and atmospherical in red). The two vertical dashed lines represent the acquisition dates of the two Sentinel-2 images (01-09-2021 and 06-09-2021); D: Sentinel-2 reflectance of seagrass leaves before the HW for both categories of QPs; E: Sentinel-2 reflectance of seagrass leaves during the HW for both categories of QPs. Average spectral signatures were obtained in areas where QPs corresponded to green and brown seagrasses leaves (green and orange circles, respectively) as identified during the field survey. The shaded areas around the reflectance spectra represent the standard deviation.

4.2.2 SHSI metric applied to Sentinel-2

Using Sentinel-2 data and the QPs, SHSI of green seagrass areas that appeared unimpacted by the HW (Unimpacted QPs ; Figure 9 C), showed minimal change between 1st and 6th of September (3 % ; Figure 10). In contrast, seagrass impacted by the HW that turned brown (Impacted QPs ; Figure 9 C) exhibited significant SHSI changes, showing an increase of 97 % during the HW exposure (Figure 10). One month after the event, on the 8th of October 2021, the SHSI of unimpacted seagrass had increased by 14 % compared to the 1st of September. Regarding impacted seagrass, one month after the event, the SHSI decreased to values comparable to those of unimpacted seagrass. This change reflects an increase of 15 % compared to values recorded on the 1st of September.

Figure 10: Changes in the relative SHSI estimated from Sentinel-2, before (1st of September 2021), during (6th of September 2021) and after (8th of October 2021) a HW in the seagrass meadow of Quiberon (South Brittany, France). The relative SHSI difference was calulated using the 1st of September as a reference. SHSI was calculated for two categories of Quadrat Points (QPs; Figure 9): unimpacted seagrass (green) and impacted seagrass (orange). Points represent the estimated value of the SHSI using a GLM, while the error bar represents the 89% confidence interval.

Using the SHSI (Equation 3), we detected large darkened seagrass patches in the meadow on 6th of September (Figure 11). A total of 26.9 hectares of seagrass turned brown between 1st and 6th of September. The largest brown patch covered nearly 8 hectares. Overall, 24 % of the total seagrass meadow area showed signs of darkening between 1st and 6th of September. Comparing the spatial distribution of darkened patches with the site’s topography revealed that 94.6 % of darkened areas were located above a bathymetric level of 3.9 meters (Figure 11, A and B). One month later, on 8th of October, the previously darkened areas appeared to have regained their green coloration (Figure 11 C).

Figure 11: Sentinel-2 color composition of the seagrass meadow of Quiberon, South Brittany, France, Before (A ; 1st of September 2021), During (B ; 6th of September 2021) and After (C ; 8th of October 2021) the HW and SHSI applied to the same Sentinel-2 images Before (D), During (E) and After (F) the HW.

Additionally, there was a clear relationship between seagrass emersion time and seagrass darkening (Figure 12 ; Section 7.3). Seagrass emerged during less than 13 hours a day were not impacted, whereas seagrass emerged during more than 13 hours turned brown with the maximum darkening occuring on seagrass emerged during more than 14.5 hours daily.

Figure 12: Relative change of the SHSI before and during the HW events as a function of the daily emersion time of seagrass. The line represents a GAM estimate, and the shaded area indicates the standard error. Shaded points represent raw data, each corresponding to a single pixel of the meadow.

5 Discussion

5.1 Effect of heatwaves on Zostera noltei reflectance

This study explored the dual effect of MHWs and AHWs on spectral reflectance of intertidal seagrasses. Significant changes in the reflectance of seagrasses exposed to HWs were observed both in field measurements and the laboratory experiment. Effects were drops of reflectance around 560 nm and around 740 nm (Figure 7) and of the second derivative at 665 nm (R’’665), NDVI, and GLI (Figure 7 A, B & C) indexes. By contrast, SHSI showed a marked increase, reaching up to 600% for some samples by day 3 (Figure 7 D), which demonstrated its sensitivity at quantifying seagrass browning. These changes suggest a progressive reduction in photosynthetic activity as well as potential structural or physiological changes to the leaves, such as degraded pigmentation or altered light absorption. Moreover, we observed a considerable decrease in the reflectance in the NIR region (700 - 1300 nm) in stressed plants, which is mainly related with the destruction of the membranes (chloroplasts, thylakoids, cell walls) due to thermal oxidative stress, and the consequent decrease of light scattering through the leaf surface (Knipling, 1970). The leaf darkening observed in plants subjected to thermal stress was mainly explained by the oxidation of phenolic compounds by the enzyme polyphenol oxidase and with a damage of the photosynthetic apparatus (Allakhverdiev et al., 2008). High temperatures destabilize chloroplasts membranes, disrupting photosystem II and impairing recovery of photosynthetic function. As chlorophyll-a degrades, the ratio of pigments shifts, with pigments like carotenoids becoming more prominent (Dascaliuc et al., 2007; Jones and Clayton-Greene, 1992). These processes can eventually result in growth impairment and plant death.

The Seagrass Heat Shock Index (Equation 3) developed in this study makes it possible to detect reductions in the green and the Red-Edge regions caused by thermal stress using only three reflectance bands (560 740 and 840 nm). These spectral bands are available in current (i.e. Sentinel-2, Pleiades-Neo, WorldView-3, SkySat, and GeoSat-2) and future multi-spectral missions (i.e. Sentinel-2 Next Generation and Landsat Next), opening new perspectives for the detection and monitoring of heatwave effects on intertidal seagrass. As shown in the next section, SHSI can be used to assess the extent and severity of darkening events across intertidal seagrass meadows from high resolution satellite images.

5.2 Satellite observations of HW effects on intertidal seagrass

Using our new Seagrass Heat Shock Index, Sentinel-2 images acquired before, during, and after a HW event made it possible to detect seagrass leaf darkening at high spatial resolution, as well as to describe the spatial distribution of the heatwave impact over the whole meadow. The analyse of satellite images in Quiberon Bay during the 2021 HW revealed that most of the meadow darkening took place in the upper intertidal area, where plants were exposed to air for more than 13 hours per day, while no significant browning was observed in the lower intertidal zone. This contrasted with the experiment, in which seagrasses were exposed to a symmetric tidal cycle, spending equal amounts of time submerged and exposed to air. However, in the field, this configuration occurs only in meadows located at hydrographic zero height (mean sea level). Thus, under the same atmospheric temperatures, bathymetry, and therefore, exposure times will modulate the stress effects observed, as seagrasses located lower in the intertidal may be more resistant to HWs than those in upper fringes. Field observations further revealed that leaf detachment began after the heatwave, causing an apparent decline in seagrass density at the upper intertidal. Zostera noltei, as a species inhabiting the intertidal zone and regularly exposed to air, has developed adaptations to minimize hydric stress, such as its narrow leaves, which help reduce water loss during air exposure periods (Cabaço et al., 2009). Nonetheless, if plants at the upper subtidal are exposed to intense heating and irradiance for prolonged periods, these mechanisms may not be able to prevent degradation of photosynthetic pigments, decreases in photosynthetic capacity and tissue degradation.

During the HW, Sentinel-2 imagery indicated that impacted seagrass patches (i.e., areas exhibiting seagrass darkening) experienced a sharp decline in NDVI. Moreover, the Sentinel-2 image acquired one month after the HW showed that the NDVI of impacted seagrass patches had not recovered to its initial level (Figure 10). NDVI is commonly used as a proxy for seagrass cover (Davies et al., 2024a; Zoffoli et al., 2020), with high NDVI associated with greater seagrass coverage and low NDVI indicating reduced cover. Our findings demonstrate that heatwaves altered the spectral signature of seagrasses, causing a decrease in NDVI—and consequently, in the satellite-derived estimate of seagrass cover—while the actual coverage observed in the field remained stable (see Figure 5). This discrepancy underscores a potential limitation of remote sensing in accurately capturing true seagrass cover during stress events. We demonstrate that such bias could be corrected by the combined analysis of both NDVI and SHSI indices, as proxies of seagrass cover and heatwave impact.

One month after the heatwave, the SHSI in impacted quadrats recovered to similar values to those of the unimpacted ones. Moreover, Sentinel-2 images revealed a partial recovery in areas of the mid to low intertidal, as shown by their return to green color. These results suggested a certain the resilience capacity of Z. noltei patches after thermal disturbance. Zostera noltei is a dynamic, fast-growing colonizer species (Roca et al., 2016), with fast leaf turnover times and growth rates (Borum et al., 2004) and fast rhizome elongation rates (Duarte, 1991). The return to green color in some parts of the meadow after recovery might be due to new leaf growth, which was likely sustained by the energy reserves of the surviving below ground biomass (Hemminga, 1998).

5.3 Ecological implications of heatwaves impact on Zostera noltei

In the present study we observed that heatwave-driven thermal stress changed the spectral reflectance of Z. noltei leaves, indicating photosynthetic pigments degradation, with potential negative consequences on health, growth and survival of the above ground biomass. In this sense, previous experiments on the effects of warming on Zostera noltei have shown contrasting results, depending on the duration and intensity of heating. In some cases, there were no significant changes in photosynthetic performance or survival under short-term moderate stress (Franssen et al., 2014), or under gradual temperature increases (Román et al., 2022). By contrast, studies focused on extreme temperature events have evidenced decreases of photosynthetic performance and leaf tissue integrity (Massa et al., 2009; Román et al., 2023). Therefore, our findings, jointly with previous research, suggest that Z. noltei meadows at the upper intertidal could be a vulnerable species under the future global warming scenarios, in which heatwaves will be more frequent and intense (Climate Change (IPCC), 2023; Stillman, 2019).

The observed tissue degradation, both in the laboratory experiment and the field, can eventually result in decreases of meadows density and cover under recurrent heatwaves. This decline could have cascading effect on species that rely on the provision of refuge and/or food by the seagrass canopy (Zoffoli et al., 2023). Commercially important fish and shellfish populations that rely on seagrass for sustenance and refuge could also be impacted, with negative consequences on local fisheries’ productivity and the livelihoods of coastal communities (Unsworth and Cullen-Unsworth, 2014). Furthermore, as seagrass meadows decline, their capacity to act as a blue carbon sink—critical for climate mitigation—could diminish, contributing to increased atmospheric carbon levels (Armitage and Fourqurean, 2016; Samper-Villarreal et al., 2020). Moreover, the sediment stabilization and wave attenuation could be hampered, ultimately increasing the risk of coastal erosion (Calleja et al., 2007; Folmer et al., 2012; Gacia et al., 1999).

Given the combined effects of temperature extremes, eutrophication, and other anthropogenic pressures, targeted management strategies are essential for enhancing seagrass resilience (Loarie et al., 2009). Approaches such as reducing local stressors, cultivating heat-tolerant genotypes, and investing in restoration initiatives are vital to supporting these ecosystems in a warming climate. Although challenges remain, the adaptability and potential resilience of certain seagrass species offer hope for their persistence amid accelerating ecological shifts.

6 Conclusion

This research investigated the effects of both marine and atmospheric heatwaves on the intertidal seagrass Zostera noltei, a critical component of coastal ecosystems facing increased thermal stress due to climate change. Heat stress effects on seagrass reflectance were documented in the laboratory using a controlled experiment in order to understand how extreme heat events affect seagrass health and to assess the potential of remote sensing for HW impacts monitoring. Our findings revealed that heatwaves lead to substantial declines in seagrass reflectance, particularly in the green and near-infrared regions, likely driven by pigment degradation and structural damage. This change was reflected in significant reductions in key vegetation indices such as NDVI and GLI. The Seagrass Heat Shock Index (SHSI), developed in this study, successfully detected seagrass darkening, a visible symptom of heatwave stress, demonstrating the viability of spectral monitoring to capture early-stage impacts of heat events on intertidal ecosystems. We applied and validated the newly developed Seagrass Heat Shock Index to Sentinel-2 satellite images acquired during an heatwave event, and documented the broad spatial impact of heatwaves on seagrass meadows in Quiberon Bay, France. The correlation between heatwave exposure and darkening of seagrass suggests that remote sensing, combined with targeted field observations, can enhance our understanding of ecosystem responses to climate-driven thermal events. These results advocate for integrating regular spectral monitoring into conservation strategies, as it can help predict seagrass resilience and guide adaptive management practices. Future work should focus on refining remote sensing tools and examining the cumulative effects of repeated heatwave events to support the conservation of intertidal seagrass meadows in a warming world.

7 Annexes

7.1 Annexes A - Temperatures of the experiment

Annexe A1: Comparison of daily maximum temperatures in August measured using an in-situ sensor (blue) and retrieved from Meteo France (orange). The solid line in the middle of the boxplot represents the median, the two ends of the box represent the 25th and 75th percentiles, and the whiskers represent values that are no more than 1.5 times the interquartile range.

Annexe A1: Comparison of daily maximum temperatures in August measured using an in-situ sensor (blue) and retrieved from Meteo France (orange). The solid line in the middle of the boxplot represents the median, the two ends of the box represent the 25th and 75th percentiles, and the whiskers represent values that are no more than 1.5 times the interquartile range.

On average, in situ temperatures were 3 ± 3.2°C higher than those recorded by Meteo France. Additionally, temperatures recorded by Meteo France were more stable than those from the in situ sensors, likely due to the sheltered and shaded location of the Meteo France equipment. This difference was used to adjust heatwave temperatures measured by Meteo France to better reflect the conditions experienced by the seagrasses.

7.2 Annexes B - Outputs of GLMM

Table 1: Outputs of the Generalised Linear model assessing the relationship between time (days) and the relative change to the control of different vegetation indices (Figure 7, Second Derivative, NDVI, GLI, and SHSI). Each index has two terms: Intercept and Slope. The Intercept represents the expected value of the index at the reference time point (e.g., day 0), while the Slope quantifies the effect of time, indicating how the index changes per unit increase in days. The Estimate represents the posterior mean of the regression coefficient, while Std. Error indicates its standard deviation. The Lower 95% CI and Upper 95% CI define the 95% credible interval, showing the range within which the true parameter is expected to fall with 95% probability. The Rhat statistic assesses model convergence (values close to 1 indicate proper convergence). Bulk_ESS and Tail_ESS denote the effective sample sizes for bulk estimation and tail uncertainty, respectively, reflecting the reliability of the posterior estimates.

Index

Term

Estimate

Std.Error

Lower 95% CI

Upper 95% CI

Rhat

Bulk_ESS

Tail_ESS

Second Derivative

Intercept

0.54

0.35

-0.17

1.29

1

5,296.50

3,668.84

Slope

-0.41

0.04

-0.48

-0.34

1

31,580.27

12,321.97

NDVI

Intercept

0.11

0.12

-0.11

0.35

1

3,436.59

2,239.77

Slope

-0.14

0.02

-0.17

-0.10

1

32,903.97

13,585.03

GLI

Intercept

0.20

0.33

-0.45

0.77

1

3,953.87

2,318.63

Slope

-0.26

0.03

-0.32

-0.20

1

27,154.17

13,247.23

SHSI

Intercept

-1.65

1.03

-3.84

0.37

1

7,862.03

8,468.20

Slope

1.82

0.20

1.42

2.22

1

24,434.98

13,093.41

7.3 Annexes C - Outputs of GAM

Table 2: Outputs of the Generalized Additive Model (GAM) assessing the relationship between emersion time per day (hours) and the relative change in SHSI (Figure Figure 12). The Intercept represents the expected SHSI value at the reference level of emersion time. The Estimate indicates the mean effect, with its Standard Error quantifying variability. The t-value and p-value assess statistical significance. The smooth term captures the nonlinear effect of emersion time, where the effective degrees of freedom (edf) reflect the function’s complexity. The F-value and p-value evaluate the significance of this effect. Model fit is described by the adjusted R², deviance explained, and Generalized Cross-Validation (GCV) score, while N represents the sample size.

Component

Term

Estimate

Std Error

t-value

p-value

A. parametric coefficients

(Intercept)

0.488

0.006

81.574

<0.0001

***

Component

Term

edf

Ref. df

F-value

p-value

B. smooth terms

s(emersion_per_day)

4.953

4.999

973.282

<0.0001

***

Signif. codes: 0 <= '***' < 0.001 < '**' < 0.01 < '*' < 0.05

Adjusted R-squared: 0.357, Deviance explained 0.357

GCV : 0.313, Scale est: 0.313, N: 8768

8 References

Allakhverdiev, S.I., Kreslavski, V.D., Klimov, V.V., Los, D.A., Carpentier, R., Mohanty, P., 2008. Heat stress: An overview of molecular responses in photosynthesis. Photosynthesis research 98, 541–550.
Arias-Ortiz, A., Serrano, O., Masqué, P., Lavery, P.S., Mueller, U., Kendrick, G.A., Rozaimi, M., Esteban, A., Fourqurean, J.W., Marbà, N., others, 2018. A marine heatwave drives massive losses from the world’s largest seagrass carbon stocks. Nature Climate Change 8, 338–344.
Armitage, A., Fourqurean, J.W., 2016. Carbon storage in seagrass soils: Long-term nutrient history exceeds the effects of near-term nutrient enrichment. Biogeosciences 13, 313–321.
Bastos, M., Roebeling, P., Alves, F.L., Villasante, S., Magalhães Filho, L., 2023. High risk water pollution hazards affecting aveiro coastal lagoon (portugal)–a habitat risk assessment using InVEST. Ecological Informatics 76, 102144.
Borum, J., Duarte, C.M., Krause-Jensen, D., Greve, T.M., 2004. European seagrasses: An introduction to monitoring and management. Monitoring; Managing of European Seagrasses Project.
Boyer, M., Miller, J., Belanger, M., Hare, E., Wu, J., 1988. Senescence and spectral reflectance in leaves of northern pin oak (quercus palustris muenchh.). Remote Sensing of Environment 25, 71–87.
Brun, F.G., Vergara, J.J., Navarro, G., Hernández, I., Pérez-Lloréns, J.L., 2003. Effect of shading by ulva rigida canopies on growth and carbon balance of the seagrass zostera noltii. Marine Ecology Progress Series 265, 85–96.
Bürkner, P.-C., 2021. Bayesian item response modeling in r with brms and stan. Journal of Statistical Software 100, 1–54. https://doi.org/10.18637/jss.v100.i05
Cabaço, S., Machás, R., Santos, R., 2009. Individual and population plasticity of the seagrass zostera noltii along a vertical intertidal gradient. Estuarine, Coastal and Shelf Science 82, 301–308.
Calleja, M.L., Marbà, N., Duarte, C.M., 2007. The relationship between seagrass (posidonia oceanica) decline and sulfide porewater concentration in carbonate sediments. Estuarine, Coastal and Shelf Science 73, 583–588.
Carlson, D.F., Yarbro, L.A., Scolaro, S., Poniatowski, M., McGee-Absten, V., Carlson Jr, P.R., 2018. Sea surface temperatures and seagrass mortality in florida bay: Spatial and temporal patterns discerned from MODIS and AVHRR data. Remote Sensing of Environment 208, 171–188.
Carpenter, B., Gelman, A., Hoffman, M.D., Lee, D., Goodrich, B., Betancourt, M., Brubaker, M., Guo, J., Li, P., Riddell, A., 2017. Stan: A probabilistic programming language. Journal of Statistical Software 76, 1–32. https://doi.org/10.18637/jss.v076.i01
Climate Change (IPCC), I.P. on, 2023. Summary for policymakers. IPCC, Geneva, Switzerland. https://doi.org/10.59327/IPCC/AR6-9789291691647.001
CMEMS, 2024. European north west shelf/iberia biscay irish seas – high resolution ODYSSEA sea surface temperature multi-sensor L3 observations reprocessed, e.u. Copernicus marine service information (CMEMS). Marine data store (MDS). (Accessed on 17-10-2024). https://doi.org/10.48670/moi-00311
Dascaliuc, A., Ralea, T., Cuza, P., 2007. Influence of heat shock on chlorophyll fluorescence of white oak (quercus pubescens willd.) leaves. Photosynthetica 45, 469–471. https://doi.org/10.1007/s11099-007-0084-3
Davies, B.F.R., Gernez, P., Geraud, A., Oiry, S., Rosa, P., Zoffoli, M.L., Barillé, L., 2023. Multi-and hyperspectral classification of soft-bottom intertidal vegetation using a spectral library for coastal biodiversity remote sensing. Remote Sensing of Environment 290, 113554.
Davies, B.F.R., Oiry, S., Rosa, P., Zoffoli, M.L., Sousa, A.I., Thomas, O.R., Smale, D.A., Austen, M.C., Biermann, L., Attrill, M.J., others, 2024a. A sentinel watching over inter-tidal seagrass phenology across western europe and north africa. Communications Earth & Environment 5, 382. https://doi.org/10.1038/s43247-024-01543-z
Davies, B.F.R., Oiry, S., Rosa, P., Zoffoli, M.L., Sousa, A.I., Thomas, O.R., Smale, D.A., Austen, M.C., Biermann, L., Attrill, M.J., others, 2024b. Intertidal seagrass extent from sentinel-2 time-series show distinct trajectories in western europe. Remote Sensing of Environment 312, 114340. https://doi.org/10.1016/j.rse.2024.114340
Deguette, A., Barrote, I., Silva, J., 2022. Physiological and morphological effects of a marine heatwave on the seagrass cymodocea nodosa. Scientific Reports 12, 7950.
Douay, F., Verpoorter, C., Duong, G., Spilmont, N., Gevaert, F., 2022. New hyperspectral procedure to discriminate intertidal macroalgae. Remote Sensing 14, 346.
Duarte, C.M., 1991. Allometric scaling of seagrass form and productivity. Marine ecology progress series. Oldendorf 77, 289–300.
Electric Blue, 2023. Intertidal chamber.
European Space Agency, 2024b. Sen2Cor: Sentinel-2 atmospheric correction processor.
European Space Agency, 2024a. Copernicus open access hub.
Folmer, E.O., Geest, M. van der, Jansen, E., Olff, H., Michael Anderson, T., Piersma, T., Gils, J.A. van, 2012. Seagrass–sediment feedback: An exploration using a non-recursive structural equation model. Ecosystems 15, 1380–1393.
Franssen, S.U., Gu, J., Winters, G., Huylmans, A.-K., Wienpahl, I., Sparwel, M., Coyer, J.A., Olsen, J.L., Reusch, T.B., Bornberg-Bauer, E., 2014. Genome-wide transcriptomic responses of the seagrasses zostera marina and nanozostera noltii under a simulated heatwave confirm functional types. Marine Genomics 15, 65–73.
Gacia, E., Granata, T., Duarte, C., 1999. An approach to measurement of particle flux and sediment retention within seagrass (posidonia oceanica) meadows. Aquatic Botany 65, 255–268.
Garbulsky, M.F., Peñuelas, J., Gamon, J., Inoue, Y., Filella, I., 2011. The photochemical reflectance index (PRI) and the remote sensing of leaf, canopy and ecosystem radiation use efficiencies: A review and meta-analysis. Remote sensing of environment 115, 281–297.
Gardner, R.C., Finlayson, C., 2018. Global wetland outlook: State of the world’s wetlands and their services to people, in: Ramsar Convention Secretariat. pp. 2020–5.
Gladstone-Gallagher, R.V., Hughes, R.W., Douglas, E.J., Pilditch, C.A., 2018. Biomass-dependent seagrass resilience to sediment eutrophication. Journal of Experimental Marine Biology and Ecology 501, 54–64.
Green, A.E., Unsworth, R.K., Chadwick, M.A., Jones, P.J., 2021. Historical analysis exposes catastrophic seagrass loss for the united kingdom. Frontiers in plant science 12, 629962.
Guerrero-Meseguer, L., Marı́n, A., Sanz-Lázaro, C., 2020. Heat wave intensity can vary the cumulative effects of multiple environmental stressors on posidonia oceanica seedlings. Marine Environmental Research 159, 105001.
Hemminga, M., 1998. The root/rhizome system of seagrasses: An asset and a burden. Journal of Sea Research 39, 183–196.
Hobday, A.J., Alexander, L.V., Perkins, S.E., Smale, D.A., Straub, S.C., Oliver, E.C., Benthuysen, J.A., Burrows, M.T., Donat, M.G., Feng, M., others, 2016. A hierarchical approach to defining marine heatwaves. Progress in oceanography 141, 227–238.
Hobday, A.J., Oliver, E.C., Gupta, A.S., Benthuysen, J.A., Burrows, M.T., Donat, M.G., Holbrook, N.J., Moore, P.J., Thomsen, M.S., Wernberg, T., others, 2018. Categorizing and naming marine heatwaves. Oceanography 31, 162–173.
Infantes, E., Hoeks, S., Adams, M.P., Heide, T. van der, Katwijk, M.M. van, Bouma, T.J., 2022. Seagrass roots strongly reduce cliff erosion rates in sandy sediments. Marine Ecology Progress Series 700, 1–12.
IOC, 2024. Intergovernmental oceanographic commission ; sea level monitoring station - le conquet, france (LECY).
Jankowska, E., Michel, L.N., Lepoint, G., Włodarska-Kowalczuk, M., 2019. Stabilizing effects of seagrass meadows on coastal water benthic food webs. Journal of Experimental Marine Biology and Ecology 510, 54–63.
Jones, R.B., Clayton-Greene, K.A., 1992. The role of photosynthesis and oxidative reactions in leaf blackening of protea neriifolia r. Br. leaves. Scientia Horticulturae 50, 137–145. https://doi.org/10.1016/S0304-4238(05)80016-0
Knipling, E.B., 1970. Physical and physiological basis for the reflectance of visible and near-infrared radiation from vegetation. Remote Sensing of Environment 1, 155–159. https://doi.org/10.1016/S0034-4257(70)80021-9
Lehnert, L.W., Meyer, H., Bendix, J., 2017. Hsdar: Manage, analyse and simulate hyperspectral data in r.
Loarie, S.R., Duffy, P.B., Hamilton, H., Asner, G.P., Field, C.B., Ackerly, D.D., 2009. The velocity of climate change. Nature 462, 1052–1055.
Louhaichi, M., Borman, M.M., Johnson, D.E., 2001. Spatially located platform and aerial photography for documentation of grazing impacts on wheat. Geocarto International 16, 65–70.
Malvern Panalytical, 2023. RS3 software.
Marbà, N., Duarte, C.M., 2010. Mediterranean warming triggers seagrass (posidonia oceanica) shoot mortality. Global change biology 16, 2366–2375.
Mariën, B., Balzarolo, M., Dox, I., Leys, S., Lorène, M.J., Geron, C., Portillo-Estrada, M., AbdElgawad, H., Asard, H., Campioli, M., 2019. Detecting the onset of autumn leaf senescence in deciduous forest trees of the temperate zone. New Phytologist 224, 166–176.
Massa, S., Arnaud-Haond, S., Pearson, G., Serrão, E., 2009. Temperature tolerance and survival of intertidal populations of the seagrass zostera noltii (hornemann) in southern europe (ria formosa, portugal). Hydrobiologia 619, 195–201.
Moore, K.A., Jarvis, J.C., 2008. Environmental factors affecting recent summertime eelgrass diebacks in the lower chesapeake bay: Implications for long-term persistence. Journal of Coastal Research 135–147.
Oiry, S., Davies, B.F.R., Sousa, A.I., Rosa, P., Zoffoli, M.L., Brunier, G., Gernez, P., Barillé, L., 2024. Discriminating seagrasses from green macroalgae in european intertidal areas using high-resolution multispectral drone imagery. Remote Sensing 16. https://doi.org/10.3390/rs16234383
Olmedo-Masat, O.M., Raffo, M.P., Rodrı́guez-Pérez, D., Arijón, M., Sánchez-Carnero, N., 2020. How far can we classify macroalgae remotely? An example using a new spectral library of species from the south west atlantic (argentine patagonia). Remote Sensing 12, 3870.
Peñuelas, J., Munné-Bosch, S., Llusià, J., Filella, I., 2004. Leaf reflectance and photo-and antioxidant protection in field-grown summer-stressed phillyrea angustifolia. Optical signals of oxidative stress? New Phytologist 162, 115–124.
Perkins, S.E., Alexander, L.V., 2013. On the measurement of heat waves. Journal of climate 26, 4500–4517.
R Core Team, 2023. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria.
Ramesh, C., Mohanraju, R., 2020. Seagrass ecosystems of andaman and nicobar islands: Status and future perspective. Environmental & Earth Sciences Research Journal 7.
Roca, G., Alcoverro, T., Krause-Jensen, D., Balsby, T.J.S., Katwijk, M.M. van, Marbà, N., Santos, R., Arthur, R., Mascaró, O., Fernández-Torquemada, Y., others, 2016. Response of seagrass indicators to shifts in environmental stressors: A global review and management synthesis. Ecological Indicators 63, 310–323.
Román, A., Tovar-Sánchez, A., Olivé, I., Navarro, G., 2021. Using a UAV-mounted multispectral camera for the monitoring of marine macrophytes. Frontiers in Marine Science 8, 722698.
Román, M., Gilbert, F., Viejo, R.M., Román, S., Troncoso, J.S., Vázquez, E., Olabarria, C., 2023. Are clam-seagrass interactions affected by heatwaves during emersion? Marine Environmental Research 186, 105906.
Román, S., Vázquez, E., Román, M., Viejo, R.M., Woodin, S.A., Wethey, D.S., Troncoso, J.S., Olabarria, C., 2022. Effects of warming on biological interactions between clams and the seagrass zostera noltei: A case study using open top chambers. Estuarine, Coastal and Shelf Science 276, 108027.
Rouse, J.W., Haas, R.H., Schell, J.A., Deering, D.W., others, 1974. Monitoring vegetation systems in the great plains with ERTS. NASA Spec. Publ 351, 309.
Samper-Villarreal, J., Bolaños, R.C., Heidemeyer, M., Vargas, M.M., Vargas, R.M., 2020. Characterization of seagrasses at two new locations in the eastern tropical pacific (el jobo and matapalito, costa rica). Aquatic botany 165, 103237.
Sawall, Y., Ito, M., Pansch, C., 2021. Chronically elevated sea surface temperatures revealed high susceptibility of the eelgrass zostera marina to winter and spring warming. Limnology and Oceanography 66, 4112–4124.
Schlegel, R.W., Smit, A.J., 2018. heatwaveR: A central algorithm for the detection of heatwaves and cold-spells. Journal of Open Source Software 3, 821. https://doi.org/10.21105/joss.00821
Scott, A.L., York, P.H., Duncan, C., Macreadie, P.I., Connolly, R.M., Ellis, M.T., Jarvis, J.C., Jinks, K.I., Marsh, H., Rasheed, M.A., 2018. The role of herbivory in structuring tropical seagrass ecosystem service delivery. Frontiers in Plant Science 9, 127.
SHOM, 2022. Service hydrographique et océanographique de la marine ; références altimétriques maritimes: Ports de france métropolitaine et d’outre-mer, cotes du zéro hydrographique et niveaux caractéristiques de la marée. Shom, Brest, France.
SHOM, 2021. Service hydrographique et océanographique de la marine ; bathymétrie Litto3D® bretagne 2018-2021.
Skendzic, S., 2023. Drought stress in winter wheat-physiological responses and detection using remote and proximal sensing techniques, in: Smart Life Sciences and Technology for Sustainable Development. pp. 39–39.
Sousa, A.I., Silva, J.F. da, Azevedo, A., Lillebø, A.I., 2019. Blue carbon stock in zostera noltei meadows at ria de aveiro coastal lagoon (portugal) over a decade. Scientific reports 9, 14387.
Stan Development Team, C., others, 2020. RStan: The r interface to stan. R package version 2.21. 2.
Stillman, J.H., 2019. Heat waves, the new normal: Summertime temperature extremes will impact animals, ecosystems, and human communities. Physiology 34, 86–100.
Strydom, S., Murray, K., Wilson, S., Huntley, B., Rule, M., Heithaus, M., Bessey, C., Kendrick, G.A., Burkholder, D., Fraser, M.W., others, 2020. Too hot to handle: Unprecedented seagrass death driven by marine heatwave in a world heritage area. Global change biology 26, 3525–3538.
Thomsen, E., Herbeck, L.S., Viana, I.G., Jennerjahn, T.C., 2023. Meadow trophic status regulates the nitrogen filter function of tropical seagrasses in seasonally eutrophic coastal waters. Limnology and Oceanography 68, 1906–1919.
Thomson, J.A., Burkholder, D.A., Heithaus, M.R., Fourqurean, J.W., Fraser, M.W., Statton, J., Kendrick, G.A., 2015. Extreme temperatures, foundation species, and abrupt ecosystem change: An example from an iconic seagrass ecosystem. Global change biology 21, 1463–1474.
Unsworth, R., Cullen-Unsworth, L.C., 2014. Biodiversity, ecosystem services, and the conservation of seagrass meadows. Coast. Conserv 19, 95.
Unsworth, R.K., Butterworth, E.G., 2021. Seagrass meadows provide a significant resource in support of avifauna. Diversity 13, 363.
Unsworth, R.K., Cullen-Unsworth, L.C., Jones, B.L., Lilley, R.J., 2022. The planetary role of seagrass conservation. Science 377, 609–613.
Ustin, S.L., Jacquemoud, S., 2020. How the optical properties of leaves modify the absorption and scattering of energy and enhance leaf functionality. Remote sensing of plant biodiversity 349–384.
Waycott, M., Duarte, C.M., Carruthers, T.J., Orth, R.J., Dennison, W.C., Olyarnik, S., Calladine, A., Fourqurean, J.W., Heck Jr, K.L., Hughes, A.R., others, 2009. Accelerating loss of seagrasses across the globe threatens coastal ecosystems. Proceedings of the national academy of sciences 106, 12377–12381.
Winters, G., Nelle, P., Fricke, B., Rauch, G., Reusch, T.B., 2011. Effects of a simulated heat wave on photophysiology and gene expression of high-and low-latitude populations of zostera marina. Marine Ecology Progress Series 435, 83–95.
Zahoor, I., Mushtaq, A., 2023. Water pollution from agricultural activities: A critical global review. Int. J. Chem. Biochem. Sci 23, 164–176.
Zoffoli, M.L., Gernez, P., Godet, L., Peters, S., Oiry, S., Barillé, L., 2021. Decadal increase in the ecological status of a north-atlantic intertidal seagrass meadow observed with multi-mission satellite time-series. Ecological Indicators 130, 108033.
Zoffoli, M.L., Gernez, P., Oiry, S., Godet, L., Dalloyau, S., Davies, B.F.R., Barillé, L., 2023. Remote sensing in seagrass ecology: Coupled dynamics between migratory herbivorous birds and intertidal meadows observed by satellite during four decades. Remote Sensing in Ecology and Conservation 9, 420–433. https://doi.org/10.1002/rse2.319
Zoffoli, M.L., Gernez, P., Rosa, P., Le Bris, A., Brando, V.E., Barillé, A.-L., Harin, N., Peters, S., Poser, K., Spaias, L., others, 2020. Sentinel-2 remote sensing of zostera noltei-dominated intertidal seagrass meadows. Remote Sensing of Environment 251, 112020.

Footnotes

  1. Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France↩︎

  2. Consiglio Nazionale delle Ricerche, Istituto di Scienze Marine (CNR-ISMAR), 00133 Rome, Italy↩︎

  3. Bio-littoral, Immeuble Le Nevada, 2 Rue du Château de l’Eraudière, 44300 Nantes, France↩︎

  4. Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France↩︎

  5. Consiglio Nazionale delle Ricerche, Istituto di Scienze Marine (CNR-ISMAR), 00133 Rome, Italy↩︎

  6. Bio-littoral, Immeuble Le Nevada, 2 Rue du Château de l’Eraudière, 44300 Nantes, France↩︎

  7. Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France↩︎

  8. Consiglio Nazionale delle Ricerche, Istituto di Scienze Marine (CNR-ISMAR), 00133 Rome, Italy↩︎

  9. Bio-littoral, Immeuble Le Nevada, 2 Rue du Château de l’Eraudière, 44300 Nantes, France↩︎

  10. Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France↩︎

  11. Consiglio Nazionale delle Ricerche, Istituto di Scienze Marine (CNR-ISMAR), 00133 Rome, Italy↩︎

  12. Bio-littoral, Immeuble Le Nevada, 2 Rue du Château de l’Eraudière, 44300 Nantes, France↩︎

Citation

BibTeX citation:
@online{oiry^[institut_des_substances_et_organismes_de_la_mer,_isomer,_nantes_université,_ur_2160,_f-44000_nantes,_france]2025,
  author = {Oiry\^{}{[}Institut des Substances et Organismes de la Mer,
    ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France{]}, Simon
    and Ffinian Rowe Davies¹, Bede and Rosa¹, Philippe and Debly¹,
    Augustin and Laura Zoffoli\^{}{[}Consiglio Nazionale delle Ricerche,
    Istituto di Scienze Marine (CNR-ISMAR), 00133 Rome, Italy{]}, Maria
    and Barillé\^{}{[}Bio-littoral, Immeuble Le Nevada, 2 Rue du Château
    de l’Eraudière, 44300 Nantes, France{]}, Anne-Laure and Harin³,
    Nicolas and Román¹⁴, Marta and Gernez¹, Pierre and Barillé¹,
    Laurent},
  title = {Heatwave Impacts on Intertidal Seagrass Reflectance: From
    Laboratory Experiment to Satellite Mapping of Seagrass Heat Shock
    Index},
  date = {2025-02-26},
  langid = {en},
  abstract = {Seagrasses play a vital role in coastal ecosystems,
    providing habitat, stabilizing sediments, and contributing to carbon
    sequestration. However, climate change has increased the frequency
    and intensity of heatwaves, posing a significant threat to seagrass
    health. This study investigates the effects of marine and
    atmospheric heatwaves on the spectral reflectance of the intertidal
    seagrass *Zostera noltei*. Laboratory experiments were conducted
    under controlled heatwave conditions, where hyperspectral
    reflectance measurements were taken to assess the impacts over time.
    Heatwaves caused a substantial decline in seagrass reflectance,
    particularly in the green and near-infrared regions, corresponding
    to the browning of green leaves. Key vegetation indices, including
    the Normalized Difference Vegetation Index (NDVI) and Green Leaf
    Index (GLI), showed pronounced reductions under heatwave stress,
    with NDVI values decreasing by up to 34\% and GLI by 57\%. A novel
    metric, the Seagrass Heat Shock Index (SHSI), was developed to
    quantify the transition of seagrass leaves from green to brown,
    demonstrating a strong ability to capture the effects of heatwave
    exposure on seagrass coloration. Multispectral satellite
    observations corroborated the laboratory results, revealing
    widespread browning of seagrass leaves during marine and atmospheric
    heatwave events in South Brittany, France. Notably, darkened
    seagrass patches were observed in intertidal areas exposed to
    temperatures exceeding 32°C for over 13.5 hours per day. These
    findings highlight the potential of spectral reflectance as a tool
    for detecting early signs of heatwave-induced stress in seagrasses,
    offering a valuable method for remote sensing-based habitat
    assessment. The present study underscores the potential of remote
    sensing to capture rapid environmental changes in intertidal zones,
    enabling for continuous monitoring of seagrass meadows under the
    current and future climate regimes.}
}
For attribution, please cite this work as:
Oiry^[Institut des Substances et Organismes de la Mer, ISOMer, Nantes Université, UR 2160, F-44000 Nantes, France], S., Ffinian Rowe Davies¹, B., Rosa¹, P., Debly¹, A., Laura Zoffoli^[Consiglio Nazionale delle Ricerche, Istituto di Scienze Marine (CNR-ISMAR), 00133 Rome, Italy], M., Barillé^[Bio-littoral, Immeuble Le Nevada, 2 Rue du Château de l’Eraudière, 44300 Nantes, France], A.-L., Harin³, N., Román¹⁴, M., Gernez¹, P., Barillé¹, L., 2025. Heatwave impacts on intertidal seagrass reflectance: from laboratory experiment to satellite mapping of seagrass heat shock index [WWW Document]. Remote Sensing of Environment.