A basal slip model for Lagrangian finite element simulations of 3D landslides

A Lagrangian numerical approach for the simulation of rapid landslide runouts is presented and discussed. The simulation approach is based on the so‐called Particle Finite Element Method. The moving soil mass is assumed to obey a rigid‐viscoplastic, non‐dilatant Drucker–Prager constitutive law, which is cast in the form of a regularized, pressure‐sensitive Bingham model. Unlike in classical formulations of computational fluid mechanics, where no‐slip boundary conditions are assumed, basal slip boundary conditions are introduced to account for the specific nature of the landslide‐basal surface interface. The basal slip conditions are formulated in the form of modified Navier boundary conditions, with a pressure‐sensitive threshold. A special mixed Eulerian–Lagrangian formulation is used for the elements on the basal interface to accommodate the new slip conditions into the Particle Finite Element Method framework. To avoid inconsistencies in the presence of complex shapes of the basal surface, the no‐flux condition through the basal surface is relaxed using a penalty approach. The proposed model is validated by simulating both laboratory tests and a real large‐scale problem, and the critical role of the basal slip is elucidated. Copyright © 2016 John Wiley & Sons, Ltd.     

Time-dependent radiation transfer and a possible explanation of the interpulse in CP 0950

The time-dependent equation of radiative transfer is solved exactly and in then-th Gaussian approximation.The atmosphere is plane-parallel and semi-infinite; isotropic scattering is assumed, but the boundary condition at =0 is arbitrary.The results are used to investigate a suggested mechanism for the origin of the secondary pulses in CP 0950; it is found that a binary system of neutron stars can indeed explain formation, time delay and intensity of the observed interpulse.     

Fe<Superscript>3+</Superscript> reduction during biotite melting in graphitic metapelites: another origin of CO<Subscript>2</Subscript> in granulites

The Fe³⁺/Fe_tot of all Fe-bearing minerals has been analysed by Mössbauer spectroscopy in a suite of biotite-rich to biotite-free graphitic metapelite xenoliths, proxies of an amphibolite-granulite transition through progressive biotite melting. Biotite contains 9 to 16% Fe³⁺/Fe_tot, whereas garnet, cordierite and ilmenite are virtually Fe³⁺ -free (0–1% Fe³⁺/Fe_tot) in all samples, regardless of biotite presence. Under relatively reducing conditions (graphite-bearing assemblages), biotite is the only carrier of Fe³⁺ during high-temperature metamorphism; therefore, its disappearance by melting represents an important event of iron reduction during granulite formation, because haplogranitic melts usually incorporate small amounts of ferric iron. Iron reduction is accompanied by the oxidation of carbon and the production of CO₂, according to the redox reaction:

Far-field Gravity and Tilt Signals by Large Earthquakes: Real or Instrumental Effects?

A wide set of dynamics phenomena (i.e., Geodynamics, Post Glacial Rebound, seismicity and volcanic activity) can produce time gravity changes, which spectrum varies from short to long (more than 1 year) periods. The amplitude of the gravity variations is generally in the order of consequently their detection requires instruments with high sensitivity and stability: then, high quality experimental data. Spring and superconducting gravimeters are intensively used with this target and they are frequently jointed with tiltmeters recording stations in order to measure the elasto-gravitational perturbation of the Earth. The far-field effects produced by large earthquakes on records collected by spring gravimeters and tiltmeters are investigated here. Gravity and tilt records were analyzed on time windows spanning the occurrence of large worldwide earthquakes; the gravity records have been collected on two stations approximately 600 km distant. The background noise level at the stations was characterized, in each season, in order to detect a possible seasonal dependence and the presence of spectral components which could hide or mask other geophysical signals, such as, for instance, the highest mode of the Seismic Free Oscillation (SFO) of the Earth. Some spectral components (6.5; 8; 9; 14, 20, 51) have been detected in gravity and tilt records on the occasion of large earthquakes and the effect of the SFO has been hypothesized. A quite different spectral content of the EW and NS tiltmeter components has been detected and interpreted as a consequence of the radiation pattern of the disturbances due to the earthquakes. Through the analysis of the instrumental sensitivity, instrumental effects have been detected for gravity meters at very low frequency.F.S.E. (Fondo Sociale Europeo -European Community -)     

Natural and Anthropic Middle-term Variations in the Delta System at the Mouth of the Goro River (Po Delta, Italy)

