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.     

WEKU – a GIS-Supported stochastic model of groundwater residence times in upper aquifers for the supraregional groundwater management

 The supraregional GIS-supported stochastical model, WEKU, for the determination of groundwater residence times in the upper aquifers of large groundwater provinces is presented. Using a two-dimensional analytical model of groundwater flow, groundwater residence times are determined within two extreme cases. In the first case, maximal groundwater residence times are calculated, representing the part of groundwater, that is drained by the main surface water of a groundwater catchment area. In the second case, minimal groundwater residence times for drainage into the nearest surface water are determined. Using explicit distribution functions of the input parameters, mean values as well as potential ranges of variations of the groundwater residence times are derived. The WEKU model has been used for the determination of groundwater residence times throughout Germany. The model results – mean values and deviations of the groundwater velocity and the maximal and minimal groundwater residence times in the upper aquifers – are presented by general maps and discussed in detail. It is shown that the groundwater residence times in the upper aquifer vary regionally, differentiated between less than 1 year and more than 2000 years. Using this information, the time scales can be specified, until measures to remediate polluted groundwater resources may lead to a substantial groundwater quality improvement in the different groundwater provinces of Germany. With respect to its supraregional scale of application, the WEKU model may serve as a useful tool for the supraregional groundwater management on a state, federal or international level. Received: 15 August 1995 · Accepted: 15 October 1995     

Studying Hydrodynamic Instability Using Shock-Tube Experiments

The hydrodynamic instability, which develops on the contact surface between two fluids, has great importance in astrophysical phenomena such as the inhomogeneous density distribution following a supernova event. In this event acceleration waves pass across a material interface and initiate and enhance unstable conditions in which small perturbations grow dramatically. In the present study, an experimental technique aimed at investigating the above-mentioned hydrodynamic instability is presented. The experimental investigation is based on a shock-tube apparatus by which a shock wave is generated and initiates the instability that develops on the contact surface between two gases. The flexibility of the system enables one to vary the initial shape of the contact surface, the shock-wave Mach number, and the density ratio across the contact surface. Three selected sets of shock-tube experiments are presented in order to demonstrate the system capabilities: (1) large-initial amplitudes with low-Mach-number incident shock waves; (2) small-initial amplitudes with moderate-Mach-number incident shock waves; and (3) shock bubble interaction. In the large-amplitude experiments a reduction of the initial velocity with respect to the linear growth prediction was measured. The results were compared to those predicted by a vorticity-deposition model and to previous experiments with moderate- and high-Mach number incident shock waves that were conducted by others. In this case, a reduction of the initial velocity was noted. However, at late times the growth rate had a 1/t behavior as in the small-amplitude low-Mach number case. In the small-amplitude moderate-Mach number shock experiments a reduction from the impulsive theory was noted at the late stages. The passage of a shock wave through a spherical bubble results in the formation of a vortex ring. Simple dimensional analysis shows that the circulation depends linearly on the speed of sound of the surrounding material and on the initial bubble radius.     

RELAXING ASSUMPTIONS OF PERFECT INFORMATION IN PARK VISITATION MODELS*

Variables related to urban park awareness are identified and methods for relaxing assumptions of perfect information in park use models are discussed. Park awareness is related to park characteristics (age and degree of development of the park), population characteristics (race, age, length of residence, recreation participation), and distance. Park attributes are stronger predictors of both park awareness and use than is distance. These findings parallel similar research on the cognitive aspects of shopping decisions.     

Crustal and upper mantle structures of the brazilian coast

The Rayleigh wave phase and group velocities in the period range of 24–39 sec, obtained from two earthquakes which occurred in northeastern brazil and which were recorded by the Brazilian seismological station RDJ (Rio de Janeiro), have been used to study crustal and upper mantle structures of the Brazilian coastal region. Three crustal and upper mantle models have been tried out to explain crustal and upper mantle structures of the region. The upper crust has not been resolved, due basically to the narrow period range of the phase and group velocities data. The phase velocity inversions have exhibited good resolutions for both lower crust and upper mantle, with shear wave velocities characteristic of these regions. The group velocity data inversions for these models have showed good results only for the lower crust. The shear wave velocities of the lower crust (3.86 and 3.89 km/sec), obtained with phase velocity inversions, are similar to that (=3.89 km/sec) found byHwang (1985) to the eastern South American region, while group velocity inversions have presented shear velocity (=3.75 km/sec) similar to that (=3.78 km/sec) found byLazcano (1972) to the Brazilian shield. It was not possible to define sharply the crust-mantle transition, but an analysis of the phase and group velocity inversions results has indicated that the total thickness of the crust should be between 30 and 39 km. The crustal and upper mantle model, obtained with phase velocity inversion, can be used as a preliminary model for the Brazilian coast.     

Directions and possibilities of mathematical geology

This paper considers the present state of mathematical geology. Three directions are recognized: applied, theoretical, and mathematical. Applied mathematical geology includes formal use of mathematics to solve problems and computer processing of data. Success is achieved by a correspondence of mathematical methods used to the nature of geological data. This correspondence can be demonstrated by purely mathematical means. Theoretical mathematical geology uses mathematics as a language of geology; however, a number of methodological problems must be solved: formalization of initial geological concepts and creation of a strict conceptual basis, substantiation of initial principles of mathematical simulation, creation of theoretical geological models, problems of elementary and coincidence in geology, and methodological substantiations of possibilities of any mathematical model to approximate geological models. The essense and significance of these problems are considered. The main task of mathematical geology is to prove its correspondence to the nature of the geological objects studied, geological data obtained, and geological problems solvable. Finally, the main problems of mathematical geology are not so much mathematical as geological and methodological.     

岩体三维节点偶对分析及其应用

本文讨论了节点偶对分析的三维扩展。节点偶对分析是针对摩擦滑动节理单元的约束边界进行的序列矢量判定。在三维节理单元的分析中,几何约束和相应的力学判据极为复杂,采用节点偶对等效结合力方法才能迎刃而解。文中对三维转向节理模型的分析得到合理的结果,它表明本文所提出的方法可成功地应用于节理岩体的有限元分析。