Modelling raindrop impact and splash erosion processes within a spatial cell: a stochastic approach

A new approach is proposed to simulate splash erosion on local soil surfaces. Without the effect of wind and other raindrops, the impact of free-falling raindrops was considered as an independent event from the stochastic viewpoint. The erosivity of a single raindrop depending on its kinetic energy was computed by an empirical relationship in which the kinetic energy was expressed as a power function of the equivalent diameter of the raindrop. An empirical linear function combining the kinetic energy and soil shear strength was used to estimate the impacted amount of soil particles by a single raindrop. Considering an ideal local soil surface with size of I m x I m, the expected number of received free-failing raindrops with different diameters per unit time was described by the combination of the raindrop size distribution function and the terminal velocity of raindrops. The total splash amount was seen as the sum of the impact amount by all raindrops in the rainfall event. The total splash amount per unit time was subdivided into three different components, including net splash amount, single impact amount and re-detachment amount. The re-detachment amount was obtained by a spatial geometric probability derived using the Poisson function in which overlapped impacted areas were considered. The net splash amount was defined as the mass of soil particles collected outside the splash dish. It was estimated by another spatial geometric probability in which the average splashed distance related to the median grain size of soil and effects of other impacted soil particles and other free-falling raindrops were considered. Splash experiments in artificial rainfall were carried out to validate the availability and accuracy of the model. Our simulated results suggested that the net splash amount and re-detachment amount were small parts of the total splash amount. Their proportions were 0.15% and 2.6%, respectively. The comparison of simulated data with measured data showed that this model could be applied to simulate the soil-splash process successfully and needed information of the rainfall intensity and original soil properties including initial bulk intensity, water content, median grain size and some empirical constants related to the soil surface shear strength, the raindrop size distribution function and the average splashed distance. Copyright (c) 2007 John Wiley & Sons, Ltd.

Explicit representation of subgrid heterogeneity in a GCM land-surface scheme

Essery, R L H, Best, M J, Betts, R A, Cox, P M & Taylor, C M, Explicit representation of subgrid heterogeneity in a GCM land-surface scheme. Journal of Hydrometeorology 4, pp 530-543 (2003).

Isotopic Assessment of Sources of Surface Water Nitrate within the Oldman River Basin, Southern Alberta, Canada

