Distribution of gold and silver and its relation with hypogene ore minerals in the Esperanza porphyry deposit,Antofagasta Region,Chile

Modes of occurrence of Au‐ and Ag‐bearing phases and their relation with associated hypogene ore minerals were examined with the objective to elucidate Au‐Ag distribution at the Esperanza porphyry deposit in the Eocene Centinela copper belt, using ore‐microscope modal analysis, semi‐quantitative analyses by automated mineralogy, electron probe microanalysis, and secondary ion mass spectrometer. The Esperanza hypogene mineralization is characterized by early‐stage chalcopyrite‐rich veinlets in the potassic alteration zone and later polymetallic stage with tennantite and galena in the chlorite‐sericitic alteration zone. Only the early‐stage chalcopyrite contains fine‐grained electrum (Au₆₈Ag₃₂ ‐ Au₈₁Ag₁₉) and hessite (Ag₂Te), and thus yields positive correlations in Cu vs. Au and Cu vs. Ag grades that are clearly recognized in the hypogene sulfide zone. The early‐stage chalcopyrite grains frequently exhibit polysynthetic twinning suggestive of inversion from intermediate solid solution. These features suggest that the fine‐grained electrum and hessite are products exsolved in the cooling process with the intermediate solid solution to chalcopyrite inversion. In contrast, tennantite and galena of the later‐stage mineralization contain no detectable Ag, and it is thus proposed that the early‐stage inverted chalcopyrite is the principal storage of economically important precious metals.     

On the accuracy of the ground force estimated in vibroseis acquisition

For a linear elastic Earth the time derivative of the ground force is considered proportional to the far-field wavelet. Under the assumption that the baseplate is stiff and the bending forces of the baseplate are negligible, the ground force is also approximated by the sum of the accelerations of the baseplate and the reaction mass weighted by the respective masses. Combining these two assumptions, the time derivative of the weighted sum is considered proportional to the far-field wavelet. This result, often referred to as the far-field wavelet assumption, although convenient and most often employed is not always valid. We explore its validity using the spectral harmonic ratios of recorded data, which are used extensively in data filtering and analysis of vibratory data. We show that the far-field wavelet assumption fails particularly for harmonic components of even order. More compact soil after repeated shots further invalidates this assumption. Non-linear modelling of the ground under the vibrator point may provide a direction towards solving this discrepancy. Finally, we describe a method for the estimation of the harmonic spectral ratios.     

A unified framework for earthquake risk assessment of transportation networks and gross regional product

Natural threats like earthquakes, hurricanes or tsunamis have had serious impacts on communities. In the past, major earthquakes in the United States like Loma Prieta 1989, Northridge 1994, or recent events in Italy like L’Aquila 2009 or Emilia 2012 emphasized the importance of preparedness and awareness to reduce social impacts. In addition to that, earthquake damaged businesses dramatically reduced the gross regional product. Generating scenario earthquakes in a proper way is important to suitably assess the risk in bridge networks and social losses in terms of gross regional product reduction. Seismic hazard is traditionally assessed by means of probabilistic seismic hazard analysis (PSHA). Although PSHA well represents the hazard at a specific location it is not suitable for spatially distributed systems. Scenario earthquakes can overcome this problem; they represent the actual distribution of ground shaking for a spatially distributed system while being hazard consistent. In this work a methodology to generate scenario earthquakes has been proposed using a novel approach with the aim of being the basic step for investigating possible earthquake consequences in seismic areas and contributing to reduce losses.     

A comparison of two physics-based numerical models for simulating surface water–groundwater interactions

