Numerical study of the seasonal salinity budget of the upper ocean in the Bay of Bengal in 2014*

Qingyue WANG , Changming DONG , Junde LI , Jingsong YANG ,3, Qing TIAN

1 School of Marine Sciences, Nanjing University of Information Science and Technology, Nanjing 210044, China

2 State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography, Ministry of Natural Resources, Hangzhou 310012, China

3 Southern Marine Science and Engineering Guangdong Laboratory(Zhuhai), Zhuhai 519080, China

Abstract Impact factors on the salinity budget, especially the eddy salt fluxes and smaller-scale diff usive salt fluxes for the upper 50 m of the Bay of Bengal (BoB) in 2014 are investigated using a box model based on the Regional Ocean Modeling System (ROMS) daily outputs. The model results reproduce that the precipitation and river runoff s are the dominant factors modulating the sharp salinity decrease during the summer monsoon season. The analysis shows that the salinity increase after the summer monsoon is mostly due to the meridional advective and diff usive salt fluxes. The vertical advective salt flux, which is sensitive to the diff erent signals of the wind stress curl, plays an important role in balancing the salinity change induced by the meridional advective salt flux during both the summer and winter monsoon seasons. Distinctive spatial mesoscale structures are presented in the eddy salt flux throughout the year, and their contributions are sizeable(over 30% in the meridional direction and about 10%-30% in the vertical direction). The meridional eddy salt flux is larger in the monsoon seasons than that in the inter-monsoon seasons, and in a positive pattern near the western boundary during the winter monsoon and autumn inter-monsoon. The vertical eddy salt flux makes an important contribution to the salinity budget, especially along the coastal area and around the Andaman and Nicobar Islands. The vertical eddy salt flux becomes large when a tropical cyclone passes the area.

Keyword: salinity budget; Bay of Bengal; Regional Ocean Modeling System (ROMS); monsoon; eddy flux

1 INTRODUCTION

Ocean salinity is a key variable for understanding the eff ect of freshwater fluxes on ocean dynamics.Salinity aff ects the ocean density, which in turn modulates the transport of tracers. The Bay of Bengal(BoB), located between 5°N-22°N and 80°E-95°E(Fig.1), is one of the least salty seas in the world. The salinity budget in the upper BoB is ofimportance to understand the oceanic dynamics (e.g., Schott and McCreary, 2001; Thompson et al., 2006), ecology(e.g., Narvekar and Prasanna Kumar, 2014; Sherin et al., 2020) and climate change of the region (e.g., Trott et al., 2019b).

The salinity budget in the upper BoB is mainly controlled by the sea surface net freshwater flux(evaporation minus precipitation, E-P), river discharge, horizontal advection, vertical advection,and turbulent diff usion. According to diff erent temporal scales, the advection can be further decomposed into the mean flow and eddy flux(Meijers et al., 2007; Köhler et al., 2018).

Fig.1 Bathymetry of the Bay of Bengal from ETOPO2 (2 arc-minute global relief model of Earth’s surface that integrates land topography and ocean bathymetry)

The freshwater flux from surplus precipitation over evaporation is the dominant factor freshening the BoB during the summer monsoon (e.g., Rao and Sivakumar, 2003; Akhil et al., 2014; Pant et al., 2015;Wilson and Riser, 2016). Both observational (e.g.,Pant et al., 2015; Lekha et al., 2018) and modeling studies (e.g., Han and McCreary, 2001; Jana et al.,2015; Sandeep et al., 2018; Trott et al., 2019a) show that a large freshwater discharge from river inputs,such as the Ganges and Brahmaputra at the northern boundary of the BoB (Fig.1), significantly modulates the surface salinity during the summer period.Sengupta et al. (2006) emphasize that the net freshwater input from the precipitation minus evaporation and river discharge is over 4 000 km3/a(i.e., 0.13 Sv, where 1 Sv=106m3/s) in the BoB north of 6°N.

Horizontal advection is vital for the salinity distribution in the upper BoB, and shows remarkable seasonal variability. During the summer monsoon, the surface circulation in the northern BoB (north of 15°N) is a cyclonic gyre (Pant et al., 2015), with the East Indian Coastal Current (EICC) flowing southwestward and transporting low-salinity water southward along the western boundary of the northern bay (Shetye et al., 1996; Schott and McCreary, 2001;Sandeep et al., 2018). Meanwhile, the high-salinity water from the Arabian Sea intrudes into the southern BoB during the summer monsoon (Pant et al., 2015;Hormann et al., 2019). Currents along the eastern boundary of the BoB flow southward during this period (Hormann et al., 2019). During the late fall and early winter monsoon (November-January), the EICC still flows southwards, transporting low-salinity water along the western boundary of the bay into the whole basin (Shetye et al., 1996). During late winter monsoon, the basin-wide anticyclonic gyre controls the surface BoB circulation, coupled with the reversing northeastward EICC (Shetye et al., 1996;Schott and McCreary, 2001) and the freshwater plume along the eastern boundary of the bay (Han and McCreary, 2001; Jana et al., 2015). During most of the year, the eastern rim of the BoB, including the Andaman Sea (AS), is another freshwater outflow pathway toward the equator (Han and McCreary,2001; Akhil et al., 2014; Wilson and Riser, 2016).

