Variability of Heat and Water Fluxes in the Red Sea Using ERA5 Data (1981–2020)

: The study of heat and water ﬂuxes is one of the most essential components for understand-ing the interactions and exchanges between the ocean and atmosphere. Heat transfer across the air–sea interface is an important process in ocean–atmosphere dynamics. In this study, a 40-year (1981–2020) high-resolution (0.25 ◦ × 0.25 ◦ ) ERA-5 reanalysis dataset from the European Centre for Medium-Range Weather Forecasts (ECMWF) is used to estimate the variability and trends of heat and water ﬂux components in the Red Sea. The results show that the surface net heat ﬂux is negative (loss) in the Northern Red Sea (NRS) and positive (gain) in the Southern Red Sea (SRS). The highest seasonal surface net heat ﬂux is observed in the spring and early summer, while the lowest is reported in the winter. A signiﬁcant linear trend is found in the surface net heat ﬂux over the NRS and SRS, with values of about − 0.12 ± 0.052 (W/m 2 )/yr and +0.20 ± 0.021 (W/m 2 )/yr, respectively. The annual mean surface net water ﬂux loss to the atmosphere over the entire Red Sea is +1.46 ± 0.23 m/yr. The seasonal surface net water ﬂux peak occurs in winter as a result of the northeast monsoon wind, which increases evaporation rate over the whole length of the Red Sea. The highest surface net water ﬂux (+2.1 m/yr) is detected during 2020, while the lowest value (+1.3 m/yr) is observed during 1985.


Introduction
The Red Sea is a narrow, elongated basin that connects northeastern Africa to the Arabian Peninsula [1], as shown in Figure 1. The surface buoyancy fluxes caused by the exchange of heat and freshwater between the atmosphere and the ocean are thought to drive the Red Sea circulation [2]. It is completely bordered by dry and semiarid regions and receives very little precipitation and very little seasonal freshwater along the Eritrean coast [3,4]. At the same time, the Red Sea has one of the highest evaporation rates in the world, reaching 2 m/yr [5] and sometimes exceeding 3 m/yr in the Northern Red Sea (NRS) according to [6]. Because exchanges between the Red Sea and the Mediterranean via the Suez Canal are negligible, heat and freshwater exchanges at Bab Al Mandab act as a strong constraint on Red Sea heat and freshwater surface fluxes [5]. Furthermore, it is one of the enclosed seas, with warm and highly saline outflow [2,7], and it contributes significantly to the circulation of the Gulf of Aden and the Arabian Sea.
Ref. [9] calculated the net heat flux from seven stations on the Red Sea's east coast using meteorological data from 1970 to 1984. They discovered that the Red Sea's net heat flux was −22 W/m 2 . Ref. [17] mentioned that the Red Sea's surface deficit (7 W/m 2 ) was compensated by the advective heat surplus into the Red Sea (+6 W/m 2 ) through Bab AI-Mandab Strait, resulting in a total net heat flux of −1 w/m 2 . They used climatological monthly mean data at four coastal stations along the Red Sea's eastern coast. One station covered from 1961 to 1990, and the other three were from 1985 to 1991. Several studies have calculated the Red Sea's net heat and freshwater fluxes [1,2,5,6,[8][9][10][11][12][13][14][15][16].
Ref. [9] calculated the net heat flux from seven stations on the Red Sea's east coast using meteorological data from 1970 to 1984. They discovered that the Red Sea's net heat flux was −22 W/m 2 . Ref. [17] mentioned that the Red Sea's surface deficit (7 W/m 2 ) was compensated by the advective heat surplus into the Red Sea (+6 W/m 2 ) through Bab AI-Mandab Strait, resulting in a total net heat flux of −1 w/m 2 . They used climatological monthly mean data at four coastal stations along the Red Sea's eastern coast. One station covered from 1961 to 1990, and the other three were from 1985 to 1991.
Ref. [5] estimated the Red Sea's heat and freshwater budgets using historical air-sea flux data from a revised version of the Comprehensive Ocean-Atmosphere Data Set (UWM/COADS). This dataset was based on the initial COADS released by [18] and was objectively analyzed on a 1° × 1° global grid. It covered the 45-yr period from January 1945 Ref. [5] estimated the Red Sea's heat and freshwater budgets using historical airsea flux data from a revised version of the Comprehensive Ocean-Atmosphere Data Set (UWM/COADS). This dataset was based on the initial COADS released by [18] and was objectively analyzed on a 1 • × 1 • global grid. It covered the 45-yr period from January 1945 to December 1989. The study discovered that the Red Sea annual surface net heat flux was −8 ± 2 W/m 2 based on the UWM/COADS dataset. The annual net water flux rate evaporation (E)-precipitation (P) was found to be 1.60 ± 0.35 m/yr, the UWM/COADS estimated P was 0.15 m/yr, and evaporation rate was E = 1.75 m/yr.
Ref. [10] investigated the Red Sea's heat and freshwater budgets using direct observations at Bab Al-Mandab. The observations were collected between June 1995 and November 1996. The heat and water fluxes information came from transport resolving current-meter arrays and temperature-salinity chain moorings. According to the study the annual mean water flux over the Red Sea was estimated to be 2.06 ± 0.22 m/yr, while the annual mean heat flux loss was approximately 11 ± 5 W/m 2 .
Ref. [1] used the NASA Modern-Era Retrospective Analysis for Research and Applications (MERRA) reanalysis with a resolution of ( 2 /3 • × 1 /2 • ) and the objectively analyzed air-sea fluxes (OAFlux) satellite reanalysis datasets from 1979 to 2013 to investigate the effects of climate modes on air-sea heat exchange in the Red Sea. They noticed different water and atmospheric characteristics from north to south Red Sea according to the seasonal monsoon wind system [1,19]. The study defined the geographical latitude of 20 • N as the separating imaginary line between the northern and southern basins (see Figure 1). They discovered that the mean net heat flux in the Northern Red Sea during the winter season was −200 W/m 2 , whereas the mean net heat gain went up to 60-80 W/m 2 in the southernmost basin. During the summer season, the net heat gain across the whole basin was greater than 50 W/m 2 , with the southern area having higher values than the northern region by 100 W/m 2 .
Using a high-resolution, nonhydrostatic, hindcast MIT general circulation model, Ref.
[12] estimated the net water flow for the Red Sea from 1993 to 2013. They concluded that the net water flux for summer showed a significant increasing trend of~1.45 m/yr, while it dropped to~1.22 m/yr during winter over the 19-yr period.
The main objectives of this study are to investigate the variability and trends in Red Sea heat and water fluxes, to estimate the relative importance of different heat flux components, to assess the spatial distribution pattern of each component using ERA5 hourly high-resolution datasets (0.25 • × 0.25 • ) from 1981 to 2020, and to determine and estimate the main factors affecting the surface net water flux over the Red Sea.

