Measurement results of d18O, Mg/Ca, and Ti/Ca of sediment cores GeoB10042-1 and GeoB10043-3

The advection of relatively fresh Java Sea water through the Sunda Strait is presently responsible for the low-salinity "tongue" in the eastern tropical Indian Ocean with salinities as low as 32 per mil. The evolution of the hydrologic conditions in the eastern tropical Indian Ocean since the last glacial period, when the Sunda shelf was exposed and any advection via the Sunda Strait was cutoff, and the degree to which these conditions were affected by the Sunda Strait opening are not known. Here we have analyzed two sediment cores (GeoB 10042-1 and GeoB 10043-3) collected from the eastern tropical Indian Ocean off the Sunda Strait that cover the past ~40,000?years. We investigate the magnitude of terrigenous supply, sea surface temperature (SST), and seawater d18O (d18Osw) changes related to the sea level-driven opening of the Sunda Strait. Our new spliced records off the Sunda Strait show that during the last glacial, average SST was cooler and d18Osw was higher than elsewhere in the eastern tropical Indian Ocean. Seawater d18O decreased ~0.5 per mil after the opening of the Sunda Strait at ~10 kyr B.P. accompanied by an SST increase of 1.7°C. We suggest that fresher sea surface conditions have persisted ever since due to a continuous transport of low-salinity Java Sea water into the eastern tropical Indian Ocean via the Sunda Strait that additionally increased marine productivity through the concomitant increase in terrigenous supply.

(Table 3) Interannual variability of total mass flux at the deep sediment trap site CBmeso between 1988 and 2006

This article will review major features of the 'giant' Cape Blanc filament off Mauritania with regard to the transport of chlorophyll and organic carbon from the shelf to the open ocean. Within the filament, chlorophyll is transported about 400 km offshore. Modelled particle distributions along a zonal transect at 21°N showed that particles with a sinking velocity of 5 m d**-1 are advected offshore by up to 600 km in subsurface particle clouds generally located between 400 m and 800 m water depth, forming an Intermediate Nepheloid Layer (INL). It corresponds to the depth of the oxygen minimum zone. Heavier particles with a sinking velocity of 30 m d**-1 are transported from the shelf within the Bottom Layer (BL) of more than 1000 m thickness, largely following the topography of the bottom slope. The particles advected within the BL contribute to the enhanced winter-spring mass fluxes collected at the open-ocean mesotrophic sediment trap site CB-13 (200 nm offshore), due to a long distance advection in deeper waters. The lateral contribution to the deep sediment trap in winter-spring is estimated to be 63% and 72% for organic carbon and total mass, respectively, whereas the lateral input for both components on an annual basis is estimated to be in the order of 15%. Biogenic opal increases almost fivefold from the upper to the lower mesotrophic CB-13 trap, also pointing to an additional source for biogenic silica from eutrophic coastal waters. Blooms obviously sink in smaller, probably mesoscale-sized patches with variable settling rates, depending on the type of aggregated particles and their ballast content. Generally, particle sinking rates are exceptionally high off NW Africa. Very high chlorophyll values and a large size of the Cape Blanc filament in 1998-1999 are also documented in enhanced total mass and organic carbon fluxes. An increasing trend in satellite chlorophyll concentrations and the size of the Cape Blanc filament between 1997 and 2008 as observed for other coastal upwelling areas is not documented.

The significance of nitrogen regeneration for new production within a filament of the Mauritanian upwelling system

The Lagrangian progression of a biological community was followed in a filament of the Mauritanian upwelling system, north-west Africa, during offshore advection. The inert dual tracers sulfur hexafluoride and helium-3 labelled a freshly upwelled patch of water that was mapped for 8 days. Changes in biological, physical, and chemical characteristics were measured, including phytoplankton productivity, nitrogen assimilation, and regeneration. Freshly upwelled water contained high nutrient concentrations but was depleted in N compared to Redfield stoichiometry. The highest rate of primary productivity was measured on the continental shelf, associated with high rates of nitrogen assimilation and a phytoplankton community dominated by diatoms and flagellates. Indicators of phytoplankton abundance and activity decreased as the labelled water mass transited the continental shelf slope into deeper water, possibly linked to the mixed layer depth exceeding the light penetration depth. By the end of the study, the primary productivity rate decreased and was associated with lower rates of nitrogen assimilation and lower nutrient concentrations. Nitrogen regeneration and assimilation took place simultaneously. Results highlighted the importance of regenerated NHC 4 in sustaining phytoplankton productivity and indicate that the upwelled NO3 pool contained an increasing fraction of regenerated NO3 as it advected offshore. By calculating this fraction and incorporating it into an f ratio formulation, we estimated that of the 12:38Tg C of annual regional production, 4:73Tg C was exportable.

