Magnitude–frequency relations in debris flow

Debris flow occurs frequently in mountainous regions in China. Because of the difficulties involved in predicting and catching live debris flows, an assessment of the potential for debris flow is crucial in hazard mitigation. Magnitude–frequency (MF) relations are of special significance in such assessments. In previous studies, MF relations have been inferred by analyzing environmental factors and historical records and using empirical relations. This paper is concerned with the derivation of MF relations at regional and valley scales, using a large database of statistics. At the regional scale, it is represented by the distribution of the valley area, because the area is often taken to indicate the potential magnitude of debris flow. Statistics from over 5,000 debris flow valleys in various provinces in China show that a power law holds for the distribution, i.e., p(A) ∼ A ⁻ⁿ , where p(A) is the percentage of valleys with area A and n varies with region and thus describes regional differences. At the valley scale, a case study focusing on Jiangjia Gully (JJG) was conducted, and the MF relations derived from it were expressed by the distributions of discharge and runoff (i.e., the total volume) of living debris flows observed over the last 40 years. The distributions can be expressed as exponential functions where the exponents vary with the events. These MF relations provide not only a potential quantitative reference for practical purposes but also hint at the intrinsic properties of the debris flow.     

A compilation of Australian heat‐flow measurements

A map is presented based on all known Australian heat‐flow estimates, including five new ones. A second map, based on the first, also is presented, but excludes any determinations judged for any reason to be unreliable. The data show that heat flow over large areas of the continent is effectively uniform to within 0.5 heat flow units (i.e. to within 20 mW m^‐2), with perhaps three major regional heat‐flow provinces being defined in western, central, and eastern Australia.     

Heat flow–heat production relationship not found: what drives heat flow variability of the Western Canadian foreland basin?

International Journal of Earth Sciences - Heat flow high −80 ± 10 mW/m2 in the northern western parts of the Western Canadian foreland basin is in large...     

The NGC 1023 galaxy group: An anti-hubble flow?

We discuss recently published data indicating that the nearby galaxy group NGC 1023 includes an inner, virialized, quasi-stationary component together with an outer component comprising a flow of dwarf galaxies falling toward the center of the system. The inner component is similar to the Local Group of galaxies, but the Local Group is surrounded by a receding set of dwarf galaxies forming the local Hubble flow, rather than a system of approaching dwarfs. This clear difference in the structures of these two systems, which are very similar in other respects, may be associated with the dark energy in which they are immersed. Self-gravity dominates in the inner component of the Local Group, while the anti-gravity created by the cosmic dark-energy background dominates in the surrounding Hubble flow. In contrast, self-gravity likewise dominates throughout the NGC 1023 Group, both in its central component and in the surrounding “anti-Hubble” flow. NGC 1023 as a whole is apparently in an ongoing state of formation and virialization. We expect that there exists a receding flow similar to the local Hubble flow at distances of 1.4–3 Mpc from the center of the group, where anti-gravity should become stronger than the gravity of the system.     

A dynamic network model for two‐phase immiscible flow

A dynamic pore‐scale network model is formulated for two‐phase immiscible flow. Interfaces are tracked through the pore throats using a modified Poiseuille equation, whereas special displacement rules are used at the pore bodies. The model allows interfaces to move over several pore‐lengths within a time step. Initial computational results are presented for a drainage experiment to demonstrate some of the features of the model. This revised version was published online in July 2006 with corrections to the Cover Date.     

Spherical harmonic analysis of earth’s conductive heat flow

A reappraisal of the international heat flow database has been carried out and the corrected data set was employed in spherical harmonic analysis of the conductive component of global heat flow. Procedures used prior to harmonic analysis include analysis of the heat flow data and determination of representative mean values for a set of discretized area elements of the surface of the earth. Estimated heat flow values were assigned to area elements for which experimental data are not available. However, no corrections were made to account for the hypothetical effects of regional-scale convection heat transfer in areas of oceanic crust. New sets of coefficients for 12° spherical harmonic expansion were calculated on the basis of the revised and homogenized data set. Maps derived on the basis of these coefficients reveal several new features in the global heat flow distribution. The magnitudes of heat flow anomalies of the ocean ridge segments are found to have mean values of less than 150 mW/m². Also, the mean global heat flow values for the raw and binned data are found to fall in the range of 56–67 mW/m², down by nearly 25% compared to the previous estimate of 1993, but similar to earlier assessments based on raw data alone. To improve the spatial resolution of the heat flow anomalies, the spherical harmonic expansions have been extended to higher degrees. Maps derived using coefficients for 36° harmonic expansion have allowed identification of new features in regional heat flow fields of several oceanic and continental segments. For example, lateral extensions of heat flow anomalies of active spreading centers have been outlined with better resolution than was possible in earlier studies. Also, the characteristics of heat flow variations in oceanic crust away from ridge systems are found to be typical of conductive cooling of the lithosphere, there being little need to invoke the hypothesis of unconfined hydrothermal circulation on regional scales. Calculations of global conductive heat loss, compatible with the observational data set, are found to fall in the range of 29–34 TW, nearly 25% less than the 1993 estimate, which rely on one-dimensional conductive cooling models.     

