Efficiency of a two-step upscaling method for permeability evaluation at Darcy and pore scales

This work presents a new subdivision method to upscale absolute permeability fields. This process, called two-step method, consists in (i) solving micro-scale equations on subdomains obtained from the full domain regular decomposition and (ii) solve a second upscaling with Darcy’s law on the permeability fields obtained in the first step. The micro-scale equations used depend on the case studied. The two-step upscaling process is validated on randomly generated Darcy-scale permeability fields by measuring the numerical error induced by upscaling. The method is then applied to real domains obtained from sandstone micro-tomographic images. The method specificities due to pore-space structure are discussed. The main advantage of the two-step upscaling method resides in the drastic reduction of computational costs (CPU time and memory usage) while maintaining a numerical error similar to that of other upscaling procedures. This new upscaling method may improve permeability predictions by the use of finer meshes or larger sample volumes.     

Twenty years of evolution for the DORIS permanent network: from its initial deployment to its renovation

The ground network is one of the major components of the DORIS system. Its deployment, managed by the French national mapping agency [Institut Géographique National, (IGN)], started in 1986 at a sustained pace that allowed it to reach 32 stations upon the launch of the first DORIS-equipped satellite (SPOT-2) in 1990. For the first generation of transmitting antennas, the installation procedures were adapted to the decimetre performance objective for the DORIS system. During the second era of the deployment of an even denser network, the antenna support layouts gradually evolved towards a better quality, thus improving the long-term stability of the antenna reference point, and a new antenna model allowed a more accurate survey. As the positioning accuracy of the DORIS system improved, it was necessary to review the antenna stability for the whole network. A first stability estimation, using criteria like antenna model and support design, was followed by a major renovation effort which started in 2000 and is now almost complete. In 6 years, through the renovation or installation of 43 stations and the implementation of new installation procedures to meet more stringent stability requirements, significant improvement in network quality was achieved. Later a more analytical approach, taking into account the characteristics of each element that support the antenna, has been taken to assess the potential stability of all DORIS occupations. IGN is also in charge of its operational maintenance, an intensive activity on account of the significant failure rate of the successive generations of equipment. Nevertheless, thanks to its unique density and homogeneity, DORIS has maintained a very good coverage rate of the satellite orbits. Through 38 well-distributed current co-locations with the Global Positioning System, Satellite Laser Ranging and Very Long Baseline Interferometry techniques in its current 56-station network, DORIS contributes significantly to the realisation of the International Terrestrial Reference System. DORIS stations in areas where no other space geodesy technique is available provide a significant contribution to the study of plate tectonics. Many stations co-located with tide gauges contribute to the monitoring of sea level changes. Although it has several advantages over similar techniques, there is still room for improvement in the DORIS network.     

Monte-Carlo Simulations of Cosmic-ray and Internal Radiation Effects on ISOCAM on Board ISO

One of the main limitations to the sensitivity of the infrared camera ISOCAM on-board the Infrared Space Observatory (ISO) comes from responsivity variations and glitches caused by the impacts of charged particles in photo-detectors. Glitch rate measurements, glitch properties and removal methods have already been addressed during the first ISO detector workshop(Madrid, 1998) and published in a special issue of Experimental Astronomy. It appeared that glitch rate and most of glitch properties could be reproduced by Monte-Carlo simulations. This is very interesting in order to predict before launch the effect of charged particles in photo-detectors operated in space. This paper presents results of Monte-Carlo simulations of radiation effects on ISOCAM detectors. Glitch rates, spatial and energetic properties of glitches have been computed and are compared with measured values. This revised version was published online in July 2006 with corrections to the Cover Date.     

Rb-Sr and Sm-Nd ages and isotopic geochemistry of Archaean granodioritic gneisses from Eastern Finland

