Dynamical evolution of NEAs: Close encounters, secular perturbations and resonances

We discuss the main mechanisms affecting the dynamical evolution of Near-Earth Asteroids (NEAs) by analyzing the results of three numerical integrations over 1 Myr of the NEA (4179) Toutatis. In the first integration the only perturbing planet is the Earth. So the evolution is dominated by close encounters and looks like a random walk in semimajor axis and a correlated random walk in eccentricity, keeping almost constant the perihelion distance and the Tisserand invariant. In the second integration Jupiter and Saturn are present instead of the Earth, and the 3/1 (mean motion) and v ₆ (secular) resonances substantially change the eccentricity but not the semimajor axis. The third, most realistic, integration including all the three planets together shows a complex interplay of effects, with close encounters switching the orbit between different resonant states and no approximate conservation of the Tisserand invariant. This shows that simplified 3-body or 4-body models cannot be used to predict the typical evolution patterns and time scales of NEAs, and in particular that resonances provide some fast-track dynamical routes from low-eccentricity to very eccentric, planet-crossing orbits.On leave from the Department of Mathematics, University of Pisa, Via Buonarroti 2, 56127 Pisa, Italy, thanks to the G. Colombo fellowships of the European Space Agency.     

3D numerical modelling of fluid flow in the Val-d’Or orogenic gold district: major crustal shear zones drain fluids from overpressured vein fields

Fluid flow patterns have been determined using oxygen isotope isopleths in the Val-d’Or orogenic gold district. 3D numerical modelling of fluid flow and oxygen isotope exchange in the vein field shows that the fluid flow patterns can be reproduced if the lower boundary of the model is permeable, which represents middle or lower crustal rocks that are infiltrated by a metamorphic fluid generated at deeper levels. This boundary condition implies that the major crustal faults so conspicuous in vein fields do not act as the only major channel for upward fluid flow. The upper model boundary is impermeable except along the trace of major crustal faults where fluids are allowed to drain out of the vein field. This upper impermeable boundary condition represents a low-permeability layer in the crust that separates the overpressured fluid from the overlying hydrostatic fluid pressure regime. We propose that the role of major crustal faults in overpressured vein fields, independent of tectonic setting, is to drain hydrothermal fluids out of the vein field along a breach across an impermeable layer higher in the crust and above the vein field. This breach is crucial to allow flow out of the vein field and accumulation of metals in the fractures, and this breach has major implications for exploration for mineral resources. We propose that tectonic events that cause episodic metamorphic dehydration create a short-lived pulse of metamorphic fluid to rise along zones of transient permeability. This results in a fluid wave that propagates upward carrying metals to the mineralized area. Earthquakes along crustal shear zones cause dilation near jogs that draw fluids and deposit metals in an interconnected network of subsidiary shear zones. Fluid flow is arrested by an impermeable barrier separating the hydrostatic and lithostatic fluid pressure regimes. Fluids flow through the evolving and interconnected network of shear zones and by advection through the rock matrix. Episodic breaches in the impermeable barrier along the crustal shear zones allow fluid flow out of the vein field.     

Multibeam sonar backscatter data processing

Multibeam sonar systems now routinely record seafloor backscatter data, which are processed into backscatter mosaics and angular responses, both of which can assist in identifying seafloor types and morphology. Those data products are obtained from the multibeam sonar raw data files through a sequence of data processing stages that follows a basic plan, but the implementation of which varies greatly between sonar systems and software. In this article, we provide a comprehensive review of this backscatter data processing chain, with a focus on the variability in the possible implementation of each processing stage. Our objective for undertaking this task is twofold: (1) to provide an overview of backscatter data processing for the consideration of the general user and (2) to provide suggestions to multibeam sonar manufacturers, software providers and the operators of these systems and software for eventually reducing the lack of control, uncertainty and variability associated with current data processing implementations and the resulting backscatter data products. One such suggestion is the adoption of a nomenclature for increasingly refined levels of processing, akin to the nomenclature adopted for satellite remote-sensing data deliverables.     

Wind-forced circulation model and water exchanges through the channel in the Bay of Toulon

A hydrodynamic model of the Bay of Toulon has been developed for use as a post-accident radionuclide dispersion simulation tool. Located in a Mediterranean urban area, the Bay of Toulon is separated into two basins by a 1.4-km long seawall. The Little Bay is semi-enclosed and connected to the Large Bay by a fairway channel. This channel is the site of significant water mass exchange as a result of both wind-driven currents and bathymetry. It is therefore a focal point for marine contamination. As part of the model calibration and validation process, the first step consisted of studying the water mass exchange between the two basins. An Acoustic Doppler Current Profiler was moored in the channel for 1 year. The present study analyses in situ data to determine the current intensity and direction, and also to better understand the vertical current profile, which is highly correlated with meteorological forcing. Comparisons of model-generated and measured data are presented, and various atmospheric forcing datasets are used to enhance computed results. It appears that accurate meteorological forcing data is needed to enhance the accuracy of the hydrodynamic model. This channel is an important location for water mass renewal in the Bay of Toulon, and model results are used to quantify these exchanges. The mean calculated annual water exchange time is approximately 3.4 days. However, this duration is strongly wind dependent and shortens during windy winter months. It ranges from 1.5 days during strong wind periods to 7.5 days during calm weather. Residence time values calculated through tracer dispersion modelling after release at the back of the Little Bay are found to be comparable to the mean exchange time values, especially for windy conditions.     

