A derived probability distribution approach to stormwater quality modeling

The closed-form analytical stormwater quality models are developed for simulating urban catchment pollutant buildup and washoff processes. By integrating the rainfall–runoff transformation with pollutant buildup and washoff functions, stormwater quality measures, such as the cumulative distribution functions (CDFs) of pollutant loads, the expected value of pollutant event mean concentrations (EMCs) and the average annual pollutant load can be derived. This paper presents methodologies and major procedures for the development of urban stormwater quality models based on derived probability distribution theory. In order to investigate the spatial variation in model parameters and its impact on stormwater pollutant buildup and washoff processes as well as pollutant loads to receiving waters, an extended form of the original rainfall–runoff transformation which is based on lumped runoff coefficient approach is proposed to differentiate runoff generation mechanisms between the impervious and pervious areas of the catchment. In addition, as a contrast to the aggregated pollutant buildup models formulated with a single lumped buildup parameter, the disaggregated form of the pollutant buildup model is proposed by introducing a number of physically-based parameters associated with pollutant buildup and washoff processes into the pollutant load models. The results from the case study indicate that analytical urban stormwater management model are capable of providing results in good agreement with the field measurements, and can be employed as alternatives to continuous simulation models in the evaluation of long-term stormwater quality measures.     

A discontinuous finite element baroclinic marine model on unstructured prismatic meshes

We describe the time discretization of a three-dimensional baroclinic finite element model for the hydrostatic Boussinesq equations based upon a discontinuous Galerkin finite element method. On one hand, the time marching algorithm is based on an efficient mode splitting. To ensure compatibility between the barotropic and baroclinic modes in the splitting algorithm, we introduce Lagrange multipliers in the discrete formulation. On the other hand, the use of implicit–explicit Runge–Kutta methods enables us to treat stiff linear operators implicitly, while the rest of the nonlinear dynamics is treated explicitly. By way of illustration, the time evolution of the flow over a tall isolated seamount on the sphere is simulated. The seamount height is 90% of the mean sea depth. Vortex shedding and Taylor caps are observed. The simulation compares well with results published by other authors.     

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.     

Isolation of 12 microsatellite markers following a pyrosequencing procedure and cross-priming in two invasive cryptic species, Alexandrium catenella (group IV) and A. tamarense (group III) (Dinophyceae)

Alexandrium catenella (group IV) and Alexandrium tamarense (group III) (Dinophyceae) are two cryptic invasive phytoplankton species belonging to the A. tamarense species complex. Their worldwide spread is favored by the human activities, transportation and climate change. In order to describe their diversity in the Mediterranean Sea and understand their settlements and maintenances in this area, new microsatellite markers were developed based on Thau lagoon (France) samples of A. catenella and A. tamarense strains. In this study twelve new microsatellite markers are proposed. Five of these microsatellite markers show amplifications on A. tamarense and ten on A. catenella. Three of these 12 microsatellite markers allowed amplifications on both cryptic species. Finally, the haplotypic diversity ranged from 0.000 to 0.791 and 0.000 to 0.942 for A. catenella and A. tamarense respectively.     

Contrasting origins of Cenozoic silicic volcanic rocks from the western Cordillera of the United States

Two fundamentally different types of silicic volcanic rocks formed during the Cenozoic of the western Cordillera of the United States. Large volumes of dacite and rhyolite, mostly ignimbrites, erupted in the Oligocene in what is now the Great Basin and contrast with rhyolites erupted along the Snake River Plain during the Late Cenozoic. The Great Basin dacites and rhyolites are generally calc-alkaline, magnesian, oxidized, wet, cool (<850°C), Sr-and Al-rich, and Fe-poor. These silicic rocks are interpreted to have been derived from mafic parent magmas generated by dehydration of oceanic lithosphere and melting in the mantle wedge above a subduction zone. Plagioclase fractionation was minimized by the high water fugacity and oxide precipitation was enhanced by high oxygen fugacity. This resulted in the formation of Si-, Al-, and Sr-rich differentiates with low Fe/Mg ratios, relatively low temperatures, and declining densities. Magma mixing, large proportions of crustal assimilation, and polybaric crystal fractionation were all important processes in generating this Oligocene suite. In contrast, most of the rhyolites of the Snake River Plain are alkaline to calc-alkaline, ferroan, reduced, dry, hot (830–1,050°C), Sr-and Al-poor, and Nb-and Fe-rich. They are part of a distinctly bimodal sequence with tholeiitic basalt. These characteristics were largely imposed by their derivation from parental basalt (with low fH₂O and low fO₂) which formed by partial melting in or above a mantle plume. The differences in intensive parameters caused early precipitation of plagioclase and retarded crystallization of Fe–Ti oxides. Fractionation led to higher density magmas and mid-crustal entrapment. Renewed intrusion of mafic magma caused partial melting of the intrusive complex. Varying degrees of partial melting, fractionation, and minor assimilation of older crust led to the array of rhyolite compositions. Only very small volumes of distinctive rhyolite were derived by fractional crystallization of Fe-rich intermediate magmas like those of the Craters of the Moon-Cedar Butte trend. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The physical volcanology of the 1600 eruption of Huaynaputina, southern Peru

