In this paper, we identify a crucial numerical problem in sigma coordinate models, leading to unacceptable spurious diapycnal mixing. This error is a by-product of recent advances in numerical methods, namely the implementation of high-order diffusive advection schemes. In the case of ROMS, spurious mixing is produced by its third-order upwind advection scheme, but our analysis suggests that all diffusive advection schemes would behave similarly in all sigma models. We show that the common idea that spurious mixing decreases with resolution is generally false. In a coarse-resolution regime, spurious mixing increases as resolution is refined, and may reach its peak value when eddy-driven lateral mixing becomes explicitly resolved. At finer resolution, diffusivities are expected to decrease but with values that only become acceptable at resolutions finer than the kilometer. The solution to this problem requires a specifically designed advection scheme. We propose and validate the RSUP3 scheme, where diffusion is split from advection and is represented by a rotated biharmonic diffusion scheme with flow-dependent hyperdiffusivity satisfying the Peclet constraint. The rotated diffusion operator is designed for numerical stability, which includes improvements of linear stability limits and a clipping method adapted to the sigma-coordinate. Realistic model experiments in a southwest Pacific configuration show that RSUP3 is able to preserve low dispersion and diffusion capabilities of the original third-order upwind scheme, while preserving water mass characteristics. There are residual errors from the rotated diffusion operator, but they remain acceptable. The use of a constant diffusivity rather than the Peclet hyperdiffusivity tends to increase these residual errors which become unacceptable with Laplacian diffusion. Finally, we have left some options open concerning the use of time filters as an alternative to spatial diffusion. A temporal discretization approach to the present problem (including implicit discretization) will be reported in a following paper. 相似文献
Based upon a re-interpretation of previous data and a new field campaign, a structural evolution is proposed for the early history of Piton de la Fournaise volcano from 500,000 to 50,000 years. Conceptually, it is shown that the formation of a caldera in which lava flows are contained inside the caldera depression, gives time for erosion to excavate deep canyons on the external slopes of the volcano, for example, the Rivière des Remparts, the Rivière Langevin and the Rivière de l'Est canyons on Piton de la Fournaise volcano. These canyons are infilled when lavas, filling the caldera and overflowing its rim, are able again to flow on the external slopes of the volcano. In the past, this excavating/infilling process has occurred twice following the formation of the Rivière des Remparts and Morne Langevin calderas. The formation of the third caldera, the Plaine des Sables caldera, was followed by the excavation of the current canyons. In addition to this process, two large landslides have been documented in the field. The first, which happened about 300,000 years ago, is apparently the first episode of the break up of Piton de la Fournaise volcano, predating the formation of the four large calderas. The second landslide, which occurred 150,000 years ago and is considered to be less extensive, has carried away the entire southern flank of the Rivière des Remparts caldera. 相似文献
Trace element and isotopic characteristics of late Carboniferous to early Permian minettes and kersantites have been determined. These lamprophyres have been sampled throughout the Western European Hercynian orogen, from Brittany to the west to Schwarzwald to the east. In spite of sharp petrological differences reflected by mineralogy and major element geochemistry, minettes and kersantites exhibit close identity with respect to trace element and isotopic features. These features comprise enrichment in incompatible elements, highCs/Rb and lowCe/Pb ratios, Ta and Ti relative depletion, high abundance in transition elements and highNi/Mg ratios. Pb isotope ratios are undistinguishable from those measured on Hercynian continental crust. Initial143Nd/144Nd ratios are between0.5120 (εi −5) and0.5122 (εi −1) for minettes and kersantites whereas initial87Sr/86Sr ratios vary between 0.7055–0.710 for minettes and 0.707–0.708 for kersantites. No simple mixing relations are visible on RbSr and SmNd isochron diagrams. The exceptional homogeneity of these geochemical characteristics along a 1000 km traverse does not allow for an hypothesis of enrichment through upper level assimilation and thus leads to propose that these rocks originated through melting of a mantle enriched by recycling of crustal material. 相似文献
During the last decade, it has been shown that accurate analyses of calcite twinning in thin section may allow computation of the paleo-deviatoric stress tensor associated with this deformation. Following from the work of Larroque and Laurent (1988), and Lacombe et al. (1989), this is a new application of the Etchecopar inverse method. Weakly deformed limestones from about 150 to 200 km north of the Pyrenean chain (South France) were studied. The Etchecopar inverse method, which can be used on a microcomputer, is comparable to the method of analysis of striated fault planes proposed by Etchecopar et al. in 1981, as it determines the principal stress orientations and the stress ellipsoid shape ratio φ = (σ2 − σ3)/(σ1 − σ3). The Etchecopar method using calcite twins has been selected for this study because it is suitable for the determination of polyphase deformation in coarse-grained limestones (Lacombe et al. 1989).
Ten samples of oolitic limestones, cemented with sparry calcite, from the Bajocian (170 Ma) and the Oligocene (30 Ma) were studied. For each of them, only one or two stress tensors were determined: the first is characterized by a principal compressive stress axis σ1 which is horizontal and approximately N-S. Local departures from this orientation have been observed and are interpreted as local stress perturbations. For the second, which was defined in the Oligocene, but occurred only locally in the Jurassic limestones, the position of σ1 is nearly vertical; the direction of extension σ3 is either NE-SW or undefined into the horizontal plane (σ2 = σ3); this radial extension is interpreted as the result of two or more superposed extensive deformation phases with vertical σ1. The N-S compression was never found in the Oligocene limestones, and from its orientation, it is assumed to correspond to the Pyrenean deformation. These results are compared to the more complex deformation pattern obtained from classical microtectonical studies (Arthaud and Choukroune, 1972; Bonijoly and Blès, 1983; Granier and Blès, 1988; Blès et al., 1989). Only the major tectonic events are recorded by calcite twinning. Finally, the suitability of calcite twin analyses for the determination of complex polyphase deformation is discussed. 相似文献
The Moselle river flows in the north-east of France, from Vosges Mountains to neighboring countries Luxembourg and Germany.
One of its tributaries, the Fensch river, drains a highly industrialized watershed, strongly impacted by mining, smelting
and surfacing activities. The objective of this work, part of a general research program on Moselle watershed (Zone Atelier
Moselle) was to assess the impact of the polluted Fensch river on the global quality of the Moselle river. For that purpose,
water, sediments and suspended particulate matter were sampled in both rivers, upstream and downstream the junction. Four
main sampling campaigns were carried out, in winter during a flood event and in spring at low water level. On a first step,
mineralogical analyses (XRD and FTIR) and chemical analyses (ICP-MS, ICP-AES), were performed on sediments, suspended particulate
matters and filtered waters. Major and trace elements concentrations were obtained on two different granulometric fractions
(0–2 mm and 0–50 μm) revealing the enrichment in heavy metals of fine particles. From one collecting campaign to another,
seasonal variations could be evidenced on suspended matter composition even though major minerals (calcite, quartz and kaolinite)
were always present. Furthermore, spatial variations were evidenced for Fensch and Moselle downstream sediments. Thus, very
fine-grained sediments, poorly crystallized, displaying at the same time higher metal concentrations and higher organic matter
content than in Fensch river material, were collected downstream, in a low hydrodynamic conditions zone, assumed as a preferential
sedimentary zone or even as a placer. Strong correlations could be revealed between iron content and contaminant concentrations,
confirming the origin of polluted material. 相似文献
Caledonian eclogite-facies metamorphism partially reworking Grenvillian granulite-facies anorthosite allows us to study the processes of garnet reequilibration at high pressure and to reconstruct the evolution of the unit near metamorphic peak conditions. Our results indicate that eclogite-facies metamorphism happened in two successive phases: first, inherited granulitic garnet was fractured and reequilibrated from their boundaries (crystal or fracture rims); then eclogite-facies minerals were crystallised in the fractures as overgrowths on inherited garnets. The reequilibration of inherited garnets is achieved through Fe2+Mg−1 exchange, whereas eclogite-facies garnets crystallised during the subsequent phase are notably richer in Ca than un- and re-equilibrated granulitic garnet. Pseudosection construction shows that this lack in Ca reequilibration cannot be related to variations in thermodynamic conditions (aH2O, reacting system composition) between the two phases. From the compilation of the available data, the reequilibration of granulitic garnet seems to be controlled by the inefficient intra- and inter-granular transport properties of Ca compared to Fe2+ and Mg. While these kinetic factors confine garnet reequilibration to Fe2+Mg−1 exchange, the extent of reequilibration along this exchange vector is controlled by partitioning with adjacent omphacite. On the contrary to the diffusional reequilibration of granulitic garnet that lasted for several My according to our modelling of the diffusional relaxation, the strong compositional gradients between eclogite-facies and reequilibrated garnets, which are almost unaffected by diffusional reequilibration, provide evidence that rapid exhumation followed the crystallisation of eclogite-facies minerals. We propose that the movement reversal itself, from burial to exhumation, and associated deformation and fluid flow, triggered this crystallisation event. The resulting evolution near metamorphic peak conditions is therefore strongly asymmetrical: on the one hand, the prograde diffusional relaxation profiles indicate slow movement during the last stages of burial, whereas the unaffected retrograde overgrowth indicates fast exhumation rates. 相似文献
This study focuses on the response of lacustrine littoral chironomid communities to late glacial changes in limnological,
environmental and climate conditions in the Mediterranean context. Late glacial chironomid (Diptera: Chironomidae) assemblages,
organic petrography and geochemistry were analysed in a sediment core from the littoral zone of Lago dell’Accesa (Tuscany,
Italy), where the lake-level fluctuations and the vegetation history have been previously reconstructed. Comparison of the
chironomid stratigraphy to other proxies (pollen assemblages, organic petrography and geochemistry, lake-level) and regional
climate reconstruction suggested the predominant influence of lake-level changes on the littoral chironomid fauna. The main
lowering events that occurred during the Oldest and the Younger Dryas were followed by higher proportions of taxa typical
of littoral habitats. A complementary study of organic matter suggested the indirect impact of lake-level on the chironomids
through changes in humic status and habitat characteristics, such as the type of substrate and aquatic macrophyte development.
Several chironomid taxa, such as Glyptotendipes, Microtendipes and Cricotopus type patens, were identified as possible indicators of low lake-level in the late glacial records. Nevertheless, this study suggested
that parallel analyses of organic matter and chironomid assemblages may be needed to circumvent misinterpretation of littoral
chironomid assemblage stratigraphy. There was a weak response of the chironomid assemblages to small lake-level lowerings
that corresponded to the Older Dryas and Preboreal oscillations. A higher level of determination, e.g. to the species group
level, may be necessary to increase the sensibility of the indicators to lake-level changes. 相似文献
The second field campaign of the Cloud Ice Mountain Experiment (CIME) project took place in February 1998 on the mountain Puy de Dôme in the centre of France. The content of residual aerosol particles, of H2O2 and NH3 in cloud droplets was evaluated by evaporating the drops larger than 5 μm in a Counterflow Virtual Impactor (CVI) and by measuring the residual particle concentration and the released gas content. The same trace species were studied behind a round jet impactor for the complementary interstitial aerosol particles smaller than 5 μm diameter. In a second step of experiments, the ambient supercooled cloud was converted to a mixed phase cloud by seeding the cloud with ice particles by the gas release from pressurised gas bottles. A comparison between the physical and chemical characteristics of liquid drops and ice particles allows a study of the fate of the trace constituents during the presence of ice crystals in the cloud.In the present paper, an overview is given of the CIME 98 experiment and the instrumentation deployed. The meteorological situation during the experiment was analysed with the help of a cloud scale model. The microphysics processes and the behaviour of the scavenged aerosol particles before and during seeding are analysed with the detailed microphysical model ExMix. The simulation results agreed well with the observations and confirmed the assumption that the Bergeron–Findeisen process was dominating during seeding and was influencing the partitioning of aerosol particles between drops and ice crystals. The results of the CIME 98 experiment give an insight on microphysical changes, redistribution of aerosol particles and cloud chemistry during the Bergeron–Findeisen process when acting also in natural clouds. 相似文献
Localized shear zones along low-angle normal faults have been identified in regions of extension at the brittle-ductile transition of the continental crust. The possibility of the strain localizing at a depth of 10 km is interpreted here as a consequence of an increase in the equivalent shear stress applied to the flow of the lower crust. This enhancement of the flow stress is seen as a prerequisite for the triggering of brittle deformation mechanisms leading to strain localization. The lower crust rheology used to examine this stress increase is strain-rate, temperature and grain-size dependent, due to the coupling of dislocation and diffusion creep. The model structure proposed consists of a top layer, the upper crust, gliding rigidly above a bottom layer, the lower crust, which deforms in simple shear. During a short time interval (1400 years), the equivalent shear stress is found to increase by a factor of up to 3 (67 MPa for anorthite and 17 MPa for quartz). For anorthite, this stress could explain the activation of a Mohr-Coulomb failure with a friction coefficient of 0.2, which is reasonable at the depth of 10 km. Dislocation creep is activated during a rapid change in the prescribed velocity, whereas diffusion creep dominates if the velocity is held constant, highlighting the importance of grain-size sensitivity for lower crustal rheology. 相似文献