Study Area
The Red Sea is connected to the Gulf of Aden-Indian Ocean via the Strait of Bab Al-Mandab, which has a sill depth of 160 m and a maximum width of 25 km [21]. The study area extends from 10 • N to 31 • N and from 32 • E to 46 • E. In our study, we divided the Red Sea into the North Red Sea (NRS) and South Red Sea (SRS) at the geographical latitude line depicted in Figure 1. The separation was made due to differences in water and atmospheric characteristics from north to south in the Red Sea as a result of the seasonal monsoon wind system [1,19]. The wind system is strong in the SRS, and the wind direction changes seasonally in conjunction with the Indian monsoon system [10]. In particular, the reversal of the monsoons in the SRS modifies the exchange between the Red Sea and Gulf of Aden flow through the Strait of Bab Al-Mandab [5,22]. During the summer, from April to October, a northwest (NW) wind blows along the entire length of the sea, with speeds of 10 m/s approaching and frequently exceeding 15 m/s [23,24]. During the winter, the same northerly wind dominates over the northern part of the basin, while southeastern (SE) winds associated with the Indian Ocean's northeastern (NE) monsoon prevail over the SRS. The coexistence of northerly and southerly winds creates a convergence zone near the Red Sea's center, approximately at a latitude of 20 • N [24,25].

