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.     

Screening of polymers on selected Hawaii soils for erosion reduction and particle settling

In recent years, high‐molecular‐weight anionic polyacrylamides (PAMs) have been tested on a variety of soils, primarily in temperate climates. However, little information is available regarding the effectiveness of PAM for preventing soil loss through runoff in tropical settings. Screening tests were performed using three negatively charged PAMs and one positively charged PAM on five Hawaii soils (two Oxisols, one Vertisol, and two Aridisols) to determine erosion loss, sediment settling, and aggregate stability. A laboratory‐scale rainfall simulator was used to apply erosive rainfall at intensities from 5 to 8·5 cm h^?1 at various PAM doses applied in both dry and solution forms. Soil detachment due to splash and runoff, as well as the runoff and percolate water volumes, were measured for initial and successive storms. The impact of PAM on particle settling and aggregate stability was also evaluated for selected soil‐treatment combinations. Among the PAMs, Superfloc A‐836 was most effective, and significantly reduced runoff and splash sediment loss for the Wahiawa Oxisol and Pakini Andisol at rates varying between 10 and 50 kg ha^?1. Reduced runoff and splash sediment loss were also noted for PAM Aerotil‐D when applied in solution form to the Wahiawa Oxisol. Significant reductions in soil loss were not noted for either the Lualualei Vertisol or the Holomua Oxisol. It is believed that the high montmorillonite content of the Lualualei Vertisol and the low cation‐exchange capacity of the Holomua Oxisol diminished the effectiveness of the various PAMs tested. The polymers were also found to enhance sediment settling of all soils and helped improve their aggregate stability. This screening study shows the potential use of PAM for tropical soils for applications such as infiltration enhancement, runoff reduction, and enhanced sedimentation of detention ponds. Copyright © 2005 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.     

Spine–sheath polarization structures in four active galactic nuclei jets

We present the results of multifrequency (15 + 8 + 5 GHz) polarization Very Long Baseline Array (VLBA) observations of the three BL Lacertae objects 0745+241, 1418+546 and 1652+398 together with 5-GHz VLBI Space Observatory Programme (VSOP) observations of 1418+546 and 1.6- and 5-GHz VSOP observations of the blazar 1055+018. The jets of all these sources have polarization structure transverse to the jet axis, with the polarization E vectors aligned with the jet along the jet spine and 'sheaths' of orthogonal E vectors at one or both edges of the jet. The presence of polarization aligned with the jet near the 'spine' may indicate that the jets are associated with helical B fields that propogate outward with the jet flow; the presence of orthogonal polarization near the edges of the jet may likewise be a consequence of a helical jet B field, or may be owing to an interaction with the ambient medium on parsec scales. We have tentatively detected interknot polarization in 1055+018 with E aligned with the local jet direction, consistent with the possibility that the jet of this source is associated with a helical B field.     

A kinetic model of gas flow in a porous cometary mantle

A numerical study of gas flow through a porous cometary mantle is presented. A kinetic model based on the well-known Test Particle Monte Carlo Method for the solution of rarefied gas dynamics problems is proposed. The physical model consists of two spatial plane regions: the condensed ice phase and a porous dust mantle. The structure of the porous dust layer is described as a bundle of cylindrical inclined channels not crossing each other. A vertical temperature gradient may exist across the dust mantle. The aim is to investigate how the characteristics of molecular flow depend on the capillary length, inclination angle, and temperature gradient. Examples illustrating a significant deviation of some results from equilibrium values are shown. In particular, the gas velocity distribution at both ends of the pore is strongly non-Maxwellian if there is an important temperature contrast across the pore. The emergent gas flow rate is found to vary with the pore length/radius ratio in excellent agreement with Clausing's empirical formula. The degree of collimation of the flow is quantitatively studied as a function of the length/radius ratio, and consequences for the jet force of outgassing through a dust mantle or, indeed, a rough surface are estimated.     

The production of light elements in unconventional cosmologies with VMOs

The abundances of light elements (deuterium, helium and lithium) are calculated in unconventional cosmologies in which these elements are produced by interaction of fast -particles ejected by massive pregalactic stars with a purely hydrogenous cosmological substrate. The number density and energetics of such stars needed for suitable production of light elements are estimated. It is shown that the models always lead to overproduction of lithium. The intensity of ultraviolet, X-ray, and -ray emission caused by explosions of pregalactic massive stars is calculated. The distortions of cosmic background radiation by Compton losses of electrons heated by explosive stars are also calculated.     

Energetic solar electrons in the interplanetary medium

ISEE-3 measurements extending down to 2 keV energy have provided a new perspective on energetic solar electrons in the interplanetary medium. Impulsive solar electron events are observed, on average, several times a day near solar maximum, with 40% detected only below 15 keV. The electron energy spectra have a nearly power-law shape extending smoothly down to 2 keV, indicating that the origin of these events is high in the corona. These coronal flare-like events often produced ³He-rich particle events.In large solar flares which accelerate electrons and ions to relativistic energies, the electron spectrum appears to be modified by a second acceleration which results in a double power-law shape above 10 keV with a break near 100 keV and flattening from 10–100 keV. Large flares result in long-lived (many days) streams of outflowing electrons which dominate the interplanetary fluxes at low energies. Even in the absence of solar activity, significant fluxes of low energy electrons flow out from the Sun.Solar type-III radio bursts are produced by the escaping 2–10² keV electrons through a beam-plasma instability. The detailed ISEE-3 measurements show that electron plasma waves are generated by the bump-on-tail distribution created by the faster electrons running ahead of the slower ones. These plasma waves appear to be converted into radio emission by nonlinear wave-wave interactions.     

Hydrodynamical models of type II supernovae

A series of hydrodynamical models of type-II supernova outbursts (SNII) has been calculated. Approximate relations connecting the total outburst energy ε, the mass of envelope ejectedM, the presupernova radiusR, and the amount of ionizing quanta radiated by the supernovaeN _H with such values as the duration of the light curve plateau Δt, and absolute magnitude in the wavelength bandV and photospheric velocityU _PH observed near the middle of the plateau have been established. Advantage has been taken of the relations to obtain a preliminary evaluation for the characteristics of the average SN II: ε=7×10⁵⁰ erg,M=6M _⊙,R=500R _⊙,N _H=2×10⁵⁸. The SNIIs with plateau-like light curves seem to be accounted for by thermonuclear explosions of degenerate cores of red giant stars and result in a total disruption of the star without any stellar remnant. To the contrary, SNIIs with linear light curves have substantially different properties (in particular, they throw considerably less massive envelopes off). These SNII must signify the birth of collapsed objects—neutron stars (pulsars) or black holes.