Post-eruptive lahars at Kali Putih following the 2010 eruption of Merapi volcano,Indonesia: occurrences and impacts

Following the 2010 VEI 4 eruption of Merapi volcano, more than 250 lahars were triggered during two rainy seasons from October 2010 to March 2012. This high number of post-eruption lahars mainly occurred in the Kali (valley) Putih watershed and was mostly associated with high-magnitude rainstorms. A lahar occurring on January 8, 2011, caused significant damage to homes in several communities, bridges, sabo dams, and agricultural crops. The aims of this contribution are to document the impacts of lahars on the Kali Putih watershed and specifically (1) to analyze the lahar frequency during the period of 1969–2012 on an inter-annual and intra-annual basis and to determine the link between the volume of tephra and the frequency of lahars; (2) to detail the lahar trajectory and channel evolution following the January 8th lahar; (3) to map the spatial distribution of the thickness and geomorphic effects of the lahar deposit; and (4) to determine the impacts of the lahar on the infrastructure (sabo dams and roads) and settlements in the distal area of the volcano. The Kali Putih watershed has experienced 62 lahars, which represent 22% of all lahars triggered on 17 rivers at Merapi between 2010 and 2012. The main geomorphic impacts are: (1) excessive sedimentation in valleys, settlements and agricultural areas; (2) undercutting of the river banks by as much as 50 m, accompanied by channel widening; and (3) abrupt changes in the river channel direction in the distal area (15–20 km downstream of the volcano). About 19 sabo dams were damaged, and 3 were totally destroyed. Over 307 houses were damaged, and the National Road Yogyakarta–Semarang was regularly cut (18 times during approximately 25 days). Although the sabo dams on Kali Putih were originally constructed to protect distal areas from lahar damage, they had little effect on the 2010–2012 rain-triggered lahars. The underlying design of those dams along this river is one of the main reasons for the major destruction in this sector of the volcano’s lower slope. The catch basin capacity of the sabo dam was only 1.75?×?10⁶ m³, whereas the total volume of the 2010–2011 lahars exceeded 5?×?10⁶ m³. In order to prepare for future lahars, the government has invested in significant mitigation measures, ranging from structural approaches (e.g., building new sabo dams and developing an early warning system) to non-structural approaches (e.g., contingency and preparedness planning and hazard education).

     

A four-parameter,two degree-of-freedom block-spring model: Effect of the driver velocity

We analyze the effect of tectonic plate velocities in the earthquake pattern using a simple mass-spring model of the Burridge and Knopoff type with two blocks and a velocity-weakening friction law. Previous versions of the two-block model assume a steady driver during slip events (limit of zero driver velocity), which, in some cases makes necessary the introduction of artificial parameters to start the numerical integration of the equations of motion at impending slip of any block. Still maintaining the condition of zero driver velocity during slip, we shall introduce a procedure to start the numerical integration without introducing artificial parameters and this will be done by using a linearized version of the equations of motion valid for small velocities and considering nonzero driver velocity. We also introduce a four parameter model in which the driver velocity enters the equations during the whole simulation, and analyze the effect of the new parameter, the driver velocity, in the displacement and time patterns of blocks motion, directly related to earthquake statistics such as coseismic slips and average repeat times.     

Temporal variations in bedload transport rates associated with the migration of bedforms

Temporal variations in bedload transport rates that occur at a variety of timescales, even under steady flow conditions, are accepted as an inherent component of the bedload transport process. Rarely, however, has the cause of such variations been explained clearly. We consider three data sets, obtained from laboratory experiments, that refer to measurements of bedload transport made with continuously recording bedload traps. Each data set is characterized by a predominant low-frequency oscillation, on which additional higher-frequency oscillations generally are superimposed. The period of these oscillations, as isolated through the use of spectral analysts, ranged between 0·47 and 168 minutes, and was associated unequivocally with the migration of bedforms such as ripples, dunes, and bars. The extent to which such oscillatory behaviour may be recognized in a data set depends on the duration of sampling and the length of the sampling time, with respect to the period of a given bedform. Several theoretical probability distribution functions have been developed to describe the frequency distributions of (relative) bedload transport rates that are associated with the migration of bedforms (Einstein, 1937b; Hamamori, 1962; Carey and Hubbell, 1986). These distribution functions were derived without reference to a sampling interval. We present a modification of Hamamori's (1962) probability distribution function, generated by Monte Carlo simulation, which permits one to specify the sampling interval, in relation to the length of a bedform. Comparisons between the simulated and observed frequency distributions, that were undertaken on the basis of the data described herein, are good (significant at the 90 per cent confidence level). Finally, the implications that temporal variability, which is associated with the migration of bedforms, have for the accurate determination of bedload transport rates are considered.     

Geochemical Indicators of Intrinsic Bioremediation