Data Sets
The atmospheric data for the computation of the surface net heat flux were taken from the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 from 1 January 1981 to 31 December 2020 [26] (https://cds.climate.copernicus.eu/cdsapp#!/ dataset/reanalysis-era5-single-levels (accessed on 20 August 2021)). The data set has a 0.25 • × 0.25 • spatial resolution with an hourly temporal step. The atmospheric field used is the hourly air temperature at a 2 m elevation (T2m), surface net longwave backscatter radiation (Q b ), surface net shortwave solar radiation (Q s ), surface sensible heat flux (Q h ), surface latent heat flux (Q e ), evaporation (E), total runoff (R) and total precipitation (P). Daily means were calculated by averaging ERA5 hourly data using the same period from 1981 to 2020.

Surface Net Heat and Water Fluxes Mean and Trend Estimations
The surface net heat flux Q T (W/m 2 ) was computed using ERA5 data sets from 1 January 1981 to 31 December 2020 according to [27] and described in [28,29] as follow: where Q b and Q s represent the radiative terms [30], while Q h and Q e represent the turbulent terms [31]. The word radiative refers to the transmission of heat from the sun to the sea (Q s ), which is the primary source of energy for the worldwide ocean. The infrared longwave radiation Q b returns some of the absorbed solar energy to the atmosphere, but a portion of this radiation is dispersed by clouds and atmospheric moisture and returns to the sea. Cloud cover has the greatest influence on the radiative terms (Q b and Q s ), whereas specific humidity, air and sea surface temperature have less influence on their fluctuation [32][33][34].
The turbulent terms (Q h and Q e ) are affected by specific humidity, difference between air and sea surface temperature and wind speed [1,5,16,35,36]. The surface net water flux is calculated from the hourly ERA5 reanalysis data as follows: The surface net water flux = (E − P − R) where E is evaporation in m/yr, P is total precipitation over the Red Sea in m/yr, and R is the total runoff in m/yr. We included in this study the very little seasonal freshwater runoff along the Eritrean coast [3,4] which is around 0.01 m/yr, to improve the precision of water flux calculations across the Red Sea. Furthermore, the spatial and temporal annual linear trends of surface net heat and water fluxes were calculated using the least squares method [37]. The original two-tailed modified Mann-Kendall test with a 95% confidence interval was used to assess the statistical significance of these trends [38,39].

Results
In this section, we show the spatiotemporal variability of heat and water fluxes in the Red Sea using ERA5 data (1981-2020).

Mean and Trend of the Surface Net Heat Flux Components
In this part, we discuss the mean spatial distribution of the heat flux components, as well as the total surface net heat flux for the Red Sea, over the study period (1981-2020), based on ERA5 hourly datasets. We illustrate the spatial and temporal trends for the surface net heat flux over the Red Sea and, furthermore, the annual and seasonal surface net heat flux over the Red Sea.

Spatial Distribution of the Mean Surface Heat Flux Components
Figure 2a-d depicts the spatial distribution pattern for the surface heat flux components Q b , Q s , Q h and Q e for the Red Sea from 1981 to 2020 based on ERA5 hourly data. Table 1 provides further information on surface net heat flux components and 2 m air temperature basin averages with two times standard deviation (± 2σ) and trend values for the NRS and SRS, as well as the whole Red Sea, from 1981 to 2020. Figure 2a shows the spatial distribution of the surface net longwave radiation (Q b ), and we can see different patterns (i.e., values) in the NRS (<−80 W/m 2 ; mean: −104 ± 5.00; trend: −0.03 ± 0.00) compared to the SRS (>−80 W/m 2 ; mean: −76 ± 7.00; trend: −0.08 ± 0.007) over the study period . The difference between the two patterns can be seen at the Lat 20 • N line, which separates the Red Sea into north and south basins. Higher absolute values of Q b appeared in the NRS, including both Suez and Aqaba Gulfs. These findings show that NRS loses more heat than the SRS. Figure 2b exhibits that the highest net shortwave radiation values (Q s ) (>235 W/m 2 ) are in the NRS (i.e.,~Lat: >20 • N; mean: 233 ± 5.00; trend: 0.05 ± 0.005), except in both Aqaba and Suez Gulfs. Q s represents the transfer of heat via solar radiation to the sea. Lower values (<235 W/m 2 ) are found in the SRS (i.e.,~Lat: <20 • N; mean: 230 ± 3.00; trend: 0.12 ± 0.001). In Figure 2c, the highest absolute surface sensible heat flux values (Q h ) (>20 W/m 2 ) are found along the coastal Red Sea and in the NRS (i.e.,~Lat: >20 • N; mean: −20 ± 0.65; trend: −0.01 ± 0.001). The lowest absolute Q h are detected in the SRS (<20 W/m 2 ; mean: −16 ± 0.35; trend: −0.2 ± 0.002). Figure 2d shows that the SRS (i.e.,~Lat: <20 • N; mean: −111 ± 2.00; trend: −123 ± 2.00) has the lowest absolute values (<120 W/m 2 ) for the surface latent heat flux (Q e ), with the exception of Bab Al-Mandab Strait. The highest absolute surface latent heat flux values (>130 W/m 2 ) are recorded in the NRS (mean: −131 ± 3.00; trend: −0.13 ± 0.010) south of the Sinai Peninsula. Figure 3a,b demonstrates the mean and spatial trend of the air temperature during the entire study period . The NRS has the greatest air temperature mean (28.60 ± 0.72 • C), while the SRS has the lowest mean (25.5 ± 0.39 • C), with a total Red Sea mean of 27 ± 0.54 • C (see Figure 3a and Table 1). The SRS, off the southern coasts of Yemen and Eritrea, has the highest air temperature trend (trend: >0.5 • C/Decade), as shown in Figure 3b. The whole Red Sea air temperature trend is 0.30 ± 0.09 • C/Decade, compared to 0.38 ± 0.11 • C/Decade for the SRS and 0.19 ± 0.06 • C/Decade) for the NRS (see Table 1).
J. Mar. Sci. Eng. 2021, 9, x FOR PEER REVIEW 6 of 17        The difference between the two patterns is clearly seen at the Lat 20 • N line for the NRS (i.e.,~Lat: >20 • N; heat loss; -Q T ) and SRS (i.e.,~Lat: <20 • N; heat gain; + Q T ). The mean surface net heat flux across the Red Sea is −1 ± 0.62 W/m 2 . There is an interesting similarity between the surface latent heat flux (Q e ) and surface net heat flux (Q T ) fields across the NRS (see Figures 2d and 4a). This similarity between Q e and Q T might explain the north-south gradient in the heat exchange, as mentioned in [1]. whereas all positive trends were determined to be significant (p < 0.05) over the whole Red Sea. Over the whole study period (1981-2020), the mean surface net heat flux trend is +0.21 ± 0.020 W/m 2 for the entire Red Sea.  The maximum significant surface net heat flux trend (p < 0.05; trend: >0.5 W/m 2 /yr) was found in the SRS, off the southern coast of Yemen and Eritrea (see Figure 4b). Negative surface net heat flow trends (<0 W/m 2 /yr) were found to be nonsignificant (p > 0.05), whereas all positive trends were determined to be significant (p < 0.05) over the whole Red Sea. Over the whole study period (1981-2020), the mean surface net heat flux trend is +0.21 ± 0.020 W/m 2 for the entire Red Sea.

