Hybrid seismic modeling based on discrete-wave number and finite-difference methods

Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources.     

Reconstructing deep‐time histories from integrated thermochronology: An example from southern Baffin Island,Canada

We present a multi‐chronometric approach for reconstructing deep‐time thermal histories using southern Baffin Island as a case study. This continuous thermal history begins with the Palaeoproterozoic Trans‐Hudson Orogeny and is derived from inverse and forward models that integrate thermochronometers spanning some 500°C: new apatite U–Pb ages and K‐feldspar ⁴⁰Ar/³⁹Ar multi‐diffusion domain data, published (U–Th)/He zircon ages and new multi‐kinetic fission‐track results. Integration of data from a wider temperature range reduces ambiguities in thermal‐history modelling and permits us to constrain the timing of geological processes including, extended post‐orogenic cooling, enhanced later Proterozoic cooling, and then episodic burial and exhumation in the Palaeozoic–Mesozoic.     

Relative roles of climate sensitivity and forcing in defining the ocean circulation response to climate change

The response of the ocean’s meridional overturning circulation (MOC) to increased greenhouse gas forcing is examined using a coupled model of intermediate complexity, including a dynamic 3-D ocean subcomponent. Parameters are the increase in CO₂ forcing (with stabilization after a specified time interval) and the model’s climate sensitivity. In this model, the cessation of deep sinking in the north “Atlantic” (hereinafter, a “collapse”), as indicated by changes in the MOC, behaves like a simple bifurcation. The final surface air temperature (SAT) change, which is closely predicted by the product of the radiative forcing and the climate sensitivity, determines whether a collapse occurs. The initial transient response in SAT is largely a function of the forcing increase, with higher sensitivity runs exhibiting delayed behavior; accordingly, high CO₂-low sensitivity scenarios can be assessed as a recovering or collapsing circulation shortly after stabilization, whereas low CO₂-high sensitivity scenarios require several hundred additional years to make such a determination. We also systemically examine how the rate of forcing, for a given CO₂ stabilization, affects the ocean response. In contrast with previous studies based on results using simpler ocean models, we find that except for a narrow range of marginally stable to marginally unstable scenarios, the forcing rate has little impact on whether the run collapses or recovers. In this narrow range, however, forcing increases on a time scale of slow ocean advective processes results in weaker declines in overturning strength and can permit a run to recover that would otherwise collapse.     

Correlation of the major late Jurassic —early Tertiary low- and highstand cycles of south-west Egypt and north-west Sudan

The mainly continental deposits of northwest Sudan and south-west Egypt have been correlated with coeval shallow marine and marine deposits in northern Egypt along a north-south running cross-section, based on surface and subsurface data. The palaeodepth curve of northern Egypt illustrates the gradual seal-level rise, reaching its maximum during the Late Cretaceous with conspicuous advances during the Aptian and late Cenomanian. A general highstand is also recorded during the Campanian-Maastrichtian in north-west Sudan. A detailed facies correlation is given for the Aptian and late Cenomanian highstand in western Egypt. The correlation of the Cenomanian Bahariya and Maghrabi formations displays short-term relative sealevel fluctuations. The interpretation illustrates the extensiveness of related erosional processes in the hinterland, partly intensified by temporarily uplift of the Uweinat-Aswan High in the south. Regional uplift and constant erosion took place in south-west Egypt during Coniacian and Santonian times. The regional stratigraphic gaps and uncertain interpretation of the Bahariya Uplift are induced by the influence of the Trans-African Lineament, especially during the Late Cretaceous. Low-stand fluvial sheet sandstones characterized by non-cyclic sequence development and high facies stability occur, especially in the Neocomian and early Turonian. During the Barremian and Albian, fluvial architecture changes to more cyclic fluvial sequences and increasing soil formation, due to increasing subsidence, more humid climatic conditions and the generally rising sea level, culminating in the extensive shallow marine Abu Ballas and Maghrabi formations.     

The influence of bedrock orientation on the landscape evolution,surface morphology and denudation (10Be) at the Niesen,Switzerland

