The shape of the small satellites of Saturn: Gravitational equilibrium vs solid-state strength

In the set of small satellites of Saturn recently imaged by the Voyager probes, we can observe the transition from irregularly-shaped, strength-dominated objects to larger, gravity-dominated bodies with shapes roughly fitting the theoretical equilibrium figures. The transition occurs for a radius of 100±50 km, corresponding to a typical material strength of the order of 10⁷ dynes cm^?2. We discuss briefly the cases of Mimas, Enceladus, Hyperion, Phoebe and the small coorbital and F-ring shepherding moons, showing that an analysis of the shape data can often provide interesting results on the physical properties, origin and collisional history of these objects.     

The solution of Kepler's equation,I

Methods of iteration are discussed in relation to Kepler's equation, and various initial guesses are considered, with possible strategies for choosing them. Several of these are compared; the method of iteration used in the comparisons has local convergence of the fourth order.WANG Laboratories, Inc.     

Crustal motion and deformation in Greece from a decade of GPS measurements, 1993–2003

The Hellenic plate boundary region, located in the collision zone between the Nubian/Arabian and Eurasian lithospheric plates, is one of the seismo-tectonically most active areas of Europe. During the last 15 years, GPS measurements have been used to determine the crustal motion in the area of Greece with the aim to better understand the geodynamical processes of this region. An extended reoccupation network covering whole Greece has been measured periodically in numerous GPS campaigns since the late eighties, and a continuous GPS network has been operated in the region of the Ionian Sea since 1995. In this paper, we present a new detailed high-quality solution of continuous and campaign-type measurements acquired between 1993 and 2003. During the GPS processing, a special effort was made to obtain consistent results with highest possible accuracies and reliabilities. Data of 54 mainly European IGS and EUREF sites were included in the GPS processing in order to obtain results which are internally consistent with the European kinematic field and order to allow for a regional interpretation. After an overview of the results of the IGS/EUREF sites, the results from more than 80 stations in Greece are presented in terms of velocities, time series, trajectories and strain rates. Previous geodetic, geological and seismological findings are generally confirmed and substantially refined. New important results include the observation of deformation zones to the north and to the south of the North Aegean Trough and in the West Hellenic arc region, arc-parallel extension of about 19 mm/yr along the Hellenic arc, and compression between the Ionian islands and the Greek mainland. Due to continuous long-term observations of 4–8 years, it was possible to extract height changes from the GPS time series. In Greece, we observe a differential subsidence of the order of 2 mm/yr between the northern and central Ionian islands across the Kefalonia fault zone. The differential subsidence of the central Ionian islands with respect to the northwestern Greek mainland amounts to 4 mm/yr.     

Quaternary fault segmentation and interaction in the epicentral area of the 1561 earthquake (Mw = 6.4), Vallo di Diano, southern Apennines, Italy

The main structural characteristics of the Caggiano and Polla faults, exposed in the epicentral area of the 1561 earthquake (Mw = 6.4), southern Italy, have been investigated in detail to assess their spatial and temporal properties, and to evaluate their seismogenic potential. These right stepping normal faults show an overlap of about 7 km and an across strike separation of about 4 km. The geometric relationships between the Caggiano and Polla faults, but also the displacement distribution along each fault, demonstrate that they have been strongly interacting throughout the Pleistocene. Nevertheless, geological evidence of Holocene tectonic activity was mainly recognized along the Caggiano Fault (faulted late glacial deposits) and in the southernmost part of the Polla Fault (faulted deposits of probably Late Pleistocene age). This suggests that the Caggiano Fault can be considered as the most tectonically active fault in the Vallo di Diano Fault System. By calculating Coulomb stress changes, we have constrained modes of mechanical interactions between the two faults in a scenario compatible with the 1561 earthquake. This approach allows us to argue that both the Caggiano and the Polla Faults are probably linked at depth, and part of the same seismogenic structure which may be potentially responsible for composite ruptures with magnitude ≥ 6.5.     

Alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands)