Concentrations and isotopic compositions of NO-3 from the Oldman River (OMR) and some of its tributaries (Alberta, Canada) have been determined on a monthly basis since December 2000 to assess temporal and spatial variations of riverine NO-3 sources within the OMR basin. For the OMR sites, NO-3 -N concentrations reached up to 0.34 mg L-1, d15N-NO-3 values varied between –0.3 and +13.8‰, and d18O-NO-3 values ranged from –10.0 to +5.7‰. For the tributary sites, NO-3 -N concentrations were as high as 8.81 mg L-1, d15N-NO-3 values varied between –2.5 and +23.4‰, and d18O-NO-3 values ranged from –15.2 to +3.4‰. Tributaries in the western, relatively pristine forested part of the watershed add predominantly NO-3 to the OMR with d15N-NO-3 indicative of soil nitrification. In contrast, tributaries in the eastern agriculturally-urban-industrially-used part of the basin contribute NO-3 with d15N-NO-3 values of about +16‰ indicative of manure and/or sewage derived NO-3. This difference in d15N-NO-3 values of tributaries was found to be independent of the season, but rather indicates a spatial change in the NO-3 source, which correlates with land use changes within the OMR basin. As a consequence of tributary influx, d15N-NO-3 values in the Oldman River increased from +6‰ in the downstream direction (W to E), although [NO-3 -N] increased only moderately (generally

Response of surface horizons in an oak forest to prescribed burning

The reduction of forest floor ground cover and litter layers by prescribed fires may alter the morphology (field and micro) and physical properties of surface horizons. This study determined long-term (35 yr) changes in surface horizon bulk density, organic matter concentration and content, and morphology in response to periodic (5 yr) and annual (1 yr) prescribed fires. Soils were fine-silty, siliceous, thermic Glossic Fragiuldults, supporting mixed oak vegetation in middle Tennessee. Upper mineral soils (0- to 2-cm and 0- to 7.6-cm depths) were sampled and detailed field descriptions made. Periodic and control plots had a thin layer of Oi, Oe, and Oa horizons 5 yr after the 1993 burn, whereas on annual burn plots a 1- to 2-cm charred layer was present. Significant reductions in organic matter concentration and mean thickness of the A horizon were found from burning (A horizons thicknesses were 6.4, 4.6, and 2.9 cm in control, periodic, and annual plots, respectively). Periodic burns did not significantly alter the organic matter and bulk density of the upper 7.6 cm of mineral soil; however, annual burns did result in significantly higher bulk densities (1.01, 1.07, and 1.29 Mg m-3 in control, periodic, and annual plots, respectively) and lower organic matter concentrations and contents. Microscopic investigations confirmed that compaction was increased from annual burning. Thin sections also revealed that the granular structure of the A horizons in control and periodic plots resulted from bioterbation of macro and mesofauna, fungi, and roots. Long-term annual burning greatly affected surface soil properties, whereas periodic burning on a 5-yr cycle had only limited effects.

Dynamical instability in surface permeability characteristics of building sandstones in response to salt accumulation over time

This paper explores how the surface permeability of sandstone blocks changes over time in response to repeated salt weathering cycles. Surface permeability controls the amount of moisture and dissolved salt that can penetrate in and facilitate decay. Connected pores permit the movement of moisture (and hence soluble salts) into the stone interior, and where areas are more or less permeable soluble salts may migrate along preferred pathways at differential rates. Previous research has shown that salts can accumulate in the near-surface zone and lead to partial pore blocking which influences subsequent moisture ingress and causes rapid salt accumulation in the near-surface zone.

Two parallel salt weathering simulations were carried out on blocks of Peakmoor Sandstone of different volumes. Blocks were removed from simulations after 2, 5, 10, 20 and 60 cycles. Permeability measurements were taken for these blocks at a resolution of 20 mm, providing a grid of 100 permeability values for each surface. The geostatistical technique of ordinary kriging was applied to the data to produce a smoothed interpolation of permeability for these surfaces, and hence improve understanding of the evolution of permeability over time in response to repeated salt weathering cycles.

Results illustrate the different responses of the sandstone blocks of different volumes to repeated salt weathering cycles. In both cases, after an initial subtle decline in the permeability (reflecting pore blocking), the permeability starts to increase — reflected in a rise in mean, maximum and minimum values. However, between 10 and 20 cycles, there is a jump in the mean and range permeability of the group A block surfaces coinciding with the onset of meaningful debris release. After 60 cycles, the range of permeability in the group A block surface had increased markedly, suggesting the development of a secondary permeability. The concept of dynamic instability and divergent behaviour is applied at the scale of a single block surface, with initial small-scale differences across a surface having larger scale consequences as weathering progresses.

After cycle 10, group B blocks show a much smaller increase in mean permeability, and the range stays relatively steady — this may be explained by the capillary conditions set up by the smaller volume of the stone, allowing salts to migrate to the ‘back’ of the blocks and effectively relieving stress at the ‘front’ face.

Dinoflagellate cyst-based reconstructions of mid to late Holocene winter sea-surface temperature and productivity from an anoxic fjord in the NE Pacific Ocean

Published contemporary dinoflagellate distributional data from the NE Pacific margin and estuarine environments (n = 136) were re-analyzed using Canonical Correspondence Analysis (CCA) and partial Canonical Correspondence Analysis (pCCA). These analyses illustrated the dominant controls of winter temperature and productivity on the distribution of dinoflagellate cysts in this region. Dinoflagellate cyst-based predictive models for winter temperature and productivity were developed from the contemporary distributional data using the modern analogue technique and applied to subfossil data from two mid to late Holocene (~5500 calendar years before present–present) cores; TUL99B03 and TUL99B11, collected from Effingham Inlet, a 15 km long anoxic fjord located on the southwest coast of Vancouver Island that directly opens to the Pacific Ocean through Barkley Sound. Sedimentation within these basins largely comprises annually deposited laminated couplets, each made up of a winter deposited terrigenous layer and spring to fall deposited diatomaceous layer. The Effingham Inlet dinoflagellate cyst record provides evidence of a mid-Holocene gradual decline in winter SST, ending with the initiation of neoglacial advances in the region by ~3500 cal BP. A reconstructed Late Holocene increase in winter SST was initiated by a weakening of the California Current, which would have resulted in a warmer central gyre and more El Niño-like conditions.

Near-surface temperature cycling of stone and its implications for scales of surface deterioration

Many previous studies into internal temperature gradients within stone have assumed smooth, exponential increases and decreases in sub-surface temperatures in response, for example, to diurnal patterns of heating and cooling and these have been used to explain phenomena such as large-scale contour scaling. This high-resolution experimental study, in which a porous limestone block was subjected to alternate surface heating and cooling using an infrared lamp, demonstrates that internal temperature gradients in response to short-term environmental cycles (measured in minutes) can in fact be complex and inconsistent. Results confirm the significance of very steep temperature/stress gradients within the outer 10 mm or less of exposed stone. Below this the data indicate complex patterns of temperature reversals, the amplitudes of which are attenuated with depth and which are influenced in their intensity and location by variations in the relative duration of heating and cooling phases. It is suggested that the reversals might represent ‘interference patterns’ between incoming and outgoing thermal waves, but whatever their origin they are potentially important because they occur within the zone in which many stone decay processes, especially salt weathering, operate. These processes invariably respond to temperature and moisture fluctuations, and short-term interruptions to insolation could, for example, trigger these fluctuations on numerous occasions over a day. In particular, the reversals occur at a scale that is commensurate with decay by multiple flaking and could indicate an underlying control on this previously little-researched pattern of weathering. In the context of this publication, however, the main lesson to be learned from this study is that differing scales of behaviour require different scales of enquiry.

Detection of soil compaction using seismic surface waves

Seismic geophysical methods have rarely been used in precision agriculture, predominantly due to the perception that they are slow and results require a complex evaluation. This paper explores the possibility of using a recently developed surface wave seismic geophysical approach, the multichannel analysis of surface waves (MASW) method, for assessment of agricultural compaction. This approach has the advantage of being non-intrusive, rapid and is able to produce 2D ground models with a relatively high density of spatial sampling points. The method, which was tested on a research site in Oakpark, Ireland, detected a significant difference in shear wave velocity between a heavily compacted headland and an uncompacted location. The results from this approach compared favourably with those obtained
from measurements of bulk density and penetrometer resistance and demonstrate that the MASW approach can distinguish between the extreme states of heavily compacted and uncompacted soil.

Building sandstone surface modification by biofilm and iron precipitation: emerging block-scale heterogeneity and system response

Stone surfaces are sensitive to their environment. This means that they will often respond to exposure conditions by manifesting a change in surface characteristics. Such changes can be more than simply aesthetic, creating surface/subsurface heterogeneity in stone at the block scale, promoting stress gradients to be set up as surface response to, for example, temperature fluctuations, can diverge from subsurface response. This paper reports preliminary experiments investigating the potential of biofilms and iron precipitation as surface-modifiers on stone, exploring the idea of block-scale surface-to-depth heterogeneity, and investigating how physical alteration in the surface and near-surface zone can have implications for subsurface response and potentially for long-term decay patterns. Salt weathering simulations on fresh and surface-modified stone suggest that even subtle surface modification can have significant implications for moisture uptake and retention, salt concentration and distribution from surface to depth, over the period of the experimental run. The accumulation of salt may increase the retention of moisture, by modifying vapour pressure differentials and the rate of evaporation.
Temperature fluctuation experiments suggest that the presence of a biofilm can have an impact on energy transfer processes that occur at the stone surface (for example, buffering against temperature fluctuation), affecting surface-to-depth stress gradients. Ultimately, fresh and surface-modified blocks mask different kinds of system, which respond to inputs differently because of different storage mechanisms, encouraging divergent behaviour between fresh and surface modified stone over time.

A Coupled 2 × 2D Babcock–Leighton Solar Dynamo Model. I. Surface Magnetic Flux Evolution

The need for reliable predictions of the solar activity cycle motivates the development of dynamo models incorporating a representation of surface processes sufficiently detailed to allow assimilation of magnetographic data. In this series of papers we present one such dynamo model, and document its behavior and properties. This first paper focuses on one of the model's key components, namely surface magnetic flux evolution. Using a genetic algorithm, we obtain best-fit parameters of the transport model by least-squares minimization of the differences between the associated synthetic synoptic magnetogram and real magnetographic data for activity cycle 21. Our fitting procedure also returns Monte Carlo-like error estimates. We show that the range of acceptable surface meridional flow profiles is in good agreement with Doppler measurements, even though the latter are not used in the fitting process. Using a synthetic database of bipolar magnetic region (BMR) emergences reproducing the statistical properties of observed emergences, we also ascertain the sensitivity of global cycle properties, such as the strength of the dipole moment and timing of polarity reversal, to distinct realizations of BMR emergence, and on this basis argue that this stochasticity represents a primary source of uncertainty for predicting solar cycle characteristics.

A geomatics-based approach for the derivation of the spatial distribution of sediment transport processes in periglacial mountain environments

Within this paper modern techniques such as satellite image analysis and tools provided by geographic information systems (GIS.) are exploited in order to extend and improve existing techniques for mapping the spatial distribution of sediment transport processes. The processes of interest comprise mass movements such as solifluction, slope wash, dirty avalanches and rock- and boulder falls. They differ considerably in nature and therefore different approaches for the derivation of their spatial extent are required. A major challenge is addressing the differences between the comparably coarse resolution of the available satellite data (Landsat TM/ETM+, 30 in x 30 m) and the actual scale of sediment transport in this environment. A three-stepped approach has been developed which is based on the concept of Geomorphic Process Units (GPUs): parameterization, process area delineation and combination. Parameters include land cover from satellite data and digital elevation model derivatives. Process areas are identified using a hierarchical classification scheme utilizing thresholds and definition of topology. The approach has been developed for the Karkevagge in Sweden and could be successfully transferred to the Rabotsbekken catchment at Okstindan, Norway using similar input data. Copyright (C) 2008 John Wiley & Sons, Ltd.

A step toward the wet surface chemistry of glycine and alanine on Cu{110}: destabilization and decomposition in the presence of near-ambient water vapor

The coadsorption of water with organic molecules under near-ambient pressure and temperature conditions opens up new reaction pathways on model catalyst surfaces that are not accessible in conventional ultrahigh-vacuum surfacescience experiments. The surface chemistry of glycine and alanine at the water-exposed Cu{110} interface was studied in situ using ambient-pressure photoemission and X-ray absorption spectroscopy techniques. At water pressures above 10-5 Torr a significant pressure-dependent decrease in the temperature for dissociative desorption was observed for both amino acids, accompanied by the appearance of a newCN intermediate, which is not observed for lower pressures. The most likely reaction mechanisms involve dehydrogenation induced by O and/or OH surface species resulting from the dissociative adsorption of water. The linear relationship between the inverse decomposition temperature and the logarithm of water pressure enables determination of the activation energy for the surface reaction, between 213 and 232 kJ/mol, and a prediction of the decomposition temperature at the solidliquid interface by extrapolating toward the equilibrium vapor pressure. Such experiments near the equilibrium vapor pressure provide important information about elementary surface processes at the solidliquid interface, which can be retrieved neither under ultrahigh vacuum conditions nor from interfaces immersed in a solution.

Simulating dynamic crop growth with an adapted land surface model – JULES-SUCROS: model development and validation

The increasing demand for ecosystem services, in conjunction with climate change, is expected to signif- icantly alter terrestrial ecosystems. In order to evaluate the sustainability of land and water resources, there is a need for a better understanding of the relationships between crop production, land surface characteristics and the energy and water cycles. These relationships are analysed using the Joint UK Land Environment Simulator (JULES). JULES includes the full hydrological cycle and vegetation effects on the energy, water, and carbon fluxes. However, this model currently only simulates land surface processes in natural ecosystems. An adapted version of JULES for agricultural ecosystems, called JULES-SUCROS has therefore been developed. In addition to overall model improvements, JULES-SUCROS includes a dynamic crop growth structure that fully fits within and builds upon the biogeochemical modelling framework for natural vegetation. Specific agro-ecosystem features such as the development of yield-bearing organs and the phenological cycle from sowing till harvest have been included in the model. This paper describes the structure of JULES-SUCROS and evaluates the fluxes simulated with this model against FLUXNET measurements at 6 European sites. We show that JULES-SUCROS significantly improves the correlation between simulated and observed fluxes over cropland and captures well the spatial and temporal vari- ability of the growth conditions in Europe. Simulations with JULES-SUCROS highlight the importance of vegetation structure and phenology, and the impact they have on land–atmosphere interactions.