Numerical Homogenization of Jointed Rock Masses Using Wave Propagation Simulation

Homogenization in fractured rock analyses is essentially based on the calculation of equivalent elastic parameters. In this paper, a new numerical homogenization method that was programmed by means of a MATLAB code, called HLA-Dissim, is presented. The developed approach simulates a discontinuity network of real rock masses based on the International Society of Rock Mechanics (ISRM) scanline field mapping methodology. Then, it evaluates a series of classic joint parameters to characterize density (RQD, specific length of discontinuities). A pulse wave, characterized by its amplitude, central frequency, and duration, is propagated from a source point to a receiver point of the simulated jointed rock mass using a complex recursive method for evaluating the transmission and reflection coefficient for each simulated discontinuity. The seismic parameters, such as delay, velocity, and attenuation, are then calculated. Finally, the equivalent medium model parameters of the rock mass are computed numerically while taking into account the natural discontinuity distribution. This methodology was applied to 17 bench fronts from six aggregate quarries located in Tunisia, Spain, Austria, and Sweden. It allowed characterizing the rock mass discontinuity network, the resulting seismic performance, and the equivalent medium stiffness. The relationship between the equivalent Young’s modulus and rock discontinuity parameters was also analyzed. For these different bench fronts, the proposed numerical approach was also compared to several empirical formulas, based on RQD and fracture density values, published in previous research studies, showing its usefulness and efficiency in estimating rapidly the Young’s modulus of equivalent medium for wave propagation analysis.     

Along-strike variations of P–T conditions in accretionary wedges and syn-orogenic extension,the HP–LT Phyllite–Quartzite Nappe in Crete and the Peloponnese

Syn-orogenic detachments in accretionary wedges make the exhumation of high-pressure and low-temperature metamorphic rocks possible with little erosion. The velocity of exhumation within the subduction channel or the accretionary complex, and thus the shape of P–T paths, depend upon the kinematic boundary conditions. A component of slab retreat tends to open the channel and facilitates the exhumation. We document the effect of slab retreat on the shape of P–T paths using the example of the Phyllite–Quartzite Nappe that has been exhumed below the Cretan syn-orogenic detachment during the Miocene in Crete and the Peloponnese. Data show a clear tendency toward colder conditions at peak pressure and during exhumation where the intensity of slab retreat is larger. This spatial evolution of P–T gradient is accompanied with an evolution from a partly coaxial regime below the Peloponnese section of the detachment toward a clearly non-coaxial regime in Crete.     

Duration Magnitude Scale and Site Residuals for Northern Morocco

— The first empirical duration magnitude (M_D) formula is developed and tested for the Northern Morocco Seismic Network (NMSNET). This relationship is obtained by relating the IGN (Instituto Geografico National, Madrid) body-waves mbLg^IGN to the duration (), and the epicentral distance (), at 25 analogue stations of the NMSNET for 479 earthquakes with 2.5 mb 5.4, from March 1992 to February 2001. M_D estimates are significantly more precise while introducing a correction term for each of these stations, cSta_j. The magnitude for the i^th event (M_D)_i is the mean value of individual M_Dij=–0.14+1.63log₁₀(_ij)+0.031(_ij)+cSta_j. The cSta_j corrections reduce considerably the local site effects which influence the recorded durations and cause stations to either overestimate, or underestimate M_D up to 0.5 magnitude units. Average station M_D residuals (–cSta_j) are found to be independent of the distance from the epicenter to at least 10 degrees. It seems evident that regional geological features in the immediate behavior of stations have a systematic effect on the corresponding obtained residuals: older well-consolidated Precambrian crystalline rocks produce high negative residuals (shorter durations), younger unconsolidated sediments produce high positive residuals (longer durations), whereas, intermediate M_D site residuals appear to be the result of the effect of various factors, principally age and state of consolidation of the bedrock, combined with the local tectonic.     

Time series data for Canadian arctic vertebrates: IPY contributions to science, management, and policy