While previous studies propose that the eff ect of vertical processes is not comparable to the freshwater flux and horizontal advection (e.g., Rao and Sivakumar, 2003), recent studies demonstrate the significant eff ects of vertical processes on modulating upper-ocean salinity (Vinayachandran et al., 2013;Akhil et al., 2014; Wilson and Riser, 2016). Both observational data (Vinayachandran et al., 2013) and model results (Akhil et al., 2014) show the significant mechanism of saltier water upward pumping to keep the bay salty during the summer monsoon. Akhil et al.(2014) note that the vertical process (i.e., the sum of vertical advection, vertical mixing, and entrainment)counterbalance the sea surface salinity (SSS)anomalies along the eastern Indian coast during winter. Wilson and Riser (2016) argue that vertical salt fluxes (including the vertical advection and mixing) balance the salinity decrease within the upper 30 m of the ocean during and immediately following the summer monsoon.

Most previous studies consider that the horizontal turbulent diff usion is negligible (Akhil et al., 2014;Wilson and Riser, 2016; Köhler et al., 2018). In addition, the vertical turbulent diff usion is often analyzed with vertical advection and/or entrainment and collectively referred to as vertical processes(Akhil et al., 2014, 2016; Wilson and Riser, 2016).However, one recent study (Köhler et al., 2018) uses model outputs, as well as satellite data to conclude that the vertical diff usion is significant in counteracting the freshening due to precipitation and river runoff s,especially along the western boundary of the BoB in September.

The above reviews reveal that the seasonal monsoon significantly modulates the salinity budget,but diff erent factors contributing to the salinity budget in diff erent seasons are not well understood, especially the role of eddy flux remains unknown. These questions motivate us to conduct the present study by applying the eddy-resolved model outputs to revisit the seasonal salinity budget by examining the contributions of the above factors.

The rest of the paper is organized as follows. In Section 2, we describe the Regional Ocean Modeling System (ROMS) model configuration, observational data and introduce a box model for the salinity budget equation. The calculation of the eddy flux is also introduced in this section. Section 3 validates the model outputs, presents the salinity budget in the upper BoB in 2014 and discusses all significant factors aff ecting salinity changes. Eddy salt fluxes in each season are also analyzed. Sections 4 & 5 give the discussions and conclusions of the study, respectively.

2 MATERIAL AND METHOD

2.1 ROMS model configuration

This study uses the ROMS model, which is a threedimensional ocean model, solving the primitive Navier-Stokes equations under the Boussinesq approximation and hydrostatic balance assumption.The vertical coordinate system in the ROMS is a generalized sigma-coordinate system, and the horizontal coordinate is a coastline- and terrainfollowing curvilinear coordinate. There are other terrain-following coordinate models, such as Finite-Volume Coastal Ocean Model (FVCOM) (e.g., Chen et al., 2007; Xue et al., 2009), Semi-implicit Eulerian-Lagrangian Finite-element Model / Semi-implicit Cross-scale Hydroscience Integrated System Model(SELFE/SCHISM) (e.g., Zhang and Baptista, 2008;Zhang et al., 2016) and Delft3D-FLOW model (e.g.,Bouma et al., 2007; Hu et al., 2009), which are also widely used in the community. The ROMS is a freesurface split-explicit ocean model, with short time steps for barotropic dynamic (e.g., sea surface height and two-dimensional momentum) and larger time steps for baroclinic dynamic (e.g., temperature,salinity, and three-dimensional momentum). It has been widely used in regional ocean simulation, airsea interaction and ocean dynamics studies (e.g.,Dong et al., 2007; Jacox et al., 2016; Li et al., 2018;Zhang et al., 2018a).

The model domain spans 30°E-110°E and 30°S-30°N, with a 1/8°×1/8° horizontal grid resolution and 32 vertical levels. The vertical resolution ranges from about 2 m at the surface to about 5 m at the bottom in the regions shallower than 100 m, with over 20 levels in the upper 50 m; and it ranges from about 7 m at the surface to about 300 m at the bottom in the regions deeper than 3 500 m, with 5 levels in the upper 50 m. The vertical turbulent diff usivity follows the K-profile parameterization(KPP) scheme, which is described in detail by Large et al. (1994). The initial condition for December 1978 is obtained from the Simple Ocean Data Assimilation(SODA) reanalysis release 2.2.4 (Carton and Giese,2008). The model is forced with 6-hourly heat,freshwater and momentum fluxes at the sea surface from the Climate Forecast System Reanalysis (CFSR)data (Saha et al., 2010) from January 1979 to March 2011, and the Climate Forecast System Version 2(CFSv2) data (Saha et al., 2014) from April 2011 to December 2014. Monthly-averaged temperature,salinity, velocity, and sea surface height SODA data from January 1979 to December 2010 and the Estimating the Circulation & Climate of the Ocean(ECCO) reanalysis data from January 2011 to December 2014 are used at the lateral boundaries.The data for river runoff s are from the Climate and Global Dynamics Laboratory (CGDL). Tidal eff ects are not considered in this model. The bathymetry is from ETOPO2, a 2-arc-minute global relief model of Earth’s surface that integrates land topography and ocean bathymetry. The model is first spun up for 3 years (1979-1981). From January 1982 to December 2014, monthly-averaged data are outputted, but dailyaveraged data are outputted for 2014 only. The detailed settings of the model applied in this study can be also found in Li et al. (2016).

