Spatiotemporal Variability and Trends of Marine Heat Waves in the Red Sea over 38 Years

: Marine heat waves (MHWs) can have catastrophic consequences for the socio-environmental system. Especially in the Red Sea, which has the world’s second longest coral reef system. Here, we investigate the sea surface temperature (SST) variability and trends, as well as the spatiotemporal characteristics of marine heat waves (MHWs) in the Red Sea, using high resolution daily gridded (1/20 ◦ ) SST data obtained from the Copernicus Marine Environment Monitoring Service (CMEMS) for the period 1982–2019. Results show that the average warming rate was about 0.342 ± 0.047 ◦ C/decade over the entire Red Sea over the whole study period. The Empirical Orthogonal Function (EOF) analysis reveals that the maximum variability is over the central part of the Red Sea, while the minimum variability is in the southernmost part of the Red Sea. Over the last two decades (2000–2019), we have discovered that the average MHW frequency and duration increased by 35% and 67%, respectively. The results illustrate that the MHW frequency and duration trends have increased by 1.17 counts/decade and 1.79 days/decade, respectively, over the study period. The highest annual MHW frequencies were detected in the years 2018, 2019, 2010, and 2017. A strong correlation (R = 0.89) was found between the annual MHW frequency and the annual mean SST. We found that all MHW characteristic (frequency, duration, and cumulative intensity (i cum )) trends increased over the study period. The highest annual MHW frequencies were detected in the years 2018, 2019, 2010, and 2017. There were no MHWs found during the year 1984. A strong correlation was found between the annual MHW frequency and the annual mean SST. Over the period from 2000 to 2019, we discovered that the average MHW frequency and duration increased by 35% and 67%, respectively. Further research is needed to investigate the causes of MHWs and their effects on marine life in the Red Sea. We believe that our study is just the beginning of a series of studies on this hot topic.


Introduction
Heat waves are characterized as extended periods of unusually hot weather that have become more frequent and longer in recent decades, posing a threat to human health and ecosystems. A similar phenomenon known as a 'Marine Heat Wave' (MHW) has been observed in the oceans, endangering marine ecosystems and productivity [1][2][3][4]. The most recent definition of MHWs is discrete periods of abnormally warm sea surface temperatures that can last from five days to months above the 90th percentile threshold of SST climatology [5][6][7]. Despite extensive knowledge of global SST changes, the study and detection of MHWs and associated climate processes remains limited [2].
The Red Sea has a rich and diverse ecosystem [8]. It has the world's second longest coral reef system [9]. The importance of the Red Sea is highlighted by its high tourism volume and aquaculture, and also its relation to the gas and oil and fishing industries. Furthermore, the Red Sea is a hub for shipping activity, connecting European ports with China and Eastern Asia through Suez Canal.
The Red Sea has a negative water balance, which means that evaporation exceeds precipitation and combined river runoff [10,11]. It is considered to be an arid climate region with little precipitation and runoff, as there is a lack of any permanent rivers that flow into it [12][13][14].
The impact and consequences of MHW on Red Sea marine life are now critical for decision makers and stakeholders. There are no published MHW studies of the Red Sea. However, some studies in the Rea Sea have discovered a link between extreme heat events such as ocean heat waves and the spread of the El Nino-Southern Oscillation (ENSO) [15]. In addition, the authors have found several SST studies for the Red Sea region, addressing trends and changes in ocean surface temperature [16][17][18][19][20][21][22][23][24][25][26][27].
According to [21], the annual mean SST has changed from 27.4 • C from 1985 to 1993 to 28.1 • C from 1994 to 2007, which was accompanied with significant warming in the mid-1990s. The authors of [15] discovered that the Red Sea maximum SST is warming at a rate of 0.17 • C/decade, with the northern Red Sea alone warming from 0.4 to 0.45 • C/decade from 1982 to 2015, which exceeds the global trend of 0.11 • C/decade [15]. The amplification of recent summer warming and decadal variability of the Red Sea SST was studied in [25]. They found that the significant warming which began in the mid-1990s precisely corresponded with the onset of a positive phase of the Atlantic Multi-Decadal Oscillation (AMO) and a negative phase of the Silk Road Pattern (SRP) in 1997.
In this study, we have investigated the spatiotemporal characteristics and occurrence of MHWs in the Red Sea from 1982 to 2019 using high resolution SST OSTIA (Operational SST and Ice Analysis) data (0.05 • × 0.05 • ). In addition, the SST spatial trend and variability were examined. Section 2 describes the data and methodology of SST analysis. Section 3 reports the SST climatology, trend, and interannual variability. Section 4 shows the MHW results over the study period. Section 5 presents the discussion and Section 6 provides the conclusions.