The significance of nitrogen regeneration for new production within a filament of the Mauritanian upwelling system

The Lagrangian progression of a biological community was followed in a filament of the Mauritanian upwelling system, north-west Africa, during offshore advection. The inert dual tracers sulfur hexafluoride and helium-3 labelled a freshly upwelled patch of water that was mapped for 8 days. Changes in biological, physical, and chemical characteristics were measured, including phytoplankton productivity, nitrogen assimilation, and regeneration. Freshly upwelled water contained high nutrient concentrations but was depleted in N compared to Redfield stoichiometry. The highest rate of primary productivity was measured on the continental shelf, associated with high rates of nitrogen assimilation and a phytoplankton community dominated by diatoms and flagellates. Indicators of phytoplankton abundance and activity decreased as the labelled water mass transited the continental shelf slope into deeper water, possibly linked to the mixed layer depth exceeding the light penetration depth. By the end of the study, the primary productivity rate decreased and was associated with lower rates of nitrogen assimilation and lower nutrient concentrations. Nitrogen regeneration and assimilation took place simultaneously. Results highlighted the importance of regenerated NHC 4 in sustaining phytoplankton productivity and indicate that the upwelled NO3 pool contained an increasing fraction of regenerated NO3 as it advected offshore. By calculating this fraction and incorporating it into an f ratio formulation, we estimated that of the 12:38Tg C of annual regional production, 4:73Tg C was exportable.

Individual-based modelling of the development and transport of a Karenia mikimotoi bloom on the North-West European continental shelf

In 2006, a large and prolonged bloom of the dinoflagellate Karenia mikimotoi occurred in Scottish coastal waters, causing extensive mortalities of benthic organisms including annelids and molluscs and some species of fish ( Davidson et al., 2009). A coupled hydrodynamic-algal transport model was developed to track the progression of the bloom around the Scottish coast during June–September 2006 and hence investigate the processes controlling the bloom dynamics. Within this individual-based model, cells were capable of growth, mortality and phototaxis and were transported by physical processes of advection and turbulent diffusion, using current velocities extracted from operational simulations of the MRCS ocean circulation model of the North-west European continental shelf. Vertical and horizontal turbulent diffusion of cells are treated using a random walk approach. Comparison of model output with remotely sensed chlorophyll concentrations and cell counts from coastal monitoring stations indicated that it was necessary to include multiple spatially distinct seed populations of K. mikimotoi at separate locations on the shelf edge to capture the qualitative pattern of bloom transport and development. We interpret this as indicating that the source population was being transported northwards by the Hebridean slope current from where colonies of K. mikimotoi were injected onto the continental shelf by eddies or other transient exchange processes. The model was used to investigate the effects on simulated K. mikimotoi transport and dispersal of: (1) the distribution of the initial seed population; (2) algal growth and mortality; (3) water temperature; (4) the vertical movement of particles by diurnal migration and eddy diffusion; (5) the relative role of the shelf edge and coastal currents; (6) the role of wind forcing. The numerical experiments emphasized the requirement for a physiologically based biological model and indicated that improved modelling of future blooms will potentially benefit from better parameterisation of temperature dependence of both growth and mortality and finer spatial and temporal hydrodynamic resolution.

Individual-based modelling of the development and transport of a Karenia mikimotoi bloom on the North-West European continental shelf

