Immiscible two-phase Darcy flow model accounting for vanishing and discontinuous capillary pressures: application to the flow in fractured porous media

Fully implicit time-space discretizations applied to the two-phase Darcy flow problem leads to the systems of nonlinear equations, which are traditionally solved by some variant of Newton’s method. The efficiency of the resulting algorithms heavily depends on the choice of the primary unknowns since Newton’s method is not invariant with respect to a nonlinear change of variable. In this regard, the role of capillary pressure/saturation relation is paramount because the choice of primary unknowns is restricted by its shape. We propose an elegant mathematical framework for two-phase flow in heterogeneous porous media resulting in a family of formulations, which apply to general monotone capillary pressure/saturation relations and handle the saturation jumps at rocktype interfaces. The presented approach is applied to the hybrid dimensional model of two-phase water-gas Darcy flow in fractured porous media for which the fractures are modelled as interfaces of co-dimension one. The problem is discretized using an extension of vertex approximate gradient scheme. As for the phase pressure formulation, the discrete model requires only two unknowns by degree of freedom.     

Sedimentological and geochemical records of past trophic state and hypolimnetic anoxia in large,hard-water Lake Bourget,French Alps

Sedimentological, geochemical and particle-size analyses were used to reconstruct the evolution of both trophic state and hypolimnetic anoxia in Lake Bourget (French Alps) during the last century. Radionuclide dating (²¹⁰Pb, ¹³⁷Cs and ²⁴¹Am) confirmed the annual rhythm of laminations in the upper sediment profile. In Lake Bourget, biochemical varves are triplets composed of a diatom layer (spring lamina), a bio-precipitated calcite-rich layer (spring/summer lamina), and a layer rich in organic matter and detrital particles (winter lamina). The onset of eutrophication and the first appearance of an anoxic facies occurred simultaneously and were dated by laminae counting to AD 1943±1 year. Persistent anoxic conditions began in AD 1960. Eutrophication is characterised by drastic increases in the flux of biogenic silica (mostly diatoms), lacustrine organic matter, and larger calcite crystals (15–30 μm). The increase of organic matter also represents a marker of the onset of anoxic conditions in the hypolimnion. Our results show that eutrophication was the main factor controlling anoxia in the hypolimnion. This eutrophication was caused mostly by the inflow of untreated sewage effluents, and to a lesser extent, by input of fertilizer-derived phosphorus during floods of the Rhone River and run-off from the lake catchment. The Rhone River, however, can also be a source of re-oxygenation via underflows that originate during flood events. Oxygenation of the hypolimnion is also controlled by low winter temperatures, which enable turnover of the lake. Thus, global warming, associated with a forecasted reduction in precipitation, might reduce the efficiency of hypolimnetic re-oxygenation in Lake Bourget.     

An attempt to separate anthropic and natural fire signals in an archaeological context-The case of the Mousterian site Grotta Reali (Rocchetta a Volturno Molise, Central Italy)

This paper presents an image analysis method for automated quantification of charcoal total area, focusing on the charcoal fractions less than 160 μm observed on 14 pollen slides from Grotta Reali samples. Four fire signals were recorded with 4 high values in the microcharcoal concentration curve. On the basis of modern microcharcoal study, mean length/width ratio of the microcharcoal particles was applied in an archeological context as an indicator of vegetal type (wood, grass or leaves). Therefore, the 4 fire signals were separated into two types: high concentration with high length/width ratio, and high concentration with low length/width ratio. Two fire signals might be interpreted as anthropic fire based on hearth or combustion areas, as inferred from archeological remains.     

Does Manila clam cultivation affect habitats of the engineer species Lanice conchilega (Pallas, 1766)?

The major French site of Manila clam Ruditapes philippinarum (Adams and Reeves, 1850) cultivation is located in the Chausey Archipelago where the associated practices are highly mechanized: every steps of production are made with tractor-driven machinery. The Manila clam concessions are concentrated on Lanice conchilega (Pallas, 1766) bioherms, which are known to increase alpha-diversity and to locally modify sediment dynamics. This study focus on the impacts of Manila clam cultivation on (i) the natural populations of L. conchilega and on (ii) the structure of the associated benthic assemblages during the different steps of the farming production cycle. We found that the L. conchilega populations are significantly affected within the concessions where their total abundances drastically decrease, their spatial patterns are modified and the associated benthic assemblages are significantly altered. Our results are discussed in a context of a sustainable management of the Manila clam cultivation in coastal areas.     

