A new theoretical method for analyzing confined dry granular flows

A granular body may deform in a continuous fashion such that the solid particles remain in close contact. Previous research works have always used the frictional Coulomb-like continuum treatment for analyzing granular bodies. However, this approach is only applicable for quasi-static conditions and cannot capture the complicated granular contact behavior of solid particles inside a failing granular body. This paper applies a revised Savage–Hutter equation to model granular flows moving down a confined, sloping channel. The Coulomb contact friction law is modified to consider the effect of the shear rate inside a granular body. This new method also considers the confinement effect of a sloping channel on granular flow mobility. The derived depth-averaged equations of motion bear a resemblance to nonlinear shallow-water wave equations. Results computed using the derived equations are compared with measurements from flume model tests, and consistency is found between the two.     

Interpretative matrices approach to ranking lake sub-basin pollution potential: an applied study in Brazil

In 2005 a Cylindrospermopsis raciborskii bloom occurred in the Rio Verde Lake Basin (Brazil). To address this concern, a field analysis was performed to measure physicochemical variables and flows in 14 sub-basins, between 2008 and 2009. Measurements of mean total P (0.039 mg/L ± 0.018 mean SD), mean total Kjeldahl N (0.260 mg/L ± 0.226 mean SD), and mean BOD (1.2 mg/L ± 0.4 mean SD) concentrations were low in most streams, while COD reached a high of 27.1 mg/L (±4.9 mean SD). One tributary was responsible for 85 % of TP load, 77.1 % of TKN load, 78 % of t-BOD load, and 79 % of t-COD load. These concentrations and loads were used to develop the pollution potential assessment matrix (2PAM), which considered three different perspectives: stream water quality, reservoir ecosystem equilibrium, and sub-basin management. Each factor (TP, TKN, BOD and COD) was weighted based on concentration, total load and unit-area load. Pollution potential differed depending on which perspective was considered. The matrix developed, 2PAM, provides a new way to analyze concentrations and loads, enabling basin managers to prioritize action plans according to desired use within the basin.     

Zircon evidence for a ~200 k.y. supereruption-related thermal flare-up in the Miocene southern Black Mountains,western Arizona,USA

The Silver Creek caldera (southern Black Mountains, western Arizona) is the source of the 18.8 Ma, >700 km³ Peach Spring Tuff (PST) supereruption, the largest eruption generated in the Colorado River Extensional Corridor (CREC) of the southwestern United States. Within and immediately surrounding the caldera is a sequence of volcanics and intrusions ranging in age from ~19 to 17 Ma. These units offer a record of magmatic processes prior to, during, and immediately following the PST eruption. To investigate the thermal evolution of the magmatic center that produced the PST, we applied a combination of Ti-in-zircon thermometry, zircon saturation thermometry, and high-precision U–Pb CA–TIMS zircon dating to representative pre- and post-supereruption volcanic and intrusive units from the caldera and its environs. Similar to intracaldera PST zircons, zircons from a pre-PST trachytic lava (19 Ma) and a post-PST caldera intrusion (18.8 Ma) yield exceptionally high-Ti concentrations (most >20 ppm, some up to nearly 60 ppm), corresponding to calculated temperatures that exceed 900 °C. In these units, Ti-in-zircon temperatures typically surpass zircon saturation temperatures (ZSTs), suggesting the entrainment of zircon that had grown in hotter environments within the magmatic system. Titanium concentrations in younger volcanic and intrusive units (~18.7–17.5 Ma) decline through time, corresponding to an average cooling rate of 10^?3.5 °C/year. The ~200 k.y. thermal peak evident at Silver Creek caldera is spatially limited: elsewhere in the Miocene record of the northern CREC, Ti-in-zircon concentrations and ZSTs are much lower, suggesting that felsic magmas were generally substantially cooler.     

Building zoned ignimbrites by recycling silicic cumulates: insight from the 1,000 km<Superscript>3</Superscript> Carpenter Ridge Tuff,CO