Surface Net Water Flux (E-P-R)
In this section, the authors present the spatial distribution patterns for the surface water flux components (evaporation (E) and precipitation (P)), the surface net water flux for the Red Sea during the entire study period (1981-2020) and the temporal and seasonal distribution for the surface net water flux.

Surface Net Water Flux (E-P-R)
In this section, the authors present the spatial distribution patterns for the surface water flux components (evaporation (E) and precipitation (P)), the surface net water flux for the Red Sea during the entire study period (1981-2020) and the temporal and seasonal distribution for the surface net water flux.

Spatial Distribution for the Evaporation, Precipitation and Surface Net Water Flux
Figure 6a-c describes the spatial distribution patterns for the surface water flux components E, P and the surface net water flux for the Red Sea during the entire study period (1981-2020). Table 2 shows the surface net water flux component basin averages (E, P and R) and the surface net water flux for the NRS and SRS, as well as the whole Red Sea, during the aforementioned time period. Figure 6a demonstrates that the NRS (i.e.,~Lat: >20 • N; mean: 1.68 ± 0.33 m/yr) has the largest evaporation rate (>1.6 m/yr). The SRS (i.e.,~Lat: <20 • N; mean: 1.40 ± 0.22 m/yr) has the lowest values (<1.5 m/yr). The spatial mean evaporation over the whole Red Sea is 1.54 ± 0.26 m/yr, which is consistent with the findings of [5]. Figure 6b shows that the SRS (i.e.,~Lat: <20 • N; mean: 0.12 ± 0.03) has more precipitation than the NRS (i.e.,~Lat: >20 • N; mean: 0.04 ± 0.01). The average annual precipitation over the entire Red Sea is 0.08 ± 0.02 m/yr. Figure 6c indicates that the annual mean water loss to the atmosphere over the Red Sea is greater in the NRS (i.e.,~Lat: >20 • N; mean: 1.64 ± 0.31) than the SRS (i.e.,~Lat: <20 • N; mean: 1.27 ± 0.18). This might be due to the high precipitation rate and the runoff from the Eritrean coast (0.01 m/yr) in the SRS. Figure 6c shows the average surface net water flux of 1.46 ± 0.23 m/yr for the entire Red Sea over the study period 1981-2020. It is obvious that the rate of evaporation surpasses the rate of precipitation over the Red Sea. The evaporation and surface net water flux fields are quite similar. Overall, evaporation rate is a significant factor in estimating surface net water flux over the Red Sea. . The spatial mean evaporation over the whole Red Sea is 1.54 ± 0.26 m/yr, which is consistent with the findings of [5]. Figure 6b shows that the SRS (i.e., ~Lat: <20° N; mean: 0.12 ± 0.03) has more precipitation than the NRS (i.e., ~Lat: >20° N; mean: 0.04 ± 0.01). The average annual precipitation over the entire Red Sea is 0.08 ± 0.02 m/yr. Figure 6c indicates that the annual mean water loss to the atmosphere over the Red Sea is greater in the NRS (i.e., ~Lat: >20° N; mean: 1.64 ± 0.31) than the SRS (i.e., ~Lat: <20° N; mean: 1.27 ± 0.18). This might be due to the high precipitation rate and the runoff from the Eritrean coast (0.01 m/yr) in the SRS. Figure 6c shows the average surface net water flux of 1.46 ± 0.23 m/yr for the entire Red Sea over the study period 1981-2020. It is obvious that the rate of evaporation surpasses the rate of precipitation over the Red Sea. The evaporation and surface net water flux fields are quite similar. Overall, evaporation rate is a significant factor in estimating surface net water flux over the Red Sea.

