Dispersive Particle Transport: Identification of Macroscale Behavior in Heterogeneous Stratified Subsurface Flows

Dispersive mass transport processes in naturally heterogeneous geological formations (porous media) are investigated based on a particle approach to mass transport and on its numerical implementation using LPT3D, a Lagrangian Particle Tracking 3D code. We are currently using this approach for studying microscale and macroscale space–time behavior (advection, diffusion, dispersion) of tracer plumes, solutes, or miscible fluids, in 1,2,3-dimensional heterogeneous and anisotropic subsurface formations (aquifers, petroleum reservoirs). Our analyses are based on a general advection-diffusion model and numerical scheme where concentrations and fluxes are discretized in terms of particles. The advection-diffusion theory is presented in a probabilistic framework, and in particular, a numerical analysis is developed for the case of advective transport and rotational flows (numerical stability of the explicit Euler scheme). The remainder of the paper is devoted to the behavior of concentration, mass flux density, and statistical moments of the transported tracer plume in the case of heterogeneous steady flow fields, where macroscale dispersion occurs due to geologic heterogeneity and stratification. We focus on the case of perfectly stratified or multilayered media, obtained by generating many horizontal layers with a purely random transverse distribution of permeability and horizontal velocity. In this case, we calculate explicitly the exact mass concentration field C(x, t), mass flux density field f(x, t), and moments. This includes spatial moments and dispersion variance ² _x (t) on a finite domain L, and temporal moments on a finite time scale T, e.g., the mass variance of arrival times ² _T (x). The moments are related to flux concentrations in a way that takes explicitly into account finite space–time scales of analysis (time-dependent tracer mass; spatially variable flow through mass). The multilayered model problem is then used in numerical experiments for testing different ways of recovering information on tracer plume migration, dispersion, concentration and flux fields. Our analyses rely on a probabilistic interpretation that emerges naturally from the particle approach; it is based on spatial moments (particle positions), temporal moments (mass weighted arrival times), and probability densities (both concentrations and fluxes). Finally, as an alternative to direct estimations of the flux and concentration fields, we formulate and study the Moment Inverse Problem. Solving the MIP yields an indirect method for estimating the space–time distribution of flux concentrations based on observed or estimated moments of the plume. The moments may be estimated from field measurements, or numerically computed by particle tracking as we do here.     

神农架大九湖四万年以来的植被与气候变化

通过对大九湖6 m长连续沉积岩芯剖面(DJH-1 孔) 7 个样品的AMS¹⁴C年龄测定和151 块孢粉样品的鉴定分析, 揭示了神农架区域4 万以年来植被和气候演变。末次冰期阶段大九湖附近发育森林草地或草地-草甸植被。MIS 3 晚期, ~42-39 cal ka BP之间, 气候相对干冷, 发育森林草地;~39-31 cal ka BP 之间, 气候较为湿润, 草甸扩张并伴随低海拔阔叶树种的发育。MIS 2 阶段, 草甸组分中蒿属显著增加, 高海拔可能分布有荒漠草地, 气候极端干冷;该时期植被带垂直下降达到1000 余m, 按垂直温度递减率推算, 冰盛期阶段该区域温度下降约7℃左右。从冰消期开始, 森林植被开始扩张, 北温带、暖温带和亚热带乔木组分依次增加。约在9.4-4 cal ka BP之间, 演变为亚热带常绿阔叶落叶林, 属全新世适宜期;从约4 cal ka BP以来, 北温带阔叶和针叶树开始增加, 气候趋于凉干。通过对比区域高分辨率的洞穴石笋及高纬冰芯氧同位素记录, 表明神农架区域植被环境变化对气候变化敏感, 并记录了H1, YD气候突变事件;进一步体现出该区域气候环境演变主要与北半球太阳辐射控制的东亚夏季风强度变化有关, 且与北半球高纬气候变化一致。     

Towards a microzonation of the Greater Beirut area: an instrumental approach combining earthquake and ambient vibration recordings

