Increase in Dissolved Silica of Rivers Due to a Volcanic Eruption in an Estuarine Bay (Sorsogon Bay,Philippines)

Mount Bulusan, the Philippines’ fourth most active volcano, erupted in February 21, 2011, sending volcanic ash and pyroclastic materials to its surrounding rivers. The waters drained into the estuary of harmful algal blooms plagued Sorsogon Bay. We aim to determine the impact of the 2011 volcanic eruption and the preceding volcanic ash emissions to the dissolved silica concentration of rivers draining the flanks of Mt. Bulusan and its possible implications to the phytoplankton assemblage of the bay. Six river water sampling periods from August 2010 to October 2012 overlapped with Mt. Bulusan’s active phase of volcanism. Our data shows that mean river silica from pre-eruption levels of ~?500 μM increased by more than 200% during and post-eruption. Highest Si concentration of 2270 μM was measured from Cadacan River in August 2011. Here, we argue that the sustained general increase of dissolved silica is due to the silica-containing materials from Mt. Bulusan’s eruption and that their concentration in river waters is also a function of watershed lithology and precipitation. Increase in dissolved silica and other nutrients caused a shift to diatom domination and, possibly, termination of Pyrodinium bahamense var. compressum blooms. Silica load increase in embayments is a natural process that controls the dominance of algae. Our study also highlights the importance of Philippine rivers to the global ocean silica budget as a function of high precipitation, tectonics in general, and volcanism in particular.     

Potential Salinity and Temperature Futures for the Chesapeake Bay Using a Statistical Downscaling Spatial Disaggregation Framework

Estuaries are productive and ecologically important ecosystems, incorporating environmental drivers from watersheds, rivers, and the coastal ocean. Climate change has potential to modify the physical properties of estuaries, with impacts on resident organisms. However, projections from general circulation models (GCMs) are generally too coarse to resolve important estuarine processes. Here, we statistically downscaled near-surface air temperature and precipitation projections to the scale of the Chesapeake Bay watershed and estuary. These variables were linked to Susquehanna River streamflow using a water balance model and finally to spatially resolved Chesapeake Bay surface temperature and salinity using statistical model trees. The low computational cost of this approach allowed rapid assessment of projected changes from four GCMs spanning a range of potential futures under a high CO₂ emission scenario, for four different downscaling methods. Choice of GCM contributed strongly to the spread in projections, but choice of downscaling method was also influential in the warmest models. Models projected a ~2–5.5 °C increase in surface water temperatures in the Chesapeake Bay by the end of the century. Projections of salinity were more uncertain and spatially complex. Models showing increases in winter-spring streamflow generated freshening in the Upper Bay and tributaries, while models with decreased streamflow produced salinity increases. Changes to the Chesapeake Bay environment have implications for fish and invertebrate habitats, as well as migration, spawning phenology, recruitment, and occurrence of pathogens. Our results underline a potentially expanded role of statistical downscaling to complement dynamical approaches in assessing climate change impacts in dynamically challenging estuaries.     

Revisit the calibration errors on experimental slant total electron content (TEC) determined with GPS

The calibration errors on experimental slant total electron content (TEC) determined with global positioning system (GPS) observations is revisited. Instead of the analysis of the calibration errors on the carrier phase leveled to code ionospheric observable, we focus on the accuracy analysis of the undifferenced ambiguity-fixed carrier phase ionospheric observable determined from a global distribution of permanent receivers. The results achieved are: (1) using data from an entire month within the last solar cycle maximum, the undifferenced ambiguity-fixed carrier phase ionospheric observable is found to be over one order of magnitude more accurate than the carrier phase leveled to code ionospheric observable and the raw code ionospheric observable. The observation error of the undifferenced ambiguity-fixed carrier phase ionospheric observable ranges from 0.05 to 0.11 total electron content unit (TECU) while that of the carrier phase leveled to code and the raw code ionospheric observable is from 0.65 to 1.65 and 3.14 to 7.48 TECU, respectively. (2) The time-varying receiver differential code bias (DCB), which presents clear day boundary discontinuity and intra-day variability pattern, contributes the most part of the observation error. This contribution is assessed by the short-term stability of the between-receiver DCB, which ranges from 0.06 to 0.17 TECU in a single day. (3) The remaining part of the observation errors presents a sidereal time cycle pattern, indicating the effects of the multipath. Further, the magnitude of the remaining part implies that the code multipath effects are much reduced. (4) The intra-day variation of the between-receiver DCB of the collocated stations suggests that estimating DCBs as a daily constant can have a mis-modeling error of at least several tenths of 1 TECU.     

The Effect of Key Parameters on the Strength of a Dispersive Soil Stabilized with Sustainable Binders

Geotechnical and Geological Engineering - Dispersive soils are susceptible to phenomena of internal and external erosion when in contact with relatively pure water due to its particle’s...     

The Colider and Roosevelt Volcanic rocks (sw amazonian craton): geochemistry and sm-nd isotope characteristics of a silicic large igneous province

The volcanic rocks of the Colíder and Roosevelt formations are extensively exposed in the south-central portion of the Amazonian Craton where effusive and pyroclastic rocks have been mapped. Both units, topped by chemical sediments and oceanic facies as rhyolite and andesite lavas, rhyodacite, and porphyritic dacite, with frequent intercalations of pyroclastic and epiclastic deposits. Whole-rock geochemistry for 55 samples of rhyolitic to andesitic composition suggests the involvement of fertile mantle-derived components with E-MORB to OIB compositions. The analyzed rocks display calc-alkaline to shoshonitic affinity consistent with generation related to an active continental margin. The whole-rock Sm-Nd isotope data from selected felsic volcanic rocks of the Colíder and Roosevelt formations yield negative initial ε_Nd values between –3 and –9, indicating the predominantly crustal nature of the parental magmas with early Archean to late Paleoproterozoic (ca. 2.5–2.0 Ga) depleted mantle model ages.

     

