Mineral weathering and element cycling in soil-microorganism-plant system

Soil is an essential part of the critical zone,and soil-microbe-plant system serves as a key link among lithosphere,biosphere,atmosphere and hydrosphere.As one of the habitats with the richest biodiversity,soil plays a critical role in element biogeochemistry on the earth surface(weathered crust).Here we review the soil biological processes that are relevant to mineral weathering,element cycling,and transformation,with an emphasis on rock weathering mediated by soil microbes,plant root and the rhizosphere.     

Soil formation and weathering in a permafrost environment of the Swiss Alps: a multi‐parameter and non‐steady‐state approach

Spatially discontinuous permafrost conditions frequently occur in the European Alps. How soils under such conditions have evolved and how they may react to climate warming is largely unknown. This study focuses on the comparison of nearby soils that are characterised by the presence or absence of permafrost (active‐layer thickness: 2–3 m) in the alpine (tundra) and subalpine (forest) range of the Eastern Swiss Alps using a multi‐method (geochemical and mineralogical) approach. Moreover, a new non‐steady‐state concept was applied to determine rates of chemical weathering, soil erosion, soil formation, soil denudation, and soil production. Long‐term chemical weathering rates, soil formation and erosion rates were assessed by using immobile elements, fine‐earth stocks and meteoric ¹⁰Be. In addition, the weathering index (K + Ca)/Ti, the amount of Fe‐ and Al‐oxyhydroxides and clay minerals characteristics were considered. All methods indicated that the differences between permafrost‐affected and non‐permafrost‐affected soils were small. Furthermore, the soils did not uniformly differ in their weathering behaviour. A tendency towards less intense weathering in soils that were affected by permafrost was noted: at most sites, weathering rates, the proportion of oxyhydroxides and the weathering stage of clay minerals were lower in permafrost soils. In part, erosion rates were higher at the permafrost sites and accounted for 79–97% of the denudation rates. In general, soil formation rates (8.8–86.7 t/km²/yr) were in the expected range for Alpine soils. Independent of permafrost conditions, it seems that the local microenvironment (particularly vegetation and subsequently soil organic matter) has strongly influenced denudation rates. As the climate has varied since the beginning of soil evolution, the conditions for soil formation and weathering were not stable over time. Soil evolution in high Alpine settings is complex owing to, among others, spatio‐temporal variations of permafrost conditions and thus climate. This makes predictions of future behaviour very difficult. Copyright © 2016 John Wiley & Sons, Ltd.     

Effect of basin physical characteristics on solute fluxes in nine alpine/subalpine basins,Colorado, USA

Alpine/subalpine basins may exhibit substantial variability in solute fluxes despite many apparent similarities in basin characteristics. An evaluation of controls on spatial patterns in solute fluxes may allow development of predictive tools for assessing basin sensitivity to outside perturbations such as climate change or deposition of atmospheric pollutants. Relationships between basin physical characteristics, determined from geographical information system (GIS) tools, and solute fluxes and mineral weathering rates were explored for nine alpine/subalpine basins in Rocky Mountain National Park, Colorado, using correlation analyses for 1993 and 1994 data. Stream‐water nitrate fluxes were correlated positively with basin characteristics associated with the talus environment; i.e., the fractional amounts of steep slopes (≥ 30°), unvegetated terrain and young debris (primarily Holocene till) in the basins, and were correlated negatively with fractional amounts of subalpine meadow terrain. Correlations with nitrate indicate the importance of the talus environment in promoting nitrate flux and the mitigating effect of areas with established vegetation, such as subalpine meadows. Total mineral weathering rates for the basins ranged from about 300 to 600 mol ha^?1 year^?1. Oligoclase weathering accounted for 30 to 73% of the total mineral weathering flux, and was positively correlated with the amount of old debris (primarily Pleistocene glacial till) in the basins. Although calcite is found in trace amounts in bedrock, calcite weathering accounted for up to 44% of the total mineral weathering flux. Calcite was strongly correlated with steep slope, unvegetated terrain, and young debris—probably because physical weathering in steep‐gradient areas exposes fresh mineral surfaces that contain calcite for chemical weathering. Oligoclase and calcite weathering are the dominant sources of alkalinity in the basins. However, atmospherically deposited acids consume much of the alkalinity generated by weathering of calcite and other minerals in the talus environment. Published in 2001 by John Wiley & Sons, Ltd.     