In 2006, a large and prolonged bloom of the dinoflagellate Karenia mikimotoi occurred in Scottish coastal waters, causing extensive mortalities of benthic organisms including annelids and molluscs and some species of fish ( Davidson et al., 2009). A coupled hydrodynamic-algal transport model was developed to track the progression of the bloom around the Scottish coast during June–September 2006 and hence investigate the processes controlling the bloom dynamics. Within this individual-based model, cells were capable of growth, mortality and phototaxis and were transported by physical processes of advection and turbulent diffusion, using current velocities extracted from operational simulations of the MRCS ocean circulation model of the North-west European continental shelf. Vertical and horizontal turbulent diffusion of cells are treated using a random walk approach. Comparison of model output with remotely sensed chlorophyll concentrations and cell counts from coastal monitoring stations indicated that it was necessary to include multiple spatially distinct seed populations of K. mikimotoi at separate locations on the shelf edge to capture the qualitative pattern of bloom transport and development. We interpret this as indicating that the source population was being transported northwards by the Hebridean slope current from where colonies of K. mikimotoi were injected onto the continental shelf by eddies or other transient exchange processes. The model was used to investigate the effects on simulated K. mikimotoi transport and dispersal of: (1) the distribution of the initial seed population; (2) algal growth and mortality; (3) water temperature; (4) the vertical movement of particles by diurnal migration and eddy diffusion; (5) the relative role of the shelf edge and coastal currents; (6) the role of wind forcing. The numerical experiments emphasized the requirement for a physiologically based biological model and indicated that improved modelling of future blooms will potentially benefit from better parameterisation of temperature dependence of both growth and mortality and finer spatial and temporal hydrodynamic resolution.

Surface water subduction during a downwelling event in a semi-enclosed bay

The Ría de Vigo is a bay strongly influenced by upwelling-downwelling cycles along the adjacent coast of NW Iberia. Moored and ship-board observations during September 2006 showed that subduction, initially associated with an estuarine circulation, strengthened when a strong downwelling circulation, resulting from northward wind over the coastal ocean, was generated in the outer ría causing ambient waters to be advected outward in the lower layer. Incoming surface waters confined the estuarine circulation to the shallow interior and displaced isopleths downward through the water column at ∼10 m d−1. As the estuarine circulation retreated inward, strong flow convergence developed between middle and inner ria in the layer above 15 m, while divergence developed beneath. The convergence increased through the period of downwelling-favorable wind at a rate consistent with the observed isopleth displacement velocities. The coefficient of turbulent diffusion Kt, from a microstructure profiler, indicated that mixing was strong in the estuarine circulation and subsequently in the downwelling zone, where localized instabilities and temperature-salinity inversions were observed. During the downwelling, concentrations of phytoplankton, including potentially harmful species, increased, especially in the middle and inner ria, as a result of inward advection, subduction and the ability of the dinoflagellates to maintain their position in the water column by swimming. In the course of the 5 day event, the water mass of all but the innermost ría was flushed completely and replaced by waters originating in the coastally-trapped poleward flow along the Atlantic coastline.

Surface water subduction during a downwelling event in a semi-enclosed bay

The Ría de Vigo is a bay strongly influenced by upwelling-downwelling cycles along the adjacent coast of NW Iberia. Moored and ship-board observations during September 2006 showed that subduction, initially associated with an estuarine circulation, strengthened when a strong downwelling circulation, resulting from northward wind over the coastal ocean, was generated in the outer ría causing ambient waters to be advected outward in the lower layer. Incoming surface waters confined the estuarine circulation to the shallow interior and displaced isopleths downward through the water column at ∼10 m d−1. As the estuarine circulation retreated inward, strong flow convergence developed between middle and inner ria in the layer above 15 m, while divergence developed beneath. The convergence increased through the period of downwelling-favorable wind at a rate consistent with the observed isopleth displacement velocities. The coefficient of turbulent diffusion Kt, from a microstructure profiler, indicated that mixing was strong in the estuarine circulation and subsequently in the downwelling zone, where localized instabilities and temperature-salinity inversions were observed. During the downwelling, concentrations of phytoplankton, including potentially harmful species, increased, especially in the middle and inner ria, as a result of inward advection, subduction and the ability of the dinoflagellates to maintain their position in the water column by swimming. In the course of the 5 day event, the water mass of all but the innermost ría was flushed completely and replaced by waters originating in the coastally-trapped poleward flow along the Atlantic coastline.

Progenitor-Dependent Explosion Dynamics in Self-Consistent, Axisymmetric Simulations of Neutrino-Driven Core-Collapse Supernovae