Temporal Evolution of the Surface Net Water Flux (E-P-R)
To understand the seasonal behavior of surface water flux, we first focus on the evaporation and precipitation seasonal cycles (main water flux components). Figure 7a depicts a time series of the surface seasonal averaged daily mean evaporation rate (m/yr) against a time series of seasonal daily mean average precipitation. This double-scale graphic shows that the evaporation rate of the Red Sea is the most important factor influencing the surface net water flux. The highest evaporation rate (~2 m/yr) is recorded during the winter season (i.e., December and January), and this finding is similar to [17]. The minimum evaporation rate (~1 m/yr) is detected during spring and early summer (i.e., May and June). The precipitation rate is fluctuating between ~0 m/yr (summer: June and July) and 0.25 m/yr (late autumn: mid-November). Figure 7b demonstrates the seasonal daily climatology (1981-2020) mean cycle of the surface net water flux (m/yr). The surface water flux seasonal cycle is following exactly the seasonal cycle of the evaporation rate. The highest (1.9 m/yr) surface net water flux is observed during winter (i.e., December). The lowest (<+1 m/yr) surface net water flux was found during late spring and early summer (i.e., May and June), which agrees with [1,5].

Temporal Evolution of the Surface Net Water Flux (E-P-R)
To understand the seasonal behavior of surface water flux, we first focus on the evaporation and precipitation seasonal cycles (main water flux components). Figure 7a depicts a time series of the surface seasonal averaged daily mean evaporation rate (m/yr) against a time series of seasonal daily mean average precipitation. This double-scale graphic shows that the evaporation rate of the Red Sea is the most important factor influencing the surface net water flux. The highest evaporation rate (~2 m/yr) is recorded during the winter season (i.e., December and January), and this finding is similar to [17]. The minimum evaporation rate (~1 m/yr) is detected during spring and early summer (i.e., May and June). The precipitation rate is fluctuating between~0 m/yr (summer: June and July) and 0.25 m/yr (late autumn: mid-November). Figure 7b demonstrates the seasonal daily climatology (1981-2020) mean cycle of the surface net water flux (m/yr). The surface water flux seasonal cycle is following exactly the seasonal cycle of the evaporation rate. The highest (1.9 m/yr) surface net water flux is observed during winter (i.e., December). The lowest (<+1 m/yr) surface net water flux was found during late spring and early summer (i.e., May and June), which agrees with [1,5].