Debris flow affecting the Cubatão Oil Refinery,Brazil

On February 6, 1994, a large debris flow developed because of intense rains in a 800-m-high mountain range called Serra do Cubatão, the local name for the Serra do Mar, located along the coast of the state of São Paulo, Brazil. It affected the Presidente Bernardes Refinery, owned by Petrobrás, in Cubatão. The damages amounted to about US $40 million because of the muck cleaning, repairs, and 3-week interruption of the operations. This prompted Petrobrás to conduct studies, carried out by the authors, to develop protection works, which were done at a cost of approximately US $12 million. The paper describes the studies conducted on debris flow mechanics. A new criteria to define rainfall intensities that trigger debris flows is presented, as well as a correlation of slipped area with soil porosity and rain intensity. Also presented are (a) an actual grain size distribution of a deposited material, determined by laboratory and a large-scale field test, and (b) the size distribution of large boulders along the river bed. Based on theory, empirical experience and back-analysis of the events, the main parameters as the front velocity, the peak discharge and the volume of the transported sediments were determined in a rational basis for the design of the protection works. Finally, the paper describes the set of the protection works built, emphasizing their concept and function. They also included some low-cost innovative works.     

Block-preconditioned Newton–Krylov solvers for fully coupled flow and geomechanics

The focus of this work is efficient solution methods for mixed finite element models of variably saturated fluid flow through deformable porous media. In particular, we examine preconditioning techniques to accelerate the convergence of implicit Newton–Krylov solvers. We highlight an approach in which preconditioners are built from block-factorizations of the coupled system. The key result of the work is the identification of effective preconditioners for the various sub-problems that appear within the block decomposition. We use numerical examples drawn from both linear and nonlinear hydromechanical models to test the robustness and scalability of the proposed methods. Results demonstrate that an algebraic multigrid variant of the block preconditioner leads to mesh-independent convergence, good parallel efficiency, and insensitivity to the material parameters of the medium.     

Lower crust indentation or horizontal ductile flow during continental collision?

Conditions for indentation and channelised flow are investigated with two-dimensional thermomechanical models of Alpine-type continental collision. The models mimic the development of an orogen at an initial central portion of weakened lithosphere 150 km wide, coherent with several geological reconstructions. We study in particular the role of lower crustal strength in developing peculiar geometries after 20 Ma of shortening at 1 cm/year. Crustal layers produce geometries of imbricate layers, which result from two contrasted mechanisms of either channelised ductile lateral flow or horizontal rigid-like indentation:

– Channelised lateral flow develops when the lateral lower crust has a viscosity less than 10²¹ Pa s, exhibiting velocities opposite to the direction of convergence. This mechanism of deformation produces subhorizontal shear zones at the boundaries between the lower crust and the more competent upper crust and lithospheric mantle. It is also associated with a topographic plateau that equilibrates with a wide (about 200 km) but quasi-constant crustal root about 50 km deep.

– In contrast, indentation occurs with lateral lower crust layers that have a viscosity greater than about 10²³ Pa s, producing significant shortening and thickening of the central crust. In this case topography develops steep and narrow (around 100 km wide), associated with a thickened crust exceeding 60 km depth. A crustal-scale pop-up forms bounded by subvertical shear zones that root into the mantle lithosphere.

Compartmented flow at the Bátaapáti site in Hungary

Integration of hydrogeological and geological data into a conceptual model is critical for site investigation programs, since this model is the basis for hydrogeological modeling and for engineering. In the Hungarian Radioactive Waste Disposal Investigation Program, several methods have been used to characterize the potential host rock (granite) at the Bátaapáti site for a repository for low and intermediate level radioactive waste. Hydrogeological data acquisition revealed some characteristic aspects of the site. One of the most important is the presence of extensive rock deformation zones (faults) with low hydraulic conductivity, which strongly reduce the direct hydraulic communication between adjacent blocks of rock (compartmentalization). This characteristic of the rock mass results in a mosaic-like distribution of rock compartments, each with an almost constant hydraulic head. Within the compartments, hydraulic tests have shown that transmissivity is strongly scale dependent: the larger the scale, the higher the measured transmissivity. The extensive highly transmissive zones cause very low hydraulic gradients within each block, thus the transport processes are strongly influenced by the low or average transmissivity zones and the rock matrix.     

