Scale of pluton/wall rock interaction near May Lake,Yosemite National Park,CA, USA

Interaction of magma with wall rock is an important process in igneous petrology, but the mechanisms by which interactions occur are poorly known. The western outer granodiorite of the Cretaceous Tuolumne Intrusive Suite of Yosemite National Park, California, intruded a variety of metasedimentary and igneous wall rocks at 93.1 Ma. The May Lake metamorphic screen is a metasedimentary remnant whose contact zone exhibits a variety of interaction phenomena including xenolith incorporation, disaggregation, and partial melting. The chemical contrast of these metasedimentary rocks with the invading pluton provides an excellent measure of pluton/wall rock interactions. Wall rock xenoliths (mostly pelitic quartzite) are predominantly located in an elongate horizon surrounded by a hybridized fine-grained granodiorite. Initial Sr and Nd isotopic ratios of the hybridized granodiorite indicate significant local incorporation of crustal material. Major- and trace-element geochemical data indicate that contamination of the granodiorite occurred via selective assimilation of both high-K and low-K, high-silica partial melts derived from pelitic quartzite. Although the hybridized granodiorite shows significant amounts of contamination, adjacent to xenoliths the proportion of contamination is undetectable more than a meter away. These results indicate that the chemical and isotopic variability of the Tuolumne Intrusive Suite is not caused by magma contamination via in situ wall rock assimilation.     

Determining magma flow in sills,dykes and laccoliths and their implications for sill emplacement mechanisms

Three-dimensional seismic data from the Faeroe-Shetland Basin provides detailed information on the relationships between sills, dykes, laccoliths and contemporaneous volcanic activity. The data shows that sills are predominantly concave upwards, being complete or partial versions of radially or bilaterally symmetrical forms that possess flat inner saucers connected to a flat outer rim by a steeply inclined sheet. Such morphologies are only partially modified by pre-existing faults. Sills can be sourced from dykes or the steep climbing portions of deeper sills. Both sills and dykes can provide magma to overlying volcanic fissures and sills can be shown to feed shallow laccoliths. Magma flow patterns, as revealed by opacity rendering, suggest that sills propagate upwards and outwards away from the magma feeder. As an individual sill can consist of several leaves emplaced at different stratigraphic levels, and as a sill or dyke can provide magma to volcanic fissures, other sills and laccoliths, the data suggests that neutral buoyancy concepts may not provide a complete explanation for the mechanism and level of sill emplacement. Instead, the data suggests that the presence of lithological contrasts, particularly ductile horizons such as overpressured shales may permit sill formation at any level below the neutrally buoyant level. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Ken Thomson–deceased, April 2007     

Ecological Effects of Shoreline Armoring on Intertidal Habitats of a Puget Sound Urban Estuary

Shoreline armoring is extensive in urban areas worldwide, but the ecological consequences are poorly documented. We mapped shoreline armoring along the Duwamish River estuary (Washington State, USA) and evaluated differences in temperature, invertebrates, and juvenile salmon (Oncorhynchus spp.) diet between armored and unarmored intertidal habitats. Mean substrate temperatures were significantly warmer at armored sites, but water temperature similar to unarmored habitats. Epibenthic invertebrate densities were over tenfold greater on unarmored shorelines and taxa richness double that of armored locations. Taxa richness of neuston invertebrates was also higher at unarmored sites, but abundance similar. We did not detect differences in Chinook (O. tshawytscha) diet, but observed a higher proportion of benthic prey for chum (O. keta) from unarmored sites. Given that over 66% of the Duwamish shoreline is armored—similar to much of south and central Puget Sound—our results underscore the need for further ecological study to address the impacts of estuary armoring.     

Strength criterion for cross-anisotropic sand under general stress conditions