Data and Methods
The Red Sea is an almost enclosed marginal tropical sea characterized by an extended shape that separates the African and Asian continents and is connected to the Indian Ocean through the Bab El-Mandab Strait, as shown in Figure 1. It has a moderate size of 4.5 × 10 6 km 3 , which is roughly one fifth the size of the Mediterranean Sea. The Bab El-Mandeb strait, which connects the basin to the open ocean, is a small and shallow channel (sill depth of 160 m and the narrowest point is 25 km wide). The study area extending from 12.5 • N to 30 • N and from 32 • E to 44 • E has an average width of~220 km and a length of 2000 km, with a maximum centered depth of 2200 m ( Figure 1). However, some studies in the Rea Sea have discovered a link between extreme heat events such as ocean heat waves and the spread of the El Nino-Southern Oscillation (ENSO) [15]. In addition, the authors have found several SST studies for the Red Sea region, addressing trends and changes in ocean surface temperature [16][17][18][19][20][21][22][23][24][25][26][27]. According to [21], the annual mean SST has changed from 27.4 °C from 1985 to 1993 to 28.1 °C from 1994 to 2007, which was accompanied with significant warming in the mid-1990s. The authors of [15] discovered that the Red Sea maximum SST is warming at a rate of 0.17 °C/decade, with the northern Red Sea alone warming from 0.4 to 0.45 °C/decade from 1982 to 2015, which exceeds the global trend of 0.11 °C/decade [15]. The amplification of recent summer warming and decadal variability of the Red Sea SST was studied in [25]. They found that the significant warming which began in the mid-1990s precisely corresponded with the onset of a positive phase of the Atlantic Multi-Decadal Oscillation (AMO) and a negative phase of the Silk Road Pattern (SRP) in 1997.
In this study, we have investigated the spatiotemporal characteristics and occurrence of MHWs in the Red Sea from 1982 to 2019 using high resolution SST OSTIA (Operational SST and Ice Analysis) data (0.05° × 0.05°). In addition, the SST spatial trend and variability were examined. Section 2 describes the data and methodology of SST analysis. Section 3 reports the SST climatology, trend, and interannual variability. Section 4 shows the MHW results over the study period. Section 5 presents the discussion and Section 6 provides the conclusions.

Data and Methods
The Red Sea is an almost enclosed marginal tropical sea characterized by an extended shape that separates the African and Asian continents and is connected to the Indian Ocean through the Bab El-Mandab Strait, as shown in Figure 1. It has a moderate size of 4.5 × 10 6 km 3 , which is roughly one fifth the size of the Mediterranean Sea. The Bab El-Mandeb strait, which connects the basin to the open ocean, is a small and shallow channel (sill depth of 160 m and the narrowest point is 25 km wide). The study area extending from 12.5° N to 30° N and from 32° E to 44° E has an average width of ~220 km and a length of ~2000 km, with a maximum centered depth of 2200 m ( Figure 1).