Landscape evolution and surface morphology in mountainous settings are a function of the relative importance between sediment transport processes acting on hillslopes and in channels, modulated by climate variables. The Niesen nappe in the Swiss Penninic Prealps presents a unique setting in which opposite facing flanks host basins underlain by identical lithologies, but contrasting litho‐tectonic architectures where lithologies either dip parallel to the topographic slope or in the opposite direction (i.e. dip slope and non‐dip slope). The north‐western facing Diemtigen flank represents such a dip slope situation and is characterized by a gentle topography, low hillslope gradients, poorly dissected channels, and it hosts large landslides. In contrast, the south‐eastern facing Frutigen side can be described as non‐dip slope flank with deeply incised bedrock channels, high mean hillslope gradients and high relief topography. Results from morphometric analysis reveal that noticeable differences in morphometric parameters can be related to the contrasts in the relative importance of the internal hillslope‐channel system between both valley flanks. While the contrasting dip‐orientations of the underlying flysch bedrock has promoted hillslope and channelized processes to contrasting extents and particularly the occurrence of large landslides on the dip slope flank, the flank averaged beryllium‐10 (¹⁰Be)‐derived denudation rates are very similar and range between 0.20 and 0.26 mm yr^?1. In addition, our denudation rates offer no direct relationship to basin's slope, area, steepness or concavity index, but reveal a positive correlation to mean basin elevation that we interpret as having been controlled by climatically driven factors such as frost‐induced processes and orographic precipitation. Our findings illustrate that while the landscape properties in this part of the northern Alpine border can mainly be related to the tectonic architecture of the underlying bedrock, the denudation rates have a strong orographic control through elevation dependent mean annual temperature and precipitation. Copyright © 2013 John Wiley & Sons, Ltd.     

On the Application of the Novel Thin Spray-on Liner (TSL): A Progress Report in Mining Operations

Stability is the main problem in underground excavations that can result in rock-related accidents or closure of mines. This is often controlled by the geological condition of the rockmass and the in situ stress field where excavation is being carried out. Over the years, conventional support systems such as bolts, shotcrete and mesh have been the standard support systems used in difficult ground conditions to maintain stability of underground openings. The novel Thin Spray Liner (TSL) is a technology with an improvement in terms of support capabilities for difficult geological conditions. This technology is predicted to replace conventional support systems, as it offers a significant yielding capacity that is better than that of the traditional techniques of controlling strata. The performance of the TSL lies between that of shotcrete and wire mesh. The TSL has numerous advantages regarding rock reinforcement such as rapid curing period, fast application, high tensile strength, high adhesion properties that prevent movement of rock keyblocks, improved cycle times, reduction in material-handling compared to shotcrete, and a significant ability to penetrate joints. These advantages increase the allowance for mechanisation and improvement of safety in excavations. Despite various benefits of TSLs, there are shortcomings which make end-users to be sceptical about the technology. This study provides a detailed review of the mechanism, properties, composition and types, as well as the general and mining applications of TSLs. It also describes the challenges faced by the mining industry and other end-users in accepting the technology as a solution for stability issues in underground openings. In addition, some insightful suggestions are made as recommendations for future work on the development of suitable, effective and efficient thin surface rock support technology within limitations/constraints.

     

Presolar grains in the CM2 chondrite Sutter's Mill

The Sutter's Mill (SM) carbonaceous chondrite is a regolith breccia, composed predominantly of CM2 clasts with varying degrees of aqueous alteration and thermal metamorphism. An investigation of presolar grains in four Sutter's Mill sections, SM43, SM51, SM2‐4, and SM18, was carried out using NanoSIMS ion mapping technique. A total of 37 C‐anomalous grains and one O‐anomalous grain have been identified, indicating an abundance of 63 ppm for presolar C‐anomalous grains and 2 ppm for presolar oxides. Thirty‐one silicon carbide (SiC), five carbonaceous grains, and one Al‐oxide (Al₂O₃) were confirmed based on their elemental compositions determined by C‐N‐Si and O‐Si‐Mg‐Al isotopic measurements. The overall abundance of SiC grains in Sutter's Mill (55 ppm) is consistent with those in other CM chondrites. The absence of presolar silicates in Sutter's Mill suggests that they were destroyed by aqueous alteration on the parent asteroid. Furthermore, SM2‐4 shows heterogeneous distributions of presolar SiC grains (12–54 ppm) in different matrix areas, indicating that the fine‐grained matrix clasts come from different sources, with various thermal histories, in the solar nebula.     