The daily-averaged model outputs in 2014 are used for the analysis of the seasonal salinity variation.There are no specific reasons to choose the year 2014,but the analysis of the daily-sampled numerical outputs is good enough for the understanding of the seasonal variation.

2.2 Argo data

The Array for Real-time Geostrophic Oceanography(Argo; http://www.argodatamgt.org/) is a global array of free-drifting floats, with currently over 4 000 floats in the global ocean. Each float records one vertical profile every 10 days from the surface to about 2 000 m in depth, continuously measuring water temperature and salinity.

In this study, monthly salinity data from the Barnes Objective Analysis-Argo (BOA-Argo) dataset in 2014 are used to validate the model outputs. These data, which are provided by the China Argo Real-time Data Center (Lu et al., 2020) have 58 vertical levels in the upper 1 975 m and the grid resolution of 1°×1°.According to the Argo data center (http://www.argodatamgt.org/Access-to-data/Argo-data-selection),there are a total of 1 935 Argo vertical profiles in 2014, which is suffi cient for the generation of the gridded product.

2.3 The climatological drifter data

Drifters are oceanographic devices floating on the sea surface to record the near-surface currents and other parameters like water temperature and salinity.The Global Drifter Program (GDP) is the main component of the Global Surface Drifting Buoy Array, which aims to provide a global array of about 1 300 satellite-tracked surface drifting buoys to obtain an accurate and global in situ observations of ocean currents, sea surface temperature, salinity, atmospheric pressure, and surface winds. The sea surface currents derived from the drifter trajectory data are used to validate the model outputs.

In this study, we use “the Drifter-Derived Climatology of Global Near-Surface Currents” (ftp://ftp.aoml.noaa.gov/phod/pub/lumpkin/drifter_climatology). This monthly data derived from the GDP drifters have the horizontal grid resolution of 1/4° (Laurindo et al., 2017).

2.4 The box model for salinity budget equation

A simplified form of the averaged volumeintegration of the salinity budget equation is conducted in a fixed box. The domain is from 10°N to 23°N and 80°E to 99°E, including the BoB and the AS (Fig.1),and the bottom of the study region is set as 50 m.According to previous studies, the mixed layer thickness in the BoB is approximately 20-50 m(Narvekar and Prasanna Kumar, 2014; Wilson and Riser, 2016; Li et al., 2017; Zhang et al., 2018b; Shee et al., 2019). Therefore, the 50-m depth of the study domain includes the upper mixed layer.

There are three open boundaries: the 10°N section(hereafter referred to as the southern boundary), the sea surface and the 50-m level (hereafter referred to as the bottom boundary). The other boundaries at 80°E, 99°E, and 23°N are closed, with no advection and diff usion across these sections. The simplified salinity budget equation integrated from the free sea surface (η) to the depth of 50 m (h=50 m), is shown as follows:

whereStis the salinity variation rate.SvandSware the advective salt fluxes across the southern boundary and bottom boundary, respectively.SydiffandSzdiffare the diff usive salt fluxes across the southern boundary and bottom boundary, respectively.SEPandSRare the evaporation minus precipitation and river runoff ,respectively. [S] is time varying and volume-averaged salinity of the domain. This method was proposed by Lee et al. (2004) and has been widely used in heat and salt budget analysis (Zhang et al., 2018c; Lee et al.,2019; Soares et al., 2019). Details for the derivation of the box model for the salinity budget equation are given in Appendix A.

2.5 Calculation of the eddy flux

The eddy salt fluxes can be estimated by calculating the deviations from the time-averaged advective salt fluxes:

3 RESULT

3.1 Model validation

The above model product has been used in a series of studies: Li et al. (2016, 2018) and the present study.In Li et al. (2016), the model product has been validated with the Argo data and the SODA reanalysis dataset in terms of sea surface salinity anomalies in the tropical Indian Ocean. In Li et al. (2018), the observational and reanalysis data including the Argo data, the Research Moored Array for African-Asian-Australian monsoon analysis and prediction (RAMA),and the European Centre for Medium-Range Weather Forecasts (ECMWF) Ocean Reanalysis System 4(ORAS4) have been used in the model validation for temperature, salinity, and zonal current anomalies along the equatorial Indian Ocean.

In the present study, the modeled currents and SSS are validated against the observational data. We choose January and July as the representative months for the winter and summer monsoons, respectively, to validate the model performance in terms of the sea surface circulation seasonality (Fig.2) and the SSS(Fig.3).

Fig.2 Monthly-averaged surface horizontal currents in the BoB derived from drifter data in January (a), ROMS data in January (b), drifter data in July (c), and ROMS data in July (d)

Figure 2a-b display the monthly-averaged horizontal currents during the winter monsoon from drifter and ROMS data, respectively. It can be seen that ROMS reproduces the spatial distribution of the surface currents: a cyclonic circulation north of 15°N,the southward EICC southward of 18°N and a strong northward incoming flow at the southern boundary(about 10°N). However, the ROMS outputs underestimate the southward EICC and the circulation in the AS.

Figure 2c-d plot the monthly-averaged horizontal currents during the summer monsoon from the drifter and ROMS data, respectively. A similar general spatial distribution of the surface currents can also be seen: an anticyclonic circulation northward of 15°N,the southward EICC south of 15°N, a strong northward outflow at the southern boundary (around 85°E) and an approximately eastward flow east of 87°E. Similar to the winter monsoon comparison, the ROMS outputs also underestimate the EICC all along the western coast of the bay and the flow along the eastern coast, and overestimate the anticyclonic circulation in the northern bay area.

To further validate the model performance in terms of the climatological circulation simulation at the sea surface, the climatological model outputs (between 1982-2014) are compared with the drifter data (not shown). Even though there are diff erences in magnitude and along the west coast, the numerical results are consistent with the drifter data in terms of the spatial distribution patterns: in the northern BoB,the model reproduces the cyclonic (anticyclonic) gyre in summer (winter).

Many factors could contribute to the errors between the modeled and the drifter-derived circulations in this area. Firstly, the model is forced by 6-hourly wind field at the sea surface from the CFSv2 data in 2014,which could carry uncertainties for the oceanic circulation, especially in summer, when the local extreme weather conditions cannot be reproduced by the global model. Secondly, the drifter data used for the comparison are the climatological data, which are not a single-year oceanic circulation as presented in this study for the year 2014. Moreover, the number of the drifter records used in this set of observational data is small and unevenly distributed in the BoB,especially along the coastal region (Fig.1 in Laurindo et al., 2017). Such insuffi ciency in the statistical sampling number could severely aff ect the reliability in the circulation in the southwest of the BoB, which could be the major reason for the inconsistent results in the comparison between the modeled and the drifter-derived circulations in the area.

Fig.3 Monthly averaged SSS in 2014 in the BoB derived from Argo data in January (a), ROMS data in January (b), Argo data in July (c), and ROMS data in July (d)

Fig.4 Annually averaged SSS in 2014 in the BoB derived from Argo (a) and ROMS data (b); SSS diff erence between the ROMS and Argo data (model data minus buoy data) (c); histogram of the SSS deviation of ROMS outputs from the Argo data (d)

Compared to Argo data, the ROMS model reproduces the SSS spatial distribution, albeit a slight overestimation ofless than 0.6 over most of the region, except in the AS where a small underestimation of ～0.2 (absolute deviation value) is observed (Fig.4).Compared with previous numerical simulations of SSS in the same area (Akhil et al., 2014; Wilson and Riser, 2016), the current results exhibit smaller errors.The good agreement between the ROMS and Argo data justifies the application of the model results in the present study.

Fig.5 Comparison between the two sides of Eq.1 (blue line: left side of the equation, red line: right side of the equation) and salinity variation (black line) in the upper 50 m of the BoB in 2014

As shown in Fig.3, salinity from both the model and Argo data shows a freshwater tongue extending from the east of the bay all the way to the west of the bay in January, but the freshwater tongue disappears in July. Instead, the salinity gradually decreases northeastward in July. Although the model results are in agreement with observation, there are discrepancies between the two data, especially north of 15°N. A possible explanation is that the ROMS model is forced with climatological river runoff s from CGDL,which is diff erent from the river discharge in 2014.

To further validate the model performance in terms of the salinity simulation, we also make the comparison of the salinity vertical distribution in the domain with the BOA-Argo data (not shown), and the results show that the model reproduces the similar pattern as the observation.

3.2 Salinity variance balance

The comparison of the time series for the terms on the two sides of Eq.1 shows good agreement (Fig.5).The correlation coeffi cient between the salinitychanging-rate term (left side of the equation) and the sum of all influencing terms (right side of the equation)is 0.99, which means that the left and right sides of the equation are highly correlated. The root-mean-square error between the two sides of the equation is 0.09×10-2, which is negligible as compared to the terms of Eq.1 in magnitude of 1×10-2. AlthoughStis slightly larger than the sum of all influencing terms from February to April, it is within the acceptable error range. The change in total salt volume (Stotal)shows a sharp decrease during the summer monsoon and increase during the winter monsoon. The change in salinity is not significant during the inter-monsoon periods (Fig.5).

3.3 Salinity budget analysis

3.3.1 Eff ects of net freshwater flux

SEP, the evaporation minus precipitation term, is the dominant factor resulting in the sharp salinity decrease in the upper BoB (Fig.6b), and it presents a significant positive correlation of 0.97 (not shown) withSt(Fig.6a). During the winter monsoon and early spring inter-monsoon (January to April),SEPis positive,which means that evaporation is larger than precipitation, thus causing an increase in salinity.From May to November (mainly including the summer monsoon and autumn inter-monsoon),SEPis negative, which suggests that precipitation is more than evaporation, hence resulting in the sharp decrease in salinity.SEPreaches its minimum value in July when the region experiences heavy rain. Overall, the trend ofSEPresembles that ofSt, suggesting that it contributes significantly to the change ofSt(Fig.6a).

Fig.6 Monthly (a) and annual (b) averages of salinity budget terms

SR, the river runoff , is the other significant factor making the decrease of the salinity in the BoB(Fig.6b). It also has a similar trend to that of bothStandSEP, but the freshwater injection is smaller than inSEP(Fig.6a).Therefore, the large net freshwater flux during the summer monsoon, especially the precipitation, is one of the main factors leading to the sharp salinity decrease in the upper BoB.

3.3.2 Eff ects of advection

Sv, the meridional advective salt flux across the southern boundary, is the most important process that balances the excess precipitation and runoff over the full annual cycle (Fig.6b). It is usually positive and becomes the most considerable term among all the influencing factors during the winter monsoon,reaching its peak in December and January and playing a leading role in increasing the salinity in the bay (Fig.6a). During this period, the strong southward EICC transports a relatively salty water out of the bay.A northward flow carries salty water into the BoB through the middle area of the 10°N section between 81°E and 95°E, which is saltier than the water exported southward by the EICC. Moreover, the transport is stronger in the upper 30 m (Fig.7a). In the eastern AS, a southward current exports the relatively fresh water out of the bay.

The magnitude ofSvdecreases sharply during the spring inter-monsoon (Fig.6a). During this period, the EICC shifts to a northward current with salty water injection (Fig.7b). A negative advective salt flux appears between 81°E and 89°E because of the southward flow and positive salinity diff erence, while a positive salt flux comes up between 89°E and 92°E because the salinity diff erence turns to negative. The northward current transports fresher water into the AS, resulting in a net salinity decrease (Fig.7b).

Fig.7 Seasonally-averaged meridional velocity ( v, shaded) eff ects on the salinity diff erence ( S-[ S], lines) at the southern boundary in the winter monsoon (a), the spring inter-monsoon (b), the summer monsoon (c), and the autumn inter-monsoon (d) in 2014

Svremains negative and presents a minimal magnitude in the early summer monsoon (June to July). It changes to positive in late summer monsoon,balancing the salinity freshening caused by the huge net freshwater injection (Fig.6a). During this period,currents at the southern boundary are mostly southward except for the northward flow between 85°E and 96°E and below 20-m depth. The salinity diff erence mainly shows a positive pattern except for the AS, which in turn causes the negative salt flux basically located in the west of the southern boundary and the positive advective salt flux mostly located in the middle of the section below 20 m and the AS(Fig.7c).

Svis positive during the autumn inter-monsoon(Fig.6a). During this period, the strong southward EICC transports the low-salinity water off the western BoB. A northward flow carries salty water through the middle area of the southern boundary between 83°E and 90°E. The southward currents along with the positive salinity diff erence induce a negative advective salt flux between 90°E and 94°E and the northward currents along with the negative salinity diff erence also induce a negative advective salt flux in the upper 40 m of the AS (Fig.7d).

The annual mean ofSw, the advective salt flux across the bottom boundary, is smaller thanSv(Fig.6b), while it is important especially during the monsoon seasons (December to February and June to July) to balance the salinity change induced bySv(Fig.6a). The vertical advective salt flux is sensitive to the Ekman pumping velocity (Ashin et al., 2019) at the depth of 50 m, that is,

and then is highly related to the wind stress curl. In Eq.7,f0is the reference Coriolis parameter,ρ0is the reference density of seawater, k is the vertical unit vector, andτis the wind stress. To balance the positive(negative) wind stress curl,w|z=-50cmmust be upward(downward) according to Eq.7.

During the winter monsoon,Swis significantly negative, making the upper bay fresher and counterbalancing the net northward advective salt flux at the southern boundary (Figs.6a & 7a). The wind stress curl is mostly negative in the basin, except for the coastal region and the southern boundary of the study area (Fig.8c). The negative wind stress, in turn, causes downwelling in the center of the bay,transporting the salty water out of the upper BoB at 50 m. The upwelling is concentrated along the eastern coast of the bay, carrying the fresher water upwards at 50 m (Fig.8a). It should be noted that the salinity anomaly at the bottom boundary is derived from the time varying and volume-averaged salinity [S] of the domain. Due to this reason, even though the salinity generally increases with depth, an upwelling can transport water that is fresher than the control volume.

Fig.8 Seasonally-averaged vertical velocity ( w, shaded) eff ects on the salinity diff erence ( S-[ S], lines) at the bottom boundary(the upper row), and wind eff ects (vectors: wind stress; shading: wind stress curl) at the surface (the lower row) in winter monsoon (a, c) and summer monsoon (b, d) in 2014

During the summer monsoon,Swis significantly positive, and counterbalance the negative advective salt flux due to the large freshwater injection, thus making the upper bay saltier (Fig.6a). The wind stress curl is remarkably positive in the northern bay and along the western coast (Fig.8d). It, in turn, causes upwelling in the northern and western bay, transporting the salty water into the upper BoB at 50 m (Fig.8b).

3.3.3 Eff ects of turbulent diff usions

Turbulent diff usions also influence the salinity budget (Fig.6b).Sydiff, the meridional diff usive salt flux across the southern boundary, becomes positive from June to December (Fig.6a), balancing the sharp decrease of salinity during the summer monsoon and increasing it in the autumn inter-monsoon and early winter monsoon.Szdiff, the vertical diff usive salt flux across the bottom boundary, has a similar trend withSwbut with a smaller magnitude (Fig.6a). It is also positive during the summer monsoon.

3.4 Eddy salt fluxes

For seasonal analysis, the advection can be decomposed into mean and eddy salt fluxes by applying Eqs.4 & 5. The spatial distribution of seasonally-averaged eddy salt fluxes at the interfaces reveals the detailed mesoscale structures of strong eddy signals (Figs.9 & 10).

Sve, the eddy salt flux across the southern boundary,displays stronger eddy signals in the monsoon seasons than in the inter-monsoon seasons. Also, the positive eddy salt flux pattern exists near the western boundary throughout the year where the EICC located, and its magnitude is extremely large during the winter monsoon and the autumn inter-monsoon (Fig.9).Sveis mostly positive during the winter monsoon and shows strong positive signatures on either side of the Andaman and Nicobar Islands (Figs.1 & 9a). The absolute value ofSvedecreases in the spring intermonsoon and eddy salt flux turns negative in the middle of the section below 30 m and beside the eastern boundary (Fig.9b).Svebecomes strongly positive between 85°E and 90°E under the surface during the summer monsoon and shows abundant mesoscale signals east of 90°E (Fig.9c).Svedecreases again in the autumn inter-monsoon except for the strong positive eddy salt flux near the western boundary. Other weaker eddy-flux signals are mainly concentrated in the upper 20 m and near the eastern boundary (Fig.9d).

Fig.9 Seasonally-averaged meridional eddy flux at the southern boundary during winter monsoon (a), spring inter-monsoon(b), summer monsoon (c), and autumn inter-monsoon (d) in 2014

Fig.10 Seasonally-averaged vertical eddy flux at the bottom boundary during winter monsoon (a), spring inter-monsoon (b),summer monsoon (c), and autumn inter-monsoon (d) in 2014

The spatial distribution ofSwe(eddy salt flux across the bottom boundary) reveals abundant eddy-scale secondary salt flux circulation, with large magnitudes,except for the spring inter-monsoon (Fig.10). During the winter monsoon, strong positive eddy salt flux can be found along the western coast andSweis extremely negative around the Andaman and Nicobar Islands(Fig.10a). A positive pattern is present between 80°E and 90°E in the southwestern BoB, and a weaker positive and negative value is observed in the middle of the bay during this period (Fig.10a). During the spring inter-monsoon, the absolute value ofSwesignificantly decreases, while the eddy-scale signals are still obvious (Fig.10b). During the summer monsoon,Swemakes a significantly negative contribution, especially along the western boundary and in the middle of the bay south of 15°N, while positive eddy salt flux can still be seen in the eastern bay, and around the Andaman and Nicobar Islands(Fig.10c). During the autumn inter-monsoon,Sweplays a remarkable role in the western region of the BoB, with a negative eddy salt flux near the boundary and a positive flux in the middle of the bay (Fig.10d).This strong eddy flux signal is caused by intensive upper-ocean mixing induced by the extremely severe tropical cyclone HudHud. From 7-13 October 2014,HudHud strengthened into a tropical storm near the Andaman and Nicobar Islands around 12°N, moved northwest, strengthened into a typhoon on 10 October,and made landfall on the east coast of India near 18°N(Fig.10d), during which the maximum wind speed exceeded 50 m/s (Murty et al., 2020 and references therein).

4 DISCUSSION

In the present study, we choose 2014 as a representative year to examine the seasonal salinity variations in the BoB. Given the strong seasonal variation in the BoB (controlled by the winter and summer monsoons), whose magnitude is much larger than the inter-annual variations, the analysis of the daily-sampled numerical outputs can shed more light on the understanding of the seasonal salinity variation.

Even though the 50-m depth of the BoB includes the upper mixed layer, it is useful to examine how the depth choice aff ects the salinity budget. We have conducted the calculation of the salinity budget for diff erent bottom depths of the box (not shown).Svis the main process that balances the surface freshwater input by the precipitation and river discharge in diff erent depths in the upper BoB, although the magnitude ofSvdecreases with increasing depth.Eff ects ofSware diff erent for diff erent depths. It is the most important process that balances the excess precipitation and runoff over the full annual cycle in a shallower depth within the mixed layer, and gets smaller if the base of the control volume becomes deeper.Sydiffhas a similar trend for all the diff erent depths as tested, only the magnitude decreases as the depth increases.Szdiffis the term that has the greatest changes due to the choice of diff erent depths, which causes the decrease (increase) of salinity in the shallower (deeper) depth.

The ROMS model is forced with the climatological data from the CGDL, with a lack of accurate river runoff forcing in the specific year. The arbitrary horizontal diff usivity set in the model required for the computational stability is another limitation, resulting in an overestimation of the horizontal diff usions.

At the sea surface, it is mainly the net freshwater flux that influences salinity transport, includingSEPandSR.SEPincreases salinity in the winter monsoon and spring inter-monsoon, and plays a leading role during the summer monsoon and autumn intermonsoon, reducing the salinity by injecting a large amount of freshwater (Fig.11c & d).SRdecreases salinity in all seasons and its contribution is most significant during the summer monsoon and autumn inter-monsoon (Fig.11).

At the southern boundary, the seasonal mean meridional advective salt transport (Svm) generated by the background large-scale circulation, the meridional eddy salt flux (Sve) induced by the mesoscale signals and the meridional diff usive salt transport (Sydiff)caused by the sub-grid processes (e.g., small-scale eddies not resolved by the model) aff ect the salt transport.Svmdominates the salt transport during the winter monsoon, balancing the sharp salinity decrease induced by the large freshwater injection in the wet season (Fig.11).Svemakes the salinity increase in all seasons, and the transport is strong (over 1 Sv) during the monsoon seasons (Fig.11). In addition, the magnitude ofSveis sizeable (over 30%, Fig.11)compared with the absolute value of the total meridional advective salt transport (i.e., the sum of the absolute value ofSvmandSve).Sydiffshows the positive diff usive salt transport throughout the year except for the spring inter-monsoon, and become significant in the summer monsoon and the autumn inter-monsoon (Fig.11).

Fig.11 Seasonal salt transports during winter monsoon (a), spring inter-monsoon (b), summer monsoon (c), and autumninter-monsoon (d) in the upper 50 m BoB in 2014

At the bottom boundary, the seasonal mean vertical advective salt transport (Swm) generated by the background large-scale circulation, the vertical eddy salt flux (Swe) induced by the mesoscale signals and the vertical diff usive salt transport (Szdiff) caused by the sub-grid processes aff ect the salt transport.Swmplays a significant role during monsoon seasons,making the salinity increase (decrease) during the summer (winter) monsoon (Fig.11a & c).Swemakes the salinity increase in all seasons except for the summer monsoon. Though the total magnitude ofSweis smaller thanSve, it cannot be neglected (about 10%-30%) compared with the absolute value of the total vertical advective salt transport (i.e., the sum of the absolute value ofSwmandSwe), especially during the summer monsoon and the autumn inter-monsoon(Fig.11).Szdiffshows the positive (negative) diff usive salt transport in the spring (autumn) inter-monsoon and summer (winter) monsoon, and the magnitude ofSzdiffis comparable withSwe(Fig.11).The above results show that the horizontal eddy flux is larger than the vertical eddy flux (Fig.11). One possible reason is the low vertical (～10 m in most regions of the upper 50 m) resolutions of the model outputs, which underestimates the influence of unresolved eddy fluxes. Another possibility is that the KPP scheme cannot describe all the sub-grid processes accurately, such as the mixed layer submesoscale motions, Langmuir circulation and the mixing induced by strong tropical cyclones. Similarly, the vertical diff usive salt flux could also be underestimated due to the lack of vertical resolution of the model.

5 CONCLUSION

In this study, the factors influencing the salinity budget, including eddy salt fluxes, are analyzed in a fixed box model of the upper BoB in 2014, using daily averaged outputs from ROMS. The control volume is taken as the upper 50 m of the BoB, north of 10°N and the three open boundaries are the 10°N section, the sea surface and the 50-m level. The salinity budget indicates that the time series on the two sides of Eq.1 are balanced within the acceptable error range.

The salinity budget analysis demonstrates that the net freshwater flux is the dominant factor leading to the sharp salinity decrease during the summer monsoon. The advections, especially the meridional advective salt flux, are the main factors that cause the salinity to rise steadily after the summer monsoon,primarily aff ected by seasonal changes in circulations.The vertical advective salt flux plays an important role in balancing the salinity change induced by meridional advective salt flux during the monsoon seasons. In the winter monsoon, the downward advective salt flux at 50-m depth balances the northward high salinity flow into the bay at the southern boundary, while the opposite phenomenon exists in the summer monsoon. The downward(upward) vertical advective salt flux is caused by the negative (positive) wind stress curl during the winter(summer) monsoon. Turbulent diff usions can also influence the salinity budget. The meridional diff usive salt flux balances the sharp decrease of salinity during the summer monsoon and increases it in the autumn inter-monsoon and early winter monsoon. The vertical diff usive salt flux is also positive during the summer monsoon.

The eddy salt flux is decomposed to analyze its contribution. The spatial distributions of eddy salt fluxes show the distinctive mesoscale structures in all seasons, which magnitude is sizeable (over 30% for the meridional eddy salt flux and about 10%-30% for the vertical eddy salt flux). The eddy salt flux across the southern boundary displays stronger eddy signals during the monsoon seasons than during the intermonsoon seasons, and is strongly positive near the west boundary where the EICC is located during the winter monsoon and autumn inter-monsoon. The eddy salt flux across the bottom boundary shows abundant eddy-scale signals, and is especially strong along the coast area and around the Andaman and Nicobar Islands. Strong eddy flux signals can be seen during the autumn inter-monsoon because of the tropical cyclonic HudHud causing intensive mixing in the upper ocean.

6 DATA AVAILABILITY STATEMENT

The datasets generated and/or analyzed for the current study are available from the corresponding author.

7 ACKNOWLEDGMENT

The BOA-Argo dataset used in this study is produced by the China Argo Real-time Data Center.These data were collected and made freely available by the International Argo Program and the national programs that contribute to it (http://www.argo.ucsd.edu, http://argo.jcommops.org). The Argo Program is part of the Global Ocean Observing System. The Drifter-Derived Climatology of Global Near-Surface Currents data are from ftp://ftp.aoml.noaa.gov/phod/pub/lumpkin/drifter_climatology. The best track data are from the Joint Typhoon Warning Center (https://www.metoc.navy.mil/jtwc/jtwc.html?north-indian-ocean).

whereSis the salinity,u,vandware the eastward, northward and vertical (upward) velocities, respectively, andκSHandκSVare the horizontal and vertical turbulent diff usivities, respectively.FSdenotes the internal forcing for the salinity variation, such as external sources and sinks of salt within the ocean, which is not shown in the present study.

Turbulent diff usion terms typically refer to the parameterization of sub-grid processes in the model, which contain the stirring eff ects of small-scale eddies not resolved by the model. The horizontal turbulent diff usivity in the ROMS model is chosen as a constant background diff usivity and the numerical scheme for the advection terms carries implicit numerical diff usivity. The vertical turbulent diff usivity follows the K-profile parameterization(KPP) scheme, which is described in detail by Large et al. (1994).

For the volume-integration salinity budget calculation, the domain extends from 10°N to 23°N and 80°E to 99°E, including the BoB and the AS (Fig.1), and the bottom of the study region is set as 50 m. According to previous studies, the mixed layer thickness in the BoB is approximately in the range of 20-50 m (Narvekar and Prasanna Kumar, 2014; Wilson and Riser, 2016; Li et al., 2017; Zhang et al., 2018b; Shee et al., 2019). Therefore,the 50-m depth of the study domain includes the upper mixed layer.

An averaged volume-integration of the salinity budget equation, integrated from the free sea surface (η) to 50 m (h=50 m), is shown below:

This box model is greatly simplified because only the 10°N section, the sea surface and the 50-m depth are open boundaries, and 80°E, 99°E, and 23°N sections are closed boundaries (i.e., no advection occur across these sections). Furthermore, the surface vertical velocity is much smaller than the vertical velocity at 50 m (Fig.A1).We can, therefore, ignore the influence of surface vertical advective salt flux. To eliminate the dependence of the reference to zero salinity value, the advective salt flux into the control volume is referenced to the volumeaveraged salinity of the domain. This method was proposed by Lee et al. (2004) and has been widely used in heat and salt budget analysis (Zhang et al., 2018c; Lee et al., 2019; Soares et al., 2019). Thus, the advection terms in Eq.A2 can be rewritten as the advective salt flux through the interfaces of the box:

where [S] is time varying and volume-averaged salinity of the domain; the subscript E, W, N, S, Sur, and B indicate the eastern, western, north, south, surface (η), and bottom (h=50 m) boundaries, respectively.

Because the boundaries at 80°E, 99°E, and 23°N are closed, and the vertical boundary condition of the turbulent diff usion at the sea surface is

whereρ0is the constant reference density;Eis the evaporation;Pis the precipitation; andRis the river runoff ,the diff usion terms in Eq.A2 can be described as:

Including all the terms above, the original volume-integration salinity budget equation (Eq.A2) can be written as:

Fig.A1 Averaged vertical velocities at the surface (red line) and 50-m depth (blue line) in the BoB in 2014

The average is for the whole study area at the surface and the depth of 50 m.