The ~1,000 km³ Carpenter Ridge Tuff (CRT), erupted at 27.55 Ma during the mid-tertiary ignimbrite flare-up in the western USA, is among the largest known strongly zoned ash-flow tuffs. It consists primarily of densely welded crystal-poor rhyolite with a pronounced, highly evolved chemical signature (high Rb/Sr, low Ba, Zr, Eu), but thickly ponded intracaldera CRT is capped by a more crystal-rich, less silicic facies. In the outflow ignimbrite, this upper zone is defined mainly by densely welded crystal-rich juvenile clasts of trachydacite composition, with higher Fe–Ti oxide temperatures, and is characterized by extremely high Ba (to 7,500 ppm), Zr, Sr, and positive Eu anomalies. Rare mafic clasts (51–53 wt% SiO₂) with Ba contents to 4,000–5,000 ppm and positive Eu anomalies are also present. Much of the major and trace-element variations in the CRT juvenile clasts can be reproduced via in situ differentiation by interstitial melt extraction from a crystal-rich, upper-crustal mush zone, with the trachydacite, crystal-rich clasts representing the remobilized crystal cumulate left behind by the melt extraction process. Late recharge events, represented by the rare mafic clasts and high-Al amphiboles in some samples, mixed in with parts of the crystal cumulate and generated additional scatter in the whole-rock data. Recharge was important in thermally remobilizing the silicic crystal cumulate by partially melting the near-solidus phases, as supported by: (1) ubiquitous wormy/sieve textures and reverse zoning patterns in feldspars and biotites, (2) absence of quartz in this very silicic unit stored at depths of >4–5 km, and (3) heterogeneous melt compositions in the trachydacite fiamme and mafic clasts, particularly in Ba, indicating local enrichment of this element due mostly to sanidine and biotite melting. The injection of hot, juvenile magma into the upper-crustal cumulate also imparted the observed thermal gradient to the deposits and the mixing overprint that partly masks the in situ differentiation process. The CRT provides a particularly clear perspective on processes of in situ crystal-liquid separation into a lower crystal-rich zone and an upper eruptible cap, which appears common in incrementally built upper-crustal magma reservoirs of high-flux magmatic provinces.     

Matrix‐rich and associated matrix‐poor sandstones: Avulsion splays in slope and basin‐floor strata

A common facies observed in deep‐water slope and especially basin‐floor rocks of the Neoproterozoic Windermere Supergroup (British Columbia, Canada) is structureless, coarse‐tail graded, medium‐grained to coarse‐grained sandstone with from 30% to >50% mud matrix content (i.e. matrix‐rich). Bed contacts are commonly sharp, flat and loaded. Matrix‐rich sandstone beds typically form laterally continuous units that are up to several metres thick and several tens to hundreds of metres wide, and commonly adjacent to units of comparatively matrix‐poor, scour‐based sandstone beds with large tabular mudstone and sandstone clasts. Matrix‐rich units are common in proximal basin‐floor (Upper Kaza Group) deposits, but occur also in more distal basin‐floor (Middle Kaza Group) and slope (Isaac Formation) deposits. Regardless of stratigraphic setting, matrix‐rich units typically are directly and abruptly overlain by architectural elements comprising matrix‐poor coarse sandstone (i.e. channels and splays). Despite a number of similarities with previously described matrix‐rich beds in the literature, for example slurry beds, linked debrites and co‐genetic turbidites, a number of important differences exist, including the stratal make‐up of individual beds (for example, the lack of a clean sandstone turbidite base) and their stratigraphic occurrence (present throughout base of slope and basin‐floor strata, but most common in proximal lobe deposits) and accordingly suggest a different mode of emplacement. The matrix‐rich, poorly sorted nature of the beds and the abundance and size of tabular clasts in laterally equivalent sandstones imply intense upstream scouring, most probably related to significant erosion by an energetic plane‐wall jet or within a submerged hydraulic jump. Rapid energy loss coupled with rapid charging of the flow with fine‐grained sediment probably changed the rheology of the flow and promoted deposition along the margins of the jet. Moreover, these distinctive matrix‐rich strata are interpreted to represent the energetic initiation of the local sedimentary system, most probably caused by a local upflow avulsion.     

Sea-level and reef accretion history of Marine Oxygen Isotope Stage 7 and late Stage 5 based on age and facies of submerged late Pleistocene reefs,Oahu, Hawaii

In situ Pleistocene reefs form a gently sloping nearshore terrace around the island of Oahu. TIMS Th–U ages of in situ corals indicate that most of the terrace is composed of reefal limestones correlating to Marine Oxygen Isotope Stage 7 (MIS 7, ~ 190–245 ka). The position of the in situ MIS 7 reef complex indicates that it formed during periods when local sea level was ~ 9 to 20 m below present sea level. Its extensiveness and geomorphic prominence as well as a paucity of emergent in situ MIS 7 reef-framework deposits on Oahu suggest that much of MIS 7 was characterized by regional sea levels below present. Later accretion along the seaward front of the terrace occurred during the latter part of MIS 5 (i.e., MIS 5a–5d, ~ 76–113 ka). The position of the late MIS 5 reefal limestones is consistent with formation during a period when local sea level was below present. The extensiveness of the submerged Pleistocene reefs around Oahu compared to the relative dearth of Holocene accretion is due to the fact that Pleistocene reefs had both more time and more accommodation space available for accretion than their Holocene counterparts.     

Periodic solar EUV flux monitored near Venus