THE DISTRIBUTION OF PREFERENTIAL PATHS AND ITS RELATION TO THE SOIL CHARACTERISTICS IN THE THREE GORGES AREA, CHINA

To study the characteristics of the distribution of the preferential paths and the affecting factors in the Three Gorges area, four soil profiles were dug to observe the distribution of preferential paths in the Quxi watershed in the Yangtze River basin. The Morisita exponential test method was used to examine the distribution type of preferential paths. The physical properties and infiltration characteristics of the soil were also measured to evaluate their relationship to preferential paths. The results showed that in this area, preferential paths clustered and mainly distributed in the 80-100 cm soil layer, and along the interface between the weathered layer and semi-weathered layer. There were more non-capillary pores in the 83-110 cm layer than in the other layers. It can be derived that most non-capillary pores in this layer were preferential paths caused by geological processes and rotten plant roots. The percentage of coarse soil particles increased with the depth of the soil layer. In the deeper soil layer, the coarse soil particles helped the formation of preferential paths. The fastest steady infiltration rate was observed in the of 83-110cm layer, which is inferred to be due to the greater number of preferential paths.     

Comparison of the effects of litter covering and incorporation on infiltration and soil erosion under simulated rainfall

Plant litter can either cover on soil surface or be incorporated into top-soil layer in natural ecosystems. Their effects on infiltration and soil erosion are likely quite different. This study was performed to compare the effects of litter covering on soil surface and being incorporated into top-soil layer on infiltration and soil erosion under simulated rainfall. Four litter types (needle-leaf, broad-leaf, brush, and herb) were collected from fields and applied to cover on soil surface or to be incorporated into top-soil layer (5 cm) at the same rate (0.2 kg/m²). The simulated rainfalls (40 and 80 mm/hr) were run at two slope angles (10° and 20°). The results showed that the mean infiltration rate of litter covering treatment was 1.4 times as great as that of litter incorporated. Litter covering enhanced infiltration via protecting surface from soil sealing. Whereas, litter incorporation affected infiltration by its water repellency. Soil erosion of litter incorporated treatment was 5.4 times as large as that of litter covered treatment, which was attributed to the changes in surface litter coverage and soil erosion resistance. Litter type affected soil erosion through the variations in litter coverage and litter morphology. For litter covering treatment, litter coverage can explain the major variance of soil loss on the slopes. Whereas, for litter incorporated treatment, both the influences of litter coverage and litter length on soil erosion resistance were considered necessary to well explain the variance of soil loss. The results also showed that the benefits of litter to control soil erosion declined with rainfall intensity and slope gradient for both covering and incorporated treatments. The results of this study are helpful to understand the mechanisms of litter influencing hydrological and erosion processes on hillslopes.     

Sediment detachment in piping-prone soils: Cohesion sources and potential weakening mechanisms

Soil piping is a widespread land degradation process that may lead to gully formation. However, the processes involved in sediment detachment from soil pipe walls have not been well studied, although their recognition is a crucial step to protect soils from piping erosion. This study aims to recognize the factors affecting cohesion and to identify the mechanisms which are likely to be responsible for the disintegration of soil. The study has been conducted in mid-altitude mountains under a temperate climate (the Bieszczady Mountains, the Carpathians, SE Poland). The research was based on the detailed field and laboratory analyses of morphology, and the physical and chemical properties of soil profiles with and without soil pipes. Moreover, experiments with flooding the undisturbed soil samples using different solutions (deionized water, ammonium oxalate, dithionate citrate, 35% hydrochloric acid and 30% hydrogen peroxide) were conducted in order to check the role of air slaking, the removal of soil organic carbon (SOC), and Fe and Al oxides on sediment detachment. The obtained results have confirmed that soil pipes develop in quite cohesive soils (silt loams), which allow the formation and maintenance of pipes with a diameter up to 30 cm. Soil cohesion, and thus susceptibility to piping, are impacted by the content of major oxides, soil particle size distribution, biological activity and porosity. The tested soils affected by piping erosion have a lower content of Al₂O₃ and Fe₂O₃, and free Fe (Fe_(DCB)), lower clay content, higher biological activity (more roots and animal burrows), higher porosity, and more and larger pores than the profile without soil pipes. The experiments have indicated that especially SOC along with Fe and Al oxides are an important cohesion source in the study area. This suggests that the removal of SOC, and Fe and Al oxides may weaken and disintegrate aggregates in soil pipes. Further study of soil leaching and tensile strength will broaden understanding of which chemical processes control where pipes will develop in other cohesive piping-prone soils. © 2020 John Wiley & Sons, Ltd.     