Numerical simulation of air–soil two-phase flow based on turbulence modeling

Soil flow and induced air blasts are of great harm to humanity, and historically they have caused a lot of damage to infrastructure. However, these phenomena cannot be described by traditional analog modeling methods that limit their use in disaster prevention efforts. Computational fluid dynamics (CFD) is an applied technique commonly used in a range of fields including the chemical industry, and aircraft and automobile manufacturing, but little is reported on the use of this method to simulate flowing soil in geotechnical engineering applications. The CFD method can effectively make up for the deficiency of normal calculation methods in the analysis of soil flow and air blasts. This paper uses the FLUENT (version 6.3) CFD calculation software to simulate the processes of soil flow and induced air blast changes during soil flow with an Eulerian air–soil two-phase model included in a standard k-ε turbulence model. Velocity vectors of air blasts at different times during soil flow are obtained, and the characteristics of turbulent flow can be found based on the velocity vectors. The numerical simulation techniques adopted in this paper captured precise configurations of soil flow. The results show that the CFD method is especially suitable for simulating the process of soil flow; hazard assessments can be implemented, and the performance of structures involved with disaster prevention can be improved based on the numerical simulation of changing air blasts.     

Catastrophe revisited – disastrous flow failures of mine and municipal solid waste

Catastrophic flow failures have occurred with alarming frequency in mine tailings dams and dumps of discards and other mine waste. In recent years, catastrophic flow failures have also occurred in dumps of municipal solid waste and even in what were intended to be carefully controlled and well-engineered landfills. Apart from the environmental devastation caused by these flows, they are also dangerous to human life and society. Examples include the Buffalo Creek disaster in the USA in 1972 (Vick, Planning, Design and Analysis of Tailings Dams, Wiley, 1983) that killed 118 people, made 4000 homeless and destroyed 50 million US dollars worth of property and facilities. The flow slide that occurred in the Umraniye– Hekimbasi refuse dump in Turkey in 1993 (Kocasoy and Curi, Waste Management and Research, 13, 1995, 305) killed 39 people, destroying their homes in the process. This paper will briefly review some of the more notable flow slides in waste materials, analyzing conditions necessary for a flow failure to occur and pointing to ways of preventing this type of failure by a combination of sound design and operating procedures.     

Deformation band populations in fault damage zone—impact on fluid flow

Fault damage zones in highly porous reservoirs are dominated by deformation bands that generally have permeability-reducing properties. Due to an absence of sufficiently detailed measurements and the irregular distribution of deformation bands, a statistical approach is applied to study their influence on flow. A stochastic model of their distribution is constructed, and band density, distribution, orientation, and flow properties are chosen based on available field observations. The sensitivity of these different parameters on the upscaled flow is analyzed. The influence of a heterogeneous permeability distribution was also studied by assuming the presence of high permeability holes within bands. The fragmentation and position of these holes affect significantly the block-effective permeability. Results of local upscaling with a diagonal and full upscaled permeability tensor are compared, and qualitatively similar results for the flow characteristics are obtained. Further, the procedure of iterative local–global upscaling is applied to the problem.     

Linear flow behaviour of matched joints – case study on standard JRC profiles

Roughness on rock joints produces a variable aperture across the joints and increases the flow path length. These conditions should be taken into account for a good approximation from cubic law. In this paper, the concept of local true aperture and tortuosity is applied to assumed joints where surfaces are matched to each other and correspond with standard Joint Roughness Coefficient (JRC) profiles. Furthermore, the hydraulic behaviour of JRC profiles is studied by a new laboratory experiment setup. The analytical approach provides new insights into the effects of roughness on hydraulic properties of rock joints. The results indicate that for a constant mechanical aperture, both the minimum local aperture and hydraulic aperture decrease with increasing JRC. Furthermore, tortuosity and standard deviation of local true aperture increase with JRC increment. The trend obtained between different parameters and JRC shows an obvious fluctuation for JRC lower than 10. On one hand, the results of this study along with a critical review of previous studies demonstrate that JRC profiles cannot present a precise roughness increment when JRC is less than 10. A new laboratory setup was designed to study the flow behaviour of JRC profiles. The results obtained from laboratory experiments under linear flow conditions validate the accuracy of the applied analytical method.     

How appropriate is the Thiem equation for describing groundwater flow to actual wells?