Long-term data are critically important to science, management, and policy formation. Here we describe a number of data collections from arctic Canada that monitor vertebrate population trends of freshwater and marine fish, marine birds, marine and terrestrial mammals. These time series data cover the last ca. 30?years and capture a period from the onset of global changes affecting the Arctic up to recent years with a rapid increase in temperature. While many of these data collections were initiated through a variety of government and university programs, they also include a surge in polar research launched with the recent International Polar Year (2007?C2008). We estimated the long-term vertebrate index from our data that summarizes various taxa abundance trends within a global context and observed a continuous decline of about 30?% in population abundance since the 1990s. Though most data collections are biased towards few taxa, we conduct time-series analyses to show that the potential value of long-term data emerges as individual monitoring sites can be spread across space and time scales. Despite covering a handful of populations, the different time series data covered a large spectrum of dynamics, cyclic to non-cyclic, including coherence with the North Atlantic Oscillation, lag effects, and density dependence. We describe a synthesis framework to integrate ecological time-series research and thereby derive additional benefits to management, science, and policy. Future requirements include: (1) continuation of current observation systems; (2) expansion of current monitoring sites to include additional trophic links and taxonomic indicators; (3) expansion beyond the existing program to include greater spatial coverage into less-sampled ecosystems and key representative locations; and (4) integration of circumpolar observations and comprehensive analyses. Development of a circumpolar observation system is necessary for innovative science, large-scale adaptive management, and policy revision essential to respond to rapid global change.     

Cambrian–early Ordovician volcanism across the South Armorican and Occitan domains of the Variscan Belt in France: Continental break-up and rifting of the northern Gondwana margin

The Cambrian–lower Ordovician volcanic units of the South Armorican and Occitan domains are analysed in a tectonostratigraphic survey of the French Variscan Belt. The South Armorican lavas consist of continental tholeiites in middle Cambrian–Furongian sequences related to continental break-up. A significant volcanic activity occurred in the Tremadocian, dominated by crustal melted rhyolitic lavas and initial rifting tholeiites. The Occitan lavas are distributed into five volcanic phases: (1) basal Cambrian rhyolites, (2) upper lower Cambrian Mg-rich tholeiites close to N-MORBs but crustal contaminated, (3) upper lower–middle Cambrian continental tholeiites, (4) Tremadocian rhyolites, and (5) upper lower Ordovician initial rift tholeiites. A rifting event linked to asthenosphere upwelling took place in the late early Cambrian but did not evolve. It renewed in the Tremadocian with abundant crustal melting due to underplating of mixed asthenospheric and lithospheric magmas. This main tectono-magmatic continental rift is termed the “Tremadocian Tectonic Belt” underlined by a chain of rhyolitic volcanoes from Occitan and South Armorican domains to Central Iberia. It evolved with the setting of syn-rift coarse siliciclastic deposits overlain by post-rift deep water shales in a suite of sedimentary basins that forecasted the South Armorican–Medio-European Ocean as a part of the Palaeotethys Ocean.     

Geochemical indicators of interbasin groundwater flow within the southern Rio Grande Valley, southwestern USA

Increased groundwater withdrawals for the growing population in the Rio Grande Valley and likely alteration of recharge to local aquifers with climate change necessitates an understanding of the groundwater connection between the Jornada del Muerto Basin and the adjoining and more heavily used aquifer in the Mesilla Basin. Separating the Jornada and Mesilla aquifers is a buried bedrock high from Tertiary intrusions. This bedrock high or divide restricts and/or retards interbasin flow from the Jornada aquifer into the Mesilla aquifer. The potentiometric surface of the southern Jornada aquifer near part of the bedrock high indicates a flow direction away from the divide because of a previously identified damming effect, but a groundwater outlet from the southern Jornada aquifer is necessary to balance inputs from the overall Jornada aquifer. Differences in geochemical constituents (major ions, δD, δ¹⁸O, δ³⁴S, and ⁸⁷Sr/⁸⁶Sr) indicate a deeper connection between the two aquifers through the Tertiary intrusions where Jornada water is geochemically altered because of a geothermal influence. Jornada groundwater likely is migrating through the bedrock high in deeper pathways formed by faults of the Jornada Fault Zone, in addition to Jornada water that overtops the bedrock high as previously identified as the only connection between the two aquifers. Increased groundwater withdrawals and lowering of the potentiometric surface of the Jornada aquifer may alter this contribution ratio with less overtopping of the bedrock high and a continued deeper flowpath contribution that could potentially increase salinity values in the Mesilla Basin near the divide.     

