Phytoplankton versus macrophyte contribution to primary production and biogeochemical cycles of a coastal mesotidal system. A modelling approach

This study presents an assessment of the contributions of various primary producers to the global annual production and N/P cycles of a coastal system, namely the Arcachon Bay, by means of a numerical model. This 3D model fully couples hydrodynamic with ecological processes and simulates nitrogen, silicon and phosphorus cycles as well as phytoplankton, macroalgae and seagrasses. Total annual production rates for the different components were calculated for different years (2005, 2007 and 2009) during a time period of drastic reduction in seagrass beds since 2005. The total demand of nitrogen and phosphorus was also calculated and discussed with regards to the riverine inputs. Moreover, this study presents the first estimation of particulate organic carbon export to the adjacent open ocean. The calculated annual net production for the Arcachon Bay (except microphytobenthos, not included in the model) ranges between 22,850 and 35,300 tons of carbon. The main producers are seagrasses in all the years considered with a contribution ranging from 56% to 81% of global production. According to our model, the -30% reduction in seagrass bed surface between 2005 and 2007, led to an approximate 55% reduction in seagrass production, while during the same period of time, macroalgae and phytoplankton enhanced their productions by about +83% and +46% respectively. Nonetheless, the phytoplankton production remains about eightfold higher than the macroalgae production. Our results also highlight the importance of remineralisation inside the Bay, since riverine inputs only fulfill at maximum 73% nitrogen and 13% phosphorus demands during the years 2005, 2007 and 2009. Calculated advection allowed a rough estimate of the organic matter export: about 10% of the total production in the bay was exported, originating mainly from the seagrass compartment, since most of the labile organic matter was remineralised inside the bay.

A New Formulation of the Spectral Energy Budget of the Atmosphere, with Application to Two High-Resolution General Circulation Models

International audience

A bistable reaction-diffusion system in a stretching flow

We examine the evolution of a bistable reaction in a one-dimensional stretching flow, as a model for chaotic advection. We derive two reduced systems of ordinary differential equations (ODEs) for the dynamics of the governing advection-reaction-diffusion partial differential equations (PDE), for pulse-like and for plateau-like solutions, based on a non-perturbative approach. This reduction allows us to study the dynamics in two cases: first, close to a saddle-node bifurcation at which a pair of nontrivial steady states are born as the dimensionless reaction rate (Damkoehler number) is increased, and, second, for large Damkoehler number, far away from the bifurcation. The main aim is to investigate the initial-value problem and to determine when an initial condition subject to chaotic stirring will decay to zero and when it will give rise to a nonzero final state. Comparisons with full PDE simulations show that the reduced pulse model accurately predicts the threshold amplitude for a pulse initial condition to give rise to a nontrivial final steady state, and that the reduced plateau model gives an accurate picture of the dynamics of the system at large Damkoehler number. Published in Physica D (2006)

Dissolved inorganic carbon budgets in the eastern subpolar North Atlantic in the 2000s from in situ data

The subpolar North Atlantic (SPNA) is important in the global carbon cycle because of the deep water ventilation processes that lead to both high uptake of atmospheric CO2 and large inventories of anthropogenic CO2 (C-ant). Thus, it is crucial to understand its response to increasing anthropogenic pressures. In this work, the budgets of dissolved inorganic carbon (DIC), C-ant and natural DIC (DICnat) in the eastern SPNA in the 2000s, are jointly analyzed using in situ data. The DICnat budget is found to be in steady state, confirming a long-standing hypothesis from in situ data for the first time. The biological activity is driving the uptake of natural CO2 from the atmosphere. The C-ant increase in the ocean is solely responsible of the DIC storage rate which is explained by advection of C-ant from the subtropics (65%) and C-ant air-sea flux (35%). These results demonstrate that the C-ant is accumulating in the SPNA without affecting the natural carbon cycle.

Simulação numérica da convecção mista em cavidade preenchida com meio poroso heterogêneo e homogêneo

In this work is presented mixed convection heat transfer inside a lid-driven cavity heated from below and filled with heterogeneous and homogeneous porous medium. In the heterogeneous approach, the solid domain is represented by heat conductive equally spaced blocks; the fluid phase surrounds the blocks being limited by the cavity walls. The homogeneous or pore-continuum approach is characterized by the cavity porosity and permeability. Generalized mass, momentum and energy conservation equations are obtained in dimensionless form to represent both the continuum and the pore-continuum models. The numerical solution is obtained via the finite volume method. QUICK interpolation scheme is set for numerical treatment of the advection terms and SIMPLE algorithm is applied for pressure-velocity coupling. Aiming the laminar regime, the flow parameters are kept in the range of 102≤Re≤103 and 103≤Ra≤106 for both the heterogeneous and homogeneous approaches. In the tested configurations for the continuous model, 9, 16, 36, and 64 blocks are considered for each combination of Re and Ra being the microscopic porosity set as constant φ=0,64 . For the pore-continuum model the Darcy number (Da) is set according to the number of blocks in the heterogeneous cavity and the φ. Numerical results of the comparative study between the microscopic and macroscopic approaches are presented. As a result, average Nusselt number equations for the continuum and the pore continuum models as a function of Ra and Re are obtained.

