1. Introduction
The eastern Gulf of Thailand (eGOT), depicted in
Figure 1, spans over 25,000 km
2 and constitutes a shallow coastal water body situated at the northeastern extremity of the Gulf of Thailand. With an average depth of 35 m and the deepest area located offshore in the southwestern direction, reaching a maximum depth of 65 m, the eGOT encompasses approximately 500 km of coastline, which represents approximately 15% of Thailand’s total coastline. It serves as a crucial site for marine environmental and socio-economic services in Thailand. The sea within the eGOT borders the coastal zones of Rayong, Chanthaburi, and Trat provinces in Thailand, as well as western regions of Cambodia, hosting vital industries, tourism facilities, cities, fisheries, and farmland. Despite its significant footprint, there exists limited documentation regarding its hydrographic characteristics. A comprehensive understanding of the hydrographic properties and their variations is essential for facilitating improved sea and resource management in this region.
Previous hydrographic investigations pertaining to the eastern Gulf of Thailand (eGOT) have predominantly focused on the broader characteristics of the Gulf of Thailand on a regional scale, with the eGOT area being merely a subset of their scope. Notably, ref. [
1] utilized data from the NAGA expedition to provide an initial description of the intricate hydrographic dynamics within the South China Sea and the Gulf of Thailand, highlighting their seasonal variability linked to the Asian–Australian monsoon. Subsequent numerical studies, including those by [
2,
3,
4,
5] have further elucidated the seasonal variations in Gulf of Thailand circulation. These variations are shown to be influenced by a combination of factors, including tides, prevailing winds, river discharge, direct rainfall to the sea, and sea surface heat exchanges. Ref. [
6] conducted an analysis utilizing data obtained from High-Frequency Radar, unveiling seasonal variations in tidal residual surface current patterns. Additionally, ref. [
7] corroborated these observations, providing further insight into the intricate complexities of these variations through the utilization of remotely sensed data.
Early works that provide detailed insights into the eGOT are scarce [
8], utilizing remotely sensed observations to discern the presence of a warm water pool during winter in the northeastern region of the Gulf of Thailand, adjacent to the eGOT. This phenomenon was attributed to orographic effects from the Cardamom Mountains in Cambodia. Furthermore, ref. [
9] conducted limited field measurements using current meters at Ban Leam Sork in Trat province, revealing robust tidal currents with maximum speeds reaching 24 cm/s. They observed alternating flow directions during flood and ebb tides, along with residual currents predominantly towards the southeastern direction in December. Additionally, they documented the intricate seasonal patterns of residual currents around Trat province. These patterns were found to be influenced by strength of the monsoonal winds, which were further modified by tides and local topography.
Expanding upon previous research, field evidence and unpublished reports have revealed peculiar hydrographic features within the eastern Gulf of Thailand (eGOT). Notably, despite the absence of significant river inflows, instances of low salinity seawater (<28 ppt) have been documented, with notable salinity disparities (>2–3 ppt) between near-surface and sub-surface layers. We hypothesize that this reduced salinity stems from direct rainfall into the sea. Although lacking major river inputs, the eGOT experiences rainfall rates exceeding the national average. Direct rainfall likely diminishes coastal water salinity, thereby reducing overall water density. Consequently, a freshwater layer can form atop denser seawater, resulting in stratification that hampers vertical mixing and restricts exchanges of heat, nutrients, and dissolved gases between surface and deeper layers. Such stratification may induce buoyancy-driven flows and horizontal currents, influencing coastal hydrodynamics and distinguishing the eGOT from the broader Gulf of Thailand (GOT) system. Similar phenomena of rainfall influence on coastal dynamics have been documented both in the open ocean [
10,
11] and coastal environments [
12,
13,
14]. Moreover, studies by [
5,
15] have demonstrated the substantial impact of atmospheric freshwater fluxes of the Asian–Australian monsoon system on the Bay of Bengal and the Gulf of Thailand, respectively. Alterations in coastal hydrodynamics induced by direct rainfall can yield ecological ramifications. Changes in circulation patterns and stratification can affect nutrient availability, primary production, and the distribution of plankton and other marine organisms, thereby influencing ecosystem dynamics [
16,
17,
18].
This paper aims to offer detailed insights into the topographic, meteorological and three-dimensional hydrographic characteristics of the eGOT intricate three-dimensional hydrography of the region. Furthermore, we present influence of direct rainfall on local three-dimensional hydrography at the eGOT and the emergence of the Chanthaburi Coastal Current (CCC) through numerical experiments. This finding holds relevance for hydrodynamic modeling endeavors in tropical shallow coastal regions that receive significant freshwater inputs, not only from river systems but also from direct rainfall.
3. Results and Discussion
3.1. Field Observation Results
The seasonal spatial distributions of field data are depicted in
Figure 4,
Figure 5 and
Figure 6, while the temporal variations are illustrated in
Figure 7. The temporal data are from dataloggers that are deployed continuously and cover over a year, representing timescales not documented elsewhere. The results reveal daily, tidal, non-tidal longer-term, and seasonal variations in hydrographic conditions. There are significant seasonal fluctuations in water temperature, ranging from 24 to 33 °C. The highest water temperature is observed during Transition 1, contrasting with the lowest temperature recorded during the NEM period. In terms of salinity at Koh Monnai Island, levels range between 26 and 33 PSU, with peak salinity occurring during the NEM season and the lowest during the SWM and Transition 2. Consistent water temperatures and salinities are observed between near-surface and near-bottom layers at Koh Monnai Island and Koh Kood Island, except for the SWM period when near-surface temperatures are warmer/fresher than near-bottom values.
Water level data exemplify pronounced tidal fluctuations, with a tidal range of 1–2.5 m during neap and spring tides. Tides are mixed tide prevailing diurnal which governed by M2, S2, O1, P1, K1 tidal. The amplitudes were 0.34, 0.16, 0.53, 0.18, 0.1 m, respectively. The Formzahl number (K1 + P1 + O1)/(M2 + S2) has a value of 1.62. The annual sea level variation in this area, from measured data, is approximately 0.24 m, with the highest average sea level occurring at the start of the year and the lowest in the middle of the year. Additionally, current measurements depict significant tidal variations. Specifically, the current flow field at Koh Kood Island’s measurement station during ebb tide reveals a complex pattern, with opposite flow directions near the curved coastline compared to the main southeastward current offshore.
Field observations provide initial insights into the significance of tidal influences and complex coastal topography in the eGOT. Although not included here, spatial flow velocity measurements using ADCP at Koh Kood Island reveal complex alternating flow velocity patterns across different tidal phases. An example of the measurement during ebb tide is presented in
Figure 6.
3.2. Model Validation
Model validation against long-term measured data between October 2019 and December 2021 includes the water level at Leamsing Station, and water temperature and salinity at Koh Monnai Island. The model demonstrates good agreement with the observed water level at Leamsing, with a Root Mean Square Error (RMSE) of <14.8 cm from a tidal range during spring tide of approximately 250 cm and an R
2 value of 0.91. For water temperature, the RMSE is 0.76 °C with an R
2 value of 0.81 (see
Figure 7). Salinity levels and trends are also well reproduced by the model. The absence of reported RMSE and R
2 values for salinity is due to concerns regarding the absolute accuracy of the measured salinity data. It is anticipated that while the salinity data can capture long-term variations in this region, there may be inaccuracies present stemming from bio-fouling caused by inadequate maintenance of the HOBO-Conductivity logger deployed in remote area. Frequency analysis of water levels at Koh Monnai Station identifies M2, S2, O1, P1, and K1 as dominant tidal constituents, with amplitudes/phases of 0.34/193, 0.16/340, 0.53/328, 0.18/359, and 0.1/135 m/degrees, respectively. The numerical model similarly reproduces realistic water levels, identifying the same primary constituents with corresponding amplitudes/phases of 0.34/191, 0.15/337, 0.50/325, 0.12/329, and 0.04/113 m/degrees, respectively. Effects of tide and complex coastline topography. The complex alternating flow velocity fields governed by tide and irregular coastline are also well reproduced as seen when compare
Figure 6 and
Figure 8b.
3.3. Tide and Tidal Currents
Tidal currents play a pivotal role in mobilizing and mixing the entire water column within coastal seas. In the eGOT, tides drive water movement alternately along the coastline in a northwest to southeast axis. During the flood tide, tidal currents flow northwestward, while they shift southeastward during the ebb tide. Tidal currents ease during high and low tides. Maximum tidal current speeds during spring tides range between 0.3 and 0.5 m/s, halving during neap tides. Offshore tidal current speeds generally exceed those near the shoreline, except in channels between land and island and island and island (see
Figure 8a). Local shoreline irregularities and bathymetry intricately influence tidal flow patterns. At the lee side of headlands, currents change direction, forming clockwise and counterclockwise eddies, evident in both field measurements (see
Figure 5) and well reproduced by modeling (see
Figure 8b).
3.4. Water Temperature and Salinity
The results from field observations (
Figure 4,
Figure 5 and
Figure 7) and modeling (
Figure 9 and
Figure 10) demonstrate good agreement in depicting the spatial and seasonal variability of water temperature and salinity in the eGOT, with ranges spanning between 24 and 30 °C and 29 and 32.5 ppt, respectively. Water temperature is lowest during the NEM, rapidly rising to its annual peak during Transition 1. This temperature fluctuation is most pronounced along the shoreline, where waters are coolest during the NEM and warmest during Transition 1. In contrast, salinity exhibits marked changes during the SWM and Transition 2, with reductions observed along the shoreline, particularly along the coasts of Chanthaburi and Trat provinces. It is important to note that since the model does not account for freshwater input from rivers, the observed reduction in salinity is attributed to direct rainfall into the sea surface.
The model results reveal notable water column stratification characterized by alternating water temperature and salinity profiles. This phenomenon persists throughout the year, except during the NEM when the water column exhibits relatively homogeneous conditions. During Transition 1, a distinct temperature stratification emerges, with the near-surface layer up to 30 m depth exhibiting higher temperatures compared to the lower layers. Following Transition 1, water temperature becomes more evenly distributed throughout the water column. Conversely, salinity does not exhibit significant stratification during Transition 1 but becomes strongly stratified during the SWM and Transition 2. Both water temperature and salinity revert to a more uniform vertical distribution during the NEM.
3.5. Seasonal Water Circulation
The monthly mean simulated flow velocity fields depict the diverse circulation patterns of the eGOT influenced by monsoonal effects, as illustrated in
Figure 11. Two primary circulation patterns emerge, characterized by flow directions to the southeast during Transition 1 and the SWM, and to the northwest during Transition 2 and the NEM. Monthly average current speeds peak (>0.15 m/s) during the SWM and the NEM, with the strongest currents observed offshore during the SWM and approximately 40 km from the shoreline during the NEM. Current speeds diminish during the monsoon transitions, albeit with a notable exception during the SWM and Transition 2, where a robust coastal current (within 10 km from the shoreline) flows northeastward with speeds exceeding 15 cm/s. This coastal current, flowing along Chanthaburi Province in Thailand, is dubbed the “Chanthaburi Coastal Current (CCC)” and forms during the SWM, becoming particularly strong and stable during Transition 2.
Additionally, the model results highlight significant differences between near-surface and sub-surface circulation, which vary seasonally. Near-surface currents tend to flow towards the land during Transition 1 and the SWM, shifting to offshore flow during Transition 2 and the NEM. Sub-surface counter-currents persist during periods of water column stratification. In
Figure 11, averaged current profiles along and across section A-A, location as shown in
Figure 1, further elucidate the layered flow patterns during stratified periods, displaying a three-layer and two-layer flow system during the SWM and Transition 2, respectively. Conversely, during the NEM, when the water column is well mixed, flow tends to be more homogeneous.
3.6. Influence of the Direct Rainfall into the Sea
The substantial spatial and temporal variations observed in salinity patterns within the eGOT are expected to exert significant influence on water column dynamics and circulation regimes. Given the prominent role of rainfall in shaping salinity, this study endeavors to explore the impact of direct rainfall through the establishment of an additional model. Simulations conducted with and without direct precipitation yield insightful outcomes, as depicted in
Figure 12 and
Figure 13, respectively.
A comparative analysis of these outcomes highlights the pronounced effect of direct rainfall, particularly during the SWM and Transition 2 phases, characterized by intensified rainfall and reduced wind activity within the eGOT. While the influence of rainfall on other temporal intervals and water temperature remains relatively subdued, the marked decrease in salinity significantly contributes to water column stratification during the SWM and Transition 2. This stratification profoundly alters overarching circulation patterns, precipitating the emergence of the CCC during these intervals. Notably, the removal of rainfall from the simulations results in the cessation of the CCC and the attendant two-layer circulation.
The influence of freshwater influx from direct rainfall into the marine environment has been previously documented to affect various phenomena, such as the sea surface microlayer [
11] and the formation of freshwater lenses under calm wind conditions within the California Current [
10]. These effects subsequently exert a notable impact on the hydrodynamic characteristics of the marine domain. Our study marks the first documentation of the effects of direct rainfall on the generation of upper-layer low salinity water, thereby contributing to the genesis of the CCC. Interestingly, a similar occurrence of upper-layer low salinity coastal currents was recently observed during Transition 2 along the western coast of the Gulf of Thailand [
26]. However, in that instance, the diminished salinity levels were attributed to river discharge.
4. Conclusions
Noted from field observations, hydrographic conditions at the eGOT exhibit significant spatial and temporal variability, with pronounced variations in water temperature, salinity, and flow velocities influenced by the Asian–Australian monsoon. Offshore and deep areas typically experience lower water temperatures, reaching their lowest levels during the NEM. During Transition 1, lower layer temperatures are 2–3 degrees lower than near-surface values. The SWM brings higher freshwater input from rainfall, resulting in a significant reduction in the salinity of near-surface layers. This salinity decrease is particularly prominent during the SWM and Transition 2. Differences between near-surface and sub-surface, and near-shore and offshore salinity can be 2–3 ppt, significantly impacting seawater density. The influx of freshwater from rainfall diminishes the overall salinity of coastal waters, leading to decreased water density. During the SWM and Transition 2, strong stratification is observed in shallow areas (water depth < 15 m), influencing coastal hydrodynamics. Salinity levels are low during these periods when rainfall rates significantly increase and wind is relatively calm. Stratification affects vertical mixing and circulation patterns of the water column, altering hydrodynamics. It can hinder vertical mixing and restrict exchange of heat, nutrients, and dissolved gases between surface and deeper layers [
10,
11,
12,
13,
27].
The water circulation in the eGOT exhibits two distinct modes, flowing southeastward during Transition 1 and the SWM, and northwestward during Transition 2 and NEM, closely aligned with predominant wind directions. The strongest circulation occurs during the SWM and NEM periods. Layered flows induced by water column stratification are particularly notable during the SWM and Transition 2. During these periods, a northwesterly Chanthaburi Coastal Current (CCC) emerges. Model simulations suggest that the CCC can be sustained by direct rainfall into the sea during the SWM and Transition 2. The presence of the CCC may significantly influence the seasonal transport of sediments, nutrients, and pollutants along the eGOT coastline, with ecological ramifications. Altered circulation patterns and stratification affect nutrient availability, primary production, and the distribution of plankton and other marine organisms, ultimately impacting the entire ecosystem and fisheries dynamics [
17,
18].
Future work is essential to conduct detailed sensitivity analyses and model runs to investigate the mechanisms governing the Chanthaburi Coastal Current (CCC) and its variations. This may include considering freshwater discharges from small rivers and creeks along the eGOT coastline, as well as grid refinement and additional continuous field observations for model validation along the coastline. These efforts could enhance the model’s capability to capture the detailed behavior of the CCC. Furthermore, investigation of the effects of changing rainfall patterns in the future, resulting from inter-annual variations in climate drivers such as the El Niño-Southern Oscillation and the Indian Ocean Dipole, is warranted. Additionally, exploring the influence of extreme rainfall events, which are expected to occur more frequently under ongoing global warming, is crucial [
14,
16,
28]. Such future research endeavors will provide valuable insights into the dynamics of the eGOT’s hydrographic features and their responses to environmental changes.