Heavy metal inflow into the floodplains at the mouth of the river Weiße Elster (Central Germany)

From 1998 to 2002 investigations were carried out to estimate both the quantity and quality (heavy metal contents) of suspended matter loads discharged into a floodplain area measuring approx. 5 km² of the Central German river Weiße Elster, a tributary of the river Saale in the Elbe river basin. Flood sediments, suspended particulate matter, and floodplain soils were investigated especially for the main pollutants Cd, Zn, Cu, Pb, Cr, and Hg. Supplementary gamma spectroscopy examinations were performed to help to identify the age of deposits in order to estimate sedimentation rates for last decades. The recent flood sediments are contaminated by up to 33 times the geogenic background (Cd). Up to 55% of the total annual SPM load of the river Weiße Elster is relocated into the investigated retention area in the flood‐stricken year 2002. The sedimentation rates for the last 50 years vary between 0.5 mm/a and 1 mm/a in far away and rarely flooded parts and between 1 mm/a and 2 mm/a in the frequently flooded parts of the inflow of the retention area.     

Chemical and textural equilibration of garnet during amphibolite facies metamorphism: The influence of coupled dissolution–reprecipitation

Metamorphic equilibration requires chemical communication between minerals and may be inhibited through sluggish volume diffusion and or slow rates of dissolution in a fluid phase. Relatively slow diffusion and the perceived robust nature of chemical growth zoning may preclude garnet porphyroblasts from readily participating in low‐temperature amphibolite facies metamorphic reactions. Garnet is widely assumed to be a reactant in staurolite‐isograd reactions, and the evidence for this has been assessed in the Late Proterozoic Dalradian pelitic schists of the Scottish Highlands. The 3D imaging of garnet porphyroblasts in staurolite‐bearing schists reveals a good crystal shape and little evidence of marginal dissolution; however, there is also lack of evidence for the involvement of either chlorite or chloritoid in the reaction. Staurolite forms directly adjacent to the garnet, and its nucleation is strongly associated with deformation of the muscovite‐rich fabrics around the porphyroblasts. “Cloudy” fluid inclusion‐rich garnet forms in both marginal and internal parts of the garnet porphyroblast and is linked both to the production of staurolite and to the introduction of abundant quartz inclusions within the garnet. Such cloudy garnet typically has a Mg‐rich, Mn‐poor composition and is interpreted to have formed during a coupled dissolution–reprecipitation process, triggered by a local influx of fluid. All garnet in the muscovite‐bearing schists present in this area is potentially reactive, irrespective of the garnet composition, but very few of the schists contain staurolite. The staurolite‐producing reaction appears to be substantially overstepped during the relatively high‐pressure Barrovian regional metamorphism reflecting the limited permeability of the schists in peak metamorphic conditions. Fluid influx and hence reaction progress appear to be strongly controlled by subtle differences in deformation history. The remaining garnet fails to achieve chemical equilibrium during the reaction creating distinctive patchy compositional zoning. Such zoning in metamorphic garnet created during coupled dissolution–reprecipitation reactions may be difficult to recognize in higher grade pelites due to subsequent diffusive re‐equilibration. Fundamental assumptions about metamorphic processes are questioned by the lack of chemical equilibrium during this reaction and the restricted permeability of the regional metamorphic pelitic schists. In addition, the partial loss of prograde chemical and textural information from the garnet porphyroblasts cautions against their routine use as a reliable monitor of metamorphic history. However, the partial re‐equilibration of the porphyroblasts during coupled dissolution–reprecipitation opens possibilities of mapping reaction progress in garnet as a means of assessing fluid access during peak metamorphic conditions.     

Coloured dissolved organic matter (CDOM) in Southern North Sea waters: Optical characterization and possible origin

The variability and origin of the Coloured Dissolved Organic Matter (CDOM) were studied in the Belgian coastal and adjacent areas including offshore waters and the Scheldt estuary, through the parameters: absorption at 375 nm, a_CDOM(375), and the slope of the absorption curve, S. a_CDOM(375) varied between 0.20 and 1.31 m⁻¹ and between 0.97 and 4.30 m⁻¹ in the marine area and Scheldt estuary, respectively. S fluctuated between 0.0101 and 0.0203 nm⁻¹ in the marine area and between 0.0167 and 0.0191 nm⁻¹ in the Scheldt estuary. The comparative analysis of a_CDOM(375) and S variations evidenced different origins of CDOM in the BCZ. The Scheldt estuarine waters showed decreasing a_CDOM(375) values with increasing salinity but constant S value of ∼0.018 nm⁻¹ suggesting a dominant terrestrial origin of CDOM. On the contrary, samples collected in the marine domain showed a narrow range of a_CDOM(375) but highly variable S suggesting the additional presence of autochthonous sources of CDOM. This source was evidenced based on the sorting of the marine offshore data according to the stage of the phytoplankton bloom when they were collected. A clear distinction was made between CDOM released during the growth stage characterized by high S (∼0.017 nm⁻¹) and low a_CDOM(375) and the decay phase characterized by low S (∼0.013 nm⁻¹) and high a_CDOM(375). This observation was supported by CDOM measurements performed on pure phytoplankton cultures which showed increased CDOM release along the wax and wane of the bloom but decreasing S. We concluded that the high variability of the CDOM signature in offshore waters is explained by the local biological production and processing of CDOM.