Problems in hydrology and water management that involve both surface water and groundwater are best addressed with simulation models that can represent the interactions between these two flow regimes. In the current generation of coupled models, a variety of approaches is used to resolve surface–subsurface interactions and other key processes such as surface flow propagation. In this study we compare two physics-based numerical models that use a 3D Richards equation representation of subsurface flow. In one model, surface flow is represented by a fully 2D kinematic approximation to the Saint–Venant equations with a sheet flow conceptualization. In the second model, surface routing is performed via a quasi-2D diffusive formulation and surface runoff follows a rill flow conceptualization. The coupling between the land surface and the subsurface is handled via an explicit exchange term resolved by continuity principles in the first model (a fully-coupled approach) and by special treatment of atmospheric boundary conditions in the second (a sequential approach). Despite the significant differences in formulation between the two models, we found them to be in good agreement for the simulation experiments conducted. In these numerical tests, on a sloping plane and a tilted V-catchment, we examined saturation excess and infiltration excess runoff production under homogeneous and heterogeneous conditions, the dynamics of the return flow process, the differences in hydrologic response under rill flow and sheet flow parameterizations, and the effects of factors such as grid discretization, time step size, and slope angle. Low sensitivity to vertical discretization and time step size was found for the two models under saturation excess and homogeneous conditions. Larger sensitivity and differences in response were observed under infiltration excess and heterogeneous conditions, due to the different coupling approaches and spatial discretization schemes used in the two models. For these cases, the sensitivity to vertical and temporal resolution was greatest for processes such as reinfiltration and ponding, although the differences between the hydrographs of the two models decreased as mesh and step size were progressively refined. In return flow behavior, the models are in general agreement, with the largest discrepancies, during the recession phase, attributable to the different parameterizations of diffusion in the surface water propagation schemes. Our results also show that under equivalent parameterizations, the rill and sheet flow conceptualizations used in the two models produce very similar responses in terms of hydrograph shape and flow depth distribution.     

Model‐based attenuation for scattered dispersive waves

Coherent noise in land seismic data primarily consists in source‐generated surface‐wave modes. The component that is traditionally considered most relevant is the so‐called ground roll, consisting in surface‐wave modes propagating directly from sources to receivers. In many geological situations, near?surface heterogeneities and discontinuities, as well as topography irregularities, diffract the surface waves and generate secondary events, which can heavily contaminate records. The diffracted and converted surface waves are often called scattered noise and can be a severe problem particularly in areas with shallow or outcropping hard lithological formations. Conventional noise attenuation techniques are not effective with scattering: they can usually address the tails but not the apices of the scattered events. Large source and receiver arrays can attenuate scattering but only in exchange for a compromise to signal fidelity and resolution. We present a model?based technique for the scattering attenuation, based on the estimation of surface‐wave properties and on the prediction of surface waves with a complex path involving diffractions. The properties are estimated first, to produce surface?consistent volumes of the propagation properties. Then, for all gathers to filter, we integrate the contributions of all possible diffractors, building a scattering model. The estimated scattered wavefield is then subtracted from the data. The method can work in different domains and copes with aliased surface waves. The benefits of the method are demonstrated with synthetic and real data.     

An overview of laboratory apparatuses to measure seismic attenuation in reservoir rocks

Intrinsic wave attenuation at seismic frequencies is strongly dependent on rock permeability, fluid properties, and saturation. However, in order to use attenuation as an attribute to extract information on rock/fluid properties from seismic data, experimental studies on attenuation are necessary for a better understanding of physical mechanisms that are dominant at those frequencies. An appropriate laboratory methodology to measure attenuation at seismic frequencies is the forced oscillation method, but technical challenges kept this technique from being widely used. There is a need for the standardization of devices employing this method, and a comparison of existing setups is a step towards it. Here we summarize the apparatuses based on the forced oscillation method that were built in the last 30 years and were used to measure frequency‐dependent attenuation in fluid‐saturated and/or dry reservoir rocks under small strains (10^?8–10^?5). We list and discuss important technical aspects to be taken into account when working with these devices or in the course of designing a new one. We also present a summary of the attenuation measurements in reservoir rock samples performed with these apparatuses so far.     

DEM assessment of impact forces of dry granular masses on rigid barriers

In the design of sheltering structures/embankments for the mitigation of the risk due to rapid and long spreading landslides, a crucial role is generally played by the assessment of the impact force exerted by the flowing mass on the artificial obstacle. This paper is focused on this issue and in particular on the evaluation of the maximum impact force on the basis of the results obtained by performing an extensive numerical campaign by means of a 3D discrete element code, in which a dry granular mass is schematised as a random distribution of rigid spherical particles. The granular mass is generated just in front of the obstacle: its initial volume, velocity distribution, height, length and porosity are arbitrarily assigned, and the impact process is exclusively analysed. The initial conditions are varied to take a large variety of geometrical/mechanical factors, such as the initial front inclination, its height, the initial void ratio, the length of the impacting mass and the inter-particle friction angle, into consideration. A design formula is also proposed on the base of the obtained results and critically compared with the literature data.     

The Second World Landslide Forum, Rome, 3�C9 October 2011: state of art at May 2011

