Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Sensitivity analysis of a deterministic water temperature model to forest canopy and soil temperature in Catamaran Brook (New Brunswick,Canada)

A coupled deterministic hydrological and water temperature model, CEQUEAU, was modified to include soil temperature and crown closure in its calculation of local advective terms in the heat budget. The modified model was than tested to verify its sensitivity to these modifications. An analysis of the heat budget of a small forested catchment in eastern Canada revealed that the advective term related to interflow plays a significant role in the daily water heat budget, providing on average 28% of the local advective budget (which also includes advective heat terms from surface runoff and groundwater) and nearly 14% of the total heat budget (which includes all radiative terms at the water surface, convection and evaporation, as well as the local advective terms). Relative sensitivity indices (RSIs) were used to verify the impact of the newly introduced parameters and variables. Among them, parameters related to the forest cover (crown closure and leaf area index) have a maximum RSI of ?0·6; i.e. a 100% increase in value produces a 60% decrease in the local advective term. Parameters with the greatest influence are the volume of water contributing to interflow and the amplitude of the net radiative flux at the soil surface, which, if doubled, would double the contribution of the local interflow advective term to the heat budget. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling the effect of low soil temperatures on transpiration by Scots pine

For ecosystem modelling of the Boreal forest it is important to include processes associated with low soil temperature during spring‐early summer, as these affect the tree water uptake. The COUP model, a physically based SVAT model, was tested with 2 years of soil and snow physical measurements and sap flow measurements in a 70‐year‐old Scots pine stand in the boreal zone of northern Sweden. During the first year the extent and duration of soil frost was manipulated in the field. The model was successful in reproducing the timing of the soil warming after the snowmelt and frost thaw. A delayed soil warming, into the growing season, severely reduced the transpiration. We demonstrated the potential for considerable overestimation of transpiration by the model if the reduction of the trees' capacity to transpire due to low soil temperatures is not taken into account. We also demonstrated that the accumulated effect of aboveground conditions could be included when simulating the relationship between soil temperature and tree water uptake. This improved the estimated transpiration for the control plot and when soil warming was delayed into the growing season. The study illustrates the need of including antecedent conditions on root growth in the model in order to catch these effects on transpiration. The COUP model is a promising tool for predicting transpiration in high‐latitude stands. Copyright © 2006 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.