ALBERT

Description: The northern Arabian Sea hosts a winter chlorophyll bloom, triggered by convective overturning in response to cold and dry northeasterly monsoon winds. There is currently no consensus about the processes responsible for the interannual variations of the magnitude of this bloom. The current study aims at identifying these processes using both observations (including remotely sensed chlorophyll data and physical parameters derived from Argo data) and a coupled biophysical ocean model simulation. Six different remotely sensed chlorophyll products are compared and they show a good phase agreement at seasonal and interannual timescales, but significant discrepancies in bloom amplitude. Both model and observations indicate that the interannual fluctuations of the winter bloom amplitude are strongly tied to mixed layer depth interannual anomalies (correlation ~ 0.6 to 0.7), which are themselves controlled by the net heat flux at the air-sea interface. Our results suggest that mixed layer depth control of the bloom amplitude ensues from the modulation of nutrient entrainment into the mixed layer. In contrast, our results show insignificant correlations between the bloom amplitude and thermocline depth, which precludes a control of the bloom amplitude by daily dilution down to the thermocline depth, as suggested in a previous study.

Description: The principal objective of hydrodynamical-biological models is to provide estimates of the main carbon fluxes such as total and export oceanic production. These models are nitrogen based, that is to say that the variables are expressed in terms of their nitrogen content. Moreover models are calibrated using chlorophyll data sets. Therefore carbon to chlorophyll (C:Chl) and carbon to nitrogen (C:N) ratios have to be assumed. This paper addresses the problem of the representation of these ratios. In a 1D framework at the DYFAMED station (NW Mediterranean Sea) we propose a model which enables the estimation of the basic biogeochemical fluxes and in which the spatio-temporal variability of the C:Chl and C:N ratios is fully represented in a mechanical way. This is achieved through the introduction of new state variables coming from the embedding of a phytoplankton growth model in a more classical Redfieldian NNPZD-DOM model (in which the C:N ratio is assumed to be a constant). Following this modelling step, the parameters of the model are estimated using the adjoint data assimilation method which enables the assimilation of chlorophyll and nitrate data sets collected at DYFAMED in 1997.Comparing the predictions of the new Mechanistic model with those of the classical Redfieldian NNPZD-DOM model which was calibrated with the same data sets, we find that both models reproduce the reference data in a comparable manner. Both fluxes and stocks can be equally well predicted by either model. However if the models are coinciding on an average basis, they are diverging from a variability prediction point of view. In the Mechanistic model biology adapts much faster to its environment giving rise to higher short term variations. Moreover the seasonal variability in total production differs from the Redfieldian NNPZD-DOM model to the Mechanistic model. In summer the Mechanistic model predicts higher production values in carbon unit than the Redfieldian NNPZD-DOM model. In winter the contrary holds.

Description: We evaluate the performance of the Community Atmosphere Model's (CAM) spectral element method on variable-resolution grids using the shallow-water equations in spherical geometry. We configure the method as it is used in CAM, with dissipation of grid scale variance, implemented using hyperviscosity. Hyperviscosity is highly scale selective and grid independent, but does require a resolution-dependent coefficient. For the spectral element method with variable-resolution grids and highly distorted elements, we obtain the best results if we introduce a tensor-based hyperviscosity with tensor coefficients tied to the eigenvalues of the local element metric tensor. The tensor hyperviscosity is constructed so that, for regions of uniform resolution, it matches the traditional constant-coefficient hyperviscosity. With the tensor hyperviscosity, the large-scale solution is almost completely unaffected by the presence of grid refinement. This later point is important for climate applications in which long term climatological averages can be imprinted by stationary inhomogeneities in the truncation error. We also evaluate the robustness of the approach with respect to grid quality by considering unstructured conforming quadrilateral grids generated with a well-known grid-generating toolkit and grids generated by SQuadGen, a new open source alternative which produces lower valence nodes.

Description: In-situ observations indicate that the Arabian Sea oxygen minimum zone (OMZ) is only weakly influenced by the strong seasonal cycle of ocean dynamic and biogeochemistry forced by the asian monsoon system and it is spatially decorrelated from the coastal upwelling systems where the biological production is the strongest. In this study we examine the factors controlling the seasonality and the spatial distribution of the OMZ in the Arabian Sea using a coupled bio-physical model. We find that the oxygen concentration in the OMZ displays a seasonal cycle with an amplitude of 5–15 % of the annual mean oxygen concentration. The OMZ is ventilated by lateral ventilation along the western boundary current and in the coastal undercurrent along India during the summer monsoon and by coastal downwelling and negative Ekman pumping during the fall intermonsoon and winter monsoon. This ventilation is counterbalanced by strong coastal upwelling and positive Ekman pumping of low oxygen waters at the base of the OMZ during the spring intermonsoon. Although the factors controlling the OMZ seasonality are associated with the men circulation, we find that mesoscale dynamics modulates them by limiting the vertical ventilation during winter and enhancing it through lateral advection during the rest of the year. Processes explaining the establishment and spatial distribution of the OMZ were quantified using a perturbation experiment initialised with no OMZ. As expected, the oxygen depletion is triggered by strong biological activity in central Arabian Sea during winter and in western and eastern boundary coastal upwelling systems during summer. We find that the 3-D ocean dynamic largely controls the spatial distribution of the OMZ. The eastward shift ensues from the northward lateral transport of ventilated waters along the western and eastern coasts and the advection offshore of low oxygen waters formed in the upwelling system.