By incorporating the fabric effect and Lode’s angle dependence into the Mohr–Coulomb failure criterion, a strength criterion for cross-anisotropic sand under general stress conditions was proposed. The obtained criterion has only three material parameters which can be specified by conventional triaxial tests. The formula to calculate the friction angle under any loading direction and intermediate principal stress ratio condition was deduced, and the influence of the degree of the cross-anisotropy was quantified. The friction angles of sand in triaxial, true triaxial, and hollow cylinder torsional shear tests were obtained, and a parametric analysis was used to detect the varying characteristics. The friction angle becomes smaller when the major principal stress changes from perpendicular to parallel to the bedding plane. The loading direction and intermediate principal stress ratio are unrelated in true triaxial tests, and their influences on the friction angle can be well captured by the proposed criterion. In hollow cylinder torsional shear tests with the same internal and external pressures, the loading direction and intermediate principal stress ratio are related. This property results in a lower friction angle in the hollow cylinder torsional shear test than that in the true triaxial test under the same intermediate principal stress ratio condition. By comparing the calculated friction angle with the experimental results under various loading conditions (e.g., triaxial, true triaxial, and hollow cylinder torsional shear test), the proposed criterion was verified to be able to characterize the shear strength of cross-anisotropic sand under general stress conditions.     

Uncertainty assessment of the multilayer perceptron (MLP) neural network model with implementation of the novel hybrid MLP-FFA method for prediction of biochemical oxygen demand and dissolved oxygen: a case study of Langat River

Accurate prediction of the chemical constituents in major river systems is a necessary task for water quality management, aquatic life well-being and the overall healthcare planning of river systems. In this study, the capability of a newly proposed hybrid forecasting model based on the firefly algorithm (FFA) as a metaheuristic optimizer, integrated with the multilayer perceptron (MLP-FFA), is investigated for the prediction of monthly water quality in Langat River basin, Malaysia. The predictive ability of the MLP-FFA model is assessed against the MLP-based model. To validate the proposed MLP-FFA model, monthly water quality data over a 10-year duration (2001–2010) for two different hydrological stations (1L04 and 1L05) provided by the Irrigation and Drainage Ministry of Malaysia are used to predict the biochemical oxygen demand (BOD) and dissolved oxygen (DO). The input variables are the chemical oxygen demand (COD), total phosphate (PO₄), total solids, potassium (K), sodium (Na), chloride (Cl), electrical conductivity (EC), pH and ammonia nitrogen (NH₄-N). The proposed hybrid model is then evaluated in accordance with statistical metrics such as the correlation coefficient (r), root-mean-square error, % root-mean-square error and Willmott’s index of agreement. Analysis of the results shows that MLP-FFA outperforms the equivalent MLP model. Also, in this research, the uncertainty of a MLP neural network model is analyzed in relation to the predictive ability of the MLP model. To assess the uncertainties within the MLP model, the percentage of observed data bracketed by 95 percent predicted uncertainties (95PPU) and the band width of 95 percent confidence intervals (d-factors) are selected. The effect of input variables on BOD and DO prediction is also investigated through sensitivity analysis. The obtained values bracketed by 95PPU show about 77.7%, 72.2% of data for BOD and 72.2%, 91.6% of data for DO related to the 1L04 and 1L05 stations, respectively. The d-factors have a value of 1.648, 2.269 for BOD and 1.892, 3.480 for DO related to the 1L04 and 1L05 stations, respectively. Based on the values in both stations for the 95PPU and d-factor, it is concluded that the neural network model has an acceptably low degree of uncertainty applied for BOD and DO simulations. The findings of this study can have important implications for error assessment in artificial intelligence-based predictive models applied for water resources management and the assessment of the overall health in major river systems.     

Enhancement of removal efficiency of ammonia nitrogen in sequencing batch reactor using natural zeolite

Two sets of lab-scale sequencing batch reactors (SBR), i.e., control SBR and SBR using zeolite as carrier (zeo-SBR), were applied to assess nitrogen removal efficiency. The test results revealed that zeolite powder added in SBR could improve its performance. Due to the combination of zeolite adsorption for NH₄ ⁺–N and enhanced simultaneous nitrification and de-nitrification (SND), a higher removal ratio of ammonia nitrogen in wastewater was observed in the test reactor, and the introduction of zeolite powder was helpful to inhabit sludge bulging comparing with the control SBR, in other words, activated sludge immobilized by zeolite powder could remove NH₄ ⁺–N, COD, and PO₄ significantly in a shorter cycle time. Applied two hydraulic retention times (HRTs) showed that the nitrogen and phosphorus removal could be improved while adapting to load variations.     

Large aftershocks triggering by Coulomb failure stress following the 2001<Emphasis Type="Italic">M</Emphasis><Subscript>S</Subscript>=8.1 great Kunlun earthquake