Abstract. The modern delta system of the Po River is the result of an evolution strictly controlled by anthropic factors and mainly developed during a period of climatic deterioration (Little Ice Age 1550/1850).
At the beginning of the 17^th century the main branch of the Po River, Fornaci Po, was modified with the "Cut of Viro" in order to avoid landfill of the Venetian Lagoon. Simultaneously, the southernmost branch of Ariano, or Goro, was forced to flow southwards. Subsequent sedimentary variations, caused by changes in the climatic conditions, contributed to make the area supplied by the Goro Po River increasingly lagoonary. In the second half of the 19^th century (end of the Little Ice Age – L.I.A. –), the first coastal spits started to take shape and developed in different systems throughout the entire following century. Coastal spits turned out to produce in the original bay, a progressive closure towards the sea. The strong decrease in solid loads, characterising the second half of the 20^th century and mainly due to anthropic factors, marks the beginning of a new phase extending to the overall Po Delta system, in which the mouth layout is controlled by wave energy.     

Some limitations in using <Superscript>222</Superscript>Rn to assess river–groundwater interactions: the case of Castel di Sangro alluvial plain (central Italy)

²²²Rn was used to assess river–groundwater interactions within Castel di Sangro alluvial aquifer (Italy). The effectiveness of results obtained through this indicator was verified by also analyzing δ¹⁸O, major ions and temperature in both surface and groundwater, and carrying out piezometric head monitoring and discharge measurements. Hydrogeological investigations suggested that the river infiltrates into the aquifer in the south-eastern aquifer portion, while groundwater discharges into the river in the north-eastern portion. The latter phenomenon is supported by ²²²Rn data. Nevertheless, flow-through conditions cause the modelled discharge along this river reach, estimated by ²²²Rn data in a degassing-corrected two-component mixing model, to be greater than the measured discharge. Concerning river infiltration into the aquifer, δ¹⁸O, major ions and temperature data show that the river contribution is negligible in terms of aquifer recharge. Thus, the observed increase in ²²²Rn concentration in that portion of the aquifer is due to the enrichment process caused by infiltration of rainwater (²²²Rn free) which flows from the local divide area. Hence, in the study site, the use of only ²²²Rn to predict river–groundwater interactions causes some estimation inaccuracies and it must be coupled with other hydrochemical and hydrogeological parameters to gain a thorough understanding of such interactions.     

Performance of intertidal topography video monitoring of a meso-tidal reflective beach in South Portugal

This study discusses site-specific system optimization efforts related to the capability of a coastal video station to monitor intertidal topography. The system consists of two video cameras connected to a PC, and is operating at the meso-tidal, reflective Faro Beach (Algarve coast, S. Portugal). Measurements from the period February 4, 2009 to May 30, 2010 are discussed in this study. Shoreline detection was based on the processing of variance images, considering pixel intensity thresholds for feature extraction, provided by a specially trained artificial neural network (ANN). The obtained shoreline data return rate was 83%, with an average horizontal cross-shore root mean square error (RMSE) of 1.06 m. Several empirical parameterizations and ANN models were tested to estimate the elevations of shoreline contours, using wave and tidal data. Using a manually validated shoreline set, the lowest RMSE (0.18 m) for the vertical elevation was obtained using an ANN while empirical parameterizations based on the tidal elevation and wave run-up height resulted in an RMSE of 0.26 m. These errors were reduced to 0.22 m after applying 3-D data filtering and interpolation of the topographic information generated for each tidal cycle. Average beach-face slope tan(β) RMSE were around 0.02. Tests for a 5-month period of fully automated operation applying the ANN model resulted in an optimal, average, vertical elevation RMSE of 0.22 m, obtained using a one tidal cycle time window and a time-varying beach-face slope. The findings indicate that the use of an ANN in such systems has considerable potential, especially for sites where long-term field data allow efficient training.     

A Variable Mesh Spacing for Large-Eddy Simulation Models in the Convective Boundary Layer

A variable vertical mesh spacing for large-eddy simulation (LES) models in a convective boundary layer (CBL) is proposed. The argument is based on the fact that in the vertical direction the turbulence near the surface in a CBL is inhomogeneous and therefore the subfilter-scale effects depend on the relative location between the spectral peak of the vertical velocity and the filter cut-off wavelength. From the physical point of view, this lack of homogeneity makes the vertical mesh spacing the principal length scale and, as a consequence, the LES filter cut-off wavenumber is expressed in terms of this characteristic length scale. Assuming that the inertial subrange initial frequency is equal to the LES filter cut-off frequency and employing fitting expressions that describe the observed convective turbulent energy one-dimensional spectra, it is feasible to derive a relation to calculate the variable vertical mesh spacing. The incorporation of this variable vertical grid within a LES model shows that both the mean quantities (and their gradients) and the turbulent statistics quantities are well described near to the ground level, where the LES predictions are known to be a challenging task.     

African, southern Indian and South American cratons were not part of the Rodinia supercontinent: evidence from field relationships and geochronology

