Resurfacing Saturn's satellites: Models of partial differentiation and expansion

The faulting seen on the surfaces of Saturn's icy moons may have been caused either by external events, such as large impacts, or internal stresses caused by the expansion of the moons as long-lived radionuclides produced internal heating and phase changes. We estimate the stress as a function of radius expansion is σ = 44 (Δr/r) kbar. The extensional stress needed for fracture is probably something less than 40 bar so extensional fracture is likely to occur when Δr/r is greater than one part in a thousand. The radius change for these moons can be calculated analytically, given suitable assumptions; in addition, detailed time-dependent computer models of the thermal and physical evolution of Tethys, Dione, Rhea, and Iapetus were carried out. From these calculations we conclude that the most reasonable cause for rifting on Dione and Rhea is the refreezing of an ammonia-water eutectic melt inside these moons roughly two billion years after their formation, while the rift on Tethys was caused by a large impact, and little rifting should be expected on Iapetus.     

Type C Ca, Al-rich inclusions from Allende: Evidence for multistage formation

The coarse-grained, igneous, anorthite-rich (Type C) CAIs from Allende studied (100, 160, 6-1-72, 3529-40, CG5, ABC, TS26, and 93) have diverse textures and mineralogies, suggesting complex nebular and asteroidal formation histories. CAIs 100, 160, 6-1-72, and 3529-40 consist of Al,Ti-diopside (fassaite; 13-23 wt% Al₂O₃, 2-14 wt% TiO₂), Na-bearing åkermanitic melilite (0.1-0.4 wt% Na₂O; Åk_30-75), spinel, and fine-grained (∼5-10 μm) anorthite groundmass. Most of the fassaite and melilite grains have “lacy” textures characterized by the presence of abundant rounded and prismatic inclusions of anorthite ∼5-10 μm in size. Lacy melilite is pseudomorphed to varying degrees by grossular, monticellite, and pure forsterite or wollastonite. CAI 6-1-72 contains a relict Type B CAI-like portion composed of polycrystalline gehlenitic melilite (Åk_10-40), fassaite, spinel, perovskite, and platinum-group element nuggets; the Type B-like material is overgrown by lacy melilite and fassaite. Some melilite and fassaite grains in CAIs 100 and 160 are texturally similar to those in the Type B portion of 6-1-72. CAIs ABC and TS26 contain relict chondrule fragments composed of forsteritic olivine and low-Ca pyroxene; CAI 93 is overgrown by a coarse-grained igneous rim of pigeonite, augite, and anorthitic plagioclase. These three CAIs contain very sodium-rich åkermanitic melilite (0.4-0.6 wt% Na₂O; Åk_63-74) and Cr-bearing Al,Ti-diopside (up to 1.6 wt% Cr₂O₃, 1-23 wt% Al₂O, 0.5-7 wt% TiO₂). Melilite and anorthite in the Allende Type C CAI peripheries are replaced by nepheline and sodalite, which are crosscut by andradite-bearing veins; spinel is enriched in FeO. The CAI fragment CG5 is texturally and mineralogically distinct from other Allende Type Cs. It is anorthite-poor and very rich in spinel poikilitically enclosed by Na-free gehlenitic melilite (Åk_20-30), fassaite, and anorthite; neither melilite nor pyroxene have lacy textures; secondary minerals are absent. The Al-rich chondrules 3655b-2 and 3510-7 contain aluminum-rich and ferromagnesian portions. The Al-rich portions consist of anorthitic plagioclase, Al-rich low-Ca pyroxene, and Cr-bearing spinel; the ferromagnesium portions consist of fosteritic olivine, low-Ca pyroxene, and opaque nodules.We conclude that Type C CAIs 100, 160, 6-1-72, and 3529-40 formed by melting of coarse-grained Type B-like CAIs which experienced either extensive replacement of melilite and spinel mainly by anorthite and diopside (traces of secondary Na-bearing minerals, e.g., nepheline or sodalite, might have formed as well), or addition of silica and sodium during the melting event. CG5 could have formed by melting of fine-grained spinel-melilite CAI with melilite and spinel partially replaced anorthite and diopside. CAIs ABC, 93, and TS-26 experienced melting in the chondrule-forming regions with addition of chondrule-like material, such as forsteritic olivine, low-Ca pyroxene, and high-Ca pyroxene. Anorthite-rich chondrules formed by melting of the Al-rich (Type C CAI-like) precursors mixed with ferromagnesian, Type I chondrule-like precursors. The Allende Type C CAIs and Al-rich chondrules experienced fluid-assisted thermal metamorphism, which resulted in pseudomorphic replacement of melilite and anorthite by grossular, monticellite, and forsterite (100, 160, 6-1-72, 3592-40) or by grossular, monticellite, and wollastonite (ABC, 93, TS-26). The pseudomorphic replacement was followed or accompanied by iron-alkali metasomatic alteration resulting in replacement of melilite and anorthite by nepheline and sodalite, enrichment of spinel in FeO, and precipitation of salite-hedenbergite pyroxenes, wollastonite, and andradite in fractures and pores in and around CAIs.     