Biotic controls on shallow translational landslides

In undisturbed tropical montane rainforests massive organic layers accommodate the majority of roots and only a small fraction of roots penetrate the mineral soil. We investigated the contribution of vegetation to slope stability in such environments by modifying a standard model for slope stability to include an organic layer with distinct mechanical properties. The importance of individual model parameters was evaluated using detailed measurements of soil and vegetation properties to reproduce the observed depth of 11 shallow landslides in the Andes of southern Ecuador. By distinguishing mineral soil, organic layer and above‐ground biomass, it is shown that in this environment vegetation provides a destabilizing effect mainly due to its contribution to the mass of the organic layer (up to 973 t ha^{? 1} under wet conditions). Sensitivity analysis shows that the destabilizing effect of the mass of soil and vegetation can only be effective on slopes steeper than 37.9°. This situation applies to 36% of the study area. Thus, on the steep slopes of this megadiverse ecosystem, the mass of the growing forest promotes landsliding, which in turn promotes a new cycle of succession. This feedback mechanism is worth consideration in further investigations of the impact of landslides on plant diversity in similar environments. Copyright © 2012 John Wiley & Sons, Ltd.     

Quantifying erosional equilibrium across a slowly eroding,soil mantled landscape

The textbook concept of an equilibrium landscape, which posits that soil production and erosion are balanced and equal channel incision, is rarely quantified for natural systems. In contrast to mountainous, rapidly eroding terrain, low relief and slow-eroding landscapes are poorly studied despite being widespread and densely inhabited. We use three field sites along a climosequence in South Africa to quantify very slow (2-5 m/My) soil production rates that do not vary across hillslopes or with climate. We show these rates to be indistinguishable from spatially invariant catchment-average erosion rates while soil depth and chemical weathering increase strongly with rainfall across our sites. Our analyses imply landscape-scale equilibrium although the dominant means of denudation varies from physical weathering in dry climates to chemical weathering in wet climates. In the two wetter sites, chemical weathering is so significant that clay translocates both vertically in soil columns and horizontally down hillslope catenas, resulting in particle size variation and the accumulation of buried stone lines at the clay-rich depth. We infer hundred-thousand-year residence times of these stone lines and suggest that bioturbation by termites plays a key role in exhuming sediment into the mobile soil layer from significant depths below the clay layer. Our results suggest how tradeoffs in physical and chemical weathering, potentially modulated by biological processes, shape slowly eroding, equilibrium landscapes. © 2019 John Wiley & Sons, Ltd.     

Modelling to analyse the impacts of animal treading effects on soil infiltration

Animal treading can change soil physical properties, and thus is an important factor in hydrological modelling. We investigated the impacts of animal treading on infiltration by using a series of rainfall simulation experiments at Whatawhata Research Center, Waikato, New Zealand. The study identified significant variables for estimating soil steady‐state infiltration at a micro‐site (0·5 m²) and fitted the Green and Ampt equation by modifying or including variables for soil and water parameters and animal activities on grazing paddocks. A regression function for estimating steady‐state infiltration rate was created for each of four scenarios: between tracks (inter‐track), track, easy slope with ash soil, and easy slope with clay soil. Significant variables included the number of days after treading, antecedent soil moisture, field capacity, percentage of bare ground, bulk density, and the high degree of soil damage (damage not compacted). Regression models explained more than 71% of the variance in steady‐state infiltration for three scenarios, but only 53% for the easy slope with clay soil. The remodified Green and Ampt equation provided satisfactory estimation of infiltration for all scenarios (accuracy > 80%), and thus enables us to use the modified model for Waikato hill country pastures of different topography, soil physical condition, season and grazing management. Copyright © 2006 John Wiley & Sons, Ltd.     

Petrophysics and mineral exploration: a workflow for data analysis and a new interpretation framework