Intraseasonal variability of mixed layer depth in the tropical Indian Ocean

In this paper, we use an observational dataset built from Argo in situ profiles to describe the main large-scale patterns of intraseasonal mixed layer depth (MLD) variations in the Indian Ocean. An eddy permitting (0.25A degrees) regional ocean model that generally agrees well with those observed estimates is then used to investigate the mechanisms that drive MLD intraseasonal variations and to assess their potential impact on the related SST response. During summer, intraseasonal MLD variations in the Bay of Bengal and eastern equatorial Indian Ocean primarily respond to active/break convective phases of the summer monsoon. In the southern Arabian Sea, summer MLD variations are largely driven by seemingly-independent intraseasonal fluctuations of the Findlater jet intensity. During winter, the Madden-Julian Oscillation drives most of the intraseasonal MLD variability in the eastern equatorial Indian Ocean. Large winter MLD signals in northern Arabian Sea can, on the other hand, be related to advection of continental temperature anomalies from the northern end of the basin. In all the aforementioned regions, peak-to-peak MLD variations usually reach 10 m, but can exceed 20 m for the largest events. Buoyancy flux and wind stirring contribute to intraseasonal MLD fluctuations in roughly equal proportions, except for the Northern Arabian Sea in winter, where buoyancy fluxes dominate. A simple slab ocean analysis finally suggests that the impact of these MLD fluctuations on intraseasonal sea surface temperature variability is probably rather weak, because of the compensating effects of thermal capacity and sunlight penetration: a thin mixed-layer is more efficiently warmed at the surface by heat fluxes but loses more solar flux through its lower base.

Chaotic mixing of a competitive-consecutive reaction

The evolution of a competitive-consecutive chemical reaction is computed numerically in a two-dimensional chaotic fluid flow with initially segregated reactants. Results from numerical simulations are used to evaluate a variety of reduced models commonly adopted to model the full advection-reaction-diffusion problem. Particular emphasis is placed upon fast reactions, where the yield varies most significantly with Peclet number (the ratio of diffusive to advective time scales). When effects of the fluid mechanical mixing are strongest, we find that the yield of the reaction is underestimated by a one-dimensional lamellar model that ignores the effects of fluid mixing, but overestimated by two other lamellar models that include fluid mixing.

A four-dimensional analysis of the thermal structure in the Gulf of Lion

The Theoretical and Experimental Tomography in the Sea Experiment (THETIS 1) took place in the Gulf of Lion to observe the evolution of the temperature field and the process of deep convection during the 1991-1992 winter. The temperature measurements consist, of moored sensors, conductivity-temperature-depth and expendable bathythermograph surveys, ana acoustic tomography. Because of this diverse data set and since the field evolves rather fast, the analysis uses a unified framework, based on estimation theory and implementing a Kalman filter. The resolution and the errors associated with the model are systematically estimated. Temperature is a good tracer of water masses. The time-evolving three-dimensional view of the field resulting from the analysis shows the details of the three classical convection phases: preconditioning, vigourous convection, and relaxation. In all phases, there is strong spatial nonuniformity, with mesoscale activity, short timescales, and sporadic evidence of advective events (surface capping, intrusions of Levantine Intermediate Water (LIW)). Deep convection, reaching 1500 m, was observed in late February; by late April the field had not yet returned to its initial conditions (strong deficit of LIW). Comparison with available atmospheric flux data shows that advection acts to delay the occurence of convection and confirms the essential role of buoyancy fluxes. For this winter, the deep. mixing results in an injection of anomalously warm water (Delta T similar or equal to 0.03 degrees) to a depth of 1500 m, compatible with the deep warming previously reported.

Coastal zone oil spill modeling

This work presents the development and application of a three-dimensional oil spill model for predicting the movement of an oil slick in the coastal waters of Singapore. In the model, the oil slick is divided into a number of small elements for simulating of the oil processes of spreading, advection, turbulent diffusion. This model is capable of predicting the horizontal movement of surface oil slick. Satellite images and field observations of oil slicks on the surface in the Singapore Straits are used to validate the newly developed model. Compared with the observations, the numerical results of the oil spill model show good conformity. In this study, the 3d model was generated using the geometrical data of Singapore Straits waters by GAMBIT which is a pre-processor of FLUENT programmed.