A detailed field investigation has been completed at a gasoline-contaminated aquifer near Rocky Point, NC, to examine possible indicators of intrinsic bioremediation and identify factors that may significantly influence the rate and extent of bioremediation. The dissolved plume of benzene, toluene, ethylbenzene, and xylene (BTEX) in ground water is naturally degrading. Toluene and o-xylene are most rapidly degraded followed by m-, p-xylene, and benzene. Ethylbenzene appears to degrade very slowly under anaerobic conditions present in the center of the plume. The rate and extent of biodegradation appears to be strongly influenced by the type and quantity of electron acceptors present in the aquifer. At the upgradient edge of the plume, nitrate, ferric iron, and oxygen are used as terminal electron acceptors during hydrocarbon biodegradation. The equivalent of 40 to 50 mg/I of hydrocarbon is degraded based on the increase in dissolved CO₂ relative to background ground water. Immediately downgradient of the source area, sulfate and iron are the dominant electron acceptors. Toluene and o-xylene are rapidly removed in this region. Once the available oxygen, nitrate, and sulfate are consumed, biodegradation is limited and appears to be controlled by mixing and aerobic biodegradation at the plume fringes.     

Ocean storage of nuclear wastes? Experiences from the Russian Arctic

An international demonstration (RD&D) project for ocean storage of radioactive wastes should be proposed, to study the feasibility of the concept of ocean storage of nuclear waste. This international project should utilize the scientific, engineering and technical capabilities of selected universities, oceanographic institutions, NGOs and industries. This project would need to be an independent (non-governmental) study, utilizing the capabilities of selected universities, oceanographic institutions, environmental NGOs (Non-Governmental Organizations) and industries. Scientists and engineers first need to conduct an engineering, environmental, and economic feasibility study of the concept. The goal of the project would be to determine if ocean-based storage reduced the risks to the environment and public health to a greater degree than land-based storage. This would require comparing the risks and factors involved and making the data and information available to anyone, anywhere, anytime on the internet. The mere presence of an investigation of the ocean storage option could facilitate scientific and engineering competition between the two options, could subsequently reduce environmental and public risks and provide better protection and cost benefits in the system utilized. One of the primary concerns of the scientific community would be related to the sensitivity and precision of the monitoring of individaul containers on the ocean bottom. An advantage of the land-based option is that if there is a release, its presence could be detected at very low levels and be contained in the storage facility. On the ocean bottom, a release from a container might not be easily detected due to dispersion. Therefore the containment system would have to be a system within a system with monitoring between the two providing greater protection. Ocean storage may have greater technical and political hurdles than land-based options, but it may provide greater protection over time, because it negates the threat of terrorism, it therefore merit further study. In the future if the use of nuclear energy and nuclear wastes increases, the global society could benefit from this international project, because it could reduce environmental and public health risks and promote energy independence.     

Controls on sediment export from the Waipaoa River basin, New Zealand

A stability model of drainage basin mass balance is used to interpret historic and prehistoric patterns of sediment production, storage and output from the Waipaoa River basin, New Zealand and assess the sensitivity of basin sediment yield to land use change in the historic period. Climate and vegetation cover changed during the late Holocene, but the drainage basin mass balance system was stable before the basin was deforested by European colonists in the late 19th and early 20th centuries. In this meso‐scale dispersal system sediment sources and sinks are closely linked, and before that time there was also little variability in the rate of terrigenous mass accumulation on the adjacent continental shelf. However, despite strong first‐order geologic controls on erosion and extensive alluvial storage, sediment delivery to the continental shelf is sensitive and highly responsive to historic hillslope destabilization driven by land use change. Alluvial buffering can mask the effects of variations in sediment production within a basin on sediment yield at the outlet, but this is most likely to occur in basins where alluvial storage is large relative to yield and where the residence time of alluvial sediment is long relative to the time scale of environmental change. At present, neither situation applies to the Waipaoa River basin. Thus, the strength of the contemporary depositional signal may not only be due to the intensity of the erosion processes involved, but also to the fact that land use change in the historic period destabilized the drainage basin mass balance system.     

An experimental study of sapped drainage network development

In a basin developed on a stream table, effluent subsurface flow supported a channel network that evolved by a combination of headward growth, lateral widening and divide decay. The area occupied by the developing network increased with time. Circularity was used to characterize network evolution which occurred in three phases (initiation, extension and abstraction). Basin sediment discharge declined exponentially with time. Pronounced quasi-cyclic variability was superimposed upon this general trend. Some of the variability was directly linked to changes in the amount of sediment supplied to the channel. The variation of mean network sediment yield (mean sediment discharge scaled by network area) with time adequately described the general decline in sediment discharge as the network evolved.     

Small scale solar magnetic fields and their relation to various solar activity indices

Correlation studies between various solar activity indices and a long time series of annual sums of the maximum value of solar magnetic field intensity, observed for each group of sunspots during each passage of it over the visible solar hemisphere, have pointed out a couple of interesting points. First, the faculae have a significant contribution to the numerical representation of the small scale solar magnetic coefficients and low standard errors of estimation to the above mentioned maximum values of the solar magnetic field. These properties give to the area index an important physical meaning which is a first approximation to the small scale solar magnetic fields expressed by the above-mentioned maximum values of it. Finally, the main point which comes out is that long term studies of the solar magnetic fields, especially extrapolated studies to the past, could be supported by photospheric indices of the solar activity. This paper constitutes the expanded version of a report presented to theIAU Symposium No. 102 ‘Solar and Stellar magnetic fields: Origins and coronal effects’, held in Zürich 2–6 August, 1982.