Seasonal variability of hydrography off the east coast of Qatar, central Arabian Gulf

Seasonal variability of hydrography along a nearshore-offshore transect in the eastern part of the EEZ of Qatar has been analyzed using in situ measurements carried out during 5 different months. The study reveals distinct features in temperature, salinity, density, dissolved oxygen (DO), and chlorophyll fluorescence. The sea surface temperature (SST) varies from nearshore to offshore with a positive gradient during January, June, and August, of the order of 1.0–2.0 °C, and with a negative gradient during October and April, of the order of 1.0 °C. Thermal stratification began in June, reached a well-defined thermocline with a vertical difference in temperature of about 10 °C during August, and started to de-stratify during October. The low salinity and low-density inflow of IOSW is evident in the surface layer of the middle of the transect during August, which has enriched the DO in the surface layer up to 20 m depth, while hypoxia prevailed below 50 m depth. The lowest chlorophyll fluorescence was measured in April 2013 (~ 1.0 μg l−1), moderate in June 2013 (~ 1.5 μg l−1), and relatively high in August 2013 (~ 6.0 μg l−1).


Introduction
The Arabian/Persian Gulf (hereafter the "Gulf") is a semienclosed, shallow marginal sea connected to the Arabian Sea through the Strait of Hormuz (Fig. 1). It is a hypersaline and warmest marine area in the world, where evaporation exceeds precipitation (Al-Ansari et al. 2015). The evaporation rate estimated in the Gulf is 1.4-2.1 m year −1 (Xue and Eltahir, 2015), while the precipitation and river run-off are in the range of 0.03-0.11 and 0.1-0.46 m year −1 , respectively (Al Senafi and Anis 2015;Al-Mazroui et al. 2012).
Recent studies indicate that harmful environmental impacts of anthropogenic origin in the Gulf increase in frequency (Al-Ansari et al. 2015;Feary et al. 2011;Sheppard et al. 2010;Al-Ansi et al. 2002). Besides the increasing natural and anthropogenic stresses that affect the nearshore regions, the water mass exchanges and circulation greatly help to maintain stable conditions in most of the Gulf. For instance, Ibrahim et al. (2020) identified that the Gulf basin salinity is relatively stable and resilient to the existing brine discharges from the desalination plants around the Gulf. However, the impacts of brine charges on the marine ecosystem are crucial Responsible Editor: Zhihua Zhang Abdel-Moati has passed away during the preparation of this study.
in the nearshore waters of the Gulf (Hosseini et al. 2021;Joydas et al. 2015).
The basin-scale circulations, seawater properties, and water masses were investigated in the Gulf (Al-Shehhi et al. 2021;Campos et al. 2020;Azizpour et al. 2014;Yao and Johns, 2010a, b;Swift and Bower, 2003;Reynolds, 1993). The exchange of water masses between the Gulf and the Arabian Sea maintains a stable salt budget in the Gulf (Johns et al. 2003). The high-saline and welloxygenated Persian Gulf Water (PGW) enters the Arabian Sea through the Strait of Hormuz, ventilating the world's thickest oxygen minimum zone (OMZ), and flows further down into the Indian Ocean (Lachkar et al. 2019). The outflow water during the first half of the year mostly originates from the southern Gulf coast, while that during the second half is mainly transported from the northwestern Gulf (Lorenz et al. 2020). Traces of PGW have been identified in the northwestern Indian Ocean and Bay of Bengal, indicating that the saltier PWG is strong enough to impact the biological activity in these regions (L'Hégaret et al. 2021;Sheehan et al. 2020). Whereas the low-saline Indian Ocean Surface Water (IOSW) flows into the Gulf through the Strait of Hormuz, which is more evident during summer (Kampf and Sadrinasab, 2006). The distribution of water masses and the formation of a strong pycnocline during summer reflect the combined influence of horizontal and vertical spreading of low-salinity surface inflow and excessive atmospheric heating (Al-Ansari et al. Fig. 1 a Geographical map and C-MAP bathymetry of the Arabian Gulf and b EEZ of Qatar (dotted lines indicate the boundary). Red dots numbering 1 to 6 are the sampling stations. Black dots numbered in 5 digits are the Coriolis drifting buoy positions during the study period 2015). The presence of IOSW has been identified in the deep waters of the Exclusive Economic Zone (EEZ) of Qatar (Rakib et al. 2021;Elobaid et al. 2022). Horizontal and vertical gradients in temperature, salinity, and density are evident in the Gulf basin during summer, caused by the extreme surface heating (Pous et al. 2015;Kampf and Sadrinasab, 2006). However, vertical gradients are diminished during the winter due to well-mixing by the winter shamal winds (Al Azhar et al. 2016).
The biogeochemistry and ecosystems of the Gulf are experiencing massive changes in recent years, caused by a reduction in river runoff in the north and an increase in nutrient loading from the urban areas (Al-Yamani and Naqvi, 2019). Increasing sea surface temperature (SST) in response to global warming is another factor that influences the physical and biogeochemical distribution in the Gulf (Noori et al. 2019;Strong et al. 2011). The phytoplankton biomass (chlorophyll a) is reported to decrease from approximately 94 µg l −1 in the Shatt Al-Arab estuary in the north to about 2-4 µg l −1 around Qatar and off the Abu Dhabi coast (Al-Naimi, et al. 2017;Mezhoud et al. 2016;Quigg et al. 2013). Low biomass and productivity in the Gulf are attributed to its shallow nature, moderate to strong wind and current regimes, and well-mixed water column driven by geomorphology and physical dynamics (Al-Ansari et al. 2015;Jones, 1985). Due to excess respiration and microbial decomposition in the subsurface (depth > 45 m), dissolved oxygen (DO) concentration in the deepest part of the central Gulf and on the Iranian side of the Strait of Hormuz decreases throughout the summer and reaches hypoxic levels by early autumn (Saleh et al. 2021;Al-Ansari et al. 2015). Short episodes of summer hypoxia were identified at Dhabiya Reef in the southern Gulf (De Verneil et al. 2021) and Sulaibikhat Bay in Kuwait (Al-Mutairi et al. 2014).
The physical and biogeochemical characteristics of the Gulf are highly complex due to bathymetric-anthropogenicclimatic interactions. While climatic features are regional in nature, bathymetric and anthropogenic influences are mostly localized. Although basin-scale features have been explored to some extent, focused studies pertaining to the localized physical-biogeochemical interactions and features are very limited, especially in the EEZ of Qatar. The summer hypoxia explored in the central Gulf (Al-Ansari et al. 2015) and the studies carried out on the hydrography and biogeochemistry (Elobaid et al. 2022;Rakib et al. 2021) point out the necessity of having further investigations in the EEZ of Qatar. The main constraint on this is the lack of measurements in this region. The present study utilizes seasonal datasets measured on a transect off the east coast of Qatar and the available Coriolis Ocean Data to investigate the distribution and seasonal variability of temperature, salinity, density, DO, and chlorophyll fluorescence in the water column.

