Wavelet analysis of inter‐annual variability in the runoff regimes of glacial and nival stream catchments,Bow Lake,Alberta

Continuous wavelet analyses of hourly time series of air temperature, stream discharge, and precipitation are used to compare the seasonal and inter‐annual variability in hydrological regimes of the two principal streams feeding Bow Lake, Banff National Park, Alberta: the glacial stream draining the Wapta Icefields, and the snowmelt‐fed Bow River. The goal is to understand how water sources and flow routing differ between the two catchments. Wavelet spectra and cross‐wavelet spectra were determined for air temperature and discharge from the two streams for summers (June–September) 1997–2000, and for rainfall and discharge for the summers of 1999 and 2000. The diurnal signal of the glacial runoff was orders of magnitude higher in 1998 than in other years, indicating that significant ice exposure and the development of channelized glacial drainage occurred as a result of the 1997–98 El Niño conditions. Early retreat of the snowpack in 1997 and 1998 led to a significant summer‐long input of melt runoff from a small area of ice cover in the Bow River catchment; but such inputs were not apparent in 1999 and 2000, when snow cover was more extensive. Rainfall had a stronger influence on runoff and followed quicker flow paths in the Bow River catchment than in the glacial catchment. Snowpack thickness and catchment size were the primary controls on the phase relationship between temperature and discharge at diurnal time scales. Wavelet analysis is a fast and effective means to characterize runoff, temperature, and precipitation regimes and their interrelationships and inter‐annual variability. The technique is effective at identifying inter‐annual and seasonal changes in the relative contributions of different water sources to runoff, and changes in the time required for routing of diurnal meltwater pulses through a catchment. However, it is less effective at identifying changes/differences in the type of the flow routing (e.g. overland flow versus through flow) between or within catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

Comparison between polar regions of Mars from HEND/Odyssey data

In this paper, we have analyzed neutron spectroscopy data gathered by the High Energy Neutron Detector (HEND) instrument onboard Mars Odyssey for comparison of polar regions. It is known that observation of the neutron albedo of Mars provides important information about the distribution of water-ice in subsurface layers and about peculiarities of the CO₂ seasonal cycle. It was found that there are large water-rich permafrost areas with contents of up to ∼50% water by mass fraction at both the north and south Mars polar regions. The water-ice layers at high northern latitudes are placed close to the surface, but in the south they are covered by a dry and relatively thick (10-20 cm) layer of soil. Analysis of temporal variations of neutron flux between summer and winter seasons allowed the estimation of the masses of the CO₂ deposits which seasonally condense at the polar regions. The total mass of the southern seasonal deposition was estimated as 6.3×¹⁰¹⁵ kg, which is larger than the total mass of the seasonal deposition at the north by 40-50%. These results are in good agreement with predictions from the NASA Ames Research Center General Circulation Model (GCM). But, the dynamics of the condensation and sublimation processes are not quite as consistent with these models: the peak accumulation of the condensed mass of CO₂ occurred 10-15 degrees of Ls later than is predicted by the GCM.     

A coupled microphysical/radiative transfer model of albedo and emissivity of planetary surfaces covered by volatile ices

A new model of albedo and emissivity of planetary surfaces covered by volatile ices in the form of porous slab-like deposits is described. In the model, a radiative transfer model is coupled with a microphysical model in order to link changes in albedo and emissivity to changes in porosity caused by ice metamorphism. Preliminary results for Triton, Pluto, and Io are presented (the martian CO₂ caps will be the subject of a separate publication). The coupled model will aid in the interpretation of ground-based and spacecraft observations and should lead to advances in surface and atmospheric modeling.     

Longevity of fluorine-bearing tremolite on Venus

There is a general belief that hydrous minerals cannot exist on Venus under current surface conditions. This view was challenged when Johnson and Fegley (2000, Icarus 146, 301-306) showed that tremolite (Ca₂Mg₅Si₈O₂₂(OH)₂), a hydrous mineral, is stable against thermal decomposition at current Venus surface temperatures, e.g., 50% decomposition in 4 Ga at 740 K. To further explore hydrous mineral thermal stability on Venus, we experimentally determined the thermal decomposition kinetics of fluorine-bearing tremolite. Fluor-tremolite is thermodynamically more stable than OH-tremolite and should decompose more slowly. However how much slower was unknown. We measured the decomposition rate of fluorine-bearing tremolite and show that its decomposition is several times to greater than ten times slower than that of OH-tremolite. We also show that F-bearing tremolite is depleted in fluorine after decomposition and that fluorine is lost as a volatile species such as HF gas. If tremolite ever formed on Venus, it would probably also contain fluorine. The exceptional stability of F-bearing tremolite strengthens our conclusions that if hydrous minerals ever formed on Venus, they could still be there today.     

Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Reliability assessment of serviceability performance of braced retaining walls using a neural network approach

In urban environments, one major concern with deep excavations in soft clay is the potentially large ground deformations in and around the excavation. Excessive movements can damage adjacent buildings and utilities. There are many uncertainties associated with the calculation of the ultimate or serviceability performance of a braced excavation system. These include the variabilities of the loadings, geotechnical soil properties, and engineering and geometrical properties of the wall. A risk‐based approach to serviceability performance failure is necessary to incorporate systematically the uncertainties associated with the various design parameters. This paper demonstrates the use of an integrated neural network–reliability method to assess the risk of serviceability failure through the calculation of the reliability index. By first performing a series of parametric studies using the finite element method and then approximating the non‐linear limit state surface (the boundary separating the safe and ‘failure’ domains) through a neural network model, the reliability index can be determined with the aid of a spreadsheet. Two illustrative examples are presented to show how the serviceability performance for braced excavation problems can be assessed using the reliability index. Copyright © 2005 John Wiley & Sons, Ltd.     