Twenty granodioritic rocks and one amphibolitic enclave of the “basement” of the Suomussalmi-Kuhmo Archaean (2.65 Ga) greenstone belts (central-eastern Finland), have been chosen together with one greenstone sample for Rb-Sr and Sm-Nd geochronological and isotopic studies.The granitoïd rocks are subdivided into three groups: two generations of grey gneisses and a post-belt augen gneiss. The Rb-Sr ages of the first and second generation of grey gneisses are 2.86 ± 0.09 and 2.62 ± 0.07 Ga, respectively. These results are corroborated by Sm-Nd data. The post-belt augen gneiss gives an age of 2.51 ± 0.11 Ga. The results show that the two generations of grey gneisses, the greenstone belts and the post-greenstone augen gneiss, were developed over a period > 350 Ma. The two generations of grey gneisses show identical I_Sr values (0.7023 ± 8 and 0.7024 ± 6) which contrast with that of the augen gneiss (0.7049 ± 8). The low I_Sr and the near-chondritic ?_T^CHUR values indicate that the grey gneisses cannot derived from much older continental materials. Trace element studies suggest that these grey gneisses have had a multi-stage development. The augen gneiss with a moderately high I_Sr is likely to be derived from a granodiorite originated by partial melting of older sialic crust. The more probable parent rock seems to be the first generation grey gneisses. The I_Sr and average Rb/Sr values preclude the greenstone belt and the second generation of grey gneisses as the protolith.     

Converging P-T paths of Mesozoic HP-LT metamorphic units (Diego de Almagro Island, Southern Chile): evidence for juxtaposition during late shortening of an active continental margin

Summary On Diego de Almagro Island in Chilean Patagonia (51°30S), a convergent strike slip zone, the Seno Arcabuz shear zone, separates the Diego de Almagro Metamorphic Complex from very low grade metagreywackes in the east, which were intruded by Jurassic granitoids. The Diego de Almagro Metamorphic Complex is composed of a metapsammopelitic sequence containing blueschist intercalations in the west and (garnet) amphibolite lenses in the east. Peak metamorphic conditions (stage I) at 9.5–13.5kbar, 380–450°C in the blueschist and at 11.2–13.2, 460–565°C in the amphibolite indicate subduction and accretion at different positions within the deepest part of the accretionary wedge. A K–Ar age of 117±28Ma of amphibole approximately dates the peak of metamorphism in the amphibolite. The early retrograde stage of metamorphism occurred under static conditions and resulted in localized equilibration (stage II) at 6.3–9.6kbar, 320–385°C in the blueschist and 6.1–8.4kbar, 310–504°C in the amphibolite. Both P-T paths converge within a midcrustal level.In contrast, an orthogneiss of trondhjemitic composition occurring within the Seno Arcabuz shear zone is associated with a garnet mica-schist containing a high temperature/intermediate pressure assemblage formed at 4.9–6.5kbar, 580–690°C. A muscovite K–Ar age of 122.2±4.6Ma dates cooling after this event which is related to a concomitant magmatic arc. These rocks were overprinted by a mylonitic deformation, which is caused by convergent strike slip shearing and ends during formation of a retrograde phengite-chlorite-stilpnomelane assemblage at a minimum pressure of approximately 5.7kbar (at 300°C).Zircon fission track ages from rocks of the Seno Arcabuz shear zone are 64.9±2.7 and 64.9±2.7Ma; they record the end of shearing in the Seno Arcabuz shear zone that juxtaposed all rocks in the middle crust. Zircon fission track ages ranging from 78 to 105Ma in the South Patagonian batholith to the east indicate earlier cooling through 280°C. The rocks of the Diego de Almagro Metamorphic Complex were initially slowly exhumed and resided at a midcrustal level before being emplaced via shearing in the Seno Arcabuz shear zone. Apatite fission track ages (54±8Ma) from the Seno Arcabuz shear zone show that exhumation and cooling rates increased after this event. The incorporation of continental crust within the subduction system was a late process, which modified the Cretaceous accretionary wedge, resulting in considerable shortening of the convergent margin.     

Événements panafricains dans l'Adrar Souttouf (Sahara marocain)

The Adrar Souttouf belt (Western Sahara) lies in the northern part of the Mauritanide belt and was considered as a part of the Hercynian Mauritanian-Appalachian belt. Recently the presence of a Panafrican protolith (595 Ma) was evidenced by the age of inherited zircons from an eclogitic unit and led to consider complex Panafrican events included in the Hercynian belt, as in the central and southern Mauritanide belt. Our new field survey completed by ten KAr dating (whole rocks and separate minerals) confirmed that multiphased Panafrican formations outcrop within the center of the Adrar Souttouf massif. Diverse events ranged from ca. 1000 to 500 Ma, some of them being indicative of an oceanic type environment and metamorphism. The 1000 Ma metamorphic event is recorded for the first time in the West African craton. A large Hercynian remobilisation is also evidenced. To cite this article: M. Villeneuve et al., C. R. Geoscience 338 (2006).     