The World Landslide Forum (WLF) is a triennial mainstream conference aimed at gathering scientists, stakeholders, policy makers and industry members dealing with the management of landslide risk. The First World Landslide Forum, organised by the International Consortium of Landslides, UNESCO, WMO, Food and Agriculture Organization, UNISDR, UNU, UNEP, IBRD, UNDP, ICSU, WFEO, KU and the Japan Landslide Society, was held in 2008 at the United Nations University, Tokyo. The 1st WLF adopted the 2008 Tokyo Declaration “Strengthening the International Programme on Landslides with UNISDR”. The Second World Landslide Forum has the objective to further develop the outcomes of the First Forum in Tokyo 2008 by providing a global cross-cutting information and cooperation platform for all types of organisations representing academia, United Nations organisations, governments, private enterprises and individuals that contribute to landslide research, practice, education and decision making and are willing to strengthen landslide and other related Earth system risk reduction strategies. The emphasis of this forum will be “Putting Science into Practice” with special attention given to actual implementation of technology and research in everyday applications and procedures with the direct involvement of researchers, engineers, private enterprises, stakeholders as well as policy and decision makers. The abstract submission and pre-registration of participants was quite successful with 649 received abstracts and more than 800 participants at May 2011. The full organisation of the event is now in progress taking into consideration the new figure and interest in scientific community and stakeholders.     

Spatially distributed rainfall thresholds for the initiation of shallow landslides

The evaluation of the combined influence of rainfall patterns (in terms of mean intensity and duration) and the geomorphological and mechanical characteristics of hillslopes on their stability conditions is a major goal in the assessment of the shallow landslide triggering processes. Geographic Information Systems (GIS) represent an important tool to develop models that combine hydrological and geomechanical analyses for the evaluation of slope stability, as they allow to combine information concerning rainfall characteristics with topographic and mechanical properties of the slopes over wide areas. In this paper, a GIS-based code is developed to determine physically based intensity/duration rainfall thresholds at the local scale. Given the rainfall duration and the local geometric, hydrological and mechanical characteristics of the slopes, the code evaluates the spatial distribution of the minimum rainfall intensity that triggers shallow landslides and debris flows over a given area. The key feature of the code is the capability of evaluating the time t _p required to reach the peak pore pressure head on the failure surface and computing the corresponding critical intensity/duration thresholds based on post-event peak pore pressures. The reliability of the model is tested using a set of one-dimensional analyses, comparing the physically based thresholds obtained for three different slopes with some empirical rainfall thresholds. In a log–log scale, the thresholds provided by the model decrease linearly with increased rainfall duration and they are bracketed by the empirical thresholds considered. Finally, an example of application to a study area of the Umbria region in central Italy is presented, describing the capability of the model of providing site-specific thresholds for different rainfall scenarios.     

Soil-precipitation feedbacks during the South American Monsoon as simulated by a regional climate model

We summarize the recent progress in regional climate modeling in South America with the Rossby Centre regional atmospheric climate model (RCA3-E), with emphasis on soil moisture processes. A series of climatological integrations using a continental scale domain nested in reanalysis data were carried out for the initial and mature stages of the South American Monsoon System (SAMS) of 1993–92 and were analyzed on seasonal and monthly timescales. The role of including a spatially varying soil depth, which extends to 8 m in tropical forest, was evaluated against the standard constant soil depth of the model of about 2 m, through two five member ensemble simulations. The influence of the soil depth was relatively weak, with both beneficial and detrimental effects on the simulation of the seasonal mean rainfall. Secondly, two ensembles that differ in their initial state of soil moisture were prepared to study the influence of anomalously dry and wet soil moisture initial conditions on the intraseasonal development of the SAMS. In these simulations the austral winter soil moisture initial condition has a strong influence on wet season rainfall over feed back upon the monsoon, not only over the Amazon region but in subtropical South America as well. Finally, we calculated the soil moisture–precipitation coupling strength through comparing a ten member ensemble forced by the same space–time series of soil moisture fields with an ensemble with interactive soil moisture. Coupling strength is defined as the degree to which the prescribed boundary conditions affect some atmospheric quantity in a climate model, in this context a quantification of the fraction of atmospheric variability that can be ascribed to soil moisture anomalies. La Plata Basin appears as a region where the precipitation is partly controlled by soil moisture, especially in November and January. The continental convective monsoon regions and subtropical South America appears as a region with relatively high coupling strength during the mature phase of monsoon development.