Identification of Mineralization in Geochemistry Along a Transect Based on the Spatial Curvature of Log-Ratios

Detecting subcropping mineralizations but also deeply buried mineralizations is one important goal in geochemical exploration. The identification of useful indicators for mineralization is a difficult task, as mineralization might be influenced by many factors, including location, investigated media and depth. Here, a statistical method is proposed which indicates chemical elements related to mineralization along a transect. Moreover, the method determines the potential area of the deposit along a transect. The identification is based on general additive models (GAMs) for the element concentrations across the spatial coordinate(s). The log-ratios of the GAM fits are taken to compute the curvature, where high and narrow curvature is supposed to indicate the mineralization area. By defining a measure for the quantification of high curvature, the log-ratios can be ranked, and elements can be identified that are indicative of the anomaly patterns.

     

Mesoscale surface circulation and variability of Southern Indian Ocean derived by combining satellite altimetry and drifter observations

High resoultion Eulerian mean velocity field has been derived by combining the satellite tracked surface drifter data with satellite altimetry and ocean surface winds. The drifter data used in this study includes Argos and surface drifter data from Global Drifter Program. Maps of Sea Level Anomaly(MSLA) weekly files with a resolution of(1/3)° in both Latitude and Longitude for the period 1993–2012 have been used. The Ekman current is computed using ocean surface mean wind fields from scatterometers onboard ERS 1/2,Quikscat and ASCAT. The derived mean velocity field exhibits the broad flow of Antarctic Circumpolar Current with speeds up to 0.6 m/s.Anomalous field is quite significant in the western part between 20° and 40°E and in the eastern part between 80°E and 100°E with velocity anomaly up to 0.3 m/s. The estimated mean flow pattern well agrees with the dynamic topography derived from in-situ observations. Also,the derived velocity field is consistent with the in-situ ADCP current measurements. Eddy kinetic energy illustrates an increasing trend during 1993–2008 and is in phase coherence with the Southern Annular Mode by three month lag. Periodic modulations are found in the eddy kinetic energy due the low frequency Antarctic Circumpolar Wave propagation.     

Evaluation of Inference Adequacy in Cumulative Logistic Regression Models: An Empirical Validation of ISWRidge Relationships

Internal solitary wave propagation over a submarine ridge results in energy dissipation, in which the hydrodynamic interaction between a wave and ridge affects marine environment. This study analyzes the effects of ridge height and potential energy during wave-ridge interaction with a binary and cumulative logistic regression model. In testing the Global Null Hypothesis, all values are p 〈0.001, with three statistical methods, such as Likelihood Ratio, Score, and Wald. While comparing with two kinds of models, tests values obtained by cumulative logistic regression models are better than those by binary logistic regression models. Although this study employed cumulative logistic regression model, three probability functions p^1, p^2 and p^3, are utilized for investigating the weighted influence of factors on wave reflection. Deviance and Pearson tests are applied to cheek the goodness-of-fit of the proposed model. The analytical results demonstrated that both ridge height （X1 ） and potential energy （X2 ） significantly impact （p 〈 0. 0001 ） the amplitude-based refleeted rate; the P-values for the deviance and Pearson are all 〉 0.05 （0.2839, 0.3438, respectively）. That is, the goodness-of-fit between ridge height （ X1 ） and potential energy （X2） can further predict parameters under the scenario of the best parsimonious model. Investigation of 6 predictive powers （ R2, Max-rescaled R^2, Sorners＇ D, Gamma, Tau-a, and c, respectively） indicate that these predictive estimates of the proposed model have better predictive ability than ridge height alone, and are very similar to the interaction of ridge height and potential energy. It can be concluded that the goodness-of-fit and prediction ability of the cumulative logistic regression model are better than that of the binary logistic regression model.