We present self-consistent, axisymmetric core-collapse supernova simulations performed with the Prometheus-Vertex code for 18 pre-supernova models in the range of 11–28 M ⊙, including progenitors recently investigated by other groups. All models develop explosions, but depending on the progenitor structure, they can be divided into two classes. With a steep density decline at the Si/Si–O interface, the arrival of this interface at the shock front leads to a sudden drop of the mass-accretion rate, triggering a rapid approach to explosion. With a more gradually decreasing accretion rate, it takes longer for the neutrino heating to overcome the accretion ram pressure and explosions set in later. Early explosions are facilitated by high mass-accretion rates after bounce and correspondingly high neutrino luminosities combined with a pronounced drop of the accretion rate and ram pressure at the Si/Si–O interface. Because of rapidly shrinking neutron star radii and receding shock fronts after the passage through their maxima, our models exhibit short advection timescales, which favor the efficient growth of the standing accretion-shock instability. The latter plays a supportive role at least for the initiation of the re-expansion of the stalled shock before runaway. Taking into account the effects of turbulent pressure in the gain layer, we derive a generalized condition for the critical neutrino luminosity that captures the explosion behavior of all models very well. We validate the robustness of our findings by testing the influence of stochasticity, numerical resolution, and approximations in some aspects of the microphysics.

The Roles of Land and Orography on Precipitation and Ocean Circulation in Global Climate Models

Thesis (Ph.D.)--University of Washington, 2016-08

Investigating Marine Boundary Layer Aerosol Budgets and Variability