As mineral exploration seeks deeper targets, there will be a greater reliance on geophysical data and a better understanding of the geological meaning of the responses will be required, and this must be achieved with less geological control from drilling. Also, exploring based on the mineral system concept requires particular understanding of geophysical responses associated with altered rocks. Where petrophysical datasets of adequate sample size and measurement quality are available, physical properties show complex variations, reflecting the combined effects of various geological processes. Large datasets, analysed as populations, are required to understand the variations. We recommend the display of petrophysical data as frequency histograms because the nature of the data distribution is easily seen with this form of display. A petrophysical dataset commonly contains a combination of overlapping sub-populations, influenced by different geological factors. To understand the geological controls on physical properties in hard rock environments, it is necessary to analyse the petrophysical data not only in terms of the properties of different rock types. It is also necessary to consider the effects of processes such as alteration, weathering, metamorphism and strain, and variables such as porosity and stratigraphy. To address this complexity requires that much more supporting geological information be acquired than in current practice. The widespread availability of field portable instruments means quantitative geochemical and mineralogical data can now be readily acquired, making it unnecessary to rely primarily on categorical rock classification schemes. The petrophysical data can be combined with geochemical, petrological and mineralogical data to derive explanations for observed physical property variations based not only on rigorous rock classification methods, but also in combination with quantitative estimates of alteration and weathering. To understand how geological processes will affect different physical properties, it is useful to define three end-member forms of behaviour. Bulk behaviour depends on the physical properties of the dominant mineral components. Density and, to a lesser extent, seismic velocity show such behaviour. Grain and texture behaviour occur when minor components of the rock are the dominate controls on its physical properties. Grain size and shape control grain properties, and for texture properties the relative positions of these grains are also important. Magnetic and electrical properties behave in this fashion. Thinking in terms of how geological processes change the key characteristics of the major and minor mineralogical components allows the resulting changes in physical properties to be understood and anticipated.     

Determining short-term changes in the hydraulic properties of a sandy-loam soil by a three-run infiltration experiment

ABSTRACT

Soil structure-dependent parameters can vary rapidly as a consequence of perturbing events such as intense rainfall. Investigating their short-term changes is therefore essential to understand the general behaviour of a porous medium. The aim of this study is to gain insight into the effects of wetting, perturbation and recovery processes through different sequences of Beerkan infiltration experiments performed on a sandy-loam soil. Two different three-run infiltration experiments (LHL and LLL) were carried out by pouring water at low (L, non-perturbing) and high (H, perturbing) heights above the soil surface and at short time intervals (hours, days). The results demonstrate that the proposed method allows one to capture short-term variations in soil structure-dependent parameters. The developed methodology is expected to simplify the parameterization of hydrological models with temporally variable soil hydraulic properties.     

An examination of the steady‐state assumption in soil development models with application to landscape evolution

Soil depth and soil production are highly complicated phenomena, generated from a complex interaction of physical, biological and chemical processes. It has, nevertheless, become increasingly clear that soil formation rates are closely related to chemical weathering rates. Somewhat paradoxically, it is likewise becoming apparent that such biogeochemical reactions as slowly transform rock to soil are limited by physical processes, such as flowing water and the formation of fractures. We have formulated a theoretical approach that relates soil formation rates to chemical weathering rates, and those, likewise, to solute transport rates. For such a theoretical framework to be relevant, the solute transport rates cannot equal those of the flowing water, as is the case in Gaussian solute transport. Rather, solute transport must be slowed in accordance with heavy‐tailed solute arrival time distributions. The inference is that the traditional advection–dispersion equation formulation for solute transport is inadequate in the typically heterogeneous geological media that weather to form soils. Here we examine the implications of this soil production model on the assumption of the approach to steady state. Particularly at slow erosion rates we find that many soil columns are not in equilibrium. This tendency may be accentuated in dry climates. Copyright © 2017 John Wiley & Sons, Ltd.     

Element geochemistry of weathering profile of dolomitite and its implications for the average chemical composition of the upper-continental crust--Case studies from the Xinpu profile,northern Guizhou Province,China