New insights into the transport processes controlling the sulfate-methane-transition-zone near methane vents

Over the past years, several studies have raised concerns about the possible interactions between methane hydrate decomposition and external change. To carry out such an investigation, it is essential to characterize the baseline dynamics of gas hydrate systems related to natural geological and sedimentary processes. This is usually treated through the analysis of sulfate-reduction coupled to anaerobic oxidation of methane (AOM). Here, we model sulfate reduction coupled with AOM as a two-dimensional (2D) problem including, advective and diffusive transport. This is applied to a case study from a deep-water site off Nigeria’s coast where lateral methane advection through turbidite layers was suspected. We show by analyzing the acquired data in combination with computational modeling that a two-dimensional approach is able to accurately describe the recent past dynamics of such a complex natural system. Our results show that the sulfate-methane-transition-zone (SMTZ) is not a vertical barrier for dissolved sulfate and methane. We also show that such a modeling is able to assess short timescale variations in the order of decades to centuries.

THE IMPACT OF UPPER-LEVEL PROCESSES ON THE INTENSITY AND STRUCTURAL CHANGES OF HURRICANE SANDY (2012)

The first part of this study examines the relative roles of frontogenesis and tropopause undulation in determining the intensity and structural changes of Hurricane Sandy (2012) using a high-resolution cloud-resolving model. A 138-h simulation reproduces Sandy’s four distinct development stages: (i) rapid intensification, (ii) weakening, (iii) steady maximum surface wind but with large continued sea-level pressure (SLP) falls, and (iv) re-intensification. Results show typical correlations between intensity changes, sea-surface temperature and vertical wind shear during the first two stages. The large SLP falls during the last two stages are mostly caused by Sandy’s moving northward into lower-tropopause regions associated with an eastward-propagating midlatitude trough, where the associated lower-stratospheric warm air wraps into the storm and its surrounding areas. The steady maximum surface wind occurs because of the widespread SLP falls with weak pressure gradients lacking significant inward advection of absolute angular momentum (AAM). Meanwhile, there is a continuous frontogenesis in the outer region during the last three stages. Cyclonic inward advection of AAM along each frontal rainband accounts for the continued expansion of the tropical-storm-force wind and structural changes, while deep convection in the eyewall and merging of the final two survived frontal rainbands generate a spiraling jet in Sandy’s northwestern quadrant, leading to its re-intensification prior to landfall. The physical, kinematic and dynamic aspects of an upper-level outflow layer and its possible impact on the re-intensification of Sandy are examined in the second part of this study. Above the outflow layer isentropes are tilted downward with radius as a result of the development of deep convection and an approaching upper-level trough, causing weak subsidence. Its maximum outward radial velocity is located above the cloud top, so the outflow channel experiences cloud-induced long-wave cooling. Because Sandy has two distinct convective regions (an eyewall and a frontal rainband), it has multiple outflow layers, with the eyewall’s outflow layer located above that of the frontal rainband. During the re-intensification stage, the eyewall’s outflow layer interacts with a jet stream ahead of the upper-level trough axis. Because of the presence of inertial instability on the anticyclonic side of the jet stream and symmetric instability in the inner region of the outflow layer, Sandy’s secondary circulation intensifies. Its re-intensification ceases when these instabilities disappear. The relationship between the intensity of the secondary circulation and dynamic instabilities of the outflow layer suggests that the re-intensification occurs in response to these instabilities. Additionally, it is verified that the long-wave cooling in the outflow layer helps induce symmetric instability by reducing static stability.

Reduced models of chemical reaction in chaotic flows

We describe and evaluate two reduced models for nonlinear chemical reactions in a chaotic laminar flow. Each model involves two separate steps to compute the chemical composition at a given location and time. The “manifold tracking model” first tracks backwards in time a segment of the stable manifold of the requisite point. This then provides a sample of the initial conditions appropriate for the second step, which requires solving one-dimensional problems for the reaction in Lagrangian coordinates. By contrast, the first step of the “branching trajectories model” simulates both the advection and diffusion of fluid particles that terminate at the appropriate point; the chemical reaction equations are then solved along each of the branched trajectories in a second step. Results from each model are compared with full numerical simulations of the reaction processes in a chaotic laminar flow.

Changes in the storage and sink of carbon dioxide in subsurface atmospheres controlled by climate-driven processes: the case of the Ojo Guareña karst system