A detector sharing the orbital rate of Venus has a unique perspective on solar periodicities. Fourier analysis of the 8.6 year record of solar EUV output gathered by the Langmuir probe on Pioneer Venus Orbiter shows the influences of global oscillation modes located in the convective envelope and in the radiative interior. Seven of the eight lowest angular harmonic r-mode families are detected by their rotation rates which differ almost unmeasurably from ideal theoretical values. This determines a mean sidereal rotation rate for the envelope of 457.9 ± 2.0 nHz which corresponds to a period of 25.3 days. Many frequencies are aliased at ± 106 nHz by modulation from the lowest angular harmonic r-mode in the envelope. The rotation of this mode seems slightly retrograde, -1.5 ± 2.0 nHz, but small positive values are not excluded. We confirm that the rotation of the radiative interior, 381 nHz, is slower than the envelope by detecting g-mode frequencies for angular harmonics, 2 l 6, and a possible first detection of the rotation rate for the l = 1 case. Solar EUV lacks the sudden darkenings (dips) shown by visible irradiance; vortex cores in the photosphere and below are again suggested as a possible explanation.     

Implementing quality control on a random number stream to improve a stochastic weather generator

For decades, stochastic modellers have used computerized random number generators to produce random numeric sequences fitting a specified statistical distribution. Unfortunately, none of the random number generators we tested satisfactorily produced the target distribution. The result is generated distributions whose mean even diverges from the mean used to generate them, regardless of the length of run. Non‐uniform distributions from short sequences of random numbers are a major problem in stochastic climate generation, because truly uniform distributions are required to produce the intended climate parameter distributions. In order to ensure generation of a representative climate with the stochastic weather generator CLIGEN within a 30‐year run, we tested the climate output resulting from various random number generators. The resulting distributions of climate parameters showed significant departures from the target distributions in all cases. We traced this failure back to the uniform random number generators themselves. This paper proposes a quality control approach to select only those numbers that conform to the expected distribution being retained for subsequent use. The approach is based on goodness‐of‐fit analysis applied to the random numbers generated. Normally distributed deviates are further tested with confidence interval tests on their means and standard deviations. The positive effect of the new approach on the climate characteristics generated and the subsequent deterministic process‐based hydrology and soil erosion modelling are illustrated for four climatologically diverse sites. Copyright © 2007 John Wiley & Sons, Ltd.     

Improvement of distributed snowmelt energy balance modeling with MODIS‐based NDSI‐derived fractional snow‐covered area data

Describing the spatial variability of heterogeneous snowpacks at a watershed or mountain‐front scale is important for improvements in large‐scale snowmelt modelling. Snowmelt depletion curves, which relate fractional decreases in snow‐covered area (SCA) against normalized decreases in snow water equivalent (SWE), are a common approach to scale‐up snowmelt models. Unfortunately, the kinds of ground‐based observations that are used to develop depletion curves are expensive to gather and impractical for large areas. We describe an approach incorporating remotely sensed fractional SCA (FSCA) data with coinciding daily snowmelt SWE outputs during ablation to quantify the shape of a depletion curve. We joined melt estimates from the Utah Energy Balance Snow Accumulation and Melt Model (UEB) with FSCA data calculated from a normalized difference snow index snow algorithm using NASA's moderate resolution imaging spectroradiometer (MODIS) visible (0·545–0·565 µm) and shortwave infrared (1·628–1·652 µm) reflectance data. We tested the approach at three 500 m² study sites, one in central Idaho and the other two on the North Slope in the Alaskan arctic. The UEB‐MODIS‐derived depletion curves were evaluated against depletion curves derived from ground‐based snow surveys. Comparisons showed strong agreement between the independent estimates. Copyright © 2010 John Wiley & Sons, Ltd.     

Understanding geomorphic responses to environmental change: a 19 000‐year case study from semi‐arid central Texas,USA

Understanding the impact of past climatic changes on landscape stability is crucial in order to predict and mitigate the effects of future changes. However, in arid and semi‐arid environments, reconstructions are often hampered by a poor understanding of the relationship between sediment deposition and climate. We present here data from central Texas, a region that is sensitive to environmental change but has received relatively little attention. The study integrates a chronology of 29 optically stimulated luminescence (OSL) ages from six sedimentary sites in a range of depositional contexts with a 19 000 a climate record derived from pollen extracted from the adjacent Boriack Bog. By comparing the two records, we aimed to assess the relationship between climate change and geomorphic activity. Data show that extensive aeolian and colluvial deposition occurred during the mid to late Holocene, with sedimentation generally increasing during more arid phases. However, a number of depositional events on slopes were associated with moister episodes, and sediment was also deposited in a summit setting immediately after phases of increased precipitation. Linkages between climate change and geomorphic response were therefore complex, being controlled by both sediment supply and transport energy. The climatic fluctuations identified in the Boriack Bog record highlight the sensitivity of central Texas to environmental change, while disparities in conditions recorded between it and other palaeoenvironmental sites in the southern USA emphasise the need for further work in order to enhance understanding of landscape sedimentary response to climate. Copyright © 2010 John Wiley & Sons, Ltd.