Ranking of Placer Gold Prospects in Chile Through Analytic Hierarchy Process

Chile has a rich, but poorly known history of placer gold mining. At present, this sector is almost nonexistent and there are some restrictions for its revival: disperse and partial information on existing resources and limited technical expertise to assess the potential of placer gold mine sites. This paper presents the background, methodology and results of the prioritization process of known prospects of this kind in Chile. This research was part of a publicly funded project aimed to incentivize the development of this industry. The ranking was carried out using the analytic hierarchy process, which allowed to include different quantitative and qualitative variables related to the economic potential, technical aspects, contextual viability and socioeconomic factors in the analysis. The results show that, despite the increasing relevance of environmental and community issues in mining development, the business potential and the economic/technical aspects are the main factors in the early selection of a site to advance in exploration and development activities. Both variables represented around 40% and 37% of weights in the final selection, respectively. In contrast, contextual viability and local socioeconomic impacts only accounted for the remaining 23%. This study also shows that the inclusion of experts with different backgrounds in the process enriches the analysis and does not significantly distort the final outcome of the prioritization. Finally, the relevance of using MCDM tools when assessing the attractiveness of mine sites for their development is highlighted, particularly when public funds for subsequent exploration activities are committed.

     

Late Miocene evaporite geochemistry of Lorca and Fortuna basins (Eastern Betics,SE Spain): Evidence of restriction and continentalization

The Lorca and Fortuna basins are two intramontane Neogene basins located in the eastern Betic Cordillera (SE Spain). During the Late Tortonian—Early Messinian, marine and continental evaporites precipitated in these basins as a consequence of increased marine restriction and isolation. Here we show a stratigraphic correlation between the evaporite records of these basins based on geochemical indicators. We use SO₄ isotope compositions and Sr isotopic ratios in gypsum, and halite Br contents to characterize these units and to identify the marine or continental source of the waters feeding the evaporite basins. In addition, we review the available chronological information used to date these evaporites in Lorca (La Serrata Fm), including a thick saline deposit, that we correlate with the First Evaporitic Group in Fortuna (Los Baños Fm). This correlation is also supported by micropalaeontological data, giving a Late Tortonian age for this sequence. The Second Evaporitic Group, (Chicamo Fm), and the Third Evaporitic Group (Rambla Salada Fm) developed only in Fortuna during the Messinian. According to the palaeogeographical scheme presented here, the evaporites of the Lorca and Fortuna basins were formed during the Late Tortonian—Early Messinian, close to the Betic Seaway closure. Sulphate isotope compositions and Sr isotopic ratios of the Ribera Gypsum Mb, at the base of the Rambla Salada Fm (Fortuna basin), match those of the Late Messinian selenite gypsum beds in San Miguel de Salinas, in the near Bajo Segura basin (40 km to the East), and other Messinian Salinity Crisis gypsum deposits in the Mediterranean. According to these geochemical indicators and the uncertainty of the chronology of this unit, the assignment of the Rambla Salada Fm to the MSC cannot be ruled out.     

Limpets as prey: Lipophrys pholis choice for size

Patellid limpets are key species on rocky shores, structuring intertidal assemblages through their grazing. Their role as prey for fish is, however, often overlooked in studies of intertidal ecosystem functioning. The shanny Lipophrys pholis is a common predator of limpets on rocky shores in Northern Europe and, in this study, we examined the characteristics of its predation on limpets based on prey size. The limpet size range that adult L. pholis (10–12.5 cm) can eat was examined in the laboratory using Patella depressa in 0.2 cm classes, varying between 0.5 cm and 3.0 cm in maximum shell length. There was a limpet size refuge above 1.8 cm, while all smaller sizes were readily consumed by the shanny. The predator attacking behaviour was also examined and found to vary with prey size. Limpets up to 0.8 cm were crushed by the jaw‐constricting force and eaten whole, whereas larger prey were prised from the substratum and the body subsequently separated from the shell. To examine whether there was a L. pholis preference for P. depressa size, a two‐stage laboratory experiment was done using two size classes defined as small (0.7–0.8 cm maximum shell length) and large (1.5–1.6 cm maximum shell length). In the first stage, the predator was given each limpet size class separately. In the second stage, the fish was given a choice between the two classes. Lipophrys pholis was shown to have a preference for the large size class (1.5–1.6 cm). The average number of limpets consumed by the shanny was examined for the duration of one high‐tide typical of the central region of Portugal (≈ 4 h). On average, approximately five limpets (up to an eight limpet maximum) were consumed. The present study shows that L. pholis has a feeding preference based on limpet size, which suggests that the limpet population structure and intra‐ and inter‐specific interactions may be influenced by the shanny predation.     

