Estimation of regional material yield from coastal landslides based on historical digital terrain modelling

High‐resolution historical (1942) and recent (1994) digital terrain models were derived from aerial photographs along the Big Sur coastline in central California to measure the long‐term volume of material that enters the nearshore environment. During the 52‐year measurement time period, an average of 21 000 ± 3100 m³ km^?1 a^?1 of material was eroded from nine study sections distributed along the coast, with a low yield of 1000 ± 240 m³ km^?1 a^?1 and a high of 46 700 ± 7300 m³ km^?1 a^?1. The results compare well with known volumes from several deep‐seated landslides in the area and suggest that the processes by which material is delivered to the coast are episodic in nature. In addition, a number of parameters are investigated to determine what influences the substantial variation in yield along the coast. It is found that the magnitude of regional coastal landslide sediment yield is primarily related to the physical strength of the slope‐forming material. Coastal Highway 1 runs along the lower portion of the slope along this stretch of coastline, and winter storms frequently damage the highway. The California Department of Transportation is responsible for maintaining this scenic highway while minimizing the impacts to the coastal ecosystems that are part of the Monterey Bay National Marine Sanctuary. This study provides environmental managers with critical background data on the volumes of material that historically enter the nearshore from landslides, as well as demonstrating the application of deriving historical digital terrain data to model landscape evolution. Published in 2005 by John Wiley & Sons, Ltd.     

A simulation model of morphological,vegetation and sediment changes in ephemeral streams

A model to simulate channel changes in ephemeral river channels and to test the effects of hydrological changes due to climate change and[sol ]or land use change was developed under the auspices of the EU funded MEDALUS programme (Mediterranean Desertification and Land Use). The model, CHANGISM (Channel Change GIS Simulation Model), is designed to simulate the effect of channel flow events and of climate conditions on morphology, sediment and vegetation, through sequences of events and conditions, over periods of up to several decades. The modelling is based on cellular automata but with calculations for water and sediment continuity. Process rules have both deterministic and stochastic elements. An important feature of the model is that it incorporates feedback elements between each event. The main aim of the model is to indicate the likely outcomes of events and combinations of conditions. It is linked to GIS for both input and output. The modelling is based on a channel reach and state is input as GIS layers of morphology (DEM), sediment and vegetation cover and state. Other initial conditions of soil moisture, groundwater level, and overall gradient are input. Parameters for processes are read from tables and can be easily changed for successive runs of the model. The bases for decisions on process specifications are discussed in this paper. Initial tests of the operation and sensitivity of the model were made on idealized reaches. The model was then tested using data from monitored sites in SE Spain. Simulations using clearwater flow worked well but initial simulations using events with sediment loads showed some tendency for excess deposition. Further tests and modifications are taking place. Overall, the model is one of the most sophisticated that simulates the interaction of flows with sediment and vegetation and the outcomes in terms of erosion, deposition, morphology, sediment cover, vegetation cover and plant survival over periods of up to 30 years for the scale of a channel reach. Copyright © 2005 John Wiley & Sons, Ltd.     

Channel adjustments of the lower Trinity River,Texas, downstream of Livingston Dam

Channel cross‐sectional changes since construction of Livingston Dam and Lake Livingston in 1968 were studied in the lower Trinity River, Texas, to test theoretical models of channel adjustment, and to determine controls on the spatial extent of channel response. High and average flows were not significantly modified by the dam, but sediment transport is greatly reduced. The study is treated as an opportunistic experiment to examine the effects of a reduction in sediment supply when discharge regime is unchanged. Channel scour is evident for about 60 km downstream, and the general phenomena of incision, widening, coarsening of channel sediment and a decrease in channel slope are successfully predicted, in a qualitative sense, by standard models of channel response. However, there is no consistent channel response within this reach, as various qualitatively different combinations of increases, decreases or no change in width, depth, slope and roughness occur. These multiple modes of adjustment are predicted by the unstable hydraulic geometry model. Between about 60 km and the Trinity delta 175 km downstream of the dam, no morphological response to the dam is observed. Rather than a diminution of the dam's effects on fluvial processes, this is due to a fundamental change in controls of the fluvial system. The downstream end of the scour zone corresponds to the upstream extent of channel response to Holocene sea level rise. Beyond 60 km downstream, the Trinity River is characterized by extensive sediment storage and reduced conveyance capacity, so that even after dam construction sediment supply still exceeds transport capacity. The channel bed of much of this reach is near or below sea level, so that sea level rise and backwater effects from the estuary are more important controls on the fluvial system than upstream inputs. Copyright © 2005 John Wiley & Sons, Ltd.     