The effect of land surface changes on Eemian climate

Transient experiments for the Eemian (128–113 ky BP) were performed with a complex, coupled earth system model, including atmosphere, ocean, terrestrial biosphere and marine biogeochemistry. In order to investigate the effect of land surface parameters (background albedo, vegetation and tree fraction and roughness length) on the simulated changes during the Eemian, simulations with interactive coupling between climate and vegetation were compared with additional experiments in which these feedbacks were suppressed. The experiments show that the influence of land surface on climate is mainly caused by changes in the albedo. For the northern hemisphere high latitudes, land surface albedo is changed partially due to the direct albedo effect of the conversion of grasses into forest, but the indirect effect of forests on snow albedo appears to be the major factor influencing the total absorption of solar radiation. The Western Sahara region experiences large changes in land surface albedo due to the appearance of vegetation between 128 and 120 ky BP. These local land surface albedo changes can be as much as 20%, thereby affecting the local as well as the global energy balance. On a global scale, latent heat loss over land increases more than 10% for 126 ky BP compared to present-day.     

Time constant of hydraulic-head response in aquifers subjected to sudden recharge change: application to large basins

Analytical formulae are proposed to describe the first-order temporal evolution of the head in large groundwater systems (such as those found in North Africa or eastern Australia) that are subjected to drastic modifications of their recharge conditions (such as those in Pleistocene and Holocene times). The mathematical model is based on the hydrodynamics of a mixed-aquifer system composed of a confined aquifer connected to an unconfined one with a large storage capacity. The transient behaviour of the head following a sudden change of recharge conditions is computed with Laplace transforms for linear one-dimensional and cylindrical geometries. This transient evolution closely follows an exponential trend exp(?t/τ). The time constant τ is expressed analytically as a function of the various parameters characterizing the system. In many commonly occurring situations, τ depends on only four parameters: the width a _c of the main confined aquifer, its transmissivity T _c, the integrated storage situated upstream in the unconfined aquifer M?=?S _u a _u, and a curvature parameter accounting for convergence/divergence effects. This model is applied to the natural decay of large aquifer basins of the Sahara and Australia following the end of the mid-Holocene humid period. The observed persistence of the resource is discussed on the basis of the time constant estimated with the system parameters. This comparison confirms the role of the upstream water reserve, which is modelled as an unconfined aquifer, and highlights the significant increase of the time constant in case of converging flow.     

Numerical investigation of an oceanic resonant regime induced by hurricane winds

The oceanic mixed layer (OML) response to an idealized hurricane with different propagation speeds is investigated using a two-layer reduced gravity ocean model. First, the model performances are examined with respect to available observations relative to Hurricane Frances (2004). Then, 11 idealized simulations are performed with a Holland (Mon Weather Rev 108(8):1212–1218, 1980) symmetric wind profile as surface forcing with storm propagation speeds ranging from 2 to 12 m s⁻¹. By varying this parameter, the phasing between atmospheric and oceanic scales is modified. Consequently, it leads to different momentum exchanges between the hurricane and the OML and to various oceanic responses. The present study determines how OML momentum and heat budgets depend on this parameter. The kinetic energy flux due to surface wind stress is found to strongly depend on the propagation speed and on the cross-track distance from the hurricane center. A resonant regime between surface winds and near-inertial currents is clearly identified. This regime maximizes locally the energy flux into the OML. For fast-moving hurricanes (>6 m s⁻¹), the ratio of kinetic energy converted into turbulence depends only on the wind stress energy input. For slow-moving hurricanes (<6 m s⁻¹), the upwelling induced by current divergence enhances this conversion by shallowing the OML depth. Regarding the thermodynamic response, two regimes are identified with respect to the propagation speed. For slow-moving hurricanes, the upwelling combined with a sharp temperature gradient at the OML base formed in the leading part of the storm maximizes the oceanic heat loss. For fast propagation speeds, the resonance mechanism sets up the cold wake on the right side of the hurricane track. These results suggest that the propagation speed is a parameter as important as the surface wind speed to accurately describe the oceanic response to a moving hurricane.     

