Separation of PP- and PS-wave reflected seismic data using two-dimensional finite offset common-reflection-surface traveltime approximation

Recently, the interest in PS-converted waves has increased for several applications, such as sub-basalt layer imaging, impedance estimates and amplitude-versus-offset analysis. In this study, we consider the problem of separation of PP- and PS-waves from pre-stacked multicomponent seismic data in two-dimensional isotropic medium. We aim to demonstrate that the finite-offset common-reflection-surface traveltime approximation is a good alternative for separating PP- and PS-converted waves in common-offset and common shot configurations by considering a two-dimensional isotropic medium. The five parameters of the finite-offset common-reflection-surface are firstly estimated through the inversion methodology called very fast simulated annealing, which estimates all parameters simultaneously. Next, the emergence angle, one of the inverted parameters, is used to build an analytical separation function of PP and PS reflection separation based on the wave polarization equations. Once the PP- and PS-converted waves were separated, the sections are stacked to increase the signal-to-noise ratio using the special curves derived from finite-offset common-reflection-surface approximation. We applied this methodology to a synthetic dataset from simple-layered to complex-structured media. The numerical results showed that the inverted parameters of the finite offset common-reflection-surface and the separation function yield good results for separating PP- and PS-converted waves in noisy common-offset and common shot gathers.     

The state of the environment in Mexico

The present work evaluates the state of the environment in Mexico based on indicators of the present status of the country’s natural resource management, social and economical conditions and anthropogenic modifications. The Mexican environment is interpreted as a spatially open system having a historical character that is essentially determined by the continual interaction between nature, society and economy. The landscape approach is followed, considering as units of territorial analysis each one of the 145 biophysical environmental units included in the national physiographic regionalization. The assessment of 16 indicators for each biophysical environmental unit was made considering their regional environmental integrity problems, the degree of disarticulation of their structure and function, and the alteration of their territorial structure, all of which determine whether or not they accomplish their environmental functions and achieve environmental stability. The classification of the state of the environment included 5 categories in 8 combinations represented in the map of the state of the environment in Mexico for the year 2008. The map shows that nearly 47.10% of the country’s surface has an environmental status ranging between unstable and critical, the problematic areas being mostly concentrated in the southeast and center of the national territory.     

The evolution of the hydrothermal IOCG system in the Mantoverde district,northern Chile: new evidence from microthermometry and stable isotope geochemistry

New microthermometric data combined with stable isotope geochemistry and paragenetic relationships support a previously suggested cooling–mixing model for the iron oxide–copper–gold mineralization in the Mantoverde district. Fluid inclusions show characteristics of a CO₂-bearing aqueous NaCl ± CaCl₂ salt system. The evolution of the Mantoverde hydrothermal system is characterized by (1) an early hypersaline, high to moderate temperature fluid; (2) a moderate saline, moderate temperature fluid; and (3) a low saline, moderate to low temperature fluid. Early magnetite formation took place at median temperatures of 435.0°C, whereas hematite formed at median temperatures of 334.4°C. The main sulfide mineralization texturally post-dates the iron oxides and occurred before late-stage calcite, which developed at a median temperature of 244.8°C. Boiling occurs only locally and is of no relevance for the ore formation. The microthermometric and stable isotope data are supportive for a fluid cooling and mixing model, and suggestive for a predominantly magmatic–hydrothermal fluid component during the iron oxide and main sulfide mineralization. Thereafter, the incursion of a nonmagmatic fluid of ultimately meteoric or seawater gains more importance.     

Pliocene extensional tectonics in the Eastern Central Patagonian Cordillera: geochronological constraints and new field evidence

Recent field work and review of radiometric data obtained from Neogene lavas and plutonic rocks exposed in the Eastern Central Patagonian Cordillera (46–48ºS), which overlie subducted segments of the South Chile Ridge, suggest important Late Miocene to Pleistocene morphological changes in relation to base level variations and/or tectonic events. We present new field observations from a region south of Lago General Carrera‐Buenos Aires, between the main Cordillera and the Meseta del Lago Buenos Aires, demonstrating that normal faulting controlled valley incisions and occurred during lava emplacement at 5–4 Ma and after 3 Ma. We also show that the 12 Ma basaltic flows of the Meseta del Lago Buenos Aires (∼2000 m a.s.l.) have been subjected to deep incision, with younger lavas dated at 1.2 Ma partially filling the valleys. These incisions are thought to reflect progressive eastward tilting of the entire meseta. Our new observations, together with additional features from Central Patagonia, strongly suggest that tectonic events led to a regional widespread morphological change after 5 Ma. The coincidence in time and space between the subduction of segments of the South Chile Ridge at 6 and 3 Ma causing opening of a slab window, and strong base level variations in the studied area, suggests a cause‐and‐effect relationship. In Central Patagonia, compressional tectonics ended well before extensional events reported here. Causes of uplift and further extension are probably completely disconnected. The uplift is purely tectonic in origin and occurred prior to the subduction of the South Chile Ridge. Extension should be a consequence of this subduction.     