Area of study
The Gulf, located between the latitudes 24.0°-30.5° N and the longitudes 48.0°-56.5° E, is an extension of the Arabian Sea (Indian Ocean) and covers an ocean area of about 241,000 km 2 (Fig. 1). It has a length of 990 km, a width ranging from 56 to 340 km, a maximum depth of approximately 110 m, and an average depth of 36 m (Sadrinasab and Kämpf, 2004). The deeper areas of the Gulf are located along the Iranian coast, while broad shallow regions with depths of less than 35 m are located along the Arabian coast. Winds of various scales and directions dominate the Gulf, with the northwesterly shamal winds, northeasterly/easterly nashi winds, and southeasterly/ southerly kaus winds having a significant influence on sea surface dynamics (Aboobacker et al. 2021a,b c;Veerasingam et al. 2020a, b;Perrone 1979).
The EEZ of Qatar is in the central Gulf, covering 14% of the total area of the basin (Al-Qaradawi et al. 2015). The EEZ is relatively larger and deeper off the east and north coasts of Qatar, bordering Iranian and UAE waters (Fig. 1b). The west part of the EEZ is very narrow and shallow and border on Bahrain and Saudi Arabia.

Data and methods of analysis
This study uses the hydrographic and biogeochemical data measured at 6 stations along a 105-km transect off the east coast of Qatar in 5 cruises during 2012-2013 ( Fig. 1; Table 1) utilizing the R/V Janan (owned and operated by Qatar University). Samples were collected during October 2012, January 2013, April 2013, June 2013, and August 2013. The transect, which ends at the Qatari marine border, traverses more than half of the linear distance between Doha and the Iranian coast, thus providing a reasonable representation of hydrographic distributions across this relatively wider part of the Gulf. Although limited to a single transect extending from the coast to the center of the Gulf, this dataset is the first of its kind, and it has never been utilized before to investigate the seasonal variability of water column structure. The profiling of temperature and salinity for every 1.0 m depth was carried out using the SBE-25 plus CTD system in October 2012 and using the SBE-911plus CTD for the remaining 4 cruises conducted in 2013. In October 2012 and January 2013, discrete sampling was carried out for DO using SBE-32 Carousel fitted with 12-l Niskin bottles. DO samples were analyzed for concentration and saturation onboard within a few hours of collection using the titrimetric method (Winkler 1888). This method was used earlier in the EEZ of Qatar (Rakib et al. 2021;Al-Ansari et al. 2015) and other parts of the Gulf (Naqvi 2021;Saleh et al. 2021) However, in April 2013, the CTD system was fitted with new sensors for in situ DO (SBE-43) and chlorophyll fluorescence (WET Labs ECO-FL-NTU-RT-D) measurements. These sensors were well calibrated (Janzen et al. 2008). Therefore, we relied on the pre-established accuracy of the sensors, since a comparison between the sensor and titration data has not been carried out. Recently, a luminescent DO sensor (Hach IntelliCal LDO101) data in the Iranian waters of the Gulf was verified well with the titration results (Saleh et al. 2021). Ocean Data View (ODV) was used to analyze and visualize bin-averaged downcast CTD data and laboratory DO data (Schlitzer 2015).
The in situ data in the Gulf is sparse, especially in the EEZ of Qatar. We have explored the global resources of in situ data within the study region and were able to identify a few data points of Coriolis Ocean products ( Table 2). The Coriolis Ocean Dataset provides operational oceanography data around the world oceans (Tanguy et al. 2019). This dataset contains observations assimilated from different sources such as autonomous platforms (Argo profilers, fixed moorings, gliders, drifters, sea mammals) and vessels (CTDs, XBTs, ferryboxes).
Surface wind data in the nearshore and offshore regions along the transect has been obtained from the European Centre for Medium-Range Weather Forecasts (ECMWF) product, ERA5. ERA5 provides hourly wind velocities globally on a 30-km grid resolution, in which the scatterometer and in situ wind data are assimilated to improve the accuracy of predictions (Hersbach et al. 2020). In the present study, the ERA5 winds have been analyzed to characterize the variability of hydrographic parameters.