Discussion
The spatial variability of all heat flux components (Qb, Qs, Qh and Qe) differs between the NRS and SRS (i.e., higher absolute values in NRS and lower values in SRS), as shown in Figure 2a-d.
Higher absolute values of Qb (<−80 W/m 2 ) were found in the NRS (mean: −104 ± 5.00; trend: −0.03 ± 0.005), including both Suez and Aqaba Gulfs (Figure 2a). These findings are consistent with [1,5,10], which demonstrate that the NRS loses more heat than the SRS. Figure 2b shows that the highest Qs values (>235 W/m 2 ) were found in the NRS (mean: 233 ± 5.00; trend: 0.05 ± 0.005), except in the Gulfs of Aqaba and Suez, in agreement with [1,16,20]. The highest absolute Qh flux values (5-10 W/m 2 ; Figure 2c) were found along the Red Sea coast and in the NRS (i.e., ~Lat: >20° N; mean: −20 ± 0.65; trend: −0.01 ± 0.001). This is due to the large disparity in air and sea surface temperatures, as well as the high relative humidity in these locations [16], which are regarded as the primary factors influencing Qh [35,36]. Figure 2d shows that the lowest absolute values (<120 W/m 2 ) for Qe were found in the SRS (i.e., ~Lat: <20° N; mean: −111 ± 2.00; trend: −123 ± 2.00), except at Bab Al-Mandab. This might be due persistence of the strong monsoon wind system in the SRS [10]. As previously stated, the small absolute Qe value (~170 W/m 2 ) at Bab Al-Mandab Strait may be due to the significant wind variability in this specific area [6]. The highest absolute Qe values (>130 W/m 2 ) were noticed in the NRS (mean: −131 ± 3.00; trend: −0.13 ± 0.010), particularly south of the Sinai Peninsula, in agreement with [5]. According to Table  1 and Figure 3, we found that the NRS has the highest air temperature mean (28.60 ± 0.72 °C), while the SRS has the lowest mean (25.5 ± 0.39 °C), with a total Red Sea mean of 27 ± 0.54 °C. The highest air temperature trend (trend: >0.5 °C/Decade) was found off the southern coasts of Yemen and Eritrea. The whole Red Sea air temperature trend was 0.30 ± 0.09 °C/Decade, compared to 0.38 ± 0.11 °C/Decade for the SRS and 0.19 ± 0.06 °C/Decade) for the NRS (see Table 1). Figure 4a shows a surface net heat flux gain (>50 W/m 2 ) was noticed in the SRS (mean: +27 ± 0.48; trend: 0.20 ± 0.021) in front of Yemen and Eritrea's coastal regions, in accordance with [1,5,10,16]. This may be because of reduced shortwave radiation (Qs) values (<235 W/m 2 ) over the southern part of the basin, as previously mentioned by [1]. The NRS experienced a loss in surface net heat flux (QT) (−55 W/m 2 ), particularly south of the Sinai Peninsula and the Gulf of Suez (mean: −22 ± 0.82; trend: −0.12 ± 0.052), in agreement with [17]. We discovered an interesting resemblance between the surface latent heat flux (Qe)