Lebanon is situated on the 1000 km long Levant transform fault that separates the Arabic from the African tectonic plates. In Lebanon, the Levant fault splits up into a set of ramifications that had, in the past, generated major destructive earthquakes causing a lot of destruction and thousands of casualties. The most devastating one was the 551 A.D. offshore earthquake that destroyed Beirut, the capital of Lebanon. This paper presents a site effect study in Beirut, aimed at proposing a framework for future microzonation works in the city. It includes two complementary parts. A 6-month, temporary seismological experiment was first conducted to estimate the site response at 10 sites sampling the main geological units of Beirut on the basis of local and regional earthquake recordings. This spatially sparse information was then complemented by a large number (615) of microtremor measurements covering the Beirut municipality and part of its suburbs with a 400 m dense grid. The recordings were analysed with the standard site-to-reference and horizontal-to-vertical spectral ratio methods for earthquake recordings, and the horizontal-to-vertical ratio for ambient noise recordings. Significant ground motion amplification effects (up to a factor of 8) are found in a few areas corresponding to recent deposits. The consistency between results from earthquake and microtremor recordings allows proposing a map of the resonance frequencies within the city and its suburbs, with frequencies ranging from 0.5 to 5 Hz for the deepest deposits, and 5–10 Hz for shallow areas. Finally, the results are discussed and a way to combine the results obtained from the temporary stations to the great number of recordings coming from the permanent Lebanese network is proposed.     

Fish composition and assemblage structure in three Eastern English Channel macrotidal estuaries: A comparison with other French estuaries

This study has analysed for the first time fish composition and assemblage structures of three small macrotidal estuaries of the Eastern English Channel (EEC) and has explored the influences of 19 biotic and abiotic variables on the fish assemblages. Fish from Canche, Authie and Somme estuaries were collected during spring (June 2006 and May 2007) and autumn (September 2006) along the estuarine gradients using a 1.5 m beam trawl. Using identical sampling protocols, the study also analysed and compared for the first time taxonomic and functional aspects of the fish assemblages in 15 estuaries located along the Atlantic and English Channel coasts. SIMPER analysis showed high similarities in fish assemblages in the three EEC estuaries and during either spring or autumn periods. However, intra-estuary similarities were relatively low, indicating that fish assemblage structures (species richnesses or abundances) were more variable within the estuary (salinity gradient) than between estuaries and/or seasons (spring vs autumn). Although numerous environmental variables were included in the study, only 47% of the variability observed in the fish distribution was explained. Fish spatial variations in the EEC estuaries are mostly driven by abiotic variables as opposed to biological interactions. As indicated by CCA, salinity and muddy sediments were the two most important factors structuring the fish assemblages. The macrobenthos being very abundant in the EEC estuaries (580–1121 ind. m⁻²), the availability of potential prey is probably not a limiting factor in the utilization of estuaries by fish. Contrary to the majority of French estuaries dominated by estuarine species (ES), the fish assemblages of the EEC estuaries are clearly dominated by marine migrant (MM) species (65% on average) with high abundance of juveniles (mostly young-of-the-year). Cluster and SIMPROF's analyses distinguished the functional structure of the 15 estuarine fish assemblages into different clusters. The three EEC estuaries form a similar group with the largest French estuaries (Seine, Loire and Gironde). The latters were characterized by higher MM/ES ratio compared to the other French estuaries. Freshwater flow and tidal range were the main factors associated with this group of estuaries. The results of the present study do not indicate that anthropogenic impacts (chemical contamination) on estuaries affect their ecological functioning as described by the fish ecological guilds. It is suggested that the guild approach may not be useful to provide valuable information on the ecological status of estuaries.     

Assessment of soil losses and siltation of the K’sob hydrological system (semiarid area—East Algeria)

Soil losses and siltation of the hydrological system (watershed–dam) of K’sob were obtained using direct and indirect methods. The Wadi K’sob watershed of 1,484 km², average slope of 0.14, and average elevation of 1,060 m is located in a semiarid climate. The average annual rainfall is 341 mm and the mean annual water discharge is 0.89 m³/s. Data from the Medjez gauging station located 6 km upstream of the dam, are the daily liquid flow and instantaneous concentrations of suspended sediments. Over a time period from 1973 to 2010, the relationship between water and sediment discharges is quantified by the equation: Q _s?=?5.6 Q ^1.31. Thus, in view of the availability data on a daily scale, the assessment of soil erodibility of the K’sob watershed was used to estimate specific soil losses of 203 t?km^?2?year^?1or 301,000 t eroded annually from the K’sob basin. The bathymetric measurements of the sediment volumes deposited in the K’sob dam, has quantified the annual siltation of 0.8 hm³, corresponding to an average erodibility of the K’sob watershed of 809 t?km^?2?year^?1. However, when adding the volume of sediment removed by the dredging operation and de-silting by the valves during heavy floods, the value of soil losses is 2,780 t?km^?2?year^?1. The indirect assessment of soil erodibility of the basin was obtained by applying two models: the quantitative geomorphological analysis (QGA) and PISA model (prediction of silting in the artificial reservoirs, in Italian: Previsioni dell’Interimento nei Serbatoi Artificiali) using physical and climatic factors in the watershed. The obtained results by QGA method underestimate specific soil losses of 524 t?km^?2?year^?1. The PISA model gives a value of 2,915 t?km^?2?year^?1, which is close to the value obtained by bathymetric measurements. This study concludes that PISA model is most suitable to estimate soil loss and siltation of the K’sob hydrological system.     