A comprehensive environmental monitoring program was conducted in the Ojo Guareña cave system (Spain), one of the longest cave systems in Europe, to assess the magnitude of the spatiotemporal changes in carbon dioxide gas (CO2) in the cave–soil–atmosphere profile. The key climate-driven processes involved in gas exchange, primarily gas diffusion and cave ventilation due to advective forces, were characterized. The spatial distributions of both processes were described through measurements of CO2 and its carbon isotopic signal (δ13C[CO2]) from exterior, soil and cave air samples analyzed by cavity ring-down spectroscopy (CRDS). The trigger mechanisms of air advection (temperature or air density differences or barometric imbalances) were controlled by continuous logging systems. Radon monitoring was also used to characterize the changing airflow that results in a predictable seasonal or daily pattern of CO2 concentrations and its carbon isotopic signal. Large daily oscillations of CO2 levels, ranging from 680 to 1900 ppm day−1 on average, were registered during the daily oscillations of the exterior air temperature around the cave air temperature. These daily variations in CO2 concentration were unobservable once the outside air temperature was continuously below the cave temperature and a prevailing advective-renewal of cave air was established, such that the daily-averaged concentrations of CO2 reached minimum values close to atmospheric background. The daily pulses of CO2 and other tracer gases such as radon (222Rn) were smoothed in the inner cave locations, where fluctuation of both gases was primarily correlated with medium-term changes in air pressure. A pooled analysis of these data provided evidence that atmospheric air that is inhaled into dynamically ventilated caves can then return to the lower troposphere as CO2-rich cave air.

Nocturnal surface ozone enhancement over Portugal during winter: Influence of different atmospheric conditions

our distinct nocturnal surface ozone (NSO) enhancement events were observed, with NSO concentration exceeding 80μg/m3, at multiple ozone (O3) monitoring stations (32 sites) in January, November and December between year 2000–2010, in Portugal. The reasonable explanation for the observed bimodal pattern of surface ozone with enhanced NSO concentration during nighttime has to be transport processes, as the surface ozone production ceases at nighttime. Simultaneous measurements of O3 at multiple stations during the study period in Portugal suggest that horizontal advection alone cannot explain the observed NSO enhancement. Thus, detailed analysis of the atmospheric conditions, simulated with the Weather Research and Forecasting (WRF) model, were performed to evaluate the atmospheric mechanisms responsible for NSO enhancement in the region. Simulations revealed that each event occurred as a result of one or the combination of different atmospheric processes such as, passage of a cold front followed by a subsidence zone; passage of a moving surface trough, with associated strong horizontal wind speed and vertical shear; combination of vertical and horizontal transport at the synoptic scale; formation of a low level jet with associated vertical mixing below the jet stream. The study confirmed that large-scale flow pattern resulting in enhanced vertical mixing in the nocturnal boundary layer, plays a key role in the NSO enhancement events, which frequently occur over Portugal during winter months.

Extensional orogenic collapse captured by strike-slip tectonics:

The late Paleozoic collision between Gondwana and Laurussia resulted in the polyphase deformation and magmatism that characterizes the Iberian Massif of the Variscan orogen. In the Central Iberian Zone, initial con- tinental thickening (D1; folding and thrusting) was followed by extensional orogenic collapse (D2) responsible for the exhumation of high-grade rocks coeval to the emplacement of granitoids. This study presents a tectonometamorphic analysis of the Trancoso-Pinhel region (Central Iberian Zone) to ex- plain the processes in place during the transition froman extension-dominated state (D2) to a compression-dom- inated one (D3).Wereveal the existence of low-dipping D2 extensional structures later affected by several pulses of subhorizontal shortening, each of them typified by upright folds and strike-slip shearing (D3, D4 and D5, as identified by superimposition of structures). The D2 Pinhel extensional shear zone separates a low-grade domain from an underlying high-grade domain, and it contributed to the thermal reequilibration of the orogen by facil- itating heat advection from lower parts of the crust, crustal thinning, decompression melting, and magma intru- sion. Progressive lessening of the gravitational disequilibrium carried out by this D2 shear zone led to a switch from subhorizontal extension to compression and the eventual cessation and capture of the Pinhel shear zone by strike-slip tectonics during renewed crustal shortening. High-grade domains of the Pinhel shear zone were folded together with low-grade domains to define the current upright folded structure of the Trancoso-Pinhel re- gion, the D3 Tamames-Marofa-Sátão synform. Newdating of syn-orogenic granitoids (SHRIMP U\\Pb zircon dat- ing) intruding the Pinhel shear zone, together with the already published ages of early extensional fabrics constrain the functioning of this shear zone to ca. 331–311 Ma, with maximum tectonomagmatic activity at ca. 321–317 Ma. The capture and apparent cessation of movement of the Pinhel shear zone occurred at ca. 317– 311 Ma.