Comparison with Coriolis Ocean Dataset
Close to the EEZ, Copernicus data has limited drifting buoys (Fig. 1b) during Oct 2012-Sep 2013, having only sea surface temperature (SST) as given in Table 2. These buoys are away from the transect, indicating that there is a spatial data gap in the central Gulf. In addition, the CTD data also do not coincide with that of the Coriolis drifting buoys, temporally. Therefore, a quantitative comparison between the two datasets is not feasible in the context of the spatial and temporal differences. Nonetheless, we have tabulated the SST values from both the sources with the nearest available time slots to derive some inferences (Table 2).
In the EEZ of Qatar, there is a decrease in SST of the order of 2.9-4.5 °C in the EEZ of Qatar in the period between October and November. This is a transition period from summer to winter in which the shamal wind cools the sea surface faster than the inter-seasonal shifts in temperature. In addition, spatial variability exists between the west and east parts of the EEZ of Qatar (Elobaid et al. 2022). Between December and January, there is no significant temporal variation in SST. However, its spatial variability has more relevance. For instance, the drifting buoy 23,519 is the nearest Coriolis platform of Stn. 3, where the SST is only 0.3 °C higher than that of Stn.3. Stn. 6 is further offshore, where the SST is 0.9° higher than that at Stn. 3. To the west of the cruise transect (outside the EEZ), the drifting buoy data (23,754,23,749,23,756) has lower SST compared to Stn. 3 and Stn.6, with a decrement of the order of 2.0-3.5 °C. This spatial variation in SST is a major concern, especially when several water masses are active in the central Gulf. It has been discovered that the relatively warm water mass from the Sea of Oman during winter has a significant impact in Qatar's EEZ (Rakib et al. 2021). Between March and April, the increment in SST is evident due to an increase in solar radiation. During the end of August and September, within a short time span, the CTD and buoy SSTs show small spatial variations. The co-analysis of Coriolis data with the transects indicates the necessity of having more observations in the EEZ of Qatar to improve the understanding of the physical variables.