Energetic Particles Observed in the CUSP Region During a Storm Recovery Phase

In this paper we report energetic ion behavior and its composition variations observed by the Cluster/RAPID instrument when the spacecraft was travelling in the high latitude magnetospheric boundary region on the day of the 31 March, 2001, strongest magnetic storm in the past 50 years. The Dst index reached −360 nT at about 09:00 UT. During its early recovery phase, large amounts of oxygen and helium ions were observed; the ratio of oxygen to hydrogen in the RAPID energy range reached as high as 250%, which suggests that the observed energetic particles might be of magnetospheric origin. The observations further show that enhanced energetic electron fluxes are confined in a very narrow region, while protons have occupied a larger region, and heavy ions have been observed in an even larger region. The flux of energetic electrons show a slight enhancement in a region where the magnetic field magnitude is around zero. These observed energetic ions could be quasi-trapped by the current sheet in the stagnation region of the cusp.     

Seismic signature of a phreatic explosion: hydrofracturing damage at Karthala volcano, Grande Comore Island, Indian Ocean

Karthala volcano is a basaltic shield volcano with an active hydrothermal system that forms the southern two-thirds of the Grande Comore Island, off the east coat of Africa, northwest of Madagascar. Since the start of volcano monitoring by the local volcano observatory in 1988, the July 11th, 1991 phreatic eruption was the first volcanic event seismically recorded on this volcano, and a rare example of a monitored basaltic shield. From 1991 to 1995 the VT locations, 0.5<M_l<4.3, show a crack shaped pattern (3 km long, 1 km wide) within the summit caldera extending at depth from –2 km to +2 km relative to sea level. This N-S elongated pattern coincides with the direction of the regional maximum horizontal stress as deduced from regional focal mechanism solutions. This brittle signature of the damage associated with the 1991 phreatic eruption is a typical pattern of the seismicity induced by controlled fluid injections such as those applied at geothermal fields, in oil and gas recovery, or for stress measurements. It suggests the 1991 phreatic eruption was driven by hydraulic fracturing induced by forced fluid flow. We propose that the extremely high LP and VT seismicity rates, relative to other effusive volcanoes, during the climax of the 1991 phreatic explosion, are due to the activation of the whole hydrothermal system, as roughly sized by the distribution of VT hypocenters. The seismicity rate in 1995 was still higher than the pre-eruption seismicity rate, and disagrees with the time pattern of thermo-elastic stress readjustment induced by single magma intrusions at basaltic volcanoes. We propose that it corresponds to the still ongoing relaxation of pressure heterogeneity within the hydrothermal system as suggested by the few LP events that still occurred in 1995.Editorial responsibility: H Shinohara     

Caldera formation by magma withdrawal from a reservoir beneath a volcanic edifice

Caldera formation has been explained by magma withdrawal from a crustal reservoir, but little is known about the conditions that lead to the critical reservoir pressure for collapse. During an eruption, the reservoir pressure is constrained to lie within a finite range: it cannot exceed the threshold value for eruption, and cannot decrease below another threshold value such that feeder dykes get shut by the confining pressure, which stops the eruption. For caldera collapse to occur, the critical reservoir pressure for roof failure must therefore be within this operating range. We use an analytical elastic model to evaluate the changes of reservoir pressure that are required for failure of roof rocks above the reservoir with and without a volcanic edifice at Earth's surface. With no edifice at Earth's surface, faulting in the roof region can only occur in the initial phase of reservoir inflation and affects a very small part of the focal area. Such conditions do not allow caldera collapse. With a volcanic edifice, large tensile stresses develop in the roof region, whose magnitude increase as the reservoir deflates during an eruption. The edifice size must exceed a threshold value for failure of the roof region before the end of eruption. The largest tensile stresses are reached at Earth's surface, indicating that faulting starts there. Failure affects an area whose horizontal dimensions depend on edifice and chamber dimensions. For small and deep reservoirs, failure conditions cannot be achieved even if the edifice is very large. Quantitative predictions are consistent with observations on a number of volcanoes.     

The influence of particle composition on thorium scavenging in the NE Atlantic ocean (POMME experiment)

