Implicit integration under mixed controls of a breakage model for unsaturated crushable soils

This paper discusses a series of stress point algorithms for a breakage model for unsaturated granular soils. Such model is characterized by highly nonlinear coupling terms introduced by breakage‐dependent hydro‐mechanical energy potentials. To integrate accurately and efficiently its constitutive equations, specific algorithms have been formulated using a backward Euler scheme. In particular, because implementation and verification of unsaturated soil models often require the use of mixed controls, the incorporation of various hydro‐mechanical conditions has been tackled. First, it is shown that the degree of saturation can be replaced with suction in the constitutive equations through a partial Legendre transformation of the energy potentials, thus changing the thermomechanical state variables and enabling a straightforward implementation of a different control mode. Then, to accommodate more complex control scenarios without redefining the energy potentials, a hybrid strategy has been used, combining the return mapping scheme with linearized constraints. It is shown that this linearization strategy guarantees similar levels of accuracy compared with a conventional strain–suction‐controlled implicit integration. In addition, it is shown that the use of linearized constraints offers the possibility to use the same framework to integrate a variety of control conditions (e.g., net stress and/or water‐content control). The convergence profiles indicate that both schemes preserve the advantages of implicit integration, that is, asymptotic quadratic convergence and unconditional stability. Finally, the performance of the two implicit schemes has been compared with that of an explicit algorithm with automatic sub‐stepping and error control, showing that for the selected breakage model, implicit integration leads to a significant reduction of the computational cost. Such features support the use of the proposed hybrid scheme also in other modeling contexts, especially when strongly nonlinear models have to be implemented and/or validated by using non‐standard hydro‐mechanical control conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

POLARIX: a pathfinder mission of X-ray polarimetry

Since the birth of X-ray astronomy, spectral, spatial and timing observation improved dramatically, procuring a wealth of information on the majority of the classes of the celestial sources. Polarimetry, instead, remained basically unprobed. X-ray polarimetry promises to provide additional information procuring two new observable quantities, the degree and the angle of polarization. Polarization from celestial X-ray sources may derive from emission mechanisms themselves such as cyclotron, synchrotron and non-thermal bremsstrahlung, from scattering in aspheric accreting plasmas, such as disks, blobs and columns and from the presence of extreme magnetic field by means of vacuum polarization and birefringence. Matter in strong gravity fields and Quantum Gravity effects can be studied by X-ray polarimetry, too. POLARIX is a mission dedicated to X-ray polarimetry. It exploits the polarimetric response of a Gas Pixel Detector, combined with position sensitivity, that, at the focus of a telescope, results in a huge increase of sensitivity. The heart of the detector is an Application-Specific Integrated Circuit (ASIC) chip with 105,600 pixels each one containing a full complete electronic chain to image the track produced by the photoelectron. Three Gas Pixel Detectors are coupled with three X-ray optics which are the heritage of JET-X mission. A filter wheel hosting calibration sources unpolarized and polarized is dedicated to each detector for periodic on-ground and in-flight calibration. POLARIX will measure time resolved X-ray polarization with an angular resolution of about 20 arcsec in a field of view of 15 × 15 arcmin and with an energy resolution of 20% at 6 keV. The Minimum Detectable Polarization is 12% for a source having a flux of 1 mCrab and 10⁵ s of observing time. The satellite will be placed in an equatorial orbit of 505 km of altitude by a Vega launcher. The telemetry down-link station will be Malindi. The pointing of POLARIX satellite will be gyroless and it will perform a double pointing during the earth occultation of one source, so maximizing the scientific return. POLARIX data are for 75% open to the community while 25% + SVP (Science Verification Phase, 1 month of operation) is dedicated to a core program activity open to the contribution of associated scientists. The planned duration of the mission is one year plus three months of commissioning and SVP, suitable to perform most of the basic science within the reach of this instrument. A nice to have idea is to use the same existing mandrels to build two additional telescopes of iridium with carbon coating plus two more detectors. The effective area in this case would be almost doubled.     

Wave equation datuming applied to marine OBS data and to land high resolution seismic profiling

One key step in seismic data processing flows is the computation of static corrections, which relocate shots and receivers at the same datum plane and remove near surface weathering effects. We applied a standard static correction and a wave equation datuming and compared the obtained results in two case studies: 1) a sparse ocean bottom seismometers dataset for deep crustal prospecting; 2) a high resolution land reflection dataset for hydrogeological investigation. In both cases, a detailed velocity field, obtained by tomographic inversion of the first breaks, was adopted to relocate shots and receivers to the datum plane. The results emphasize the importance of wave equation datuming to properly handle complex near surface conditions. In the first dataset, the deployed ocean bottom seismometers were relocated to the sea level (shot positions) and a standard processing sequence was subsequently applied to the output. In the second dataset, the application of wave equation datuming allowed us to remove the coherent noise, such as ground roll, and to improve the image quality with respect to the application of static correction. The comparison of the two approaches evidences that the main reflecting markers are better resolved when the wave equation datuming procedure is adopted.     

Granulite-amphibolite facies metasediments from the Serre (Calabria,Southern Italy): their protoliths and the processes controlling their chemistry