The great Kunlun earthquake occurred on Nov. 14, 2001 in Qinghai Province, China. Five large aftershocks with magnitude larger than 5.0 occurred near the Kunlun fault after main shock. Calculations of the change in Coulomb failure stress reveal that 4 of 5 large aftershocks occurred in areas with Δσ_f>0 (10^?2–10^?1 MPa) and one aftershock occurred in an area with Δσ_f=?0.56 MPa. It is concluded that the permanent fault displacement due to the main shock is the main cause of activity of large aftershocks, but not the whole cause.     

Computing the Sensitivity Kernels for 2.5-D Seismic Waveform Inversion in Heterogeneous,Anisotropic Media

2.5-D modeling and inversion techniques are much closer to reality than the simple and traditional 2-D seismic wave modeling and inversion. The sensitivity kernels required in full waveform seismic tomographic inversion are the Fréchet derivatives of the displacement vector with respect to the independent anisotropic model parameters of the subsurface. They give the sensitivity of the seismograms to changes in the model parameters. This paper applies two methods, called ‘the perturbation method’ and ‘the matrix method’, to derive the sensitivity kernels for 2.5-D seismic waveform inversion. We show that the two methods yield the same explicit expressions for the Fréchet derivatives using a constant-block model parameterization, and are available for both the line-source (2-D) and the point-source (2.5-D) cases. The method involves two Green’s function vectors and their gradients, as well as the derivatives of the elastic modulus tensor with respect to the independent model parameters. The two Green’s function vectors are the responses of the displacement vector to the two directed unit vectors located at the source and geophone positions, respectively; they can be generally obtained by numerical methods. The gradients of the Green’s function vectors may be approximated in the same manner as the differential computations in the forward modeling. The derivatives of the elastic modulus tensor with respect to the independent model parameters can be obtained analytically, dependent on the class of medium anisotropy. Explicit expressions are given for two special cases—isotropic and tilted transversely isotropic (TTI) media. Numerical examples are given for the latter case, which involves five independent elastic moduli (or Thomsen parameters) plus one angle defining the symmetry axis.     

Semi-supervised least squares support vector machine algorithm: application to offshore oil reservoir

At the early stages of deep-water oil exploration and development, fewer and further apart wells are drilled than in onshore oilfields. Supervised least squares support vector machine algorithms are used to predict the reservoir parameters but the prediction accuracy is low. We combined the least squares support vector machine (LSSVM) algorithm with semi-supervised learning and established a semi-supervised regression model, which we call the semi-supervised least squares support vector machine (SLSSVM) model. The iterative matrix inversion is also introduced to improve the training ability and training time of the model. We use the UCI data to test the generalization of a semi-supervised and a supervised LSSVM models. The test results suggest that the generalization performance of the LSSVM model greatly improves and with decreasing training samples the generalization performance is better. Moreover, for small-sample models, the SLSSVM method has higher precision than the semi-supervised K-nearest neighbor (SKNN) method. The new semisupervised LSSVM algorithm was used to predict the distribution of porosity and sandstone in the Jingzhou study area.     

The Wadden Sea in transition - consequences of sea level rise

The impact of sea level rise (SLR) on the future morphological development of the Wadden Sea (North Sea) is investigated by means of extensive process-resolving numerical simulations. A new sediment and morphodynamic module was implemented in the well-established 3D circulation model GETM. A number of different validations are presented, ranging from an idealized 1D channel over a semi-idealized 2D Wadden Sea basin to a fully coupled realistic 40-year hindcast without morphological amplification of the Sylt-Rømøbight, a semi-enclosed subsystem of the Wadden Sea. Based on the results of the hindcast, four distinct future scenarios covering the period 2010–2100 are simulated. While these scenarios differ in the strength of SLR and wind forcing, they also account for an expected increase of tidal range over the coming century. The results of the future projections indicate a transition from a tidal-flat-dominated system toward a lagoon-like system, in which large fractions of the Sylt-Rømøbight will remain permanently covered by water. This has potentially dramatic implications for the unique ecosystem of the Wadden Sea. Although the simulations also predict an increased accumulation of sediment in the back-barrier basin, this accumulation is far too weak to compensate for the rise in mean sea level.