Seasonal overturning circulation in the Red Sea : 2. Winter circulation

Author Posting. ©0American Geophysical Union, 2014. This article is posted here by permission of [American Geophysical Union for personal use, not for redistribution. The definitive version was published in Journal of Geophysical Research: Oceans 119 (2014): 2263–2289, doi:10.1002/2013JC009331.The shallow winter overturning circulation in the Red Sea is studied using a 50 year high-resolution MITgcm (MIT general circulation model) simulation with realistic atmospheric forcing. The overturning circulation for a typical year, represented by 1980, and the climatological mean are analyzed using model output to delineate the three-dimensional structure and to investigate the underlying dynamical mechanisms. The horizontal model circulation in the winter of 1980 is dominated by energetic eddies. The climatological model mean results suggest that the surface inflow intensifies in a western boundary current in the southern Red Sea that switches to an eastern boundary current north of 24°N. The overturning is accomplished through a cyclonic recirculation and a cross-basin overturning circulation in the northern Red Sea, with major sinking occurring along a narrow band of width about 20 km along the eastern boundary and weaker upwelling along the western boundary. The northward pressure gradient force, strong vertical mixing, and horizontal mixing near the boundary are the essential dynamical components in the model's winter overturning circulation. The simulated water exchange is not hydraulically controlled in the Strait of Bab el Mandeb; instead, the exchange is limited by bottom and lateral boundary friction and, to a lesser extent, by interfacial friction due to the vertical viscosity at the interface between the inflow and the outflow.Partial support for this effort was provided by the Saudi Aramco Marine Environmental Research Center at KAUST.2014-10-1

Continuity of Local Time: An applied perspective

Continuity of local time for Brownian motion ranks among the most notable mathematical results in the theory of stochastic processes. This article addresses its implications from the point of view of applications. In particular an extension of previous results on an explicit role of continuity of (natural) local time is obtained for applications to recent classes of problems in physics, biology and finance involving discontinuities in a dispersion coefficient. The main theorem and its corollary provide physical principles that relate macro scale continuity of deterministic quantities to micro scale continuity of the (stochastic) local time.Comment: To appear in: "The fascination of Probability, Statistics and Their Applications. In honour of Ole E. Barndorff-Nielsen on his 80th birthday

Phytoplankton phenology indices in coral reef ecosystems : application to ocean-color observations in the Red Sea

© The Author(s), 2015. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Remote Sensing of Environment 160 (2015): 222-234, doi:10.1016/j.rse.2015.01.019.Phytoplankton, at the base of the marine food web, represent a fundamental food source in coral reef ecosystems. The timing (phenology) and magnitude of the phytoplankton biomass are major determinants of trophic interactions. The Red Sea is one of the warmest and most saline basins in the world, characterized by an arid tropical climate regulated by the monsoon. These extreme conditions are particularly challenging for marine life. Phytoplankton phenological indices provide objective and quantitative metrics to characterize phytoplankton seasonality. The indices i.e. timings of initiation, peak, termination and duration are estimated here using 15 years (1997–2012) of remote sensing ocean-color data from the European Space Agency (ESA) Climate Change Initiative project (OC-CCI) in the entire Red Sea basin. The OC-CCI product, comprising merged and bias-corrected observations from three independent ocean-color sensors (SeaWiFS, MODIS and MERIS), and processed using the POLYMER algorithm (MERIS period), shows a significant increase in chlorophyll data coverage, especially in the southern Red Sea during the months of summer NW monsoon. In open and reef-bound coastal waters, the performance of OC-CCI chlorophyll data is shown to be comparable with the performance of other standard chlorophyll products for the global oceans. These features have permitted us to investigate phytoplankton phenology in the entire Red Sea basin, and during both winter SE monsoon and summer NW monsoon periods. The phenological indices are estimated in the four open water provinces of the basin, and further examined at six coral reef complexes of particular socio-economic importance in the Red Sea, including Siyal Islands, Sharm El Sheikh, Al Wajh bank, Thuwal reefs, Al Lith reefs and Farasan Islands. Most of the open and deeper waters of the basin show an apparent higher chlorophyll concentration and longer duration of phytoplankton growth during the winter period (relative to the summer phytoplankton growth period). In contrast, most of the reef-bound coastal waters display equal or higher peak chlorophyll concentrations and equal or longer duration of phytoplankton growth during the summer period (relative to the winter phytoplankton growth period). The ecological and biological significance of the phytoplankton seasonal characteristics are discussed in context of ecosystem state assessment, and particularly to support further understanding of the structure and functioning of coral reef ecosystems in the Red Sea

A reduced-order strategy for 4D-Var data assimilation

This paper presents a reduced-order approach for four-dimensional variational data assimilation, based on a prior EO F analysis of a model trajectory. This method implies two main advantages: a natural model-based definition of a mul tivariate background error covariance matrix $\textbf{B}_r$, and an important decrease of the computational burden o f the method, due to the drastic reduction of the dimension of the control space. % An illustration of the feasibility and the effectiveness of this method is given in the academic framework of twin experiments for a model of the equatorial Pacific ocean. It is shown that the multivariate aspect of $\textbf{B}_r$ brings additional information which substantially improves the identification procedure. Moreover the computational cost can be decreased by one order of magnitude with regard to the full-space 4D-Var method

The Gulf of Aden Intermediate Water Intrusion Regulates the Southern Red Sea Summer Phytoplankton Blooms

Knowledge on large-scale biological processes in the southern Red Sea is relatively limited, primarily due to the scarce in situ, and satellite-derived chlorophyll-a (Chl-a) datasets. During summer, adverse atmospheric conditions in the southern Red Sea (haze and clouds) have long severely limited the retrieval of satellite ocean colour observations. Recently, a new merged ocean colour product developed by the European Space Agency (ESA)Ðthe Ocean Color Climate Change Initiative (OC-CCI)Ðhas substantially improved the southern Red Sea coverage of Chl-a, allowing the discovery of unexpected intense summer blooms. Here we provide the first detailed description of their spatiotemporal distribution and report the mechanisms regulating them. During summer, the monsoon-driven wind reversal modifies the circulation dynamics at the Bab-el-Mandeb strait, leading to a subsurface influx of colder, fresher, nutrient-rich water from the Indian Ocean. Using satellite observations, model simulation outputs, and in situ datasets, we track the pathway of this intrusion into the extensive shallow areas and coral reef complexes along the basin's shores. We also provide statistical evidence that the subsurface intrusion plays a key role in the development of the southern Red Sea phytoplankton blooms

A global method for coupling transport with chemistry in heterogeneous porous media

Modeling reactive transport in porous media, using a local chemical equilibrium assumption, leads to a system of advection-diffusion PDE's coupled with algebraic equations. When solving this coupled system, the algebraic equations have to be solved at each grid point for each chemical species and at each time step. This leads to a coupled non-linear system. In this paper a global solution approach that enables to keep the software codes for transport and chemistry distinct is proposed. The method applies the Newton-Krylov framework to the formulation for reactive transport used in operator splitting. The method is formulated in terms of total mobile and total fixed concentrations and uses the chemical solver as a black box, as it only requires that on be able to solve chemical equilibrium problems (and compute derivatives), without having to know the solution method. An additional advantage of the Newton-Krylov method is that the Jacobian is only needed as an operator in a Jacobian matrix times vector product. The proposed method is tested on the MoMaS reactive transport benchmark.Comment: Computational Geosciences (2009) http://www.springerlink.com/content/933p55085742m203/?p=db14bb8c399b49979ba8389a3cae1b0f&pi=1

Nonparametric Data Assimilation Scheme for Land Hydrological Applications

Data assimilation, which relies on explicit knowledge of dynamical models, is a well-known approach that addresses models' limitations due to various reasons, such as errors in input and forcing data sets. This approach, however, requires intensive computational efforts, especially for high-dimensional systems such as distributed hydrological models. Alternatively, data-driven methods offer comparable solutions when the physics underlying the models are unknown. For the first time in a hydrological context, a nonparametric framework is implemented here to improve model estimates using available observations. This method uses Takens delay coordinate method to reconstruct the dynamics of the system within a Kalman filtering framework, called the Kalman-Takens filter. A synthetic experiment is undertaken to fully investigate the capability of the proposed method by comparing its performance with that of a standard assimilation framework based on an adaptive unscented Kalman filter (AUKF). Furthermore, using terrestrial water storage (TWS) estimates obtained from the Gravity Recovery And Climate Experiment mission, both filters are applied to a real case scenario to update different water storages over Australia. In situ groundwater and soil moisture measurements within Australia are used to further evaluate the results. The Kalman-Takens filter successfully improves the estimated water storages at levels comparable to the AUKF results, with an average root-mean-square error reduction of 37.30% for groundwater and 12.11% for soil moisture estimates. Additionally, the Kalman-Takens filter, while reducing estimation complexities, requires a fraction of the computational time, that is, ~8 times faster compared to the AUKF approach

A conceptual approach to partitioning a vertical profile of phytoplankton biomass into contributions from two communities

This is the final version. Available from American Geophysical Union / Wiley via the DOI in this record. These data were collected and made freely available by the International Argo Program and the national programs that contribute to it (https://argo.ucsd.edu, https://www.ocean-ops.org). The Argo Program is part of the Global Ocean Observing System. All data and code used in the paper are provided openly on a GitHub page (https://github.com/rjbrewin/Two-community-phyto-model). This includes an example Jupyter Notebook Python Script, processing this BGC-Argo float and tuning the models. Details of how to run it without having to install software are provided as Supplementary Material to this manuscript.We describe an approach to partition a vertical profile of chlorophyll-a concentration into contributions from two communities of phytoplankton: one (community 1) that resides principally in the turbulent mixed-layer of the upper ocean and is observable through satellite visible radiometry; the other (community 2) residing below the mixed-layer, in a stably stratified environment, hidden from the eyes of the satellite. The approach is tuned to a time-series of profiles from a Biogeochemical-Argo float in the northern Red Sea, selected as its location transitions from a deep mixed layer in winter (characteristic of vertically well-mixed systems) to a shallow mixed layer in the summer with a deep chlorophyll-a maximum (characteristic of vertically stratified systems). The approach is extended to reproduce profiles of particle backscattering, by deriving the chlorophyll-specific backscattering coefficients of the two communities and a background coefficient assumed to be dominated by non-algal particles in the region. Analysis of the float data reveals contrasting phenology of the two communities, with community 1 blooming in winter and 2 in summer, community 1 negatively correlated with epipelagic stratification, and 2 positively correlated. We observe a dynamic chlorophyll-specific backscattering coefficient for community 1 (stable for community 2), positively correlated with light in the mixed-layer, suggesting seasonal changes in photoacclimation and/or taxonomic composition within community 1. The approach has the potential for monitoring vertical changes in epipelagic biogeography and for combining satellite and ocean robotic data to yield a three-dimensional view of phytoplankton distribution.Medical Research Council (MRC)Simons Foundation (SF)King Abdullah University of Science and Technology (KAUST)European Space Agency (ESA)European Space Agency (ESA

A Cell-Centred CVD-MPFA Finite Volume Method for Two-Phase Fluid Flow Problems with Capillary Heterogeneity and Discontinuity

A novel finite-volume method is presented for porous media flow simulation that is applicable to discontinuous capillary pressure fields. The method crucially retains the optimal single of freedom per control-volume being developed within the flux-continuous control-volume distributed multi-point flux approximation (CVD-MPFA) framework (Edwards and Rogers in Comput Geosci 02(04):259–290, 1998; Friis et al. in SIAM J Sci Comput 31(02):1192–1220, 2008) . The new methods enable critical subsurface flow processes involving oil and gas trapping to be correctly resolved on structured and unstructured grids. The results demonstrate the ability of the method to resolve flow with oil/gas trapping in the presence of a discontinuous capillary pressure field for diagonal and full-tensor permeability fields. In addition to an upwind approximation for the saturation equation flux, the importance of upwinding on capillary pressure flux via a novel hybrid formulation is demonstrated

Testing variational estimation of process parameters and initial conditions of an earth system model

We present a variational assimilation system around a coarse resolution Earth System Model (ESM) and apply it for estimating initial conditions and parameters of the model. The system is based on derivative information that is efficiently provided by the ESM's adjoint, which has been generated through automatic differentiation of the model's source code. In our variational approach, the length of the feasible assimilation window is limited by the size of the domain in control space over which the approximation by the derivative is valid. This validity domain is reduced by non-smooth process representations. We show that in this respect the ocean component is less critical than the atmospheric component. We demonstrate how the feasible assimilation window can be extended to several weeks by modifying the implementation of specific process representations and by switching off processes such as precipitation