Crustal thickness estimation beneath the southern central Andes at 30°S and 36°S from S wave receiver function analysis

We have used the S wave receiver function (SRF) technique to investigate the crustal thickness beneath two seismic profiles from the CHARGE project in the southern central Andes. A previous study employing the P wave receiver function method has observed the Moho interface beneath much of the profiles. They found, however, that the amplitude of the P to S conversion was diminished in the western part of the profiles and have attributed it to a reduction of the impedance contrast at the Moho due to lower crustal ecologitization. With SRF, we have successfully detected S to P converted waves from the Moho as well as possible conversions from other lithospheric boundaries. The continental South American crust reaches its maximum thickness of ∼70 km (along 30°S between 70°W and 68.5°W) beneath the Principal Cordillera and the Famatina system and becomes thinner towards the Sierras Pampeanas with a thickness of ∼40 km. Negative phases, possibly related to the base of the continental and oceanic lithosphere, can be recognized in the summation traces at different depths. By comparing our results with data obtained from previous investigations, we are able to further constrain the thickness of the crust and lithosphere beneath the central Andes.     

Interactive 3-D computation of fault surfaces using level sets

Deformable implicit surfaces, implemented with level set methods, have demonstrated a great potential in the computational sciences for applications such as modeling, simulation, and segmentation. They allow implicit handling of complex topologies deformed by operations where large changes can occur without destroying the level set representation. The use of level set techniques for the computation and segmentation of plane-like and high positive curvature features is nontrivial. We present a technique for representing and extracting fault surfaces from 3-D seismic attribute data such that 2-manifolds can be interactively computed and segmented. We present our approach as an interactive workflow that generates highly accurate fault surfaces.     

Sulfur cycling at the Mid-Atlantic Ridge: A multiple sulfur isotope approach

The role of sulfur in two hydrothermal vent systems, the Logatchev hydrothermal field at 14°45′N/44°58′W and several different vent sites along the southern Mid-Atlantic Ridge (SMAR) between 4°48′S and 9°33′S and between 12°22′W and 13°12′W, is examined by utilizing multiple sulfur isotope and sulfur concentration data. Isotope compositions for sulfide minerals and vent H₂S from different SMAR sites range from + 1.5 to + 8.9‰ in δ³⁴S and from + 0.001 to + 0.051‰ in Δ³³S. These data indicate mixing of mantle sulfur with sulfur from seawater sulfate. Combined δ³⁴S and Δ³³S systematics reveal that vent sulfide from SMAR is characterized by a sulfur contribution from seawater sulfate between 25 and 33%. This higher contribution, compared with EPR sulfide, indicates increased seawater sulfate reduction at MAR, because of a deeper seated magma chamber and longer fluid upflow path length, and points to fundamental differences with respect to subsurface structures and fluid evolution at slow and fast spreading mid-ocean ridges.Additionally, isotope data uncover non-equilibrium isotopic exchange between dissolved sulfide and sulfate in an anhydrite bearing zone below the vent systems at fluid temperatures between 335 and 400 °C. δ³⁴S values between + 0.2 to + 8.8‰ for dissolved and precipitated sulfide from Logatchev point to the same mixing process between mantle sulfur and sulfur from seawater sulfate as at SMAR. δ³⁴S values between ? 24.5 and + 6.5‰ and Δ³³S values between + 0.001 and + 0.125‰ for sulfide-bearing sediments and mafic/ultramafic host rocks from drill cores taken in the region of Logatchev indicate a clear contribution of biogenic sulfides formed via bacterial sulfate reduction. Basalts and basaltic glass from SMAR sites with Δ³³S = ? 0.008‰ reveal lower Δ³³S lower values than suggested on the basis of previously published isotopic measurements of terrestrial materials.We conclude that the combined use of both δ³⁴S and Δ³³S provides a more detailed picture of the sulfur cycling in hydrothermal systems at the Mid-Atlantic Ridge and uncovers systematic differences to hydrothermal sites at different mid-ocean ridge sites. Multiple sulfur isotope measurements allow identification of incomplete isotope exchange in addition to isotope mixing as a second important factor influencing the isotopic composition of dissolved sulfide during fluid upflow. Furthermore, based on Δ³³S we are able to clearly distinguish biogenic from hydrothermal sulfides in sediments even when δ³⁴S were identical.     