In the northwestern Serre, a near continuous section through the Palaeozoic lower continental crust crops out. This paper deals with granulite facies metasedimentary rocks of this section. Fifty-five chemical analyses (major elements, Sc, V, Cr, Co, Ni, Rb, Sr, Y, Nb and Ba) are given. Chemical data indicate that the original sedimentary sequence was made up of arenites in the lower part and mainly Fe and Mg rich pelites in the upper part. The metapelites experienced anatexis and removal of some melt, whereas the meta-arenites suffered only incipient anatexis and became dehydration granulites. The various rock types display an overall relatively low K/Rb ratio due mainly to the stability of K- and Rb-retaining phases. A weighted mean composition of the fragment of the lower continental crust outcropping in the Serre has been established by integrating the chemical data for metasediments (55 analyses), metabasites (92 analyses) and dioritic and tonalitic gneisses (10 analyses) with the field data. This composition compared to the compositions of the bulk continental crust and of other sections considered as belonging to the lower crust, confirms the hypothesis that (1) the lower crust is heterogeneous, and (2) it is more basic than the continental crust.     

Some chemical characteristics of the brines in Bannock and Tyro Basins: salinity, sulphur compounds, Ca, F, pH, At, PO4, SiO2, NH3

Results of the chemical investigation on the Bannock and Tyro Basins are reported.Both basins were found to be hypersaline ( 10 times higher than salinity of normal seawater) and anoxic. In all investigated basins a region of transition, a few meters thick, was identified at depths > 3327 dbar. It is characterized by a sharp gradient of salinity, and all concentrations of analysed species increase strongly except for dissolved oxygen and nitrate, which immediately drop to zero. This region appears as a sharp boundary that prevents mixing. As a result, in the presence of organic matter, an anoxic condition developed with the complete depletion of dissolved oxygen. At the same time, hydrogen sulphide and ammonium accumulated within the brine. Between the Bannock and the Tyro brines differences occur in the measured concentrations of H₂S, SO²⁻₄, Ca²⁺ and NH₃. There are some differences also within the Bannock area sub-basins.The Libeccio sub-basin, in the Bannock area, contains a double-layered brine: the upper layer is 140 dbar thick and the lower layer is 300 dbar thick. A second interface between upper and lower brines develops at a depth of 3500 dbar. Nearly all of the measured concentrations vs. depth show the double layer, with the exception of ammonium, the concentration of which remains nearly constant throughout the anoxic column. Profiles of the other species analyzed show remarkable differences on passing from the upper to the lower brine. Hydrogen sulphide, sulphate and fluoride concentrations appear constant and then increase at the second interface. The calcium concentration is also constant in the upper brine, but decreases at the second interface. Total alkalinity and phosphate concentrations show a maximum peak just below the first interface. However, after passing through the second interface all the chemical parameters exhibit an almost constant behaviour down to the bottom.Hypersaline conditions are attributed to the dissolution of Messinian evaporite, and anoxia is suggested to originate from the oxidation of organic matter present in sediments and from the absence of bottom water circulation in such a deep and enclosed environment.The chemical conditions can be summarized as follows: in the Libeccio Basin the values for the species analysed have the ranges: 39–321 psu for ‘salinity’, 8.2−6.5 for pH, 2.7–4.0 mM for total alkalinity, 0.2-0 mM for dissolved oxygen, 0–1669 μM for hydrogen sulphide, 0–198 μM for thiol, 31–99 mM for sulphate, 11–21 mM for calcium, 7–100 μM for fluoride, 0.2–3080 μM for ammonium, 5.8-0 μM for nitrite, 0.2–12 μM for phosphate and 8–130 μM for silicate.     

Lava flow hazard at Nyiragongo volcano,D.R.C.

The 2002 eruption of Nyiragongo volcano constitutes the most outstanding case ever of lava flow in a big town. It also represents one of the very rare cases of direct casualties from lava flows, which had high velocities of up to tens of kilometer per hour. As in the 1977 eruption, which is the only other eccentric eruption of the volcano in more than 100 years, lava flows were emitted from several vents along a N–S system of fractures extending for more than 10 km, from which they propagated mostly towards Lake Kivu and Goma, a town of about 500,000 inhabitants. We assessed the lava flow hazard on the entire volcano and in the towns of Goma (D.R.C.) and Gisenyi (Rwanda) through numerical simulations of probable lava flow paths. Lava flow paths are computed based on the steepest descent principle, modified by stochastically perturbing the topography to take into account the capability of lava flows to override topographic obstacles, fill topographic depressions, and spread over the topography. Code calibration and the definition of the expected lava flow length and vent opening probability distributions were done based on the 1977 and 2002 eruptions. The final lava flow hazard map shows that the eastern sector of Goma devastated in 2002 represents the area of highest hazard on the flanks of the volcano. The second highest hazard sector in Goma is the area of propagation of the western lava flow in 2002. The town of Gisenyi is subject to moderate to high hazard due to its proximity to the alignment of fractures active in 1977 and 2002. In a companion paper (Chirico et al., Bull Volcanol, in this issue, 2008) we use numerical simulations to investigate the possibility of reducing lava flow hazard through the construction of protective barriers, and formulate a proposal for the future development of the town of Goma.     

The first absolute gravity measurements in Indonesia

For the purposes of the calibration of the superconducting gravimeter (SG) in Bandung and the establishment of the absolute gravity (AG) points, we carried out AG measurements for the first time in Indonesia in November 2002. The measurements in Bandung were conducted between November 15th and 20th by means of a FG5 (#210), and 14,520 effective drops were obtained. The gravity value newly determined at the AG point in Bandung is 977976701.2 μgal (1 μgal = 10⁻⁸ ms⁻²) and the scale factor for the SG is −52.22 μgal/V. We also established another AG point in Yogyakarta near Merapi volcano and carried out AG measurements in Yogyakarta between November 22nd and 26th. The gravity value determined for this station is 978203093.5 μgal.