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.     

Elevation‐dependent responses of streamflow to climate warming

Warming will affect snowline elevation, potentially altering the timing and magnitude of streamflow from mountain landscapes. Presently, the assessment of potential elevation‐dependent responses is difficult because many gauged watersheds integrate drainage areas that are both snow and rain dominated. To predict the impact of snowline rise on streamflow, we mapped the current snowline (1980 m) for the Salmon River watershed (Idaho, USA) and projected its elevation after 3 °C warming (2440 m). This increase results in a 40% reduction in snow‐covered area during winter months. We expand this analysis by collecting streamflow records from a new, elevation‐stratified gauging network of watersheds contained within high (2250–3800 m), mid (1500–2250 m) and low (300–1500 m) elevations that isolate snow, mixed and rain‐dominated precipitation regimes. Results indicate that lags between percentiles of precipitation and streamflow are much shorter in low elevations than in mid‐ and high‐elevation watersheds. Low elevation annual percentiles (Q₂₅ and Q₇₅) of streamflow occur 30–50 days earlier than in higher elevation watersheds. Extreme events in low elevations are dominated by low‐ and no‐flow events whereas mid‐ and high‐elevation extreme events are primarily large magnitude floods. Only mid‐ and high‐elevation watersheds are strongly cross correlated with catchment‐wide flow of the Salmon River, suggesting that changes in contributions from low‐elevation catchments may be poorly represented using mainstem gauges. As snowline rises, mid‐elevation watersheds will likely exhibit behaviours currently observed only at lower elevations. Streamflow monitoring networks designed for operational decision making or change detection may require modification to capture elevation‐dependent responses of streamflow to warming. Copyright © 2014 John Wiley & Sons, Ltd.     

3D seismic analysis of the gas hydrate system and the fluid migration paths in part of the Niger Delta Basin,Nigeria

Comprehensive qualitative and semi-quantitative seismic analysis was carried out on 3-dimensional seismic data acquired in the deepwater compressional and shale diapiric zone of the Niger Delta Basin using an advanced seismic imaging tool. The main aim of this work is to obtain an understanding of the forming mechanism of the gas hydrate system, and the fluid migration paths associated with this part of the basin. The results showed the presence of pockmarks on the seafloor and bottom simulating reflectors (BSRs) in the field, indicating the active fluid flux and existence of gas hydrate system in the area. In the area of approximately 195 km² occupying nearly 24% of the entire study field, three major zones with continuous or discontinuous BSRs of 3 to 7 km in length which are in the northeastern, southern and eastern part of the field respectively were delineated. The BSR is interpreted to be the transition between the free gas zone and the gas hydrate zone. The geologic structures including faults (strike-slip and normal faults), chimneys and diapirs were deduced to be the main conduits for gas migration. It is concluded that the biogenic gases generated in the basin were possibly transported via faults and chimneys by advection processes and subsequently accumulated under low temperature and high pressure conditions in the free gas zone below the BSR forming gas hydrate. A plausible explanation for the presence of the ubiquitous pockmarks of different diameters and sizes in the area is the transportation of the excessive gas to the seafloor through these mapped geologic structures.     

Hydrothermal calderas

