The biogeochemical cycling in oxygen-minimum zones (OMZs) is dominated by the interactions of microbial nitrogen transformations and, as recently observed in the Chilean upwelling system, also through the energetically less favorable remineralization of sulfate reduction. The latter process is masked, however, by rapid sulfide oxidation, most likely through nitrate reduction. Thus, the cryptic sulfur cycle links with the nitrogen cycle in OMZ settings. Here, we model the physical-chemical water column structure and the observed process rates as driven by formation and sinking of organic detritus, to quantify the nitrogen and sulfur cycles in the Chilean OMZ. A new biogeochemical submodule was developed and coupled to the Regional Ocean Model System (ROMS). The model results generally agree with the observed distribution of reactive species and the measured process rates. Modeled heterotrophic nitrate reduction and sulfate reduction are responsible for 47% and 36%, respectively, of organic remineralization in a 150 m deep zone below mixed layer. Anammox contributes to 61% of the fixed nitrogen lost to N2 gas, while the rest of the loss is through canonical denitrification as a combination of organic matter oxidation by nitrite reduction and sulfide-driven denitrification. Mineralization coupled to heterotrophic nitrate reduction supplies ∼48% of the ammonium required by anammox. Due to active sulfate reduction, model results suggest that sulfide-driven denitrification contributes to 36% of the nitrogen loss as N2 gas. Our model results highlight the importance of considering the coupled nitrogen and sulfur cycle in examining open-ocean anoxic processes under present, past, and future conditions.

In the Arabian Sea (AS), spatiotemporal nutrient limitation patterns of primary production and the possible role of nutrient inputs from the atmosphere are still not well understood. Using a biogeochemical model forced by modeled aerosol deposition, we show that without high atmospheric iron inputs through dust deposition during the summer monsoon, primary production over the AS would be reduced by half. Atmospheric iron deposition also supports most of the nitrogen fixation over the AS. However, our ocean biogeochemistry modeling results suggest that dinitrogen fixation constitutes a negligible fraction of the primary production. Finally, we show that atmospheric inputs of nitrogen, mostly from anthropogenic activities in India, have a negligible impact on primary production.

The combination of high primary productivity and weak ventilation in the Arabian Sea (AS) and Bay of Bengal (BoB) generates vast areas of depleted oxygen, known as oxygen minimum zones (OMZs). The AS OMZ is the world's thickest and hosts up to 40% of global denitrification. In contrast, the OMZ in the BoB is weaker and denitrification free. Using a series of model simulations, we show that the deeper remineralization depth (RD) in the BoB, potentially associated with organic matter aggregation with riverine mineral particles, contributes to weaken its OMZ. When the RD is set uniformly across both seas, the model fails to reproduce the observed contrast between the two OMZs, irrespective of the chosen RD. In contrast, when the RD is allowed to vary spatially, the contrasting distributions of oxygen and nitrate are correctly reproduced, and water column denitrification is simulated exclusively in the AS, in agreement with observations.

Abstract. The Arabian Sea (AS) hosts one of the most intense oxygen minimum zones (OMZs) in the world. Observations suggest a decline in O2 in the northern AS over the recent decades accompanied by an intensification of the suboxic conditions there. Over the same period, the local sea surface temperature has risen significantly, particularly over the Arabian Gulf (also known as Persian Gulf, hereafter the Gulf), while summer monsoon winds may have intensified. Here, we simulate the evolution of dissolved oxygen in the AS from 1982 through 2010 and explore its controlling factors, with a focus on changing atmospheric conditions. To this end, we use a set of eddy-resolving hindcast simulations forced with winds and heat and freshwater fluxes from an atmospheric reanalysis. We find a significant deoxygenation in the northern AS, with O2 inventories north of 20∘ N dropping by over 6 % per decade between 100 and 1000 m. These changes cause an expansion of the OMZ volume north of 20∘ N at a rate of 0.6 % per decade as well as an increase in the volume of suboxia and the rate of denitrification by 14 and 15 % per decade, respectively. We also show that strong interannual and decadal variability modulate dissolved oxygen in the northern AS, with most of the O2 decline taking place in the 1980s and 1990s. Using a set of sensitivity simulations we demonstrate that deoxygenation in the northern AS is essentially caused by reduced ventilation induced by the recent fast warming of the sea surface, including in the Gulf, with a contribution from concomitant summer monsoon wind intensification. This is because, on the one hand, surface warming enhances vertical stratification and increases Gulf water buoyancy, thus inhibiting vertical mixing and ventilation of the thermocline. On the other hand, summer monsoon wind intensification causes a rise in the thermocline depth in the northern AS that lowers O2 levels in the upper ocean. Our findings confirm that the AS OMZ is strongly sensitive to upper-ocean warming and concurrent changes in the Indian monsoon winds. Finally, our results also demonstrate that changes in the local climatic forcing play a key role in regional dissolved oxygen changes and hence need to be properly represented in global models to reduce uncertainties in future projections of deoxygenation.