Role of sea surface temperature and rainfall in determining the stock and fishery of the common octopus (Octopus vulgaris,Mollusca, Cephalopoda) in Tunisia

Octopus (Octopus vulgaris, Mollusca, Cephalopoda) is an important and valuable fishery resource on the eastern and southern coasts of Tunisia, but its landings are highly variable. This paper explores the effect of environment on octopus catch per unit effort (CPUE) during a 12‐year period, through correlation analyses and the incorporation into surplus production models of sea surface temperature (SST) and rainfall data collected during cold (January–May) and hot (August–October) seasons. CLIMPROD software was used to select the appropriate model and fit it to the fishery and environment data. In both seasons, SST significantly contributed to CPUE variability; fishery production was influenced positively by cold season SST but negatively by hot season SST. Due to a poor fit with cold season data, the impact of rainfall was analysed only for the hot season, during which it has a positive effect on production. Results are discussed in view of the life‐cycle of octopus and the dynamics of the Tunisian fishery. This first study of octopus variability in Tunisia highlights the necessity to incorporate environmental influence into stock assessment and management advice.     

The theory of asynchronous relative motion I: time transformations and nonlinear corrections

Using alternative independent variables in lieu of time has important advantages when propagating the partial derivatives of the trajectory. This paper focuses on spacecraft relative motion, but the concepts presented here can be extended to any problem involving the variational equations of orbital motion. A usual approach for modeling the relative dynamics is to evaluate how the reference orbit changes when modifying the initial conditions slightly. But when the time is a mere dependent variable, changes in the initial conditions will result in changes in time as well: a time delay between the reference and the neighbor solution will appear. The theory of asynchronous relative motion shows how the time delay can be corrected to recover the physical sense of the solution and, more importantly, how this correction can be used to improve significantly the accuracy of the linear solutions to relative motion found in the literature. As an example, an improved version of the Clohessy-Wiltshire (CW) solution is presented explicitly. The correcting terms are extremely compact, and the solution proves more accurate than the second and even third order CW equations for long propagations. The application to the elliptic case is also discussed. The theory is not restricted to Keplerian orbits, as it holds under any perturbation. To prove this statement, two examples of realistic trajectories are presented: a pair of spacecraft orbiting the Earth and perturbed by a realistic force model; and two probes describing a quasi-periodic orbit in the Jupiter-Europa system subject to third-body perturbations. The numerical examples show that the new theory yields reductions in the propagation error of several orders of magnitude, both in position and velocity, when compared to the linear approach.     

The Fourier Integral Method: An efficient spectral method for simulation of random fields

The Fourier Integral Method (FIM) of spectral simulation, adapted to generate realizations of a random function in one, two, or three dimensions, is shown to be an efficient technique of non-conditional geostatistical simulation. The main contribution is the use of the fast Fourier transform for both numerical calculus of the density spectral function and as generator of random finite multidimensional sequences with imposed covariance. Results obtained with the FIM are compared with those obtained by other classic methods: Shinozuka and Jan Method in 1D and Turning Bands Method in 2D and 3D, the points for and against different methodologies are discussed. Moreover, with the FIM the simulation of nested structures, one of which can be a nugget effect and the simulation of both zonal and geometric anisotropy is straightforward. All steps taken to implement the FIM methodology are discussed.     

Dynamic faulting under rate-dependent friction