A multiproxy evaluation of Holocene environmental change from Lake Igaliku, South Greenland

This is the first integrated multiproxy study to investigate climate, catchment evolution and lake ecology in South Greenland. A 4-m-long sedimentary sequence from Lake Igaliku (61o 00?? N, 45o 26?? W, 15?m asl) documents major environmental and climatic changes in south Greenland during the last 10?ka. The chronology is based on a ²¹⁰Pb and ¹³⁷Cs profile and 28 radiocarbon dates. The paleoenvironmental history is interpreted on the basis of magnetic susceptibility, grain size, total organic carbon, total nitrogen and sulphur, sedimentation rates, pollen, and diatom assemblages. The basal radiocarbon date at ca. 10?cal?ka BP provides a minimum age for the deglaciation of the basin, which is followed by ~500?years of high sedimentation rates in a glacio-marine environment. After the glacio-isostatic emergence of the basin ca. 9.5?cal?ka BP, limnological and terrestrial proxies suggests early warmth, which may have been interrupted by a cold, dry and windy period between 8.6 and 8.1?cal?ka BP. A dry and windy event ~5.3?C4.8?cal?ka BP preceded the Neoglacial transition at Lake Igaliku, which is characterized by a shift toward moister and perhaps cooler conditions ~4.8?cal?ka BP, causing major changes in terrestrial and aquatic ecological conditions. Significant cooling is documented after ~3?cal?ka BP. Since ~1?cal?ka BP the climatic-driven changes were overprinted by the human influence of Norse and recent agriculture.     

Relative contribution of soil moisture and snow mass to seasonal climate predictability: a pilot study

Land surface hydrology (LSH) is a potential source of long-range atmospheric predictability that has received less attention than sea surface temperature (SST). In this study, we carry out ensemble atmospheric simulations driven by observed or climatological SST in which the LSH is either interactive or nudged towards a global monthly re-analysis. The main objective is to evaluate the impact of soil moisture or snow mass anomalies on seasonal climate variability and predictability over the 1986–1995 period. We first analyse the annual cycle of zonal mean potential (perfect model approach) and effective (simulated vs. observed climate) predictability in order to identify the seasons and latitudes where land surface initialization is potentially relevant. Results highlight the influence of soil moisture boundary conditions in the summer mid-latitudes and the role of snow boundary conditions in the northern high latitudes. Then, we focus on the Eurasian continent and we contrast seasons with opposite land surface anomalies. In addition to the nudged experiments, we conduct ensembles of seasonal hindcasts in which the relaxation is switched off at the end of spring or winter in order to evaluate the impact of soil moisture or snow mass initialization. LSH appears as an effective source of surface air temperature and precipitation predictability over Eurasia (as well as North America), at least as important as SST in spring and summer. Cloud feedbacks and large-scale dynamics contribute to amplify the regional temperature response, which is however, mainly found at the lowest model levels and only represents a small fraction of the observed variability in the upper troposphere.     

Summer interactions between weather regimes and surface ocean in the North-Atlantic region

This study aims at understanding the summer ocean-atmosphere interactions in the North Atlantic European region on intraseasonal timescales. The CNRMOM1d ocean model is forced with ERA40 (ECMWF Re-Analysis) surface fluxes with a 1-h frequency in solar heat flux (6 h for the other forcing fields) over the 1959–2001 period. The model has 124 vertical levels with a vertical resolution of 1 m near the surface and 500 m at the bottom. This ocean forced experiment is used to assess the impact of the North Atlantic weather regimes on the surface ocean. Composites of sea surface temperature (SST) anomalies associated with each weather regime are computed and the mechanisms explaining these anomalies are investigated. Then, the SST anomalies related to each weather regime in the ocean-forced experiment are prescribed to the ARPEGE Atmosphere General Circulation Model. We show that the interaction with the surface ocean induces a positive feedback on the persistence of the Blocking regime, a negative feedback on the persistence of the NAO-regime and favours the transition from the Atlantic Ridge regime to the NAO-regime and from the Atlantic Low regime toward the Blocking regime.