Thesis (Master's)--University of Washington, 2016-08

Séquestration géologique du CO2 par carbonatation minérale dans les résidus miniers

La carbonatation minérale dans les résidus miniers est un moyen sûr et permanent de séquestrer le CO2 atmosphérique. C’est un processus naturel et passif qui ne nécessite aucun traitement particulier et donc avantageux d’un point de vue économique. Bien que la quantité de CO2 qu’il soit possible de séquestrer selon ce processus est faible à l’échelle globale, dans le cadre d’un marché du carbone, les entreprises minières pourraient obtenir des crédits et ainsi revaloriser leurs résidus. À l’heure actuelle, il y a peu d’informations pour quantifier le potentiel de séquestration du CO2 de façon naturelle et passive dans les piles de résidus miniers. Il est donc nécessaire d’étudier le phénomène pour comprendre comment évolue la réaction à travers le temps et estimer la quantité de CO2 qui peut être séquestrée naturellement dans les piles de résidus. Plusieurs travaux de recherche se sont intéressés aux résidus miniers de Thetford Mines (Québec, Canada), avec une approche principalement expérimentale en laboratoire. Ces travaux ont permis d’améliorer la compréhension du processus de carbonatation, mais ils nécessitent une validation à plus grande échelle sous des conditions atmosphériques réelles. L’objectif général de cette étude est de quantifier le processus de carbonatation minérale des résidus miniers sous des conditions naturelles, afin d’estimer la quantité de CO2 pouvant être piégée par ce processus. La méthodologie utilisée repose sur la construction de deux parcelles expérimentales de résidus miniers situées dans l’enceinte de la mine Black Lake (Thetford Mines). Les résidus miniers sont principalement constitués de grains et de fibres de chrysotile et lizardite mal triés, avec de petites quantités d’antigorite, de brucite et de magnétite. Des observations spatiales et temporelles ont été effectuées dans les parcelles concernant la composition et la pression des gaz, la température des résidus, la teneur en eau volumique, la composition minérale des résidus ainsi que la chimie de l’eau des précipitations et des lixiviats provenant des parcelles. Ces travaux ont permis d’observer un appauvrissement notable du CO2 dans les gaz des parcelles (< 50 ppm) ainsi que la précipitation d’hydromagnésite dans les résidus, ce qui suggère que la carbonatation minérale naturelle et passive est un processus potentiellement important dans les résidus miniers. Après 4 ans d’observations, le taux de séquestration du CO2 dans les parcelles expérimentales a été estimé entre 3,5 et 4 kg/m3/an. Ces observations ont permis de développer un modèle conceptuel de la carbonatation minérale naturelle et passive dans les parcelles expérimentales. Dans ce modèle conceptuel, le CO2 atmosphérique (~ 400 ppm) se dissout dans l’eau hygroscopique contenue dans les parcelles, où l’altération des silicates de magnésium forme des carbonates de magnésium. La saturation en eau dans les cellules est relativement stable dans le temps et varie entre 0,4 et 0,65, ce qui est plus élevé que les valeurs de saturation optimales proposées dans la littérature, réduisant ainsi le transport de CO2 dans la zone non saturée. Les concentrations de CO2 en phase gazeuse, ainsi que des mesures de la vitesse d’écoulement du gaz dans les cellules suggèrent que la réaction est plus active près de la surface et que la diffusion du CO2 est le mécanisme de transport dominant dans les résidus. Un modèle numérique a été utilisé pour simuler ces processus couplés et valider le modèle conceptuel avec les observations de terrain. Le modèle de transport réactif multiphase et multicomposant MIN3P a été utilisé pour réaliser des simulations en 1D qui comprennent l’infiltration d’eau à travers le milieu partiellement saturé, la diffusion du gaz, et le transport de masse réactif par advection et dispersion. Même si les écoulements et le contenu du lixivat simulés sont assez proches des observations de terrain, le taux de séquestration simulé est 22 fois plus faible que celui mesuré. Dans les simulations, les carbonates précipitent principalement dans la partie supérieure de la parcelle, près de la surface, alors qu’ils ont été observés dans toute la parcelle. Cette différence importante pourrait être expliquée par un apport insuffisant de CO2 dans la parcelle, qui serait le facteur limitant la carbonatation. En effet, l’advection des gaz n’a pas été considérée dans les simulations et seule la diffusion moléculaire a été simulée. En effet, la mobilité des gaz engendrée par les fluctuations de pression barométrique et l’infiltration de l’eau, ainsi que l’effet du vent doivent jouer un rôle conséquent pour alimenter les parcelles en CO2.

DIRECT NUMERICAL SIMULATION OF INCOMPRESSIBLE MULTIPHASE FLOW WITH PHASE CHANGE

Phase change problems arise in many practical applications such as air-conditioning and refrigeration, thermal energy storage systems and thermal management of electronic devices. The physical phenomenon in such applications are complex and are often difficult to be studied in detail with the help of only experimental techniques. The efforts to improve computational techniques for analyzing two-phase flow problems with phase change are therefore gaining momentum. The development of numerical methods for multiphase flow has been motivated generally by the need to account more accurately for (a) large topological changes such as phase breakup and merging, (b) sharp representation of the interface and its discontinuous properties and (c) accurate and mass conserving motion of the interface. In addition to these considerations, numerical simulation of multiphase flow with phase change introduces additional challenges related to discontinuities in the velocity and the temperature fields. Moreover, the velocity field is no longer divergence free. For phase change problems, the focus of developmental efforts has thus been on numerically attaining a proper conservation of energy across the interface in addition to the accurate treatment of fluxes of mass and momentum conservation as well as the associated interface advection. Among the initial efforts related to the simulation of bubble growth in film boiling applications the work in \cite{Welch1995} was based on the interface tracking method using a moving unstructured mesh. That study considered moderate interfacial deformations. A similar problem was subsequently studied using moving, boundary fitted grids \cite{Son1997}, again for regimes of relatively small topological changes. A hybrid interface tracking method with a moving interface grid overlapping a static Eulerian grid was developed \cite{Juric1998} for the computation of a range of phase change problems including, three-dimensional film boiling \cite{esmaeeli2004computations}, multimode two-dimensional pool boiling \cite{Esmaeeli2004} and film boiling on horizontal cylinders \cite{Esmaeeli2004a}. The handling of interface merging and pinch off however remains a challenge with methods that explicitly track the interface. As large topological changes are crucial for phase change problems, attention has turned in recent years to front capturing methods utilizing implicit interfaces that are more effective in treating complex interface deformations. The VOF (Volume of Fluid) method was adopted in \cite{Welch2000} to simulate the one-dimensional Stefan problem and the two-dimensional film boiling problem. The approach employed a specific model for mass transfer across the interface involving a mass source term within cells containing the interface. This VOF based approach was further coupled with the level set method in \cite{Son1998}, employing a smeared-out Heaviside function to avoid the numerical instability related to the source term. The coupled level set, volume of fluid method and the diffused interface approach was used for film boiling with water and R134a at the near critical pressure condition \cite{Tomar2005}. The effect of superheat and saturation pressure on the frequency of bubble formation were analyzed with this approach. The work in \cite{Gibou2007} used the ghost fluid and the level set methods for phase change simulations. A similar approach was adopted in \cite{Son2008} to study various boiling problems including three-dimensional film boiling on a horizontal cylinder, nucleate boiling in microcavity \cite{lee2010numerical} and flow boiling in a finned microchannel \cite{lee2012direct}. The work in \cite{tanguy2007level} also used the ghost fluid method and proposed an improved algorithm based on enforcing continuity and divergence-free condition for the extended velocity field. The work in \cite{sato2013sharp} employed a multiphase model based on volume fraction with interface sharpening scheme and derived a phase change model based on local interface area and mass flux. Among the front capturing methods, sharp interface methods have been found to be particularly effective both for implementing sharp jumps and for resolving the interfacial velocity field. However, sharp velocity jumps render the solution susceptible to erroneous oscillations in pressure and also lead to spurious interface velocities. To implement phase change, the work in \cite{Hardt2008} employed point mass source terms derived from a physical basis for the evaporating mass flux. To avoid numerical instability, the authors smeared the mass source by solving a pseudo time-step diffusion equation. This measure however led to mass conservation issues due to non-symmetric integration over the distributed mass source region. The problem of spurious pressure oscillations related to point mass sources was also investigated by \cite{Schlottke2008}. Although their method is based on the VOF, the large pressure peaks associated with sharp mass source was observed to be similar to that for the interface tracking method. Such spurious fluctuation in pressure are essentially undesirable because the effect is globally transmitted in incompressible flow. Hence, the pressure field formation due to phase change need to be implemented with greater accuracy than is reported in current literature. The accuracy of interface advection in the presence of interfacial mass flux (mass flux conservation) has been discussed in \cite{tanguy2007level,tanguy2014benchmarks}. The authors found that the method of extending one phase velocity to entire domain suggested by Nguyen et al. in \cite{nguyen2001boundary} suffers from a lack of mass flux conservation when the density difference is high. To improve the solution, the authors impose a divergence-free condition for the extended velocity field by solving a constant coefficient Poisson equation. The approach has shown good results with enclosed bubble or droplet but is not general for more complex flow and requires additional solution of the linear system of equations. In current thesis, an improved approach that addresses both the numerical oscillation of pressure and the spurious interface velocity field is presented by featuring (i) continuous velocity and density fields within a thin interfacial region and (ii) temporal velocity correction steps to avoid unphysical pressure source term. Also I propose a general (iii) mass flux projection correction for improved mass flux conservation. The pressure and the temperature gradient jump condition are treated sharply. A series of one-dimensional and two-dimensional problems are solved to verify the performance of the new algorithm. Two-dimensional and cylindrical film boiling problems are also demonstrated and show good qualitative agreement with the experimental observations and heat transfer correlations. Finally, a study on Taylor bubble flow with heat transfer and phase change in a small vertical tube in axisymmetric coordinates is carried out using the new multiphase, phase change method.

Um problema de dispersão de poluentes em rios e canais por meio do método de separação de variáveis

Neste trabalho obtém-se uma solução analítica para a equação de advecção-difusão aplicada a problemas de dispersão de poluentes em rios e canais. Para tanto, consideram-se os casos unidimensionais e bidimensionais em regime transiente com coeficientes de difusividade e velocidades constantes. A abordagem utilizada para a resolução deste problema é o método de Separação de Variáveis. Os modelos resolvidos foram simulados utilizando o MatLab. Apresentam-se os resultados das simulações numéricas em formato gráfico. Os resultados de algumas simulações numéricas existem na literatura e puderam ser comparados. O modelo proposto mostrou-se coerente em relação aos dados considerados. Para outras simulações não foram encontrados comparativos na literatura, todavia esses problemas governados por equações diferenciais parciais, mesmo lineares, não são de fácil solução analítica. Sendo que, muitas delas representam importantes problemas de matemática e física, com diversas aplicações na engenharia. Dessa forma, é de grande importância a disponibilidade de um maior número de problemas-teste para avaliação de desempenho de formulações numéricas, cada vez mais eficazes, já que soluções analíticas oferecem uma base mais segura para comparação de resultados.

Mixed-layer salinity budget in the tropical Indian Ocean: seasonal cycle based only on observations

The mixed-layer salinity (MLS) budget in the tropical Indian Ocean is estimated from a combination of satellite products and in situ observations over the 2004-2012 period, to investigate the mechanisms controlling the seasonal MLS variability. In contrast with previous studies in the tropical Indian Ocean, our results reveal that the coverage, resolution, and quality of available observations are now sufficient to approach a closed monthly climatology seasonal salt budget. In the South-central Arabian Sea and South-western Tropical Indian Ocean (SCAS and STIO, respectively), where seasonal variability of the MLS is pronounced, the monthly MLS tendency terms are well captured by the diagnostic. In the SCAS region, in agreement with previous results, the seasonal cycle of the MLS is mainly due to meridional advection driven by the monsoon winds. In the STIO, contrasting previous results indicating the control of the meridional advection over the seasonal MLS budget, our results reveal the leading role of the freshwater flux due to precipitation.