²³⁰Th, ²³²Th and ²³⁴Th were analyzed in sinking particles collected by moored and drifting sediment traps in the NE Atlantic Ocean (POMME experiment) in order to constrain the phase(s) carrying Th isotopes in the water column. It reveals a contrasted behaviour between ²³⁴Th and ²³⁰Th. ²³⁴Th is correlated to the particulate organic carbon suggesting that it is primarily scavenged by organic compounds in the surface waters. ²³⁰Th_xs is correlated with Mn, Ba and the lithogenic fraction that are enriched in small suspended particles and incorporated in the sinking particulate flux throughout the water column. The lack of correlation between ²³⁰Th_xs and CaCO₃ or biogenic silica (bSi) indicates that CaCO₃ and bSi are not responsible for ²³⁰Th scavenging in the deep waters of this oceanic region. ²³⁰Th is generally correlated with the lithogenic content of the trapped material but this correlation disappears in winter during strong atmospheric dust inputs suggesting that lithogenic matter is not directly responsible for ²³⁰Th scavenging in the deep waters or that sufficient time is required to achieve particle–solution equilibration. MnO₂ could be the prevalent ²³⁰Th_xs-bearing phase. The narrow range of K_{d_MnO2}^Th obtained for very contrasted oceanic environments supports a global control of ²³⁰Th_xs scavenging by MnO₂ and raises the possibility that the ²³⁰Th–²³¹Pa fractionation is controlled by the amount of colloidal MnO₂ in seawater.     

Evidence from the Northeastern Atlantic basin for variability in the rate of the meridional overturning circulation through the last deglaciation

A first study from the subtropical western Atlantic, using ²³¹Pa/²³⁰Th ratios as a kinematic proxy for deep water circulation, provided compelling evidence for a strong link between climate and the rate of Meridional Overturning Circulation (MOC) over the last deglaciation. However, these results warrant confirmation from additional locations and water depths because the interpretation of the sedimentary ²³¹Pa/²³⁰Th ratio in terms of circulation vigor can be biased by variations in particle flux and composition. We have measured ²³¹Pa/²³⁰Th in a core from the Iberian margin, in the Northeastern Atlantic basin, and have compared these new results to the data from the western Atlantic basin. We find that the reduction in the circulation during H1 and YD and the subsequent increases first recognized in the sediment deposited on Bermuda Rise are also evident in the eastern basin, in a totally different sedimentary regime, confirming that sedimentary ²³¹Pa/²³⁰Th ratios record basin-wide changes in deep water circulation. However, some differences between the eastern and western records are also recognized, providing preliminary evidence to differentiate between renewal rates in the two North Atlantic basins and between shallower and deeper overturning. Our results suggest the possible existence of two sources of Glacial North Atlantic Intermediate Deep Water (GNAIW), one in the south Labrador Sea and another west of Rockall Plateau. Both sources contributed to the meridional overturning but the two had different sensitivity to meltwater from the Laurentide and the Fennoscandian ice sheets during the deglaciation. These results indicate that additional information on the geometry and strength of the ventilation of the deep Atlantic can be obtained by contrasting the evolution of sediment ²³¹Pa/²³⁰Th in different sections of the Atlantic Ocean.     

Viscosity of peridotite liquid up to 13 GPa: Implications for magma ocean viscosities

The viscosity of synthetic peridotite liquid has been investigated at high pressures using in-situ falling sphere viscometry by combining a multi-anvil technique with synchrotron radiation. We used a newly designed capsule containing a small recessed reservoir outside of the hot spot of the heater, in which a viscosity marker sphere is embedded in a forsterite + enstatite mixture having a higher solidus temperature than the peridotite. This experimental setup prevents spheres from falling before a stable temperature above the liquidus is established and thus avoids difficulties in evaluating viscosities from velocities of spheres falling through a partially molten sample.

Experiments have been performed between 2.8 and 13 GPa at temperatures ranging from 2043 to 2523 K. Measured viscosities range from 0.019 (± 0.004) to 0.13 (± 0.02) Pa s. At constant temperature, viscosity increases with increasing pressure up to  8.5 GPa but then decreases between  8.5 and 13 GPa. The change in the pressure dependence of viscosity is likely associated with structural changes of the liquid that occur upon compression. By combining our results with recently published 0.1 MPa peridotite liquid viscosities [D.B. Dingwell, C. Courtial, D. Giordano, A. Nichols, Viscosity of peridotite liquid, Earth Planet. Sci. Lett. 226 (2004) 127–138.], the experimental data can be described by a non-Arrhenian, empirical Vogel-Fulcher-Tamman equation, which has been modified by adding a term to account for the observed pressure dependence of viscosity. This equation reproduces measured viscosities to within 0.08 log₁₀-units on average. We use this model to calculate viscosities of a peridotitic magma ocean along a liquid adiabat to a depth of  400 km and discuss possible effects on viscosity at greater pressures and temperatures than experimentally investigated.