Regolith breccia Northwest Africa 7533: Mineralogy and petrology with implications for early Mars

Northwest Africa 7533, a polymict Martian breccia, consists of fine‐grained clast‐laden melt particles and microcrystalline matrix. While both melt and matrix contain medium‐grained noritic‐monzonitic material and crystal clasts, the matrix also contains lithic clasts with zoned pigeonite and augite plus two feldspars, microbasaltic clasts, vitrophyric and microcrystalline spherules, and shards. The clast‐laden melt rocks contain clump‐like aggregates of orthopyroxene surrounded by aureoles of plagioclase. Some shards of vesicular melt rocks resemble the pyroxene‐plagioclase clump‐aureole structures. Submicron size matrix grains show some triple junctions, but most are irregular with high intergranular porosity. The noritic‐monzonitic rocks contain exsolved pyroxenes and perthitic intergrowths, and cooled more slowly than rocks with zoned‐pyroxene or fine grain size. Noritic material contains orthopyroxene or inverted pigeonite, augite, calcic to intermediate plagioclase, and chromite to Cr‐bearing magnetite; monzonitic clasts contain augite, sodic plagioclase, K feldspar, Ti‐bearing magnetite, ilmenite, chlorapatite, and zircon. These feldspathic rocks show similarities to some rocks at Gale Crater like Black Trout, Mara, and Jake M. The most magnesian orthopyroxene clasts are close to ALH 84001 orthopyroxene in composition. All these materials are enriched in siderophile elements, indicating impact melting and incorporation of a projectile component, except for Ni‐poor pyroxene clasts which are from pristine rocks. Clast‐laden melt rocks, spherules, shards, and siderophile element contents indicate formation of NWA 7533 as a regolith breccia. The zircons, mainly derived from monzonitic (melt) rocks, crystallized at 4.43 ± 0.03 Ga (Humayun et al. 2013 ) and a ¹⁴⁷Sm‐¹⁴³Nd isochron for NWA 7034 yielding 4.42 ± 0.07 Ga (Nyquist et al. 2016 ) defines the crystallization age of all its igneous portions. The zircon from the monzonitic rocks has a higher Δ¹⁷O than other Martian meteorites explained in part by assimilation of regolith materials enriched during surface alteration (Nemchin et al. 2014 ). This record of protolith interaction with atmosphere‐hydrosphere during regolith formation before melting demonstrates a thin atmosphere, a wet early surface environment on Mars, and an evolved crust likely to have contaminated younger extrusive rocks. The latest events recorded when the breccia was on Mars are resetting of apatite, much feldspar and some zircons at 1.35–1.4 Ga (Bellucci et al. 2015 ), and formation of Ni‐bearing pyrite veins during or shortly after this disturbance (Lorand et al. 2015 ).     

Flow separation and vertical motions in a tidal flow interacting with a shallow-water island

This paper reports on the case study of Rattray Island (Great Barrier Reef, northeast Australia), lying perpendicular to tidal flow in shallow waters. At ebb and flood, attached (stable) eddies develop in the wake where swirls of turbidity suggest that sediment-laden waters are brought to the surface as a result of vertical transport. Both eddy and tip upwellings are encountered in the tidal flow around Rattray Island but there is currently no clear-cut answer as to which secondary flow generates upwelling with the largest intensity. This paper addresses this specific issue through idealized and realistic high-resolution numerical experiments. The analysis is supported by physical arguments based on the theory of flow separation. Given Rattray's geometry and surrounding bathymetry, the mechanism of flow separation in shallow waters helps explain the asymmetry in size of the eddies and their intensity. The results of idealized numerical experiments also suggest that eddy and tip upwellings may be of similar intensity at Rattray Island.     

The Beijing Climate Center atmospheric general circulation model: description and its performance for the present-day climate

