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.     

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.     

New data on geology of the Shatak Complex (western slope of the Southern Urals)

The comprehensive study of sections of the Shatak Complex has revealed that conglomerates at the base of Middle Riphean rocks are not basal but intraformational rocks. Previously described angular unconformities between shales of the Sukhin Subformation (Yusha Formation, R₁) and conglomerates of the Kuz”elga Subformation (Mashak Formation, R₂) are related to late tectonic movements. Magmatic rocks developed at the base of the Middle Riphean section are represented by sheet intrusions formed in the course of emplacement of a fluid-saturated magmatic melt into partially or completely lithified terrigenous rocks at the graben formation stage during the origination of synkinematic faults that served as magma conduits. It is inferred that distribution of provenances of clastic materials and sedimentation basins in the Burzyanian and Yurmatian should be scrutinized in the study region, because the normal regressive sequence of rocks from the uppermost Yusha Formation to the lowermost Mashak Formation, which was established in the Shatak Ridge, eliminates a clear boundary distinguished between them at present. The idea about an older age of the Mashak conglomerates is substantiated.     

Deciphering polymetamorphic episodes in high-grade metamorphic orogens: Constraints from PbSL, Sm/Nd and Lu/Hf garnet dating of low- to high-grade metasedimentary rocks from the Kaoko belt (Namibia)

Three dating techniques for metamorphic minerals using the Sm–Nd, Lu–Hf and Pb isotope systems are combined and interpreted in context with detailed petrologic data from crustal segments in NW Namibia. The combination of isochron ages using these different approaches is a valuable tool to testify for the validity of metamorphic mineral dating. Here, PbSL, Lu–Hf and Sm–Nd garnet ages obtained on low- to medium-grade metasedimentary rocks from the Central Kaoko Zone of the Neoproterozoic Kaoko belt (NW Namibia) indicate that these samples were metamorphosed at around 550–560 Ma. On the other hand, granulite facies metasedimentary rocks from the Western Kaoko Zone underwent two phases of high-grade metamorphism, one at ca. 660–625 Ma and another at ca. 550 Ma providing substantial evidence that the 660–625 Ma-event was indeed a major tectonothermal episode in the Kaoko belt. Our age data suggest that interpreting metamorphic ages by applying a single dating method only is not reliable enough when studying complex metamorphic systems. However, a combination of all three dating techniques used here provides a reliable basis for geochronological age interpretation.     

Partial crustal melting beneath the Betic Cordillera (SE Spain): The case study of Mar Menor volcanic suite

The Neogene Volcanic Province (NVP) within the Betic Cordillera (SE Spain) consists of three main metapelitic enclave suites (from SW to NE: El Hoyazo, Mazarrón and Mar Menor). Since the NVP represents a singular place in the world where crustal enclaves were immediately quenched after melting, their microstructures provide a “photograph” of the conditions at depth just after the moment of the melting.

The thermobarometric information provided by the different microstructural assemblages has been integrated with the geophysical and geodynamical published data into a model of the petrologic evolution of the Mar Menor enclaves. They were equilibrated at 2–3 kbar, 850–900 °C, and followed a sequence of heating melt producing reactions. A local cooling event evidenced by minor melt crystallization preceded the eruption.

The lower crustal studies presented in this work contribute to the knowledge of: (i) the partial melting event beneath the Mar Menor volcanic suite through a petrologic detailed study of the enclaves; (ii) how the microstructures of fast cooled anatectic rocks play an important role in tracing the magma evolution in a chamber up to the eruption, and how they can be used as pseudothermobarometers; (iii) the past and current evolution of the Alborán Domain (Betic Cordillera) and Mediterranean Sea, and how the base of a metapelitic crust has melted within an active geodynamic setting.