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.     

山西省主要河流流域面雨量预报业务流程

以T213、HLAFS模式、MM5中尺度模式输出的格点资料以及日本降水量格点资料为基础，将影响山西降水的天气动力模型归纳为诊断模型，从中引出多个能够全面反映降水模型特征的综合物理因子；根据各种数值模式输出的降水量预报性能和质量优劣特点，依据数值模式的形势场预报优于要素场预报的现实，构造在不同环流形势背景下，启动不同预报方程的面雨量预报业务流程，有效地遏止了在环流形势调整时预报输出不能快速响应的弱点，提高了点和面雨量预报的准确度。     

鄂尔多斯盆地中生代晚期以来大地热流的变化及其对生态环境格局和演变的影响

鄂尔多斯盆地的西北部、东北部和南部三个区域现今大地热流平均值分别为56．3、67．3和65．3mW／m^2，对应的生态环境格局也有明显的差异。研究表明，大地热流每增加4～5mW／m^2可使年均地表温度升高约l℃，使最低月均地表温度升高2。C以上。鄂尔多斯盆地东北部的平均大地热流比西北部高出11mW／m^2，东北部年均地表温度可能比西北部高出2～3℃，其最低月均地表温度可能比西北部高出4～6℃。西北部的大地热流平均值已经低于维持地表生态系统延续所需大地热流的临界值(57mW／m^2)，其自然生态系统整体上已经处于脆弱境地；东北部和南部的大地热流均大于57mW／m^2，自然生态系统均尚较稳健。东北部的沙漠化可能是风沙侵入的结果，其生态应该是可以恢复的。整个西北部作为一个整体看，72万年以前大地热流就已衰减到临界值以下，区域生态系统渐趋脆弱，开始整体上向荒漠化演变。     

A 2D finite volume model for bebris flow and its application to events occurred in the Eastern Pyrenees

FLATModel is a 2D finite volume code that contains several original approaches to improve debris-flow simulation. Firstly, FLATModel incorporates a ＂stop-and-go＂ technique in each cell to allow continuous collapses and remobilizations of the debris-flow mass. Secondly, flow velocity and consequently yield stress is directly associated with the type of rheology to improve boundary accuracy. Thirdly, a simple approach for entrainment is also included in the model to analyse the effect of basal erosion of debris flows. FLATMODEL was tested at several events that occurred in the Eastern Pyrenees and simulation results indicated that the model can represent rather well the different characteristics observed in the field.