The sources of geogenic arsenic in aquifers at Datong basin,northern China: Constraints from isotopic and geochemical data

Mineralogical, geochemical and zircon U–Pb dating studies were carried out to identify the sources of arsenic in the shallow aquifers of Datong Basin in northern China. A sediment sample from 18 m depth containing 10.3 mg/kg arsenic showed a Zircon U–Pb concordant age of 2528 ± 20 to 271 ± 4 Ma that can be divided into two groups (2528 ± 20 to 1628 ± 21 Ma and 327 ± 4 to 271 ± 4 Ma) and is comparable to that of the sedimentary rocks of Taiyuan (upper Carboniferous) and Shanxi Formation (lower Permian) outcropping to the west of Datong Basin. In contrast, a sediment sample from 22.5 m depth containing 5.7 mg/kg arsenic displayed a Zircon U–Pb concordant age ranging from 2561 ± 21 to 1824 ± 26 Ma that is comparable to that of the Hengshan Complex (Ne-Archaean Precambrian) outcropping to the east of .     

Theoretical study of the local charge compensation and spectroscopic properties of B-type carbonate defects in apatite

The structure and spectroscopic properties of selected models of B-type carbonate defects in apatite locally compensated by fluoride or hydroxyl ions are investigated using first-principles quantum mechanical calculations. Theoretical infrared absorption spectra and ¹³C, and ¹⁹F nuclear magnetic resonance chemical shifts are determined. Among the investigated models, only the clumped (CO₃ ^2?, F^?) defect, with the carbonate group close to the sloping face of the tetrahedral site and the F^? ion at the remaining apex, corresponds to previous experimental observations performed on carbonate-fluorapatite samples. Although the substitution of hydroxyl by fluoride ions is commonly observed in minerals, the clumped (CO₃ ^2?, OH^?) defects are unlikely to occur in apatite, considering both their theoretical spectroscopic properties and relative stability. Anionic F^? for OH^? exchange between channel and B sites displays a preference of ~20 kJ/mol for the local charge compensation by fluoride ions at the B-site, pointing to a significantly different behavior of F^? and OH^? ions in the charge compensation mechanism. This difference is ascribed to the poor H-bond acceptor character of available oxygen atoms surrounding the apex of the tetrahedral site. The explicit calculation of the infrared absorption spectra of the defect models is also used to interpret the significant difference observed in the linewidth of the ν₂ and ν₃ CO₃ infrared powder absorption bands of carbonated apatite samples. It is shown that for a concentration of 4.4 wt% of CO₂, long-range electrostatic effects already significantly contribute to the broadening of the ν₃ CO₃ bands in apatite.     

A model for decoding the life cycle of granular avalanches in a rotating drum

Granular materials can behave as harmless sand dunes or as devastating landslides. A granular avalanche marks the transition between these distinct solid-like and fluid-like states. The solid-like state is typically described using plasticity models from critical state theory. In the fluid regime, granular flow is commonly captured using a visco-plastic model. However, due to our limited understanding of the mechanism governing the solid–fluid-like transition, characterizing the material behavior throughout the life cycle of an avalanche remains an open challenge. Here, we employ laboratory experiments of transient avalanches spontaneously generated by a rotating drum. We report measurements of dilatancy and grain kinematics before, during, and after each avalanche. Those measurements are directly incorporated into a rate-dependent plasticity model that quantitatively predicts the granular flow measured in experiments. Furthermore, we find that dilatancy in the solid-like state controls the triggering of granular avalanches and therefore plays a key role in the solid–fluid-like transition. With the proposed approach, we demonstrate that the life cycle of a laboratory avalanche, from triggering to run out, can be fully explained. Our results represent an important step toward a unified understanding of the physical phenomena associated with transitional behavior in granular media.