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.     

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.     

The importance of oxygen isotope provenance in relation to solute content of bulk meltwaters at Imersuaq Glacier,West Greenland

Stable oxygen isotope analysis and measurement of several dissolved cations and anions of bulk meltwater samples have provided information about the hydrochemical environment of the glacial hydrological system at Imersuaq Glacier, an outlet tongue from the Greenland ice‐sheet, West Greenland. The samples were collected at frequent intervals during the period 20–28 July 2000 in a small (<20 L s^?1) englacial meltwater outlet at the glacier margin. The results document the following findings: (i) a marked diurnal variation of δ¹⁸O is related to the composition of oxygen isotope provenances, mainly near‐marginal local superimposed ice and basal up‐sheared ice further up‐glacier; (ii) a relationship is seen between all base cations (Na⁺, K⁺, Ca²⁺, Mg²⁺), SO₄^2? and δ¹⁸O, indicating that solute acquisition is provided by solid–solution contact with the up‐sheared ice—as the relationship with Cl^? is weak the influence of seasalt‐derived solutes is small in the area; (iii) when the melt rate is high, two diurnal maxima of δ¹⁸O values and solute concentrations are measured, and it is suggested that a snow meltwater component is responsible for the second maximum of δ¹⁸O—a short residence time leads to a delayed decrease in ion concentrations. Copyright © 2003 John Wiley & Sons, Ltd.     

The Preboreal climate reversal and a subsequent solar‐forced climate shift

Accurate chronologies are essential for linking palaeoclimate archives. Carbon‐14 wiggle‐match dating was used to produce an accurate chronology for part of an early Holocene peat sequence from the Borchert (The Netherlands). Following the Younger Dryas–Preboreal transition, two climatic shifts could be inferred. Around 11 400 cal. yr BP the expansion of birch (Betula) forest was interrupted by a dry continental phase with dominantly open grassland vegetation, coeval with the PBO (Preboreal Oscillation), as observed in the GRIP ice core. At 11 250 cal. yr BP a sudden shift to a humid climate occurred. This second change appears to be contemporaneous with: (i) a sharp increase of atmospheric ¹⁴C; (ii) a temporary decline of atmospheric CO₂; and (iii) an increase in the GRIP ¹⁰Be flux. The close correspondence with excursions of cosmogenic nuclides points to a decline in solar activity, which may have forced the changes in climate and vegetation at around 11 250 cal. yr BP. Copyright © 2004 John Wiley & Sons, Ltd.     

The influence of macrophyte growth,typical of eutrophic conditions,on river flow velocities and turbulence production

The influence of emergent and submerged macrophytes on flow velocity and turbulence production is demonstrated in a 140 m reach of the River Blackwater in Farnborough, Hampshire, UK. Macrophyte growth occurs in patches and is dominated by Sparganium erectum and Sparganium emersum. In May 2001, patches of S. erectum were already established and occupied 18% of the channel area. The flow adjusted to these (predominantly lateral) patches by being channelled through a narrower cross‐section. The measured velocity profiles showed a logarithmic form, with deviations attributable to topographic control. The channel bed was the main source of turbulence. In September 2001, in‐stream macrophytes occupied 27% of the channel, and overhanging bank vegetation affected 32% of the area. Overall flow resistance, described by Manning's n, showed a threefold increase that could be attributed to the growth of S. emersum in the middle of the channel. Velocity profiles showed different characteristic forms depending on their position relative to plant stems and leaves. The overall velocity field had a three‐dimensional structure. Turbulence intensities were generally higher and turbulence profiles tended to mirror the velocity profiles. Evidence for the generation of coherent eddies was provided by ratios of the root mean square velocities. Spectral analysis identified deviations from the Kolmogorov ?5/3 power law and provided statistical evidence for a spectral short‐cut, indicative of additional turbulence production. This was most marked for the submerged vegetation and, in some instances, the overhanging bank vegetation. The long strap‐like leaves of S. emersum being aligned approximately parallel to the flow and the highly variable velocity field created by the patch arrangement of macrophytes suggest that the dominant mechanism for turbulence production is vortex shedding along shear zones. Wake production around individual stems of S. emersum close to the bed may also be important locally. Copyright © 2006 John Wiley & Sons, Ltd.     

Freshening of the Intermediate Waters in the Northern South China Sea over the PastSix Decades

The properties of salinity in the South China Sea (SCS), a significant marginal sea connecting the Pacific andIndian Oceans, are greatly influenced by the transport of fresh water flux between the two oceans. However, the long-termchanges in the intermediate water in the SCS have not been thoroughly studied due to limited data, particularly in relationto its thermodynamic variations. This study utilized reanalysis data products to identify a 60-year trend of freshening in theintermediate waters of the northern South China Sea (NSCS), accompanied by an expansion of low-salinity water. Thestudy also constructed salinity budget terms, including advection and entrainment processes, and conducted an analysis ofthe salinity budget to understand the impacts of external and internal dynamic processes on the freshening trend of theintermediate water in the NSCS. The analysis revealed that the freshening in the northwest Pacific Ocean and theintensification of intrusion through the Luzon Strait at intermediate levels are the primary drivers of the salinity changes inthe NSCS. Additionally, a weakened trend in the intensity of vertical entrainment also contributes to the freshening in theNSCS. This study offers new insights into the understanding of regional deep sea changes in response to variations in boththermodynamics and oceanic dynamic processes.     

Evapotranspiration from citrus trees growing in sandy soilunder drip irrigation with saline water

An experiment on evapotranspiration from citrus trees under irrigation with saline waterwas carried out for 4 months. Two lysimeters planted with a citrus tree in the green house wereused. One lysimeter was irrigated with saline water (NaCl and CaCl2 of 2000 mg/L equivalence,EC = 3.8 dS/m, SAR = 5.9) and the other was irrigated with freshwater using drip irrigation. Theapplied irrigation water was 1.2 times that of the evapotranspiration on the previous day.Evapotranspiration was calculated as the change in lysimeter weight recorded every 30 minutes.The lysimeters were filled with soil with 95.8% sand. The results of the experiment were as follows.(i) The evapotranspiration from citrus tree was reduced after irrigation with saline water. Theevapotranspiration returns to normal after leaching. However it takes months to exhaust the saltfrom the tree. ( ii ) To estimate the impact of irrigation with saline water on the evapotranspirationfrom citrus trees, the reduction coefficient due to salt stress (Ks) was used in this experiment.Evapotranspiration under irrigation with saline water (ETs) can be calculated from evapotranspira-tion under irrigation with freshwater (ET) by the equation ETs = Ks× ET. Ks can be expressed as afunction of ECsw. (iii) The critical soil-water electrical conductivity (ECsw) is 9.5 dS/m, beyondwhich adverse effects on evapotranspiration begin to appear. If ECsw can be controlled at below9.5 dS/m, saline water can be safely used for irrigation.