We discuss the effects of rate-dependent friction on the propagation of seismic rupture on active faults. Several physicists using Burridge and Knopoff's box and spring model of faulting have proposed that fault complexity may arise from the spontaneous development of a self-similar stress distribution on the fault plane. If this model proves to be correct, it has important consequences for the origin of the complexity of seismic sources. In order to test these ideas on a more realistic earthquake model, we developed a new boundary integral equation method for studying rupture propagation along an antiplane fault in the presence of nonlinear rate-dependent friction. We study rupture dynamics of models with single and twin asperities. In our models, asperities are places on the fault with a higher value of prestress. Othewise all fault parameters are homogeneous. We show that for models with such asperities, a slip velocity weakening friction leads to the propagation of supersonic healing phases and to the spontaneous arrest of fracture if the prestress outside the asperities is low enough. For models with asperities, we can also observe narrow slip velocity pulses, qualitatively similar to the so-called Heaton pulses observed in some earthquake accelerograms. We also observe a complex distribution of stress after the rupture that depends on details of the initial distribution of asperities and on the details of the friction law.     

Seismic activity related to the 1991 eruption of Colima Volcano,Mexico

Ten years after the last effusive eruption and at least 15 years of seismic quiescence, volcanic seismic activity started at Colima volcano on 14 February 1991, with a seismic crisis which reached counts of more than 100 per day and showed a diversity of earthquake types. Four other distinct seismic crises followed, before a mild effusive eruption in April 1991. The second crisis preceded the extrusion of an andesitic scoriaceous lava lobe, first reported on 1 March; during this crisis an interesting temporary concentration of seismic foci below the crater was observed shortly before the extrusion was detected. The third crisis was constituted by shallow seismicity, featuring possible mild degassing explosion-induced activity in the form of hiccups (episodes of simple wavelets that repeat with diminishing amplitude), and accompanied by increased fumarolic activity. The growth of the new lava dome was accompanied by changing seismicity. On 16 April during the fifth crisis which consisted of some relatively large, shallow, volcanic earthquakes and numerous avalanches of older dome material, part of the newly extruded dome, which had grown towards the edge of the old dome, collapsed, producing the largest avalanches and ash flows. Afterwards, block lava began to flow slowly along the SW flank of the volcano, generating frequent small incandescent avalanches. The seismicity associated with the stages of this eruptive activity shows some interesting features: most earthquake foci were located north of the summit, some of them relatively deep (7–11 km below the summit level), underneath the saddle between the Colima and the older Nevado volcanoes. An apparently seismic quiet region appears between 4 and 7 km below the summit level. In June, harmonic tremors were detected for the first time, but no changes in the eruptive activity could be correlated with them. After June, the seismicity decreasing trend was established, and the effusive activity stopped on September 1991.     

Dynamic displacements and stresses in the vicinity of a cylindrical cavity embedded in a half-space

The two-dimensional response of a viscoelastic half-space containing a buried, unlined, infinitely long cylindrical cavity of circular cross-section subjected to harmonic plane SH, P, SV and Rayleigh waves is obtained by an indirect boundary integral method based on the two-dimensional Green's functions for a viscoelastic half-space. An extensive critical review of the existing numerical results obtained by other techniques is presented together with some new numerical results describing the motion on the ground surface and the motion and stresses on the wall of the cavity for P, SV, SH and Rayleigh waves.     

Seismic response of a cylindrical shell embedded in a layered viscoelastic half-space. I: Formulation

A method to obtain the three-dimensional harmonic response of a infinitely long cylindrical shell of circular cross-section embedded in a layered viscoelastic half-space and subjected to harmonic plane waves impinging at an oblique angle with respect to the axis of the shell is presented. The procedure combines an indirect integral representation for the field in the exterior half-space with a model of the pipeline or tunnel based on Donnell shell theory. The integral representation for the soil is based on the use of moving Green's functions for the layered viscoelastic half-space. The accuracy of the formulation is tested by comparison of results obtained by using different discretizations. Extensive comparisons with previous two- and three-dimensional results for the case of a shell embedded in a uniform half-space and some new numerical results for a shell embedded in a multilayered half-space are presented in a companion paper.