The response of the South Asian Summer Monsoon circulation to intensified irrigation in global climate model simulations

Agricultural intensification in South Asia has resulted in the expansion and intensification of surface irrigation over the twentieth century. The resulting changes to the surface energy balance could affect the temperature contrasts between the South Asian land surface and the equatorial Indian Ocean, potentially altering the South Asian Summer Monsoon (SASM) circulation. Prior studies have noted apparent declines in the monsoon intensity over the twentieth century and have focused on how altered surface energy balances impact the SASM rainfall distribution. Here, we use the coupled Goddard Institute for Space Studies ModelE-R general circulation model to investigate the impact of intensifying irrigation on the large-scale SASM circulation over the twentieth century, including how the effect of irrigation compares to the impact of increasing greenhouse gas (GHG) forcing. We force our simulations with time-varying, historical estimates of irrigation, both alone and with twentieth century GHGs and other forcings. In the irrigation only experiment, irrigation rates correlate strongly with lower and upper level temperature contrasts between the Indian sub-continent and the Indian Ocean (Pearson’s r = ?0.66 and r = ?0.46, respectively), important quantities that control the strength of the SASM circulation. When GHG forcing is included, these correlations strengthen: r = ?0.72 and r = ?0.47 for lower and upper level temperature contrasts, respectively. Under irrigated conditions, the mean SASM intensity in the model decreases only slightly and insignificantly. However, in the simulation with irrigation and GHG forcing, inter-annual variability of the SASM circulation decreases by ~40 %, consistent with trends in the reanalysis products. This suggests that the inclusion of irrigation may be necessary to accurately simulate the historical trends and variability of the SASM system over the last 50 years. These findings suggest that intensifying irrigation, in concert with increased GHG forcing, is capable of reducing the variability of the simulated SASM circulation and altering the regional moisture transport by limiting the surface warming and reducing land–sea temperature gradients.     

Enhancing the relevance of Shared Socioeconomic Pathways for climate change impacts,adaptation and vulnerability research

This paper discusses the role and relevance of the shared socioeconomic pathways (SSPs) and the new scenarios that combine SSPs with representative concentration pathways (RCPs) for climate change impacts, adaptation, and vulnerability (IAV) research. It first provides an overview of uses of social–environmental scenarios in IAV studies and identifies the main shortcomings of earlier such scenarios. Second, the paper elaborates on two aspects of the SSPs and new scenarios that would improve their usefulness for IAV studies compared to earlier scenario sets: (i) enhancing their applicability while retaining coherence across spatial scales, and (ii) adding indicators of importance for projecting vulnerability. The paper therefore presents an agenda for future research, recommending that SSPs incorporate not only the standard variables of population and gross domestic product, but also indicators such as income distribution, spatial population, human health and governance.     

Coastal pollution hazards in southern California observed by SAR imagery: stormwater plumes, wastewater plumes, and natural hydrocarbon seeps

Stormwater runoff plumes, municipal wastewater plumes, and natural hydrocarbon seeps are important pollution hazards for the heavily populated Southern California Bight (SCB). Due to their small size, dynamic and episodic nature, these hazards are difficult to sample adequately using traditional in situ oceanographic methods. Complex coastal circulation and persistent cloud cover can further complicate detection and monitoring of these hazards. We use imagery from space-borne synthetic aperture radar (SAR), complemented by field measurements, to examine these hazards in the SCB. The hazards are detectable in SAR imagery because they deposit surfactants on the sea surface, smoothing capillary and small gravity waves to produce areas of reduced backscatter compared with the surrounding ocean. We suggest that high-resolution SAR, which obtains useful data regardless of darkness or cloud cover, could be an important observational tool for assessment and monitoring of coastal marine pollution hazards in the SCB and other urbanized coastal regions.     

First‐order exchange coefficient coupling for simulating surface water–groundwater interactions: parameter sensitivity and consistency with a physics‐based approach