The Beijing Climate Center atmospheric general circulation model version 2.0.1 (BCC_AGCM2.0.1) is described and its performance in simulating the present-day climate is assessed. BCC_AGCM2.0.1 originates from the community atmospheric model version 3 (CAM3) developed by the National Center for Atmospheric Research (NCAR). The dynamics in BCC_AGCM2.0.1 is, however, substantially different from the Eulerian spectral formulation of the dynamical equations in CAM3, and several new physical parameterizations have replaced the corresponding original ones. The major modification of the model physics in BCC_AGCM2.0.1 includes a new convection scheme, a dry adiabatic adjustment scheme in which potential temperature is conserved, a modified scheme to calculate the sensible heat and moisture fluxes over the open ocean which takes into account the effect of ocean waves on the latent and sensible heat fluxes, and an empirical equation to compute the snow cover fraction. Specially, the new convection scheme in BCC_AGCM2.0.1, which is generated from the Zhang and McFarlane’s scheme but modified, is tested to have significant improvement in tropical maximum but also the subtropical minimum precipitation, and the modified scheme for turbulent fluxes are validated using EPIC2001 in situ observations and show a large improvement than its original scheme in CAM3. BCC_AGCM2.0.1 is forced by observed monthly varying sea surface temperatures and sea ice concentrations during 1949–2000. The model climatology is compiled for the period 1971–2000 and compared with the ERA-40 reanalysis products. The model performance is evaluated in terms of energy budgets, precipitation, sea level pressure, air temperature, geopotential height, and atmospheric circulation, as well as their seasonal variations. Results show that BCC_AGCM2.0.1 reproduces fairly well the present-day climate. The combined effect of the new dynamical core and the updated physical parameterizations in BCC_AGCM2.0.1 leads to an overall improvement, compared to the original CAM3.     

Mechanisms of tropical Pacific interannual-to-decadal variability in the ARPEGE/ORCA global coupled model

Spatial and temporal structures of interannual-to-decadal variability in the tropical Pacific Ocean are investigated using results from a global atmosphere–ocean coupled general circulation model. The model produces quite realistic mean state characteristics, despite a sea surface temperature cold bias and a thermocline that is shallower than observations in the western Pacific. The periodicity and spatial patterns of the modelled El Niño Southern Oscillations (ENSO) compare well with those observed over the last 100 years, although the quasi-biennial timescale is dominant. Lag-regression analysis between the mean zonal wind stress and the 20°C isotherm depth suggests that the recently proposed recharge-oscillator paradigm is operating in the model. Decadal thermocline variability is characterized by enhanced variance over the western tropical South Pacific (~7°S). The associated subsurface temperature variability is primarily due to adiabatic displacements of the thermocline as a whole, arising from Ekman pumping anomalies located in the central Pacific, south of the equator. Related wind anomalies appear to be caused by SST anomalies in the eastern equatorial Pacific. This quasi-decadal variability has a timescale between 8 years and 20 years. The relationship between this decadal tropical mode and the low-frequency modulation of ENSO variance is also discussed. Results question the commonly accepted hypothesis that the low-frequency modulation of ENSO is due to decadal changes of the mean state characteristics.     

Arctic climate change in 21st century CMIP5 simulations with EC-Earth

The Arctic climate change is analyzed in an ensemble of future projection simulations performed with the global coupled climate model EC-Earth2.3. EC-Earth simulates the twentieth century Arctic climate relatively well but the Arctic is about 2 K too cold and the sea ice thickness and extent are overestimated. In the twenty-first century, the results show a continuation and strengthening of the Arctic trends observed over the recent decades, which leads to a dramatically changed Arctic climate, especially in the high emission scenario RCP8.5. The annually averaged Arctic mean near-surface temperature increases by 12 K in RCP8.5, with largest warming in the Barents Sea region. The warming is most pronounced in winter and autumn and in the lower atmosphere. The Arctic winter temperature inversion is reduced in all scenarios and disappears in RCP8.5. The Arctic becomes ice free in September in all RCP8.5 simulations after a rapid reduction event without recovery around year 2060. Taking into account the overestimation of ice in the twentieth century, our model results indicate a likely ice-free Arctic in September around 2040. Sea ice reductions are most pronounced in the Barents Sea in all RCPs, which lead to the most dramatic changes in this region. Here, surface heat fluxes are strongly enhanced and the cloudiness is substantially decreased. The meridional heat flux into the Arctic is reduced in the atmosphere but increases in the ocean. This oceanic increase is dominated by an enhanced heat flux into the Barents Sea, which strongly contributes to the large sea ice reduction and surface-air warming in this region. Increased precipitation and river runoff lead to more freshwater input into the Arctic Ocean. However, most of the additional freshwater is stored in the Arctic Ocean while the total Arctic freshwater export only slightly increases.