Seasonal variability in hydrography and water mass distribution
The temperature, salinity, and density (σ t ) along the track connecting the sampling stations show systematic variations (Fig. 2). Horizontal and vertical variability is evident across the transect. However, they are not consistent across the seasons. The SST is relatively higher in the nearshore region during October 2012  August 2013, the surface layer in the middle of the transect showed a lower salinity compared to the nearshore and offshore stations. This relatively fresh and low-density water inflow is caused by the IOSW. Its impact has been stretched up to the nearshore region, where the salinity has dropped to 2.0 compared to the nearshore salinity during June 2013. Elobaid et al. (2022) reported a similar nearshore-offshore gradient in SSS during late summer in a different transect beginning in the southern part of Qatar's east coast and ending in a region adjacent to Stn. 6. They reported a nearshore SSS of about 44.3, while the nearshore SSS we measured in the central nearshore region is about 42.0. Both of these locations have a distance of about 80 km apart. Thus, the dynamics that lead to the difference in nearshore salinity are quite different. The wind mixing and the proximity to the deep waters of the central Gulf are relatively higher on the central east coast compared to the southern east coast, which enables a relatively faster mixing of SSS with the adjacent regions. It is also linked with the influence of brine discharges from the desalination plants (Elobaid et al. 2022;Rakib et al. 2021). Hypersaline and relatively colder water formed in the nearshore region cascades downslope from January 2012 through April 2013. During summer and fall, the formation of high-saline deep-water masses and the intrusion of low-salinity water from the Sea of Oman into the surface layer appear as counter-acting processes. The low-saline surface layer appears only in the offshore area during summer (June 2013), which has a thickness of about 15 m. During late summer (August 2013), the layer moves shoreward and deepens down to about 20 m due to a downwelling cell, which is associated with an anti-cyclonic eddy formed in the central Gulf between two cyclonic eddies (Thoppil and Hogan 2010). During the fall (October 2012), only a modified remnant of the intrusion layer is observed offshore.
The temperature variability of the whole surface layer in the study area is minimal down to about 20-m depth in all seasons, excluding late summer when a strong thermocline is formed offshore (Fig. 3). The vertical temperature distribution in the offshore is mainly controlled by the seasonal mixing-stratification cycle. The water column is vertically homogenous during January 2013. Atmospheric cooling and wind-induced mixing during the winter collectively produce a homogenous water column in the Gulf (Thoppil and Hogan 2010;Al Azhar et al. 2016). The daily average wind speed during January 2013 goes up to 11.7 m/s in the offshore (station 6), while the monthly average wind speed is 4.6 m/s (Fig. 4). The seasonal thermocline forms during early summer (June 2013, 22 m) and strengthens throughout the season (August 2013, 15 m). From April through August, the water below 30 m depth is cooler than that during January, showing the effect of thermal isolation on the summer bottom layer. The pycnocline is stronger during late summer (August 2013), and this corresponds to a density difference of about 3.5 kg m −3 between 12 and 25 m water depths (Fig. 3c). As surface warming continues through October, the seasonal thermocline and halocline weaken and deepen (Fig. 3a, b). This is aligned with a weakening of low-salinity surface inflow. As the atmospheric heating is sustained into the fall season (October 2012), atmospheric warming extends to deeper depths; thus, the thermocline deepens below 40 m. This is consistent with the thermocline obtained in the deep water of the central transect during late September in a previous investigation (Elobaid et al. 2022;Rakib et al. 2021).
The potential temperature-salinity (Ɵ/S) diagram (Fig. 5) shows a clear seasonal variability in water mass characteristics. Salinity (37.99-41.65) and density (22.49-28.66 kg m −3 ) ranges are narrower than temperature (21.04-34.63 °C). It is inferred that salinity is the Fig. 3 Vertical profiles of a temperature (ºC), b salinity, and c density (σ t ) (kg m. −3 ) for all cruises at the deepest sampling station (station 6) main contributor to the variability in density. However, the temperature does have robust links during the transition periods. Vigorous vertical mixing due to the combined effects of decreased insolation, high wind speeds, and enhanced convection confines winter water masses to a small range of densities (σ t : 27.30-28.66 kg m −3 ). The range of densities is relatively wider in summer (σ t : 22.49-28.37 kg m −3 ), which is driven by the warming and freshening of the surface layer. Similar seasonal variations are identified in the earlier studies in the central Gulf (Al-Ansari 2006; Hassan and Mahmoud 1986;Rakib et al. 2021). The extent of low salinity (< 40) surface waters of Sea of Oman origin is particularly pronounced during the summer and to a lesser extent during the winter, consistent with earlier observations (Al-Shehhi et al. 2021;Campos et al. 2020). Offshore, deep waters below 30 m in the EEZ of Qatar remain colder and saltier throughout the summer, as vertical convection is hindered by seasonal pycnocline. By early fall (October), stratification has weakened considerably, rendering deeper waters warmer and relatively fresher than in summer.