Ancient fluvial successions often act as hydrocarbon reservoirs. Sub‐surface data on the alluvial architecture of fluvial successions are often incomplete and modelling is performed to reconstruct the stratigraphy. However, all alluvial architecture models suffer from the scarcity of field data to test and calibrate them. The purposes of this study were to quantify the alluvial architecture of the Holocene Rhine–Meuse delta (the Netherlands) and to determine spatio‐temporal trends in the architecture. Five north–south orientated cross‐sections, perpendicular to the general flow direction, were compiled for the fluvial‐dominated part of the delta. These sections were used to calculate the width/thickness ratios of fluvial sandbodies (SBW/SBT) and the proportions of channel‐belt deposits (CDP), clastic overbank deposits (ODP) and organic material (OP) in the succession. Furthermore, the connectedness ratio (CR) between channel belts was calculated for each cross‐section. Distinct spatial and temporal trends in the alluvial architecture were found. SBW/SBT ratios decrease by a factor of ca 4 in a downstream direction. CDP decreases from ca 0·7 (upstream) to ca 0·3 (downstream). OP increases from less than 0·05 in the upstream part of the delta to more than 0·25 in the downstream delta. ODP is approximately constant (0·4). CR is ca 0·25 upstream, which is approximately two times larger than in the downstream part of the delta. Furthermore, CDP in the downstream Rhine–Meuse delta increases after 3000 cal yr BP. These trends are attributed to variations in available accommodation space, floodplain geometry and channel‐belt size. For instance, channel belts tend to narrow in a downstream direction, which reduces SBW/SBT, CDP and CR. Tectonics cause local deviations in the general architectural trends. In addition, the positive correlation between avulsion frequency and the ratio of local to regional aggradation rate probably influenced alluvial architecture in the Rhine–Meuse delta. The Rhine–Meuse data set can be a great resource when developing more sophisticated models for alluvial architecture simulation, which eventually could lead to better characterizations of hydrocarbon reservoirs. To aid such usage of the Rhine–Meuse data set, constraints for relevant parameters are provided at the end of the paper.     

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.     

Discriminating between pore‐filling load and bed‐structure load: a new porosity‐based method,exemplified for the river Rhine

Sediments contained in the river bed do not necessarily contribute to morphological change. The finest part of the sediment mixture often fills the pores between the larger grains and can be removed without causing a drop in bed level. The discrimination between pore‐filling load and bed‐structure load, therefore, is of practical importance for morphological predictions. In this study, a new method is proposed to estimate the cut‐off grain size that forms the boundary between pore‐filling load and bed‐structure load. The method evaluates the pore structure of the river bed geometrically. Only detailed grain‐size distributions of the river bed are required as input to the method. A preliminary validation shows that the calculated porosity and cut‐off size values agree well with experimental data. Application of the new cut‐off size method to the river Rhine demonstrates that the estimated cut‐off size decreases in a downstream direction from about 2 to 0·05 mm, covariant with the downstream fining of bed sediments. Grain size fractions that are pore‐filling load in the upstream part of the river thus gradually become bed‐structure load in the downstream part. The estimated (mass) percentage of pore‐filling load in the river bed ranges from 0% in areas with a unimodal river bed, to about 22% in reaches with a bimodal sand‐gravel bed. The estimated bed porosity varies between 0·15 and 0·35, which is considerably less than the often‐used standard value of 0·40. The predicted cut‐off size between pore‐filling load and bed‐structure load (D_c,p) is fundamentally different from the cut‐off size between wash‐load and bed‐material load (D_c,w), irrespective of the method used to determine D_c,p or D_c,w. D_c,w values are in the order of 10^?1 mm and mainly dependent on the flow characteristics, whereas D_c,p values are generally much larger (about 10⁰ mm in gravel‐bed rivers) and dependent on the bed composition. Knowledge of D_c,w is important for the prediction of the total sediment transport in a river (including suspended fines that do not interact with the bed), whereas knowledge of D_c,p helps to improve morphological predictions, especially if spatial variations in D_c,p are taken into account. An alternative to using a spatially variable value of D_c,p in morphological models is to use a spatially variable bed porosity, which can also be predicted with the new method. In addition to the morphological benefits, the new method also has sedimentological applications. The possibility to determine quickly whether a sediment mixture is clast‐supported or matrix‐supported may help to better understand downstream fining trends, sediment entrainment thresholds and variations in hydraulic conductivity.     

Trace element distribution among rock-forming minerals from metamorphosed to partially molten basic igneous rocks in a contact aureole (Fuerteventura, Canaries)

Low pressure partial melting of basanitic and ankaramitic dykes gave rise to unusual, zebra-like migmatites, in the contact aureole of a layered pyroxenite–gabbro intrusion, in the root zone of an ocean island (Basal Complex, Fuerteventura, Canary Islands). These migmatites are characterised by a dense network of closely spaced, millimetre-wide leucocratic segregations. Their mineralogy consists of plagioclase (An_32–36), diopside, biotite, oxides (magnetite, ilmenite), +/− amphibole, dominated by plagioclase in the leucosome and diopside in the melanosome. The melanosome is almost completely recrystallised, with the preservation of large, relict igneous diopside phenocrysts in dyke centres. Comparison of whole-rock and mineral major- and trace-element data allowed us to assess the redistribution of elements between different mineral phases and generations during contact metamorphism and partial melting.

Dykes within and outside the thermal aureole behaved like closed chemical systems. Nevertheless, Zr, Hf, Y and REEs were internally redistributed, as deduced by comparing the trace element contents of the various diopside generations. Neocrystallised diopside – in the melanosome, leucosome and as epitaxial phenocryst rims – from the migmatite zone, are all enriched in Zr, Hf, Y and REEs compared to relict phenocrysts. This has been assigned to the liberation of trace elements on the breakdown of enriched primary minerals, kaersutite and sphene, on entering the thermal aureole. Major and trace element compositions of minerals in migmatite melanosomes and leucosomes are almost identical, pointing to a syn- or post-solidus reequilibration on the cooling of the migmatite terrain i.e. mineral–melt equilibria were reset to mineral–mineral equilibria.     