EVA: GPS-based extended velocity and acceleration determination

In this work, a new GPS carrier phase-based velocity and acceleration determination method is presented that extends the effective range of previous techniques. The method is named ‘EVA’, and may find applications in fields such as airborne gravimetry when rough terrain or water bodies make difficult or impractical to set up nearby GPS reference receivers. The EVA method is similar to methods such as Kennedy (Precise acceleration determination from carrier phase measurements. In: Proceedings of the 15th international technical meeting of the satellite division of the Institute of Navigation. ION GPS 2002, Portland pp 962–972, 2002b) since it uses L1 carrier phase observables for velocity and acceleration determination. However, it introduces a wide network of stations and it is independent of precise clock information because it estimates satellite clock drifts and drift rates ‘on-the-fly’, requiring only orbit data of sufficient quality. Moreover, with EVA the solution rate is only limited by data rate, and not by the available precise satellite clocks data rate. The results obtained are more robust for long baselines than the results obtained with the reference Kennedy method. An advantage of being independent of precise clock information is that, beside IGS Final products, also the Rapid, Ultra-Rapid (observed) and Ultra-Rapid (predicted) products may be used. Moreover, the EVA technique may also use the undifferenced ionosphere-free carrier phase combination (LC), overcoming baseline limitations in cases where ionosphere gradients may be an issue and very low biases are required. During the development of this work, some problems were found in the velocity estimation process of the Kennedy method. The sources of the problems were identified, and an improved version of the Kennedy method was used for this research work. An experiment was performed using a light aircraft flying over the Pyrenees, showing that both EVA and the improved Kennedy methods are able to cope with the dynamics of mountainous flight. A RTK-derived solution was also generated, and when comparing the three methods to a known zero-velocity reference the results yielded similar performance. The EVA and the improved-Kennedy methods outperformed the RTK solutions, and the EVA method provided the best results in this experiment. Finally, both the improved version of the Kennedy method and the EVA method were applied to a network in equatorial South America with baselines of more than 1,770 km, and during local noon. Under this tough scenario, the EVA method showed a clear advantage for all components of velocity and acceleration, yielding better and more robust results.     

ACE‐2 multiple angle micro‐pulse lidar observations from Las Galletas, Tenerife, Canary Islands

Multiple‐angle micro‐pulse lidar (MPL) observations were made at Las Galletas on Tenerife, Canary Islands during the Aerosol Characterization Experiment‐2 (ACE‐2) conducted June–July, 1997. A principal objective of the MPL observations was to characterize the temporal/spatial distributions of aerosols in the region, particularly to identify and profile elevated Saharan dust layers which occur intermittently during the June–July time period. Vertical and slant angle measurements taken 16 and 17 July characterize such an occurrence, providing aerosol backscatter, extinction, and optical depth profiles of the dust layer between 1 and 5 km above mean sea level (MSL). Additionally, horizontal measurements taken in Las Galletas throughout the 6‐week period provide a time profile of the varying aerosol extinction at the surface. This profile exhibits the alternating periods of clean maritime air and pollution outbreaks that typified the region. Horizontal measurements also provide some evidence suggesting the possible influx of Saharan dust from the free troposphere to the surface. This paper presents estimates of aerosol optical properties retrieved from the multi‐angle MPL measurements in addition to an outline of the methodologies employed to obtain these results.     

Time-Dependent Modulation of Cosmic Rays in the Heliosphere