Discussion
The spatial variability of all heat flux components (Q b , Q s , Q h and Q e ) differs between the NRS and SRS (i.e., higher absolute values in NRS and lower values in SRS), as shown in Figure 2a-d.
Higher absolute values of Q b (<−80 W/m 2 ) were found in the NRS (mean: −104 ± 5.00; trend: −0.03 ± 0.005), including both Suez and Aqaba Gulfs (Figure 2a). These findings are consistent with [1,5,10], which demonstrate that the NRS loses more heat than the SRS. Figure 2b shows that the highest Q s values (>235 W/m 2 ) were found in the NRS (mean: 233 ± 5.00; trend: 0.05 ± 0.005), except in the Gulfs of Aqaba and Suez, in agreement with [1,16,20]. The highest absolute Q h flux values (5-10 W/m 2 ; Figure 2c) were found along the Red Sea coast and in the NRS (i.e.,~Lat: >20 • N; mean: −20 ± 0.65; trend: −0.01 ± 0.001). This is due to the large disparity in air and sea surface temperatures, as well as the high relative humidity in these locations [16], which are regarded as the primary factors influencing Q h [35,36]. Figure 2d shows that the lowest absolute values (<120 W/m 2 ) for Q e were found in the SRS (i.e.,~Lat: <20 • N; mean: −111 ± 2.00; trend: −123 ± 2.00), except at Bab Al-Mandab. This might be due persistence of the strong monsoon wind system in the SRS [10]. As previously stated, the small absolute Q e value (~170 W/m 2 ) at Bab Al-Mandab Strait may be due to the significant wind variability in this specific area [6]. The highest absolute Q e values (>130 W/m 2 ) were noticed in the NRS (mean: −131 ± 3.00; trend: −0.13 ± 0.010), particularly south of the Sinai Peninsula, in agreement with [5]. According to Table 1 and Figure 3, we found that the NRS has the highest air temperature mean (28.60 ± 0.72 • C), while the SRS has the lowest mean (25.5 ± 0.39 • C), with a total Red Sea mean of 27 ± 0.54 • C. The highest air temperature trend (trend: >0.5 • C/Decade) was found off the southern coasts of Yemen and Eritrea. The whole Red Sea air temperature trend was 0.30 ± 0.09 • C/Decade, compared to 0.38 ± 0.11 • C/Decade for the SRS and 0.19 ± 0.06 • C/Decade) for the NRS (see Table 1). Figure 4a shows a surface net heat flux gain (>50 W/m 2 ) was noticed in the SRS (mean: +27 ± 0.48; trend: 0.20 ± 0.021) in front of Yemen and Eritrea's coastal regions, in accordance with [1,5,10,16]. This may be because of reduced shortwave radiation (Q s ) values (<235 W/m 2 ) over the southern part of the basin, as previously mentioned by [1]. The NRS experienced a loss in surface net heat flux (Q T ) (−55 W/m 2 ), particularly south of the Sinai Peninsula and the Gulf of Suez (mean: −22 ± 0.82; trend: −0.12 ± 0.052), in agreement with [17]. We discovered an interesting resemblance between the surface latent heat flux (Q e ) and surface net heat flux (Q T ) fields over the NRS. This similarity between Q e and Q T plays a crucial role in defining the north-south gradient in the heat exchange, as mentioned in [1]. The surface net heat flux fields showed a distinct difference between the NRS (heat loss) and SRS (heat gain) defined by the Lat 20 • N line, as shown in Figure 4a. We found that the mean surface net heat flux across the Res Sea was −1 ± 0.62 W/m 2 . This finding agrees with [9,10,40]. Figure 4b demonstrates that the maximum significant surface net heat flux trend (i.e., p < 0.05; trend: >0.5 W/m 2 /yr) was observed in the SRS, off the southern coast of Yemen and Eritrea. All positive trends were found to be significant (p < 0.05) over the whole Red Sea (Figure 4b). Negative surface net heat flux trends (i.e., <0 W/m 2 /yr) were found to be nonsignificant (p > 0.05). The mean surface net heat flux trend was +0.2 ± 0.020 W/m 2 for the entire Red Sea. The authors described the annual and the seasonal surface net heat flux (Figure 5a,b) for the Red Sea over the whole study period (1981-2020). We discovered a positive surface net heat flux temporal trend of +0.18 ± 0.057 (W/m 2 )/yr over the study period (i.e., 1981-2020) ( Figure 5a). As recently highlighted by [41], the rising trend of the surface net heat flux could be influenced by global warming. Additionally, according to a study by [42], the interannual variations in the heat content of the 0-100 m layer in the Red Sea are largely caused by changes in air temperature (i.e., global warming). The highest surface net heat flux was >80 W/m 2 during spring and early summer (i.e., from May to June). The lowest surface net heat flux was <−80 W/m 2 during winter (i.e., December and January), which is consistent with [1,5].
The evaporation rates were~Lat: >20 • N and mean: 1.68 ± 0.33 m/yr for the NRS,~Lat: <20 • N and mean: 1.40 ± 0.22 m/yr for the SRS and for the whole Red Sea 1.54 ± 0.26 m/yr, as shown in Figure 6a. The Red Sea evaporation rate is considered to be the greatest in the global ocean [10]. The high evaporation rate over the Red Sea might be due to the strong monsoon winds and high air temperature, as previously mentioned by [1,5,6,10,17,43]. Precipitation rates were found to be higher in the SRS (i.e.,~Lat: <20 • N; mean: 0.12 ± 0.03) than in the NRS (i.e.,~Lat: >20 • N; mean: 0.04 ± 0.01), and the average annual precipitation throughout the whole Red Sea is 0.08 ± 0.02 m/yr, as shown in Figure 6b. From the results, we deduced that the evaporation rate exceeds the precipitation rate over the entire Red Sea, in agreement with [5,6]. The annual mean water loss to the atmosphere over the Red Sea is greater in the NRS (i.e.,~Lat: >20 • N; mean: 1.64 ± 0.31) than the SRS (i.e.,~Lat: <20 • N; mean: 1.27 ± 0.18) (Figure 6c). The annual mean water loss to the atmosphere over the Red Sea was 1.46 ± 0.23 m/yr. The authors discovered that the evaporation and surface net water flux fields are quite similar (see Figure 6a,c). In summary, evaporation rate is a significant factor in estimating surface net water flux over the entire Red Sea. We created a double-scale graph (comparing seasonal daily climatological mean evaporation and precipitation rates) to demonstrate that evaporation is the primary factor influencing surface net water flux over the Red Sea. The authors found that the surface water flux seasonal cycle exactly followed the seasonal cycle of the evaporation rate. The maximum (1.9 m/yr) surface net water flux is during winter (i.e., December), while the minimum (<+1 m/yr) surface net water flux is during late spring and early summer (i.e., May and June) (Figure 7b). The surface net water flux peaks observed in winter are due to southeastern (SE) winds associated with the Indian Ocean's northeastern (NE) monsoon winds, which blow up with speeds of 10 m/s approaching and frequently exceeding 15 m/s over the whole length of the Red Sea [23,24]. The highest surface net water flux (+2.1 m/yr) was detected during 2020, while the lowest value (+1.3 m/yr) was observed during 1985 (see Figure 8).