Cooling of quark stars in the colour superconductive phase: effect of photons from glueball decay

The cooling history of a quark star in the colour superconductive phase is investigated. Here we specifically focus on the two-flavour colour (2SC) phase where the novel process of photon generation via glueball (GLB) decay has already been investigated. The picture we present here can, in principle, be generalized to quark stars entering a superconductive phase where similar photon generation mechanisms are at play. As much as 10⁴⁵–10⁴⁷ erg of energy is provided by the GLB decay in the 2SC phase. The generated photons slowly diffuse out of the quark star, keeping it hot and radiating as a blackbody (with possibly a Wien spectrum in gamma-rays) for millions of years. We discuss hot radio-quiet isolated neutron stars in our picture (such as RX J185635–3754 and RX J0720.4–3125) and argue that their nearly blackbody spectra (with a few broad features) and their remarkably tiny hydrogen atmosphere are indications that these might be quark stars in the colour superconductive phase where some sort of photon generation mechanism (reminiscent of the GLB decay) has taken place. Fits to observed data of cooling compact stars favour models with superconductive gaps of  Δ_2SC∼ 15–35 MeV  and densities  ρ_2SC= (2.5–3.0) ×ρ_N  (ρ_N being the nuclear matter saturation density) for quark matter in the 2SC phase. If correct, our model combined with more observations of isolated compact stars could provide vital information to studies of quark matter and its exotic phases.     

Numerical simulations of astrophysical jets from Keplerian discs – III. The effects of mass loading

We present 2.5D time-dependent simulations of the non-linear evolution of non-relativistic outflows from the surface of Keplerian accretion discs. The gas is accelerated from the surface of the disc (which is a fixed platform in these simulations) into a cold corona in stable hydrostatic equilibrium. We explore the dependence of the resulting jet characteristics upon the mass loading of the winds. Two initial configurations of the threading disc magnetic field are studied: a potential field and a uniform vertical field configuration.
We show that the nature of the resulting highly collimated, jet-like outflows (steady or episodic) is determined by the mass load of the disc wind. The mass load controls the interplay between the collimating effects of the toroidal field and the kinetic energy density in the outflow. In this regard, we demonstrate that the onset of episodic behaviour of jets appears to be determined by the quantity     which compares the speed for a toroidal Alfvén wave to cross the diameter of the jet, with the flow speed v _p along the jet. This quantity decreases with increasing load. For sufficiently large N (small mass loads), disturbances appear to grow leading to instabilities and shocks. Knots are then generated and the outflow becomes episodic. These effects are qualitatively independent of the initial magnetic configuration that we employed and are probably generic to a wide variety of magnetized accretion disc models.     

Eastern Mediterranean palaeoclimates from 26 to 5 ka B.P. documented by pollen and isotopic analysis of a core in the anoxic Bannock Basin

Continuous pollen and isotopic records were established for core BAN 84 09 GC retrieved from the anoxic Bannock Basin in the Eastern Mediterranean. On the basis of two ¹⁴C dates, they document the palaeoclimate between about 25.7 ka B.P. and 5.2 ka B.P. in the northern borderlands of the Ionian Basin. The upper half of the core has been redeposited.