The time-dependent modulation of galactic cosmic rays in the heliosphere is studied by computing intensities using a time-dependent modulation model. By introducing recent theoretical advances in the transport coefficients in the model, computed intensities are compared with Voyager 1, International Monitoring Platform (IMP) 8, and Ulysses proton observations in search of compatibility. The effect of different modulation parameters on computed intensities is also illustrated. It is shown that this approach produces, on a global scale, realistic cosmic-ray proton intensities along the Voyager 1 spacecraft trajectory and at Earth up to ≈?2004, whereafter the computed intensities recover much more slowly towards solar minimum than observed in the inner heliosphere. A modified time dependence in the diffusion coefficients is proposed to improve compatibility with the observations at Earth after ≈?2004. This modified time dependence led to an improved compatibility between computed intensities and the observations along the Voyager 1 trajectory and at Earth even after ≈?2004. An interesting result is that the cosmic-ray modulation during the current polarity cycle is not determined only by changes in the drift coefficient and tilt angle of the wavy current sheet, but is also largely dependent on changes in the diffusion coefficients.     

Consequences of winter upwelling events on biogeochemical and phytoplankton patterns in a western Galician ria (NW Iberian peninsula)

The consequences of two upwelling events in mid- (MW) and late (LW) winter on biogeochemical and phytoplankton patterns were studied in the Pontevedra Ria and compared with the patterns measured under typical winter conditions and under a summer upwelling event. Thermohaline patterns measured during the mid-winter upwelling event (MW-up) revealed the intrusion of saltier seawater (35.9) into the ria associated with the Iberian Poleward Current (IPC). During the late-winter upwelling event (LW-up), the seawater which had welled up into the ria showed characteristics of the Eastern North Atlantic Central Water mass (ENACW). In both cases the measured water residence time (4 days during MW-up and 10 days during LW-up) was related to both meteorological and fluvial forcing. This residence time contrasts with that of summer upwelling (7 days) and with that estimated under unfavorable upwelling atmospheric conditions (2–4 weeks). During MW-up, the ria became poor in nutrients due to continental freshwater dilution, associated with the shorter residence time of the water, and the intrusion of IPC, which is a water body poor in nutrient salts: 2.9 μM of nitrate, 0.1 μM of phosphate and 1.5 μM of silicate. During this event, the ria exported 3.4 mol_DIN s⁻¹, compared with 6.9 mol_DIN s⁻¹ in non-upwelling conditions. Phytoplankton showed a uniform distribution throughout the ria, as during unfavorable upwelling conditions, and was characterized by the dominance of diatoms, mainly Nitzschia longissima and Skeletonema costatum. During LW-up, a nutrient depletion in the photic layer also occurred, but as a result of a phytoplankton spring bloom developing at this time. The ria was a nutrient trap where 4.1 mol_DIN s⁻¹ were processed by photosynthesis. This budget is three times higher than the one under non-upwelling conditions. In contrast with the MW-up, which had no effect on primary production, during LW-up the ria became more productive, although not as productive as during a summer upwelling event (9.9 mol_DIN s⁻¹). The taxonomic composition of the phytoplankton community did not change noticeably during LW-up and the summer upwelling, with the same species present and changing only in relative proportions. Diatoms were always the dominant microphytoplankton community, with Pseudonitzschia pungens, Thalassionema nitzschioides and several species of Chaetoceros as characteristic taxons.     

Mathematical model for one-dimensional desiccation and consolidation of sedimented soils

A mathematical model has been developed to represent the physical phenomena that occur during the desiccation and one-dimensional consolidation of successive layers of hydraulically transported sediments as they are periodically deposited in a containment area. The governing boundary value problem, defined in terms of pore water pressures, consists of two field equations (one for the saturated domain and one for the unsaturated domain), a drainage boundary condition, an evapotranspiration boundary condition and a series of continuity conditions at the interfaces between different layers. A number of simplifying assumptions were made to render the field equations tractable, and a step-by-step numerical procedure was used to solve the linearized boundary value problem; at the end of each step, the errors introduced by the simplifying assumptions were corrected. Based on a thorough study of the convergence and stability conditions associated with the numerical approximation employed, a system of automatic corrections was incorporated into the computer program to reduce the time increment if stability problems originate during the solution. Based on the results of a parameter study, it was found that, although drainage conditions at the bottom of the, layer do exert some influence on the consolidation rate during deposition of the dredgings, the evapotranspiration potential renders this effect almost negligible when desiccation takes place at the surface. In contrast, transpiration plays an important, role on the consolidation rate during the early stages of desiccation, but its effect is reduced considerably as the water table approaches an equilibrium position.