Fresnel modeling of hematite crystal surfaces and application to martian hematite spherules

The gray crystalline hematite at Meridiani Planum first discovered by the Mars Global Surveyor Thermal Emission Spectrometer (MGS-TES) instrument occurs as spherules that have been interpreted as concretions. Analysis of the TES and mini-TES spectra shows that no 390 cm⁻¹ feature is present in the characteristic martian hematite spectrum. Here, we incorporate the mid-IR optical constants of hematite into a simple Fresnel reflectance model to understand the effect of emission angle and crystal morphology on the presence or absence of the 390 cm⁻¹ feature in an IR hematite spectrum. Based on the results we offer two models for the internal structure of the martian hematite spherules.     

Discriminating between pore‐filling load and bed‐structure load: a new porosity‐based method,exemplified for the river Rhine

Sediments contained in the river bed do not necessarily contribute to morphological change. The finest part of the sediment mixture often fills the pores between the larger grains and can be removed without causing a drop in bed level. The discrimination between pore‐filling load and bed‐structure load, therefore, is of practical importance for morphological predictions. In this study, a new method is proposed to estimate the cut‐off grain size that forms the boundary between pore‐filling load and bed‐structure load. The method evaluates the pore structure of the river bed geometrically. Only detailed grain‐size distributions of the river bed are required as input to the method. A preliminary validation shows that the calculated porosity and cut‐off size values agree well with experimental data. Application of the new cut‐off size method to the river Rhine demonstrates that the estimated cut‐off size decreases in a downstream direction from about 2 to 0·05 mm, covariant with the downstream fining of bed sediments. Grain size fractions that are pore‐filling load in the upstream part of the river thus gradually become bed‐structure load in the downstream part. The estimated (mass) percentage of pore‐filling load in the river bed ranges from 0% in areas with a unimodal river bed, to about 22% in reaches with a bimodal sand‐gravel bed. The estimated bed porosity varies between 0·15 and 0·35, which is considerably less than the often‐used standard value of 0·40. The predicted cut‐off size between pore‐filling load and bed‐structure load (D_c,p) is fundamentally different from the cut‐off size between wash‐load and bed‐material load (D_c,w), irrespective of the method used to determine D_c,p or D_c,w. D_c,w values are in the order of 10^?1 mm and mainly dependent on the flow characteristics, whereas D_c,p values are generally much larger (about 10⁰ mm in gravel‐bed rivers) and dependent on the bed composition. Knowledge of D_c,w is important for the prediction of the total sediment transport in a river (including suspended fines that do not interact with the bed), whereas knowledge of D_c,p helps to improve morphological predictions, especially if spatial variations in D_c,p are taken into account. An alternative to using a spatially variable value of D_c,p in morphological models is to use a spatially variable bed porosity, which can also be predicted with the new method. In addition to the morphological benefits, the new method also has sedimentological applications. The possibility to determine quickly whether a sediment mixture is clast‐supported or matrix‐supported may help to better understand downstream fining trends, sediment entrainment thresholds and variations in hydraulic conductivity.     

What controls dead-ice melting under different climate conditions? A discussion

In the geological record, hummocky dead-ice moraines represent the final product of the melt-out of dead-ice. Processes and rates of dead-ice melting in ice-cored moraines and at debris-covered glaciers are commonly believed to be governed by climate and debris-cover properties. Here, backwasting rates from 14 dead-ice areas are assessed in relation to mean annual air temperature, mean summer air temperature, mean annual precipitation, mean summer precipitation, and annual sum of positive degree days. The highest correlation was found between backwasting rate and mean annual air temperature. However, the correlation between melt rates and climate parameters is low, stressing that processes and topography play a major role in governing the rates of backwasting. The rates of backwasting from modern glacial environments should serve as input to de-icing models for ancient dead-ice areas in order to assess the mode and duration of deposition.A challenge for future explorations of dead-ice environments is to obtain long-term records of field-based monitoring of melt progression. Furthermore, many modern satellite-borne sensors have high potentials for recordings of multi-temporal Digital Elevation Models (DEMs) for detection and quantification of changes in dead-ice environments. In recent years, high-accuracy DEMs from airborne laser scanning altimetry (LiDAR) are emerging as an additional data source. However, time series of high-resolution aerial photographs remain essential for both visual inspection and high-resolution stereographic DEM production.     

Influence of Structural Non-Stationarity of Surface Roughness on Morphological Characterization and Mechanical Deformation of Rock Joints

Summary Structural non-stationarity of surface roughness affects accurate morphological characterization as well as mechanical behaviour of rock joints at the laboratory scale using samples with a size below the stationarity threshold. In this paper, the effect of structural non-stationarity of surface roughness is investigated by studying the scale dependence of surface roughness and mechanical behaviour of rock joints. The results show that the structural non-stationarity mainly affects the accurate characterization of the surface roughness of the fracture samples. It also controls the amount and location of the contact areas during shear tests, which in turn affects the mechanical properties and asperity degradation of the samples. It is concluded that for accurate determination of the morphological and mechanical properties of rock joints at laboratory and field scales, samples with size equal to or larger than the stationarity threshold are required. Author’s address: Nader Fardin, Rock Mechanics Group, Department of Mining Engineering, Faculty of Engineering, University of Tehran, P.O. Box: 11365/4563, Tehran, Iran