Distributed hydrologic models capable of simulating fully‐coupled surface water and groundwater flow are increasingly used to examine problems in the hydrologic sciences. Several techniques are currently available to couple the surface and subsurface; the two most frequently employed approaches are first‐order exchange coefficients (a.k.a., the surface conductance method) and enforced continuity of pressure and flux at the surface‐subsurface boundary condition. The effort reported here examines the parameter sensitivity of simulated hydrologic response for the first‐order exchange coefficients at a well‐characterized field site using the fully coupled Integrated Hydrology Model (InHM). This investigation demonstrates that the first‐order exchange coefficients can be selected such that the simulated hydrologic response is insensitive to the parameter choice, while simulation time is considerably reduced. Alternatively, the ability to choose a first‐order exchange coefficient that intentionally decouples the surface and subsurface facilitates concept‐development simulations to examine real‐world situations where the surface‐subsurface exchange is impaired. While the parameters comprising the first‐order exchange coefficient cannot be directly estimated or measured, the insensitivity of the simulated flow system to these parameters (when chosen appropriately) combined with the ability to mimic actual physical processes suggests that the first‐order exchange coefficient approach can be consistent with a physics‐based framework. Copyright © 2009 John Wiley & Sons, Ltd.     

Field survey of the East Java earthquake and tsunami of June 3, 1994

A field survey of the June 3, 1994 East Java earthquake tsunami was conducted within three weeks, and the distributions of the seismic intensities, tsunami heights, and human and house damages were surveyed. The seismic intensities on the south coasts of Java and Bali Islands were small for an earthquake with magnitudeM 7.6. The earthquake caused no land damage. About 40 minutes after the main shock, a huge tsunami attacked the coasts, several villages in East Java Province were damaged severely, and 223 persons perished. At Pancer Village about 70 percent of the houses were swept away and 121 persons were killed by the tsunami. The relationship between tsunami heights and distances from the source shows that the Hatori's tsunami magnitude wasm=3, which seems to be larger for the earthquake magnitude. But we should not consider this an extraordinary event because it was pointed out byHatori (1994) that the magnitudes of tsunamis in the Indonesia-Philippine region generally exceed 1–2 grade larger than those of other regions.     

A Euler number-based topological computation model for land parcel database updating

Intersection relations are important topological considerations in database update processes. The differentiation and identification of non-empty intersection relations between new updates and existing objects is one of the first steps in the automatic incremental update process for a land parcel database. The basic non-empty intersection relations are meet, overlap, cover, equal and inside, but these basic relationships cannot reflect the complex and detailed non-empty relations between a new update and the existing objects. It is therefore necessary to refine the basic non-empty topological relations to support and trigger the relevant update operations. Such relations have been refined by several researchers using topological invariants (e.g., dimension, type and sequence) to represent the intersection components. However, the intersection components often include only points and lines, and the refined types of 2-dimensional intersection components that occur between land parcels have not been defined. This study examines the refinement of non-empty relations among 2-dimensional land parcels and proposes a computation model. In this model, an entire spatial object is directly used as the operand, and two set operations (i.e., intersection (∩) and difference (\)) are applied to form the basic topological computation model. The Euler number is introduced to refine the relations with a single 2-dimensional intersection (i.e., cover, inside and overlap) and to distinguish the refined types of 2-dimensional intersection components for the relations with multiple intersections. In this study, the cover and overlap relations with single intersections between regions are refined into seven cases, and nine basic types of 2-dimensional intersection components are distinguished. A composite computation model is formed with both Euler number values and dimensional differences. In this model, the topological relations with single intersections are differentiated by the value of the dimension and the Euler number of the resulting set of the whole-object intersection and differences, whereas the relations with multiple intersections are discriminated by the value of the resulting set at a coarse level and are further differentiated by the type and sequence of the whole-object intersection component in a hierarchical manner. Based on the refined topological relations, an improved method for automatic and incremental updating of the land parcel database is presented. The effectiveness of the models and algorithms was verified by the incremental update of a land cover database. The results of this study represent a new avenue for automatic spatial data handling in incremental update processes.