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.     

Fast iterative solution of large undrained soil‐structure interaction problems

In view of rapid developments in iterative solvers, it is timely to re‐examine the merits of using mixed formulation for incompressible problems. This paper presents extensive numerical studies to compare the accuracy of undrained solutions resulting from the standard displacement formulation with a penalty term and the two‐field mixed formulation. The standard displacement and two‐field mixed formulations are solved using both direct and iterative approaches to assess if it is cost‐effective to achieve more accurate solutions. Numerical studies of a simple footing problem show that the mixed formulation is able to solve the incompressible problem ‘exactly’, does not create pressure and stress instabilities, and obviate the need for an ad hoc penalty number. In addition, for large‐scale problems where it is not possible to perform direct solutions entirely within available random access memory, it turns out that the larger system of equations from mixed formulation also can be solved much more efficiently than the smaller system of equations arising from standard formulation by using the symmetric quasi‐minimal residual (SQMR) method with the generalized Jacobi (GJ) preconditioner. Iterative solution by SQMR with GJ preconditioning also is more elegant, faster, and more accurate than the popular Uzawa method. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

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     

Sea‐surface temperature and salinity changes in the northwest Pacific since the Last Glacial Maximum

The oxygen isotope records of both benthic and planktonic Foraminifera in five piston cores, collected from the region between the Oyashio and Kuroshio Currents near Japan, clearly show the marked latitudinal shifts of these two currents during the past 25 kyr. Under the present hydrographic condition, a clear relationship between the sea‐surface temperature (SST) and oxygen isotope differences from benthic to planktonic Foraminifera is observed in this region. Using this relationship, we find decreased SSTs of 12–13°C (maximum 15°C) in the southernmost core site at the Last Glacial Maximum (LGM), indicating the Oyashio Current shifted southward. The SSTs at the southern two core sites abruptly increased more than 10°C at 10–11 ka, suggesting the Kuroshio Current shifted northward over these sites at 10–11 ka. In contrast, the northern two core sites have remained under the influence of the cold Oyashio Current for the past 25 kyr. With the reasonable estimate of bottom‐water temperature decrease of 2.5°C at the LGM, the SSTs estimated by this new method give exactly the same SST values calculated from Mg/Ca ratio of planktonic Foraminifera, allowing palaeosea‐surface salinities to be reconstructed. The result suggests that the ice volume effect was 1.0 ± 0.1‰ at the LGM. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Particle accretion and dissipation simulator: Collisional aggregation of icy particles

As a part of a project to study the cometary particle environment, a self-consistent ballistic rigid body particle colliding simulator, capable of producing fractal aggregates of icy particles modelled by hard or soft (metaball) spheres has been developed, to be used for creating test particles in future statistical studies of the cometary coma-particle interaction. Instead of a predefined sticking probability, actual surface interactions are modelled and a method for calculating internal stresses and fragmentation of an aggregate during the accretion phase is described. Simulations suggest that taking fragmentation into account has two major consequences, increase in the fractal dimension of resulting particles and a noticeable change in the mass spectrum of an ensemble.