Archean high-Mg monzodiorite–syenite,epidote skarn,and biotite–sericite gold lodes in the Granny Smith–Wallaby district,Australia: U–Pb and Re–Os chronometry of two intrusion-related hydrothermal systems

The Granny Smith (37 t Au production) and Wallaby deposits (38 t out of a 180 t Au resource) are located northeast of Kalgoorlie, in 2.7 Ga greenstones of the Eastern Goldfields Province, the youngest orogenic belt of the Yilgarn craton, Western Australia. At Granny Smith, a zoned monzodiorite–granodiorite stock, dated by a concordant titanite–zircon U–Pb age of 2,665 ± 3 Ma, cuts across east-dipping thrust faults. The stock is fractured but not displaced and sets a minimum age for large-scale (1 km) thrust faulting (D2), regional folding (D1), and dynamothermal metamorphism in the mining district. The local gold–pyrite mineralization, controlled by fractured fault zones, is younger than 2,665 ± 3 Ma. In augite–hornblende monzodiorite, alteration progressed from a hematite-stained alkali feldspar–quartz–calcite assemblage and quartz–molybdenite–pyrite veins to a late reduced sericite–dolomite–albite assemblage. Gold-related monazite and xenotime define a U–Pb age of 2,660 ± 5 Ma, and molybdenite from veins a Re–Os isochron age of 2,661 ± 6 Ma, indicating that mineralization took place shortly after the emplacement of the main stock, perhaps coincident with the intrusion of late alkali granite dikes. At Wallaby, a NE-trending swarm of porphyry dikes comprising augite monzonite, monzodiorite, and minor kersantite intrudes folded and thrust-faulted molasse. The conglomerate and the dikes are overprinted by barren (<0.01 g/t Au) anhydrite-bearing epidote–actinolite–calcite skarn, forming a 600-m-wide and >1,600-m-long replacement pipe, which is intruded by a younger ring dike of syenite porphyry pervasively altered to muscovite + calcite + pyrite. Skarn and syenite are cut by pink biotite–calcite veins, containing magnetite + pyrite and subeconomic gold–silver mineralization (Au/Ag = 0.2). The veins are associated with red biotite–sericite–calcite–albite alteration in adjacent monzonite dikes. Structural relations and the concordant titanite U–Pb age of the skarn constrain intrusion-related mineralization to 2,662 ± 3 Ma. The main-stage gold–pyrite ore (Au/Ag >10) forms hematite-stained sericite–dolomite–albite lodes in stacked D2 reverse faults, which offset skarn, syenite, and the biotite–calcite veins by up to 25 m. The molybdenite Re–Os age (2,661 ± 10 Ma) of the ore suggests a genetic link to intrusive activity but is in apparent conflict with a monazite–xenotime U–Pb age (2,651 ± 6 Ma), which differs from that of the skarn at the 95% confidence level. The time relationships at both gold deposits are inconsistent with orogenic models invoking a principal role for metamorphic fluids released during the main phase of compression in the fold belt. Instead, mineralization is related in space and time to late-orogenic, magnetite-series, high-Mg monzodiorite–syenite intrusions of mantle origin, characterized by Mg/(Mg + Fe_TOTAL) = 0.31–0.57, high Cr (34–96 ppm), Ni (22–63 ppm), Ba (1,056–2,321 ppm), Sr (1,268–2,457 ppm), Th (15–36 ppm), and rare earth elements (total REE: 343–523 ppm). At Wallaby, shared Ca–K–CO₂ metasomatism and Th-REE enrichment (in allanite) link Au–Ag mineralization in biotite–calcite veins to the formation of the giant epidote skarn, implicating a Th + REE-rich syenite pluton at depth as the source of the oxidized hydrothermal fluid. At Granny Smith, lead isotope data and the Rb–Th–U signature of early hematite-bearing wall-rock alteration point to fluid released by the source pluton of the differentiated alkali granite dikes.     

Glauconite from Paleogene volcano-terrigenous rocks in Western Kamchatka

It is shown that glauconite-bearing interbeds are widespread in the layer-by-layer studied sections on the Sea of Okhotsk coast (Mainach section) and Kheisliveem River valley (Kavran section), the volcanoterrigenous rocks of the Kovachin, Amanin, and Gakkhin formations of the Paleogene in western Kamchatka (Upper Eocene-Lower Oligocene boundary beds). Detailed mineralogical and structural-crystallochemical characteristics of glauconite from the Amanin Formation are presented. It is suggested that such glauconite should not be used for geochronological purposes.Some specific features of glauconite formation, particularly, the preservation of specific morphological forms at high accumulation rates of volcano-terrigenous rocks, are discussed. Possibility of the formation of glauconite with the active influence of bacterial metabolism is considered.