Geochemical behavior of chemical elements is studied in a dolomitite weathering profile in upland of karst terrain in northern Guizhou.Two stages can be recognized during the process of in situ weathering of dolomitite:the stage of sedentary accumulation of leaching residue of dolomitite and the stage of chemical weathering evolution of sedentary soil.Ni,Cr,Mo,W and Ti are the least mobile elements with reference to Al.The geochemical behavior of REE is similar to that observed in weathering of other types of rocks.Fractionation of REE is noticed during weathering,and the two layers of REE enrichments are thought to result from downward movement of the weathering front in response to changes in the environment.It is considered that the chemistry of the upper part of the profile,which was more intensively weathered,is representative of the mobile components of the upper curst at the time the dolomitite was formed,while the less weathered lower profile is chemically representative of the immobile constitution.Like glacial till and loess,the "insoluble" materials in carbonate rocks originating from chemical sedimentation may also provide valuable information about the average chemical composition of the upper continental crust.     

Effects of physical soil crusts on infiltration and splash erosion in three typical Chinese soils

Physical soil crusts likely have significant effects on infiltration and soil erosion, however, little is known on whether the effects of the crusts change during a rainfall event. Further, there is a lack of discussions on the differences among the crusting effects of different soil types. The objectives of this study are as follows： （i） to study the effects of soil crusts on infiltration, runoff, and splash erosion using three typical soils in China, （ii） to distinguish the different effects on hydrology and erosion of the three soils and discuss the primary reasons for these differences, and （iii） to understand the variations in real soil shear strength of the three soils during rainfall events and mathematically model the effects of the crusts on soil erosion. This study showed that the soil crusts delayed the onset of infiltration by 5 to 15 min and reduced the total amount of infiltration by 42.9 to 53.4% during rainfall events. For a purple soil and a loess soil, the initial crust increased the runoff by 2.8% and 3.4%, respectively, and reduced the splash erosion by 3.1% and 8.9%, respectively. For a black soil, the soil crust increased the runoff by 42.9% and unexpectedly increased the splash erosion by 95.2%. In general, the effects of crusts on the purple and loess soils were similar and negligible, but the effects were significant for the black soil. The soil shear strength decreased dynamically and gradually during the rainfall events, and the values of crusted soils were higher than those of incrusted soils, especially during the early stage of the rainfall. Mathematical models were developed to describe the effects of soil crusts on the splash erosion for the three soils as follows： purple soil, Fc= 0.002t- 0.384 ; black soil, Fc. =-0.022t ＋ 3.060 ; and loess soil, Fc = 0.233 In t- 1.239 . Combined with the equation Rc= Fc （Ruc - 1）, the splash erosion of the crusted soil can be predicted over time.     

Effect of a thermokarst lake on soil physical properties and infiltration processes in the permafrost region of the Qinghai-Tibet Plateau,China

Changes in the hydrological processes in alpine soil constitute one of the several key problems encountered with studying watershed hydrology and ecosystem stability against the background of global warming. A typically developing thermokarst lake was chosen as a subject for a study using model simulation based on observations of soil physical properties, infiltration processes, and soil moisture. The results showed that the selected thermokarst lake imposed certain changes on the soil infiltration processes and, with the degree of impact intensifying, the initial infiltration rate decreased. The greatest reduction was achieved in the area of moderate impact. However, the stable infiltration rate and cumulative infiltration gradually increased in the surface layer at a depth of 10 and 20 cm, both decreasing initially and then increasing, which is correlated significantly with soil textures. Moreover, the cumulative infiltration changed in line with steady infiltration rate. Based on a comparative analysis, the Horton model helps better understand the effect on the soil infiltration processes of the cold alpine meadow close to the chosen thermokarst lake. In conclusion, the formation of the thermokarst lake reduced the water holding capacity of the alpine meadow soil and caused the hydraulic conductivity to increase, resulting in the reduction of runoff capacity in the area of the thermokarst lake.     

Effects of the roots of Cynodon dactylon and Schefflera heptaphylla on water infiltration rate and soil hydraulic conductivity