The standard model of caldera formation is related to the emptying of a magma chamber and ensuing roof collapse during large eruptions or subsurface withdrawal. Although this model works well for numerous volcanoes, it is inappropriate for many basaltic volcanoes (with the notable exception of Hawaii), as these have eruptions that involve volumes of magma that are small compared to the collapse. Many arc volcanoes also have similar oversized depressions, such as Poas (Costa Rica) and Aoba (Vanuatu). In this article, we propose an alternative caldera model based on deep hydrothermal alteration of volcanic rocks in the central part of the edifice. Under certain conditions, the clay-rich altered and pressurized core may flow under its own weight, spread laterally, and trigger very large caldera-like collapse. Several specific mechanisms can generate the formation of such hydrothermal calderas. Among them, we identify two principal modes: mode 1: ripening with summit loading and flank spreading and mode II: unbuttressing with flank subsidence and flank sliding. Processes such as summit loading or flank subsidence may act simultaneously in hybrid mechanisms. Natural examples are shown to illustrate the different modes of formation. For ripening, we give Aoba (Vanuatu) as an example of probable summit loading, while Casita (Nicaragua) is the type example of flank spreading. For unbuttressing, Nuku Hiva Island (Marquesas) is our example for flank subsidence and Piton de la Fournaise (La Réunion) is our example of flank sliding. The whole process is slow and probably needs (a) at least a few tens of thousands of years to deeply alter the edifice and reach conditions suitable for ductile flow and (b) a few hundred years to achieve the caldera collapse. The size and the shape of the caldera strictly mimic that of the underlying weak core. Thus, the size of the caldera is not controlled by the dimensions of the underlying magma reservoir. A collapsing hydrothermal caldera could generate significant phreatic activity and trigger major eruptions from a coexisting magmatic complex. As the buildup to collapse is slow, such caldera-forming events could be detected long before their onset.     

Homogeneous nucleation of amorphous solid water particles in the upper mesosphere

Condensed water particles are known to exist in the high latitude upper mesosphere during the summer months. However, the mechanism or mechanisms through which they nucleate remains uncertain. It is postulated here that particles of amorphous solid water (ASW, condensed water with a non-crystalline structure) may nucleate homogeneously in the summer mesosphere. Using classical nucleation theory and a one-dimensional model, it is shown that more than 10⁵ cm^?3 amorphous solid water particles can nucleate homogeneously under mesopause conditions. Furthermore, it is shown that homogeneous nucleation competes with heterogeneous nucleation on meteoric smoke particles when the cooling rate is >0.5 K/h. The homogeneous nucleation of amorphous solid water could provide an explanation for the high density of ice particles (many thousands per cm³) thought to be required for electron depletions in the upper mesosphere. A parameterisation for homogeneous nucleation is presented which can be used in other mesospheric cloud models.     

Capturing the residence time boundary layer—application to the Scheldt Estuary

At high Peclet number, the residence time exhibits a boundary layer adjacent to incoming open boundaries. In a Eulerian model, not resolving this boundary layer can generate spurious oscillations that can propagate into the area of interest. However, resolving this boundary layer would require an unacceptably high spatial resolution. Therefore, alternative methods are needed in which no grid refinement is required to capture the key aspects of the physics of the residence time boundary layer. An extended finite element method representation and a boundary layer parameterisation are presented and tested herein. It is also explained how to preserve local consistency in reversed time simulations so as to avoid the generation of spurious residence time extrema. Finally, the boundary layer parameterisation is applied to the computation of the residence time in the Scheldt Estuary (Belgium/The Netherlands). This timescale is simulated by means of a depth-integrated, finite element, unstructured mesh model, with a high space–time resolution. It is seen that the residence time temporal variations are mainly affected by the semi-diurnal tides. However, the spring–neap variability also impacts the residence time, particularly in the sandbank and shallow areas. Seasonal variability is also observed, which is induced by the fluctuations over the year of the upstream flows. In general, the residence time is an increasing function of the distance to the mouth of the estuary. However, smaller-scale fluctuations are also present: they are caused by local bathymetric features and their impact on the hydrodynamics.     

Spatial Access to Health Care Services and Disparities in Colorectal Cancer Stage at Diagnosis in Texas

Early diagnosis of colorectal cancer (CRC) is important for improving the survival rate of the disease. Disparities in CRC diagnosis among different population groups have persisted in the United States. This study examines whether spatial access to medical services contributes to disparities of late-stage diagnosis of CRC in Texas. Analysis results suggest that there are significant disparities in late-stage CRC diagnosis by race/ethnicity, socioeconomic status, and geographic location in Texas. It is concluded that spatial access to primary care is associated with CRC stage at diagnosis, but not with racial/ethnic and socioeconomic disparities.     

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.