Volcán Huaynaputina is a group of four vents located at 16°36'S, 70°51'W in southern Peru that produced one of the largest eruptions of historical times when ~11 km³ of magma was erupted during the period 19 February to 6 March 1600. The main eruptive vents are located at 4200 m within an erosion-modified amphitheater of a significantly older stratovolcano. The eruption proceeded in three stages. Stage I was an ~20-h sustained plinian eruption on 19-20 February that produced an extensive dacite pumice fall deposit (magma volume ~2.6 km³). Throughout medial-distal and distal parts of the dispersal area, a fine-grained plinian ashfall unit overlies the pumice fall deposit. This very widespread ash (magma volume ~6.2 km³) has been recognized in Antarctic ice cores. A short period of quiescence allowed local erosion of the uppermost stage-I deposits and was followed by renewed but intermittent explosive activity between 22 and 26 February (stage II). This activity resulted in intercalated pyroclastic flow and pumice fall deposits (~1 km³). The flow deposits are valley confined, whereas associated co-ignimbrite ash fall is found overlying the plinian ash deposit. Following another period of quiescence, vulcanian-type explosions of stage III commenced on 28 February and produced crudely bedded ash, lapilli, and bombs of dense dacite (~1 km³). Activity ceased on 6 March. Compositions erupted are predominantly high-K dacites with a phenocryst assemblage of plagioclase>hornblende>biotite>Fe-Ti oxides-apatite. Major elements are broadly similar in all three stages, but there are a few important differences. Stage-I pumice has less evolved glass compositions (~73% SiO₂), lower crystal contents (17-20%), lower density (1.0-1.3 g/cm³), and phase equilibria suggest higher temperature and volatile contents. Stage-II and stage-III juvenile clasts have more evolved glass (~76% SiO₂) compositions, higher crystal contents (25-35%), higher densities (up to 2.2 g/cm³), and lower temperature and volatile contents. All juvenile clasts show mineralogical evidence for thermal disequilibrium. Inflections on a plot of log thickness vs area^1/2 for the fall deposits suggest that the pumice fall and the plinian ash fall were dispersed under different conditions and may have been derived from different parts of the eruption column system. The ash appears to have been dispersed mainly from the uppermost parts of the umbrella cloud by upper-level winds, whereas the pumice fall may have been derived from the lower parts of the umbrella cloud and vertical part of the eruption column and transported by a lower-altitude wind field. Thickness half distances and clast half distances for the pumice fall deposit suggests a column neutral buoyancy height of 24-32 km and a total column height of 34-46 km. The estimated mass discharge rate for the ~20-h-long stage-I eruption is 2.4᎒⁸ kg/s and the volumetric discharge rate is ~3.6᎒⁵ m³/s. The pumice fall deposit has a dispersal index (Hildreth and Drake 1992) of 4.4, and its index of fragmentation is at least 89%, reflecting the dominant volume of fines produced. Of the 11 km³ total volume of dacite magma erupted in 1600, approximately 85% was evacuated during stage 1. The three main vents range in size from ~70 to ~400 m. Alignment of these vents and a late-stage dyke parallel to the NNW-SSE trend defined by older volcanics suggest that the eruption initiated along a fissure that developed along pre-existing weaknesses. During stage I this fissure evolved into a large flared vent, vent 2, with a diameter of approximately 400 m. This vent was active throughout stage II, at the end of which a dome was emplaced within it. During stage III this dome was eviscerated forming the youngest vent in the group, vent 3. A minor extra-amphitheater vent was produced during the final event of the eruptive sequence. Recharge may have induced magma to rise away from a deep zone of magma generation and storage. Subsequently, vesiculation in the rising magma batch, possibly enhanced by interaction with an ancient hydrothermal system, triggered and fueled the sustained Plinian eruption of stage I. A lower volatile content in the stage-II and stage-III magma led to transitional column behavior and pyroclastic flow generation in stage II. Continued magma uprise led to emplacement of a dome which was subsequently destroyed during stage III. No caldera collapse occurred because no shallow magma chamber developed beneath this volcano.     

Seismic evidence for deep low-velocity anomalies in the transition zone beneath West Antarctica

We present a three-dimensional (3D) SV-wave velocity model of the upper mantle beneath the Antarctic plate constrained by fundamental and higher mode Rayleigh waves recorded at regional distances. The good agreement between our results and previous surface wave studies in the uppermost 200 km of the mantle confirms that despite strong differences in data processing, modern surface wave tomographic techniques allow to produce consistent velocity models, even at regional scale. At greater depths the higher mode information present in our data set allows us to improve the resolution compared to previous regional surface wave studies in Antarctica that were all restricted to the analysis of the fundamental mode. This paper is therefore mostly devoted to the discussion of the deeper part of the model. Our seismic model displays broad domains of anomalously low seismic velocities in the asthenosphere. Moreover, we show that some of these broad, low-velocity regions can be more deeply rooted. The most remarkable new features of our model are vertical low-velocity structures extending from the asthenosphere down to the transition zone beneath the volcanic region of Marie Byrd Land, West Antarctica and a portion of the Pacific-Antarctic Ridge close to the Balleny Islands hotspot. A deep low-velocity anomaly may also exist beneath the Ross Sea hotspot. These vertical structures cannot be explained by vertical smearing of shallow seismic anomalies and synthetic tests show that they are compatible with a structure narrower than 200 km which would have been horizontally smoothed by the tomographic inversion. These deep low-velocity anomalies may favor the existence of several distinct mantle plumes, instead of a large single one, as the origin of volcanism in and around West Antarctica. These hypothetical deep plumes could feed large regions of low seismic velocities in the asthenosphere.