Water infiltration rate and hydraulic conductivity in vegetated soil are two vital hydrological parameters for agriculturists to determine availability of soil moisture for assessing crop growths and yields, and also for engineers to carry out stability calculations of vegetated slopes. However, any effects of roots on these two parameters are not well‐understood. This study aims to quantify the effects of a grass species, Cynodon dactylon, and a tree species, Schefflera heptaphylla, on infiltration rate and hydraulic conductivity in relation to their root characteristics and suction responses. The two selected species are commonly used for ecological restoration and rehabilitation in many parts of the world and South China, respectively. A series of in‐situ double‐ring infiltration tests was conducted during a wet summer, while the responses of soil suction were monitored by tensiometers. When compared to bare soil, the vegetated soil has lower infiltration rate and hydraulic conductivity. This results in at least 50% higher suction retained in the vegetated soil. It is revealed that the effects of root‐water uptake by the selected species on suction were insignificant because of the small evapotranspiration (<0.2 mm) when the tests were conducted under the wet climate. There appears to have no significant difference (less than 10%) of infiltration rates, hydraulic conductivity and suction retained between the grass‐covered and the tree‐covered soil. However, the grass and tree species having deeper root depth and greater Root Area Index (RAI) retained higher suction. Copyright © 2015 John Wiley & Sons, Ltd.     

Element geochemistry of weathering profile of dolomitite and its implications for the average chemical composition of the upper-continental crust

Geochemical behavior of chemical elements is studied in a dolomitite weathering profile in upland of karst terrain in northern Guizhou. Two stages can be recognized during the process ofin situ weathering of dolomitite: the stage of sedentary accumulation of leaching residue of dolomitite and the stage of chemical weathering evolution of sedentary soil. Ni, Cr, Mo, W and Ti are the least mobile elements with reference to Al. The geochemical behavior of REE is similar to that observed in weathering of other types of rocks. Fractionation of REE is noticed during weathering, and the two layers of REE enrichments are thought to result from downward movement of the weathering front in response to changes in the environment. It is considered that the chemistry of the upper part of the profile, which was more intensively weathered, is representative of the mobile components of the upper curst at the time the dolomitite was formed, while the less weathered lower profile is chemically representative of the immobile constitution. Like glacial till and loess, the “insoluble” materials in carbonate rocks originating from chemical sedimentation may also provide valuable information about the average chemical composition of the upper continental crust.     

Water uptake of riparian plants in the lower Lhasa River Basin,South Tibetan Plateau using stable water isotopes

Riparian plants can adapt their water uptake strategies based on climatic and hydrological conditions within a river basin. The response of cold-alpine riparian trees to changes in water availability is poorly understood. The Lhasa River is a representative cold-alpine river in South Tibet and an under-studied environment. Therefore, a 96 km section of the lower Lhasa River was selected for a study on the water-use patterns of riparian plants. Plant water, soil water, groundwater and river water were measured at three sites for δ¹⁸O and δ²H values during the warm-wet and cold-dry periods in 2018. Soil profiles differed in isotope values between seasons and with the distance along the river. During the cold-dry period, the upper parts of the soil profiles were significantly affected by evaporation. During the warm-wet period, the soil profile was influenced by precipitation infiltration in the upper reaches of the study area and by various water sources in the lower reaches. Calculations using the IsoSource model indicated that the mature salix and birch trees (Salix cheilophila Schneid. and Betula platyphylla Suk.) accessed water from multiple sources during the cold-dry period, whereas they sourced more than 70% of their requirement from the upper 60–80 cm of the soil profile during the warm-wet period. The model indicated that the immature rose willow tree (Tamarix ramosissima Ledeb) accessed 66% of its water from the surface soil during the cold-dry period, but used the deeper layers during the warm-wet period. The plant type was not the dominant factor driving water uptake patterns in mature plants. Our findings can contribute to strategies for the sustainable development of cold-alpine riparian ecosystems. It is recommended that reducing plantation density and collocating plants with different rooting depths would be conducive to optimal plant growth in this environment.     

Effects of plant roots on the soil erosion rate under simulated rainfall with high kinetic energy

ABSTRACT

This study examined the effects of herbaceous plant roots on interrill erosion using two herbaceous species: clover (Trifolium repens) and oats (Avena sativa). We developed a simple rainfall simulator with relatively high normalized kinetic energy (KE; 23.2 J m^?2 mm^?1). Under simulated rainfall, we measured eroded soil for 42 boxes with various amounts of aboveground and belowground biomass. Aboveground vegetation had a significant effect on the soil erosion rate (SER). We found a clear negative relationship between the percent vegetation cover (c) and the SER. In contrast, plant roots showed no effects on the SER. The SER was not significantly different between the boxes with and without plant roots under similar c conditions. Thus, plant roots could have less of an effect on the SER under higher KE conditions.

Editor M.C. Acreman Associate editor N. Verhoest