Sea Surface Temperature Climatological Data
The daily Sea Surface Temperature (SST) data used were freely obtained from the Copernicus Marine Environment Monitoring Service website (CMEMS) (https://marine. copernicus.eu (accessed on 20 April 2021)) from 1 January 1982 to 31 December 2019. The data website is https://resources.marine.copernicus.eu/?option=com_csw&view= details&product_id=SST_GLO_SST_L4_REP_OBSERVATIONS_010_011 (accessed on 23 March 2021). The product name is 'SST_GLO_SST_L4_REP_OBSERVATIONS_010_011'. The CMEMS OSTIA (Operational SST and Ice Analysis) reprocessed analysis product is a satellite and in situ foundation SST analysis created by the OSTIA [28]. The SST dataset includes daily gap-free maps of sea surface temperature and ice concentration (known as the L4 product) with a horizontal grid resolution of 0.05 • × 0.05 • [29]. CMEMS OSTIA data for the Red Sea were extracted from the global data, providing a 16534-point regularly gridded dataset spanning 13,879 days from 1 January 1982 to 31 December 2019.

Trend and Interannual Variability of SST
The Empirical Orthogonal Functions (EOFs) analysis [30] was used to evaluate the dominant spatiotemporal characteristics of SST variability from 1982 to 2019, utilizing the monthly de-seasonalized and detrended SST dataset to focus on non-seasonal variability. Thus, before decomposing the EOF, the seasonal cycle and linear trends were removed from the monthly means of SST anomalies at each grid cell, and each point time series was divided by its standard deviation to normalize the SST [31,32]. The mean seasonal cycle was estimated by calculating the mean monthly values for each calendar month at each grid point. SST anomalies were computed by removing the historical climatological mean value at each grid point from the SST values at the same location (grid). The strong seasonal signal was removed from the datasets at each grid cell to produce de-seasonalized maps and time series [33]. This was accomplished by subtracting the climatological mean of each calendar month from the corresponding months in all years. The least squares method [34] was used to estimate linear trends in de-seasonalized monthly sea surface temperature anomalies (SSTA). The Modified Mann-Kendall test [35,36] was used to test the statistical significance of the trends. In this study, the significance level was set to p ≤ 0.05, ensuring that all given trends are statistically significant at least at the 95% confidence level. The error of linear trend at the 95% confidence level was calculated using [34] the formula which was mentioned in [37].

Marine Heat Waves (MHWs) Calculations
The authors of [5] used the "discrete prolonged anomalously warm-water events in a location" definition of MHWs, using the same daily SST dataset for the 1982-2019 period. The MATLAB toolbox [38] was used to detect periods of time where temperatures were above the 90th percentile threshold compared to the climatology (1982-2019), and with a duration of at least 5 days. The climatological mean has been calculated from the daily Sea Surface Temperature (SST) data from 1 January 1982 to 31 December 2019; the equations of MHW main characteristics are according to [5] and described in [7]. Daily SSTs were spatially averaged over the Red Sea (12.5 • -30 • N, 32 • -44 • E) in our study to provide a regional daily SST time series spanning 1982-2019. In order to identify the MHWs in the Red Sea, the MHW definition was then applied to the SST time series.
Every MHW event has a start and end date defined by a mean (i mean ) and a maximum intensity (i max ) in • C (mean and spatiotemporal maximum temperature compared to the threshold over the MHW event duration (in days)), as described in [39]. The MHW cumulative intensity (i cum ) in • C days is defined as the spatiotemporal summation of daily temperature compared to the threshold over the event duration as mentioned by [39]. Following the methodology of [4,40], the annual statistics were computed, including the frequency of events (i.e., the number of discrete events occurring each year). The MHW total number of days in each year can be used to examine how the frequency of occurrences in counts has changed over time, using the metrics for each location [5]. According to [5,7], the main characteristics of MHWs are maximum intensity (i max ) in • C, mean intensity (i mean ) in • C, and cumulative intensity (i cum ) in ( • C days): where T m is the climatological OSTIA SST mean calculated over the period from 1 January 1982 to 31 December 2019, to which all values are relative; j is the day of year; and t s and t e are the start and end dates of the MHW event, respectively. Figure 2a depicts the temporal evolution of the daily mean climatology of Red Sea SST ( • C) from 1982 to 2019. The highest temperatures (>29 • C) were in summer and early autumn (i.e., from July to October), while the lowest temperatures (<25 • C) were in winter and early spring (i.e., from January to March). The seasonal temperature difference between summer and winter is about 6 • C, which is consistent with [22]. Figure 2b shows the monthly de-seasonalized SST temporal trend over the study period. The higher deseasonalized SST values were noticed in the last 4 years of the study period (i.e., 2016-2019). The Red Sea temporal trend was about 0.32 ± 0.03 • C/decade. quency of events (i.e., the number of discrete events occurring each year). The MHW to number of days in each year can be used to examine how the frequency of occurrence counts has changed over time, using the metrics for each location [5]. According to [5 the main characteristics of MHWs are maximum intensity (imax) in °C, mean intensity (im in °C, and cumulative intensity (icum) in (°C days):

SST Climatology, Trend, and Interannual Variability
where Tm is the climatological OSTIA SST mean calculated over the period from 1 Janu 1982 to 31 December 2019, to which all values are relative; j is the day of year; and ts a te are the start and end dates of the MHW event, respectively. Figure 2a depicts the temporal evolution of the daily mean climatology of Red SST (°C) from 1982 to 2019. The highest temperatures (>29 °C) were in summer and ea autumn (i.e., from July to October), while the lowest temperatures (<25 °C) were in win and early spring (i.e., from January to March). The seasonal temperature difference tween summer and winter is about 6 °C, which is consistent with [22]. Figure 2b sho the monthly de-seasonalized SST temporal trend over the study period. The higher seasonalized SST values were noticed in the last 4 years of the study period (i.e., 20 2019). The Red Sea temporal trend was about 0.32 ± 0.03 °C/decade.  The spatial trend map of the de-seasonalized SSTA from 1982 to 2019 is shown in Figure 3. At the 95% confidence interval, a statistically significant trend was observed over the entire region. The linear trend over the Red Sea is not uniform, varying from 0.1 to 0.47 °C/decade. The basin averaged warming trend rate was about 0.342 ± 0.047 °C/decade. The highest trend was found in the deep area of the Red Sea (i.e., depth >~1200 m), while lower values were found in the Gulfs of Suez and Aqaba and the southern part of the Red Sea. The spatial trend map of the de-seasonalized SSTA from 1982 to 2019 is shown in Figure 3. At the 95% confidence interval, a statistically significant trend was observed over the entire region. The linear trend over the Red Sea is not uniform, varying from 0.1 to 0.47 • C/decade. The basin averaged warming trend rate was about 0.342 ± 0.047 • C/decade. The highest trend was found in the deep area of the Red Sea (i.e., depth >~1200 m), while lower values were found in the Gulfs of Suez and Aqaba and the southern part of the Red Sea. The first two EOF modes of SST account for 73% of the non-seasonal variance. The spatial pattern of the first EOF (EOF1) (accounting for approximately 52% of the nonseasonal variance) shows a positive anomaly over the entire Red Sea (Figure 4a), showing an in-phase oscillation. The highest variability was found over central and north part of The first two EOF modes of SST account for 73% of the non-seasonal variance. The spatial pattern of the first EOF (EOF1) (accounting for approximately 52% of the nonseasonal variance) shows a positive anomaly over the entire Red Sea (Figure 4a), showing an in-phase oscillation. The highest variability was found over central and north part of Red Sea, while the lowest variability was detected over the Gulfs of Suez and the Aqaba and the southernmost region of the Red Sea. Strong and continuous northeasterly monsoon winds from the Arabian Sea redirected by mountains into southeasterly winds along the axis of the Red Sea [41] may have a thermal effect that is responsible for the lowest variability over the southern Red Sea. The lower variability in the Gulfs of Suez and Aqaba is due to strong northwesterly winds that penetrate the western edge of the Sinai Peninsula, blow southwards, and sometimes reach the southern end of the Red Sea. Its temporal coefficient PC1 shows a strong interannual variability, reaching its peak in the summers of 1995, 1997, and 2019.

Main Characteristics and Trends of Marine Heat Waves (MHWs)
The mean and spatial trend of the main characteristics of MHWs (frequency, days, duration, and intensities) over the entire study period are shown in Figure 5. The trend maps are superimposed with non-significant values (p > 0.05). Figure 5a demonstrates that the MHW frequency varied from 1 to 2.5 counts, with maximum values of >2 counts found

Main Characteristics and Trends of Marine Heat Waves (MHWs)
The mean and spatial trend of the main characteristics of MHWs (frequency, days, duration, and intensities) over the entire study period are shown in Figure 5. The trend maps are superimposed with non-significant values (p > 0.05). Figure 5a demonstrates that the MHW frequency varied from 1 to 2.5 counts, with maximum values of >2 counts found in the coastal areas of Sudan, Eritrea, Yemen and Saudi Arabia. The MHW frequency trend values were statistically significant (p < 0.05) over the whole domain. The MHW frequency trend ranged between 0 and 2 counts/decade. The highest MHW frequency trend (>1.8 counts/decade) was detected in the deep area (i.e., depth > 1000 m) between Eritrea and south coast of Saudi Arabia, as shown in Figures 1 and 5b. The MHW total days spatial patterns in the study period (1982-2019) varied between 15 and 30 days, with maximum values of >25 days found inside the Gulfs of Suez and Aqaba (Figure 5c). This could be due to the Gulfs of Suez and Aqaba being well resolved with the high-resolution OSTIA data. In general, the northern Red Sea has higher values of MHW total days than the southern region. The MHW total days trend values were statistically significant (p < 0.05) across the whole domain (Figure 5d). We discovered a similarity among the p-values between the MHW frequency and total days trend fields (Figure 5b,d). The MHW mean duration ranged from 5 to 15 days, with maximum values of >12 days north of the Red Sea for the Gulfs of Suez and Aqaba (Figure 5e). The maximum MHW duration spatial trend significant values (p < 0.05; >8 days/decade) were found in the middle of the Red Sea (i.e., between~lat 18 • E and 21 • E) (Figure 5f). Both the MHW duration and the i cum spatial patterns revealed similarities between the two fields (Figure 5e,g). The MHW i cum values ranged from 10 to 25 ( • C days), with maximum values of >22 ( • C days) north of the Red Sea, excluding the Gulfs of Suez and Aqaba (i.e., between~lat 22 • E and 28 • E) (Figure 5g). The maximum MHW i cum significant trend values (p < 0.05; >10 • C days/decade) were discovered in the same areas as the MHW duration trend (Figure 5f ranged from 10 to 25 (°C days), with maximum values of >22 (°C days) north of the Red Sea, excluding the Gulfs of Suez and Aqaba (i.e., between ~lat 22° E and 28° E) (Figure 5g). The maximum MHW icum significant trend values (p < 0.05; >10 °C days/decade) were discovered in the same areas as the MHW duration trend (Figure 5f,h). The MHW imax showed values of >1.8 °C throughout most of the Red Sea except the southernmost part beside Bab El-Mandab (Figure 5i). The MHW maximum imax significant trend values (p < 0.05; >2.5 °C) were found in the middle region of the Red Sea (i.e., between ~lat 19° E and 21° E) and in the north of the Red Sea close to the Gulf of Suez (i.e., between ~lat 26°E and 28° E) off the Egyptian coast (Figure 5j).

Marine Heat Waves (MHWs) Temporal Variation
In this section, the authors present the temporal basin averaged annual MHW main characteristics (frequency, duration, and i cum ) from 1982 to 2019. The goal was to investigate temporal changes in Red Sea MHW characteristics such as frequency, duration, and i cum , as well as their relationship to SST ( Figure 6). The MHW matrices for the Red Sea region show a clear ascending trend over the last decade, as shown in Figure 6.
From 1982 to 1994, the annual mean MHW frequency was <1, while from 1994 to 2019, the frequency of MHWs increased significantly, reaching a peak in 2018 (7 counts) (Figure 6a). These findings are consistent with the monthly SST time series shown in Figure 2b, from 1982 to 2019. We noticed that no MHWs were detected in 1984. We discovered a strong correlation (R = 0.89) between the annual mean SST and MHW frequency, (Figure 6b).
The results show a general rise in annual MHW duration from approximately 6 days to 21 days in 2018 (Figure 6c), while i cum ranged from 8 • C days to 30 • C days (Figure 6e). The highest annual MHW duration and i cum were found during 2018 (~21 days and 30 • C days). We discovered that the average frequency and duration of MHWs increased by 35% and 67%, respectively, over the last two decades (2000-2019).
The correlation coefficient between MHW duration and annual mean SST was 0.64, while the correlation coefficient between MHW i cum and SST was 0.72.

Marine Heat Waves (MHWs) Temporal Variation
In this section, the authors present the temporal basin averaged annual MHW main characteristics (frequency, duration, and icum) from 1982 to 2019. The goal was to investigate temporal changes in Red Sea MHW characteristics such as frequency, duration, and icum, as well as their relationship to SST (Figure 6). The MHW matrices for the Red Sea region show a clear ascending trend over the last decade, as shown in Figure 6. From 1982 to 1994, the annual mean MHW frequency was <1, while from 1994 to 2019, the frequency of MHWs increased significantly, reaching a peak in 2018 (7 counts) (Figure 6a). These findings are consistent with the monthly SST time series shown in Figure 2b, from 1982 to 2019. We noticed that no MHWs were detected in 1984. We discovered a strong correlation (R = 0.89) between the annual mean SST and MHW frequency, (Figure 6b).

Examples of MHW Events
In this section, we show two example MHW events discovered in our MHW analysis between 1982 and 2019. In terms of duration and i max , we chose two MHW events in 2019. Table 1 shows examples of the Red Sea MHW event characteristics (mean duration, i max , i mean , and i cum ) over the study period (1982-2019) with the corresponding MHW onset (start date), MHW end (end date), and longitude and latitude. Table 1. Examples of the Red Sea MHW event characteristics (mean duration (days), i max ( • C), i mean ( • C), and i cum ( • C days)) over the study period (1982-2019) with the corresponding MHW onset , MHW end , and longitude and latitude (in degrees).

Discussion
Regarding to our SST climatology and trends analysis (Figures 2 and 3), the highest temperatures (>29 °C, Figure 2a) were found in summer and early autumn. The lowest temperatures (<25 °C, Figure 2a) were recorded in winter and early spring. We found the seasonal temperature difference between summer and winter to be approximately 6 °C (Figure 2a) in agreement with [22]. According to Figure 2b, the Red Sea temporal SST trend was approximately 0.32 ± 0.03 °C/decade from 1982 to 2019, which was higher than [26], possibly due to the use of longer time series and high-resolution OSTIA data.

Discussion
Regarding to our SST climatology and trends analysis (Figures 2 and 3), the highest temperatures (>29 • C, Figure 2a) were found in summer and early autumn. The lowest temperatures (<25 • C, Figure 2a) were recorded in winter and early spring. We found the seasonal temperature difference between summer and winter to be approximately 6 • C (Figure 2a) in agreement with [22]. According to Figure 2b, the Red Sea temporal SST trend was approximately 0.32 ± 0.03 • C/decade from 1982 to 2019, which was higher than [26], possibly due to the use of longer time series and high-resolution OSTIA data.
Many efforts have been made to investigate the Red Sea SST trends using NOAA OI SST coarse resolution data. The authors of [24] showed that August and February are the warmest and coolest months in the Red Sea, with a significant trend of 0.5 • C and 0.3 • C/decade, respectively. According to [15], the overall warming trend of the Red Sea maximum SST was approximately 0.17 • C/decade, while the northern Red Sea alone warmed by between 0.4 and 0.45 • C/decade from 1982 to 2015. The authors of [42] examined the warming of the Red Sea SST from 1982 to 2016 and found an intense warming trend of 0.29 • C/decade. Recently, [26] demonstrated that the SST warming trend ranged from 0.1 • C/decade in the south to 0.4 • C/decade in the north Red Sea between 1982 and 2017. Our analysis revealed a high spatial variability of SST trends over the Red Sea ( Figure 3) in the period from 1982 to 2019, varying from 0.10 to 0.47 • C/decade. The basin averaged warming rate was approximately 0.342 ± 0.047 • C/decade. The highest trend was found in the deep area of the Red Sea (i.e., depth >~1200 m), while lower values were found in the Gulf of Suez, the Gulf of Aqaba, and the southern part of the Red Sea.
The EOF analysis ( Figure 4) indicated that the highest SST variability was found over the central and northern part of the Red Sea, while the lowest variability was found over the Gulfs of Suez and Aqaba, as well as the southernmost part of the Red Sea. This could be attributed to the persistence of monsoon winds from the Arabian Sea in the southern Red Sea and northwesterly winds from the Sinai Peninsula in the northern Red Sea [41].
In Figure 5, the authors investigated the mean and spatial trend of the main characteristics of MHWs (frequency, days, duration, and intensities) for the Red Sea over the study period (1982-2019). The maximum MHW frequency values of >2 counts were found in the coastal areas of Sudan, Eritrea, Yemen, and Saudi Arabia. The MHW frequency trend values were statistically significant (p < 0.05) over the whole Red Sea. According to Figures 1 and 5b, the deep area (i.e., depth > 1000 m) between Eritrea and the south coast of Saudi Arabia has the highest MHW frequency trend (>1.8 counts/decade). The MHW total days' spatial patterns (Figure 5c) show maximum values of more than 25 days in the Gulfs of Suez and Aqaba. This may be due to the fact that the Gulfs of Suez and Aqaba are well resolved using high-resolution OSTIA data. Figure 5b,d show a similarity among the p-values between the MHW frequency and total days trend fields. We found the maximum MHW duration spatial trend significant values (p < 0.05; >8 days/decade) to be in the middle of the Red Sea (i.e., between~lat 18 • E and 21 • E), as demonstrated by Figure 5f. Maximum MHW i cum values of >22 ( • C days) were found in the northern Red Sea, excluding the Gulfs of Suez and Aqaba (i.e., between lat 22 • E and 28 • E) (Figure 5g). The MHW i max shows high values (>1.8 • C) throughout most of the Red Sea (Figure 5i), with the exception of the southernmost part near Bab El-Mandab. The MHW matrices for the Red Sea region show a clear ascending trend over the last decade ( Figure 6). The highest annual MHW duration and i cum were found during the year 2018 (~21 days and 30 • C days). The correlation coefficients between the annual mean SST, MHW duration, and MHW i cum were 0.64 and 0.72, respectively.
The MHW duration event of 236 days was detected in the Gulf of Suez in 2019 from 10 May 2019 to 31 December 2019 (Figure 7a,b). The high MHW i max event of 4.21 • C was discovered in May 2019 off the Sudanese coast (i.e., between~lat 18.5 • E and 20.2 • E), (Figure 7c,d).

Conclusions
This study provides a comprehensive analysis of SST variability, trends, and spatiotemporal patterns of marine heat waves in the Red Sea from 1982 to 2019. To the authors' knowledge, this is the first study to examine marine heat waves and SST trends in the Red Sea using high-resolution data. A statistically significant SST trend was observed over the Red Sea. The EOF analysis showed that the highest interannual variability was detected over the central and northern part of the Red Sea, while the lowest variability was found over the Gulfs of Suez and Aqaba and the southernmost part of the Red Sea.
We found that all MHW characteristic (frequency, duration, and cumulative intensity (i cum )) trends increased over the study period. The highest annual MHW frequencies were detected in the years 2018, 2019, 2010, and 2017. There were no MHWs found during the year 1984. A strong correlation was found between the annual MHW frequency and the annual mean SST. Over the period from 2000 to 2019, we discovered that the average MHW frequency and duration increased by 35% and 67%, respectively. Further research is needed to investigate the causes of MHWs and their effects on marine life in the Red Sea. We believe that our study is just the beginning of a series of studies on this hot topic.