Modelling size effect on rock aggregates strength using a DEM bonded-cell model

Empirical evidence has shown that particle breakage affects the mechanical behaviour of granular materials. The source of this mechanism takes place at the particle scale, and the main consequence on the macromechanical behaviour is increasing compressibility. Due to the inverse correlation between particle size and particle crushing strength, coarse rockfill materials are particularly vulnerable to mechanical degradation due to particle breakage. However, such coarse materials do not fit in standard laboratory devices, and the alternative of large sample testing is usually unavailable or too expensive. Alternatively, recent works have proposed multi-scale approaches using the discrete element method (DEM) to carry out numerical testing of coarse crushable materials, although few studies have focused on size effects. This article presents the application of a DEM bonded-cell model to study particle size-strength correlation on angular rock aggregates. Each particle is modelled by a cluster of perfectly rigid polyhedral cells with Mohr–Coulomb contact law. Constant cell density within particles implies that the presence of potential fragmentation planes increases with size. Therefore, particle strength decreases with size. A comprehensive sensitivity analysis was carried out through 1477 particle crushing simulations in a given particle size. Based on published experimental data on calcareous rock aggregates, part of the simulations were used for calibration, and 97 additional simulations of a coarser size fraction were performed for validation. The results show a good agreement with the empirical data in terms of size effect and data scatter through Weibull statistics.

     

Adding the missing third and fourth dimensions to trajectory analysis in carbonate systems

We developed a seismic geomorphology-based procedure to enhance traditional trajectory analysis with the ability to visualize and quantify lateral variability along carbonate prograding-margin types (ramps and rimmed shelves) in 3D and 4D. This quantitative approach analysed the shelf break geometric evolution of the Oligo-Miocene carbonate clinoform system in the Browse Basin and delineated the feedback between antecedent topography and carbonate system response as controlling factor on shelf break rugosity. Our geometrical analysis identified a systematic shift in the large-scale average shelf break strike direction over a transect of 10 km from 62° to 55° in the Oligo-Miocene interval of the Browse Basin, which is likely controlled by far-field allogenic forcing from the Timor Trough collision zone. Plotting of 3D shelf break trajectories represents a convenient way to visualize the lateral variability in shelf break evolution. Shelf break trajectories that indicate contemporaneous along-strike progradation and retrogradation correlate with phases of autogenic slope system re-organization and may be a proxy for morphological stability of the shelf break. Shelf break rugosity and shelf break trajectory rugosity are not inherited parameters and antecedent topography does not dictate long-term differential movement of the shelf margin through successive depositional sequences. The autogenic carbonate system response to antecedent topography smooths high-rugosity areas by filling accommodation and maintains a relatively constant shelf break rugosity of ~150 m. Color-coding of the vertical component in the shelf break trajectory captures the creation and filling of accommodation, and highlights areas of the transect that are likely to yield inconsistent 2D sequence stratigraphic interpretations.     

Ranking of Placer Gold Prospects in Chile Through Analytic Hierarchy Process

Chile has a rich, but poorly known history of placer gold mining. At present, this sector is almost nonexistent and there are some restrictions for its revival: disperse and partial information on existing resources and limited technical expertise to assess the potential of placer gold mine sites. This paper presents the background, methodology and results of the prioritization process of known prospects of this kind in Chile. This research was part of a publicly funded project aimed to incentivize the development of this industry. The ranking was carried out using the analytic hierarchy process, which allowed to include different quantitative and qualitative variables related to the economic potential, technical aspects, contextual viability and socioeconomic factors in the analysis. The results show that, despite the increasing relevance of environmental and community issues in mining development, the business potential and the economic/technical aspects are the main factors in the early selection of a site to advance in exploration and development activities. Both variables represented around 40% and 37% of weights in the final selection, respectively. In contrast, contextual viability and local socioeconomic impacts only accounted for the remaining 23%. This study also shows that the inclusion of experts with different backgrounds in the process enriches the analysis and does not significantly distort the final outcome of the prioritization. Finally, the relevance of using MCDM tools when assessing the attractiveness of mine sites for their development is highlighted, particularly when public funds for subsequent exploration activities are committed.

     

Variations of young oceanic intraplate seismic energy released with relation to lithosphere age. Implications in the East Pacific Rise and its convergence with the Rivera subduction zone

Marine Geophysical Research - We have estimated empirical relationships that associate the age of the Pacific oceanic lithosphere with the seismic energy released on both flanks of the East Pacific...     

Habitat continuity effects on gradients of fish biomass across marine protected area boundaries

Marine protected areas (MPAs) could be useful as fisheries management tools for the exportation of pelagic eggs, larvae and adult fish. A decreasing gradient of fish biomass across MPAs boundary may indicate export. We determine whether gradients of decreasing biomass of fish assemblage occurred in Tabarca Marine Reserve over two habitats with different continuity across the boundaries, to test if the patchy nature of the marine environment might act as a barrier for the fish export. In general, significant decreasing gradients in total fish biomass and biomass of some species were observed on P. oceanica and rocky substrates, independently of their different continuity through the reserve boundaries. Changes in the multivariate structure of the fish assemblage were correlated with the distance from integral reserve. All of these results support the hypothesis that the exportation of adult fish from Tabarca Marine Reserve occurs, and this process may influence the surrounding fished areas.