Variability of DO and chlorophyll fluorescence
Distinct seasonal variability is observed in the DO concentration along the transect (Fig. 6). In October 2012, the surface layer has a DO of above 125 µmol l −1 , while the bottom layer (below 50 m) has a lesser DO, on the border of hypoxia. On the other hand, hypoxia is evident below 40-m depth during August 2013. However, it is not evident during June 2013. Earlier, Al-Ansari et al. (2015) identified hypoxic conditions in the deeper water of the EEZ of Qatar during September. These indicate that the hypoxia in the EEZ of Qatar is intra-seasonal and occurs in the late summer, explicitly during August-September. The stratification of offshore waters by excessive surface heating (Elobaid et al. 2022;Mosaddad et al. 2012), dense and relatively cold deep water flow originating from the northern Gulf (Al Azhar et al. 2016), and relatively calm wind conditions (Fig. 4) together catalyzes the biogeochemical interactions leading to hypoxia. Weakening of stratification during October 2012 (Fig. 2) enabled the deep waters to make a slow recovery from the hypoxic conditions, and strong mixing during winter enriched the DO to a healthy level through the spring. The overall maxima in DO (> 201 µmol l −1 at about 100% saturation) are observed during April 2013 in the offshore surface layer, during which the water column in the offshore is well mixed. Although hypoxia exists in the deeper depths of the offshore region during August 2013, the surface layer in the middle of the transect is augmented with sufficient DO, of about 175 µmol l −1 . This is aligned with a relatively fresh and low-density surface water inflow (Fig. 2), caused by the IOSW.
The chlorophyll fluorescence measured during three cruises indicates that the overall phytoplankton biomass is low to moderate during April 2013 and June 2013 and relatively high during August 2013 (Fig. 7). Laterally, the surface layer up to 10 m depth has no significant variability between the stations and the seasons, where the chlorophyll fluorescence is below 1.0 μg l −1 . The subsurface layer (between 10 and 40 m depths) shows seasonal variation, with the lowest chlorophyll fluorescence (~ 1.0 mg l −1 ) in April 2013, moderate in June 2013 (~ 1.5 mg l −1 ), and relatively high in August 2013 (~ 6.0 mg l −1 ). The highest chlorophyll fluorescence was found in Stn. 2 around 20 m, followed by ~ 3.0 μg l −1 chlorophyll fluorescence at 20-m depth in station 6. Patches of moderate chlorophyll fluorescence can also be seen in the subsurface layer between the two stations mentioned above in August 2013. Rakib et al. (2021) identified a relatively higher chlorophyll concentration in the surface and sub-surface layers at an offshore station (adjacent to Stn. 6) during September 2014. However, early summer data measured during May 2016 (Elobaid et al. 2022) adjacent to Stn. 6 indicates a lower chlorophyll concentration, which lies in between the values obtained during April 2013 and June 2013 in the present study. This indicates that there is a gradual increase in chlorophyll fluorescence from spring to late summer. The highest chlorophyll fluorescence measured is aligned with the highest DO at the same depth during August 2013 (Fig. 6). This ensures sufficient biomass production during the late summer in the subsurface layer.
The chlorophyll fluorescence data during winter is very limited in the EEZ of Qatar. Rakib et al. (2021) measured a slightly higher concentration (~ 5.6 μg l −1 ) in the offshore region of Qatar's EEZ in January 2015. Together with our observations, we can infer that the winter and late summer have relatively higher chlorophyll fluorescence in the EEZ of Qatar, while the spring and early summer have lower concentrations. The central Gulf is one of the phytoplankton abundant regions in the Gulf during winter, which is related to the circulation patterns (Polikarpov et al. 2016). This seasonal abundance is also consistent with the observations in the northern Gulf, especially in the Kuwait waters, since winter is characterized by welloxygenated and nutrient-rich waters in this region (Ahmed et al. 2022). A more detailed investigation is required to investigate the temporal and spatial variability of the chlorophyll fluorescence and its impact on the productivity in the EEZ of Qatar.