Comparing object-based and pixel-based classifications for mapping savannas

The development of robust object-based classification methods suitable for medium to high resolution satellite imagery provides a valid alternative to ‘traditional’ pixel-based methods. This paper compares the results of an object-based classification to a supervised per-pixel classification for mapping land cover in the tropical north of the Northern Territory of Australia. The object-based approach involved segmentation of image data into objects at multiple scale levels. Objects were assigned classes using training objects and the Nearest Neighbour supervised and fuzzy classification algorithm. The supervised pixel-based classification involved the selection of training areas and a classification using the maximum likelihood classifier algorithm. Site-specific accuracy assessment using confusion matrices of both classifications were undertaken based on 256 reference sites. A comparison of the results shows a statistically significant higher overall accuracy of the object-based classification over the pixel-based classification. The incorporation of a digital elevation model (DEM) layer and associated class rules into the object-based classification produced slightly higher accuracies overall and for certain classes; however this was not statistically significant over the object-based using spectral information solely. The results indicate object-based analysis has good potential for extracting land cover information from satellite imagery captured over spatially heterogeneous land covers of tropical Australia.     

Effect of phosphatation and calcination on the environmental behaviour of sediments

Dredging operations produce considerable quantities of materials, to be managed and this opens an opportunity for valorization in civil engineering. However, the contamination of the dredged sediments has become a major problem to solve. The major contaminants are heavy metals and organic compounds. This study focuses on the use of phosphoric acid (H3PO4) to stabilize heavy metals from sediments and destroy organic matter by calcination at 650 °C with a goal of using sediments in roadworks. Several studies have been conducted in this eld. The stabilized materials obtained have been used in civil engineering. The main purpose of this work is to discuss the environmental behavior of marine sediment treated by phosphatation and calcination. Two types of phosphoric acids were used. The pH dependence leaching test has been used as the basic characterization to evaluate the effect of the type of phosphoric acid on the metals behavior in a valorization scenario. The standard leaching test and the Toxicity Characteristic Leaching Procedure (TCLP) were conducted as compliance tests. In regards of the obtained results, the environmental assessment has also shown a reduction in the availability of targeted heavy metals in alkaline environment whatever the type of acid used for treatment. This opens opportunities for co-valorization.     

Acoustic directional snapshot wavefield decomposition

Up–down wavefield decomposition is effectuated by a scaled addition or subtraction of the pressure and vertical particle velocity, generally on horizontal or vertical surfaces, and works well for data given on such surfaces. The method, however, is not applicable to decomposing a wavefield when it is given at one instance in time, i.e. on snapshots. Such situations occur when a wavefield is modelled with methods like finite-difference techniques, for the purpose of, for example, reverse time migration, where the entire wavefield is determined per time instance. We present an alternative decomposition method that is exact when working on snapshots of an acoustic wavefield in a homogeneous medium, but can easily be approximated to heterogeneous media, and allows the wavefield to be decomposed in arbitrary directions. Such a directional snapshot wavefield decomposition is achieved by recasting the acoustic system in terms of the time derivative of the pressure and the vertical particle velocity, as opposed to the vertical derivative in up–down decomposition for data given on a horizontal surface. As in up–down decomposition of data given at a horizontal surface, the system can be eigenvalue decomposed and the inverse of the eigenvector matrix decomposes the wavefield snapshot into fields of opposite directions, including up–down decomposition. As the vertical particle velocity can be rotated at will, this allows for decomposition of the wavefield into any spatial direction; even spatially varying directions are possible. We show the power and effectiveness of the method by synthetic examples and models of increasing complexity.