Description: The expansion of OMZs (oxygen minimum zones) due to climate change and their possible evolution and impacts on the ecosystems and the atmosphere are still debated, mostly because of the unability of global climate models to adequatly reproduce the processes governing OMZs. In this study, we examine the factors controlling the oxygen budget, i.e. the equilibrium between oxygen sources and sinks in the northern Arabian Sea OMZ using an eddy-resolving biophysical model. Our model confirms that the biological consumption of oxygen is most intense below the region of highest productivity in the western Arabian Sea. The oxygen drawdown in this region is counterbalanced by the large supply of oxygenated waters originated from the south and advected horizontally by the western boundary current. Although the biological sink and the dynamical sources of oxygen compensate on annual average, we find that the seasonality of the dynamical transport of oxygen is 3 to 5 times larger than the seasonality of the biological sink. In agreement with previous findings, the resulting seasonality of oxygen concentration in the OMZ is relatively weak, with a variability of the order of 15% of the annual mean oxygen concentration in the oxycline and 5% elsewhere. This seasonality primarily arises from the vertical displacement of the OMZ forced by the monsoonal reversal of Ekman pumping across the basin. In coastal areas, the oxygen concentration is also modulated seasonally by lateral advection. Along the western coast of the Arabian Sea, the Somali Current transports oxygen-rich waters originated from the south during summer and oxygen-poor waters from the northeast during winter. Along the eastern coast of the Arabian Sea, we find that the main contributor to lateral advection in the OMZ is the Indian coastal undercurrent that advects southern oxygenated waters during summer and northern low-oxygen waters during winter. In this region, our model indicates that oxygen concentrations are modulated seasonally by coastal Kelvin waves and westward-propagating Rossby waves. Whereas on seasonal time scales the sources and sinks of oxygen are dominated by the mean vertical and lateral advection (Ekman pumping and monsoonal currents), on annual time scales we find that the biological sink is counterbalanced by the supply of oxygen sustained by mesoscale structures (eddies and filaments). Eddy-driven advection hence promotes the vertical supply of oxygen along the western coast of the Arabian Sea and the lateral transport of ventilated waters offshore the coast of Oman and southwest India.

Description: Climate simulations codes, such as the Community Earth System Model (CESM), are especially complex and continually evolving. Their on-going state of development requires frequent software verification in the form of quality assurance to both preserve the quality of the code and instill model confidence. To formalize and simplify this previously subjective and computationally-expensive aspect of the verification process, we have developed a new tool for evaluating climate consistency. Because an ensemble of simulations allows us to gauge the natural variability of the model's climate, our new tool uses an ensemble approach for consistency testing. In particular, an ensemble of CESM climate runs is created, from which we obtain a statistical distribution that can be used to determine whether a new climate run is statistically distinguishable from the original ensemble. The CESM Ensemble Consistency Test, referred to as CESM-ECT, is objective in nature and accessible to CESM developers and users. The tool has proven its utility in detecting errors in software and hardware environments and providing rapid feedback to model developers.

Description: We evaluate the performance of the Community Atmosphere Model's (CAM) spectral element method on variable resolution grids using the shallow water equations in spherical geometry. We configure the method as it is used in CAM, with dissipation of grid scale variance implemented using hyperviscosity. Hyperviscosity is highly scale selective and grid independent, but does require a resolution dependent coefficient. For the spectral element method with variable resolution grids and highly distorted elements, we obtain the best results if we introduce a tensor-based hyperviscosity with tensor coefficients tied to the eigenvalues of the local element metric tensor. The tensor hyperviscosity is constructed so that for regions of uniform resolution it matches the traditional constant coefficient hyperviscsosity. With the tensor hyperviscosity the large scale solution is almost completely unaffected by the presence of grid refinement. This later point is important for climate applications where long term climatological averages can be imprinted by stationary inhomogeneities in the truncation error. We also evaluate the robustness of the approach with respect to grid quality by considering unstructured conforming quadrilateral grids generated with a well-known grid-generating toolkit and grids generated by SQuadGen, a new open source alternative which produces lower valence nodes.

Description: Climate simulation codes, such as the Community Earth System Model (CESM), are especially complex and continually evolving. Their ongoing state of development requires frequent software verification in the form of quality assurance to both preserve the quality of the code and instill model confidence. To formalize and simplify this previously subjective and computationally expensive aspect of the verification process, we have developed a new tool for evaluating climate consistency. Because an ensemble of simulations allows us to gauge the natural variability of the model's climate, our new tool uses an ensemble approach for consistency testing. In particular, an ensemble of CESM climate runs is created, from which we obtain a statistical distribution that can be used to determine whether a new climate run is statistically distinguishable from the original ensemble. The CESM ensemble consistency test, referred to as CESM-ECT, is objective in nature and accessible to CESM developers and users. The tool has proven its utility in detecting errors in software and hardware environments and providing rapid feedback to model developers.