Conclusions
This study provides the variability of heat and water fluxes over the Red Sea using ERA5 reanalysis dataset over 40 years (1981-2020). The authors determined the relative significance of various heat and water flux components and illustrated their geographical distribution patterns.
The highest absolute values of surface sensible heat flux (Q h ) were found in the NRS. This might be due to the large difference between air and sea surface temperatures, as well as the high relative humidity in this region. The SRS had the lowest absolute values for surface latent heat flux (Q e ), with the exception of the Bab Al-Mandab Strait. This might be due to the SRS's strong monsoon wind system, which fluctuates in conjunction with the Indian monsoon system. The spatial patterns of the latent heat flux (Q e ) and net heat flux (Q T ) coincided in the NRS. This similarity between Q e and Q T could explain the north-south gradient in heat exchange.
In general, the spatial pattern in the NRS surface net heat flux (Q T ) was negative (i.e., heat loss) and positive in the SRS (i.e., heat gain), with average values of −22 ± 0.82 W/m 2 and +27 ± 0.48 W/m 2 , respectively. We found linear trends for surface net heat flux in the NRS of −0.12 ± 0.052 (W/m 2 )/yr and +0.20 ± 0.021 (W/m 2 )/yr in the SRS, with an average trend value of +0.21 ± 0.020 (W/m 2 )/yr for the entire Red Sea. According to the temporal variability of the surface net heat flux, the Red Sea was in a different state before and after 2001. (i.e., there is an unaccelerated trend from 1981 to 2001, followed by a higher trend from 2002 to 2020).
The evaporation rate was a significant factor in estimating surface net water flow over the entire Red Sea. Over the whole study period, the results deduced that the temporal annual surface net water flux has positive values (i.e., freshwater loss to the atmosphere). The surface net water flux peak observed in winter could be attributed to the NE monsoon wind, which increases evaporation and has a direct impact on the surface net water flux along the entire length of the Red Sea.
The highest surface net water flux (+2.1 m/yr) was detected during 2020, while the lowest value (+1.3 m/yr) was observed during 1985.