We discuss the question whether the late Mesoproterozoic and early Neoproterozoic rocks of eastern, central and southern Africa, Madagascar, southern India, Sri Lanka and South America have played any role in the formation and dispersal of the supercontinent Rodinia, believed to have existed between about 1000 and 750 Ma ago. First, there is little evidence for the production of significant volumes of ˜1.4–1.0 Ga (Kibaran or Grenvillian age) continental crust in the Mozambique belt (MB) of East Africa, except, perhaps, in parts of northern Mozambique. This is also valid for most terranes related to West Gondwana, which are made up of basement rocks older than Mesoproterozoic, reworked in the Brasiliano/Pan-African orogenic cycle. This crust cannot be conclusively related to either magmatic accretion processes on the active margin of Rodinia or continental collision leading to amalgamation of the supercontinent. So far, no 1.4–1.0 Ga rocks have been identified in Madagascar. Secondly, there is no conclusive evidence for a ˜1.0 Ga high-grade metamorphic event in the MB, although such metamorphism has been recorded in the presumed continuation of the MB in East Antarctica. In South America, even the Sunsas mobile belt, which is correlated with the Grenville belt of North America, does not include high-grade metamorphic rocks. All terranes with Mesoproterozoic ages seem to have evolved within extensional, aulacogen-type structures, and their compressional deformation, where observed, is normally much younger and is related to amalgamation of Gondwana. This is also valid for the Trans-Saharan and West Congo belts of West Africa.Third, there is also no evidence for post-1000 Ma sedimentary sequences that were deposited on the passive margin(s) of Rodinia. In contrast, the MB of East Africa and Madagascar is characterized by extensive structural reworking and metamorphic overprinting of Archaean rocks, particularly in Tanzania and Madagascar, and these rocks either constitute marginal parts of cratonic domains or represent crustal blocks (terranes or microcontinents?) of unknown derivation. This is also the case for most terranes included in the Borborema/Trans-Saharan belt of northeastern Brazil and west-central Africa, as well as those of the Central Goíás Massif in central Brazil and the Mantiqueira province of eastern and southeastern Brazil.Furthermore, there is evidence for extensive granitoid magmatism in the period ˜840 to <600 Ma whose predominant calc-alkaline chemistry suggests subduction-related active margin processes during the assembly of the supercontinent Gondwana. The location of the main Neoproterozoic magmatic arcs suggests that a large oceanic domain separated the core of Rodinia, namely Laurentia plus Amazonia, Baltica and West Africa, from several continental masses and fragments now in the southern hemisphere, such as the São Francisco/Congo, Kalahari and Rio de La Plata cratons, as well as the Borborema/Trans-Saharan, Central Goiás Massif and Paraná blocks. Moreover, many extensional tectonic events detected in the southern hemisphere continental masses, but also many radiometric ages of granitois that are already associated with the process of amalgamation of Gondwana, are comprised within the 800–1000 age interval. This seems incompatible with current views on the time of disintegration of Rodinia, assumed to have occurred at around 750 Ma.     

P–T–t conditions of high-grade metamorphic rocks of the Garzon Massif, Andean basement, SE Colombia

The metamorphic evolution of the Garzón Massif, Colombia, is established on the basis of the textural, goethermobarometric, and geochronological relationships of the metamorphic minerals. The geothermobarometric data define a clockwise, nearly isothermal decompression path (ITD) for rocks from Las Margaritas migmatites, constrained by four P–T areas: 780–826 °C and 6.3–8.0 kbar, 760–820 °C and 8.0–8.8 kbar, 680–755 °C and 6.6–9.0 kbar, and 630 °C and 4 kbar. For the a garnet-bearing charnockitic gneiss from the Vergel granulites, the path is counterclockwise, constrained by geothermobarometric data of 5.3–6.2 kbar and 700–780 °C and 6.2–7.2 kbar and 685–740 °C. The clockwise ITD path represents a loop followed by the orogen during the transitional granulite–amphibolite metamorphic conditions, probably associated with a subduction process followed by a collisional tectonic event. This subduction framework produced continental crust thickening between 1148 and 1034 Ma and later collision with another continental block approximately 1000 Ma ago. The orogenic exhumation occurred with moderate uplift rate. The counterclockwise trajectory and two metamorphic events suggest a vertical displacement between the Vergel granulites and Las Margaritas migmatites units, because there is no isotopic difference that indicates the existence of different terranes. The data confirm that the metamorphic evolution for this domain was more dynamic than previously believed and includes: (1) metamorphic processes with the generation of new crust with a possible mixture of old material and (2) metamorphic recycling of continental crust. These geological processes characterize a complex Mesoproterozoic orogenic event that shares certain features with the Grenvillian basement rocks participating in the formation of Rodinia.