The Thiem equation of radial groundwater flow to a well is more than 100 years old and is still commonly used. Here, deviations caused by some of its simplifications are quantified by comparing the analytical to a numerical model that allows the implementation of more complex geometries. The assumption of horizontal flow in the Thiem equation, which necessitates uniform inflow over the entire screen length of the fully penetrating well, was found to cause deviations from actual pumping wells where the pump is placed above the screen, resulting in non-uniform inflow and additional drawdown. The same applies to partially penetrating wells, where inflow peaks and additional drawdown occur, especially when the well is screened in the lower part of the aquifer. The use of the Thiem equation in the near-field of a well should thus be restricted to situations where the screen inflow is relatively uniformly distributed, e.g. when it covers large portions of the aquifer thickness. The presence of a gravel pack and a background gradient, on the other hand, are of limited importance.     

Preparing stream flow drought severity–duration–frequency curves using threshold level method

Drought is a natural phenomenon which occurs in different climate regimes. In the present study, hydrological drought of the Roud Zard basin has been investigated based on run theory. Daily runoff data of Mashin hydrometery station during 1970 to 2012 was assessed using 70 % of mean daily runoff as threshold level. Results showed that the maximum drought duration of 309 days occurred in 1998 and 1999 and max drought deficit of 117.217 million cubic meters per second in 1983 with 275 days duration. Time series of the annual maxima values of duration and volume deficit showed similar trend of increase and decreasing. Burr statistical distribution, as the most suitable one fitted to the drought duration data, forecasted 339 days duration for drought event with 20 years return period and generalized extreme value forecasted 37.9 million cubic meters of deficit volume for this return period. Severity-duration-frequency (SDF) curves were prepared, classifying drought durations to four intervals and fitting statistical distribution to each. Resulted SDF curves showed that, in each period, increase of duration resulted in increased value of the volume deficit with a non-linear trend, though predicted drought volume with 20 years return period was 2.1 million cubic meters for 1 to 10 days duration, 6.9 for 11 to 30 days, 34.5 for 31 to 120 days, and 79.1 for more than 120 days duration drought event. Drought deficit volume increasing rate was also different in each class of duration interval. Drought SDF curves derived in this study can be used to quantify water deficit for natural stream and reservoir. SDFs could also be extended to allow for drought regional frequency analysis to be used in ungauged sites.     

High temperature centrifuge simulation — New method to study multiphase flow of partially molten zones

Present-day physical theories suggest that the mantle was partially or completely molten about 4.5 billion years ago. As is expected, the phase separation could result in the stratification of the initially homogeneous mantle into reservoirs of different chemical composition. Because of the poor knowledge of the physical phenomena connected with the multiphase system differentiation, many aspects of the evolution of partially molten zones remain disputable. In this study the simulation of liquid and crystal separation in a high- temperature centrifuge was carried out. Experimental simulation with the pyroxene + olivine + plagioclase + magnetite + melt mixture demonstrates the gravitational instability of the polycrystalline framework of partially molten zones. The melt motion in partially molten rocks is thought to be a multiphase flow involving both the melt and individual crystals. A high mobility of plagioclase crystals in the multiphase flow was ascertained. It is revealed as the formation of a monomineral plagioclase cumulate overlying the accumulated molten layer. The separation and motion of sulfides in the intercrystalline space is shown to be in an intimate relationship with the degree of melting of a silicate material.     

New radiocarbon age constraints for the 120 km-long Toomba flow,north Queensland,Australia

The Toomba flow is the youngest flow of the Nulla volcanic province, located in north Queensland. This 120?km-long flow has yielded a published ⁴⁰Ar/³⁹Ar age of 21,000?±?3000 years. In contrast, seven published conventional radiocarbon (¹⁴C) analyses of carbon-bearing material beneath the flow yielded radiocarbon ages of 16,000 to <2500?BP. These radiocarbon ages are younger than the ⁴⁰Ar/³⁹Ar age, potentially due to contamination of the charcoal by younger carbon that was not removed by the acid–base pre-treatment. We have re-examined the radiocarbon age of Toomba flow using newly sampled charcoal buried beneath the Toomba flow in combination with hydrogen pyrolysis pre-treatment and accelerated mass spectrometer (AMS) measurements. We determined a calibrated radiocarbon age of 20,815–19,726?cal BP (2σ) for the material beneath the Toomba flow. Our radiocarbon age, therefore: (1) is older than previous radiocarbon ages for the Toomba flow, (2) provides the most precise age yet available for the Toomba flow, (3) is in agreement with the ⁴⁰Ar/³⁹Ar age, and (4) validates that hydrogen pyrolysis is a robust and effective pre-treatment method, for subtropical conditions where samples are susceptible to contamination by younger carbon. The Toomba flow erupted during the Last Glacial Maximum, but the preserved surface suggests that the rate of weathering and soil formation has been almost negligible in this region, despite being situated in a subtropical climate that experiences highly variable often intense rainfall.