The isotopic record displays a correlation with pollen percentages that is strong and positive for Artemisia (sage-brush) and negative for Quercus (oak). The last glacial maximum and the deglaciation are identified by these combined taxa, together with Chenopodiaceae. The glacial maximum around 18 ka B.P. (which has elsewhere been dated from 20 to 15 ka B.P.) has pollen percentages that are high for Artemisia and low for Quercus. The climate in the pollen source area was arid, cold in winter, briefly warm in summer and sustained the vegetation of a semi-desert. The onset of deglaciation after 18 ka B.P. coincides with that of the decline in Artemisia pollen percentage. However, this decline does not indicate reduced aridity, because it is accompanied by a pollen percentage rise of the even more arid herbs Chenopodiaceae and Ephedra. Throughout the deglaciation from 18 to 11 ka B.P., the aridity progressively increases, culminating at 11 ka B.P. This trend is briefly interrupted by a more humid event, shown by a peak in Artemisia pollen percentage and a smaller peak in oak; these two peaks are coeval with the Bölling-Alleröd chronozone (13-11 ka B.P.). Maximum aridity occurs during the Younger Dryas chronozone (11-10 ka B.P.). Afterwards, the oak pollen percentage begins a steady increase, and its maximum value is coeval with the lowest isotopic value, dated at 8760 ± 170 yr B.P. This period was one of high moisture, warm summers, and, according to altitude, mild to cool winters. This climate sustained forests that were Mediterranean in the lowlands and warm temperate in the uplands. A high pollen concentration is observed during this period and reveals the presence of sapropel S1, which is otherwise unrecognizable in this entirely black core. During the following period between 8760 ± 170 and 5200 yr B.P., the δ¹⁸0 reverts to slightly higher values and the Quercus pollen percentage decreases, while the pollen percentage of the wetter Ostrya, the oriental hornbeam, increases. The high pollen concentration during the deposition of sapropel S1 cannot have been caused by increased pollen input into the sea, this pollen being wind-borne, nor by increased pollen production for all taxa, both trees and herbs. We conclude that it is entirely due to increased preservation of this allochtonous organic material by the deep anoxia of the bottom water, below a thick anoxic water column. The coincidence of sapropel deposition with warm and humid local climate as well as with the second global meltwater pulse suggests that the cessation of bottom-water ventilation was due to decreased surface water density, resulting from less saline incoming Atlantic surface water, increased local runoff, and warmer winters.     

Cristallisation fractionnée et contamination crustale dans la série magmatique jurassique transitionnelle du Haut Atlas central (Maroc)Fractional crystallisation and crustal contamination in the transitional Jurassic magmatic series of Central High Atlas (Morocco)

The Middle Jurassic plutonism of the Central High Atlas (Morocco) was emplaced in N45° trending anticlinal ridges. It is characterised by various petrographic facies including mafic rocks (troctolites), intermediate rocks (diorites, monzodiorites), and evolved rocks (syenites), together with heterogeneous facies resulting from mixing between acidic and the intermediate magmas. Mineralogical and chemical data show (i) the transitional character of the Jurassic magmatic series of the Central High Atlas and (ii) the implication of continental crust as a contaminant during fractional crystallization. To cite this article: R. Zayane et al., C. R. Geoscience 334 (2002) 97–104.     

Environmental factors structuring fish composition and assemblages in a small macrotidal estuary (eastern English Channel)

The fish assemblage structure was analyzed along an estuarine gradient of a small macrotidal estuary (the Canche, France). Fishes were collected every two months between May 2006 and July 2007 from 12 sampling stations using a 1.5-m beam trawl with a 5 mm mesh size in the cod end. To complement this information, sampling was also performed using 15-m fyke nets (8 mm mesh size in the cod end). For each sample, abiotic (temperature, salinity, pH, oxygen, turbidity, river flow, wind speed and depth) and biotic (macro crustacean species abundances) were recorded. Throughout the study, 28 fish species belonging to 20 families were collected. Fish catches were dominated by juveniles, especially Young-Of-the-Year (YOY) for the majority of the species. According to the Index of Relative Importance (IRI), common goby Pomatoschistus microps, flounder Platichtys flesus, sprat Sprattus sprattus, sea-bass Dicentrarchus labrax and plaice Pleuronectes platessa were the most abundant species, together accounting for 99.2% of the total IRI. Estuarine residents (ER = 66.2%) and marine juvenile migrants species (MJ = 31.4%) were the most important ecological guilds. The structure of the fish assemblage and its relationship to environmental variables was examined using multivariate techniques. Cluster and non-metric multidimensional scaling (nMDS) analysis defined six distinct groups in the Canche estuary, which are discriminated by specific species (SIMPER). Spatio-temporal variations in fish assemblage structure reflect the density peaks of the most abundant species. Spearman rank correlations and canonical correspondence analysis (CCA) showed that among the ten environmental variables examined, temperature, salinity and Crangon crangon (a potential predator for YOY fish or prey for older ones) are the three most important factors influencing fish species richness and abundances. Our observations reinforce the idea that certain fish species may have different life history styles in different geographic areas. The present study highlights the necessity of a better knowledge of the connectivity between estuaries and adjacent marine areas. The Canche constitutes an important ecosystem for fishes and as it is subject to little anthropogenic disturbance; its ichthyofauna can be viewed as a reference or normal assemblage for small temperate macrotidal estuaries.