Summary and conclusions
Seasonal variability in temperature, salinity, density, dissolved oxygen, and chlorophyll fluorescence along a transect from nearshore to offshore in the eastern part of the EEZ of Qatar has been investigated. The sea surface temperature varies from nearshore to offshore with a positive gradient during January, June, and August, of the order of 1.0-2.0 °C, and with a negative gradient during October and April, of the order of 1.0 °C. Thermal stratification began in June, attained a well-defined thermocline with a vertical difference in temperature of about 10 °C during August, and started to de-stratify during October.
The nearshore region has the highest salinity compared to the offshore with an increment of about 2.0 in all the seasons except August. The late summer is characterized by an inflow of a relatively fresh and low-density water mass through the surface layer of the intermediate offshore stations of the transect, which has freshened the nearshore waters with a decrement of about 2.0 compared to other seasons. This substantiates the influence of the IOSW in the EEZ of Qatar, as is also evident from recent studies (Elobaid et al. 2022;Rakib et al. 2021). This surface flow has enriched the surface layer (of about 20 m) dissolved oxygen in the middle of the transect with a concentration up to 175 µmol l −1 during August, from a concentration of about 75 µmol l −1 during June. While hypoxia is evident in deeper depths (below 50 m) of the offshore station during August, earlier studies (Elobaid et al. 2022;Al-Ansari et al. 2015) reported hypoxia during September. Our analysis does not show hypoxia during June, and there are no Fig. 7 Cross-sections of chlorophyll fluorescence (μg l −1 ) with cruise months indicated on the left datasets available during July in the EEZ of Qatar. Hence, this study infers that the hypoxia in the EEZ of Qatar evolves in the middle of summer and is prevalent during late summer (August-September), which subsequently disappears during October. The stratification of offshore waters by excessive surface heating, dense and relatively cold deep-water flow originating from the northern Gulf and relatively calm wind conditions together play a role in the evolution and existence of hypoxia in the EEZ of Qatar.
The surface layer up to 10-m depth has lower chlorophyll fluorescence (below 1.0 μg l −1 ) from April to August, whereas the sub-surface layer (between 10 and 40 m depths) has distinct seasonal variability, with the lowest chlorophyll fluorescence in April 2013 (~ 1.0 μg l −1 ), moderate in June 2013 (~ 1.5 μg l −1 ), and relatively high in August 2013 (~ 6.0 μg l −1 ). The chlorophyll fluorescence data during winter is very limited in the EEZ of Qatar. Together with earlier observations, we can deduce that the winter and late summer have relatively higher chlorophyll fluorescence in the EEZ of Qatar, while the spring and early summer have lower concentrations. Obvious in situ data gaps still exist in the EEZ of Qatar, which need to be filled with more observations in future studies. Further studies on phytoplankton species distribution, nutrient concentrations, and flow dynamics are currently underway for a more comprehensive understanding of the coupled physical-biogeochemical dynamics of the EEZ of Qatar.