The use of non-reflecting boundary conditions in soil dynamics

In this paper some aspects of the use of non-reflecting boundaries in dynamic problems, analyzed in time domain, are considered. Current trends for treating the above mentioned problems are summarized with a particular emphasis on the use of numerical techniques, such as Boundary Element Method (BEM) or mixed and hybrid formulations, Finite Element Method (FEM) plus BEM. As an alternative to these methods, an easy time domain boundary condition, obtained from the well known consistent transmitting boundary developed by Waas for frequency domain analysis, can be applied to represent the reactions of the unbounded soil on the interest zone. The behaviour of this proposed boundary condition is studied when waves of different frequency to the one used for its obtention are acting on the physical edge of the model. As an application example,an analysis is made of the soil-structure interaction of a rigid strip foundation on a horizontal non-linear elastic layer on bed rock. The results obtained suggest the need of time domain solutions for this type of problem

Modeling soil system: complexity under your feet

A clear statement in these lines textually cited (Byers et al., 1938) defines the framework of this special issue: “True soil is the product of the action of climate and living organism upon the parent material, as conditioned by the local relief. The length of time during which these forces are operative is of great importance in determining the character of the ultimate product. Drainage conditions are also important and are controlled by local relief, by the nature of the parent material or underlying rock strata, or by the amount of precipitation in relation to rate of percolation and runoff water. There are, therefore, five principal factors of soil formation: Parent material, climate, biological activity, relief and time. These soil forming factors are interdependent, each modifying the effectiveness of the others.” Due to these various processes associated to its formation and genesis soil dynamics reveals high complexity that creates several levels of structure using this term in a broad sense

The direct boundary element method: 2D site effects assessment on laterally varying layered media (methodology)

The Direct Boundary Element Method (DBEM) is presented to solve the elastodynamic field equations in 2D, and a complete comprehensive implementation is given. The DBEM is a useful approach to obtain reliable numerical estimates of site effects on seismic ground motion due to irregular geological configurations, both of layering and topography. The method is based on the discretization of the classical Somigliana's elastodynamic representation equation which stems from the reciprocity theorem. This equation is given in terms of the Green's function which is the full-space harmonic steady-state fundamental solution. The formulation permits the treatment of viscoelastic media, therefore site models with intrinsic attenuation can be examined. By means of this approach, the calculation of 2D scattering of seismic waves, due to the incidence of P and SV waves on irregular topographical profiles is performed. Sites such as, canyons, mountains and valleys in irregular multilayered media are computed to test the technique. The obtained transfer functions show excellent agreement with already published results.

Permanent deformation estimates of dynamic equipment foundations: Application to a gas turbine in granular soils

Permanent displacements of a gas turbine founded on a fine, poorly graded, and medium density sand are studied. The amplitudes and modes of vibration are computed using Barkan´s formulation, and the “High-Cycle Accumulation” (HCA) model is employed to account for accumulated deformations due to the high number of cycles. The methodology is simple: it can be easily incorporated into standard mathematical software, and HCA model parameters can be estimated based on granulometry and index properties. Special attention is devoted to ‘transient’ situations at equipment´s start-up, during which a range of frequencies – including frequencies that could be similar to the natural frequencies of the ground – is traversed. Results show that such transient situations could be more restrictive than stationary situations corresponding to normal operation. Therefore, checking the stationary situation only might not be enough, and studying the influence of transient situations on computed permanent displacements is needed to produce a proper foundation design

Effect of the Legume Plant Host on Rhizobial Soil Population Dynamics

Rhizobium leguminosarum bv viciae (Rlv) is a soil bacterium able to establish specific root-nodule symbioses with legumes of four different genera: Pisum, Vicia, Lens and Lathyrus. Rlv isolates from nodules of any of these legumes can nodulate any of them; however, it has been shown that plants select specific rhizobial genotypes from those present in the soil (1,2). We have previously shown this at the genomic level by following a population genomics approach. Pool genomic sequences from 100 isolates from each of four plant species: P. sativum, L. culinaris, V. faba and V. sativa, show different, specific profiles at the single nucleotide polymorphism (SNP) level for relevant genes. In this work, the extent of Rlv selection from a well-characterized soil population by different legume plant hosts: P. sativum, L. culinaris, V. faba and V. sativa, after a medium-term mesocosm study is described. Direct soil isolates from each of these mesocosm studies have been tested for specific rhizobial genes (glnII and fnrN) and symbiotic genes (nodC and nifH). Different populations were characterized further by Sanger sequencing of both the rpoB phylogenetic marker gene and the symbiotic genes nodC and nifH. The distribution and size of the rhizobial population for each legume host showed changes during the medium-term mesocosm study. Particularly, a non-symbiotic group of rhizobia was enriched by all four hosts, in contrast to the symbiotic rhizobia profile, which was specific for each legume plant host.

Soil compaction and vegetation cover in a Scots pine stand at the Mediterranean rangelands

Right development of ROOT SYSTEMS is essential to ensure seedling survival in the initial stages of natural regeneration processes. Soil compaction determines this development both because of its influence on soil Tª & moisture dynamics and for its direct effect on soil mechanical impedance to root growth. All this effects can be assessed as a whole through soil penetration resistance (Soil Strength) measurements. SOIL STRENGTH has been usually evaluated in forest research in connection with severe disturbances derived from heavy machinery works during forest operations. Nevertheless, undisturbed soils are also expected to show different levels of compaction for root development. Organic matter modifies soil structure and so on porosity, compaction and resultant soil resistance to penetration. Its concentration in surface layers is rather related to vegetation cover composition and density. So within forest stands, a relationship is expected to be found between VEGETATION COVER density and compaction measured as resistance to penetration (soil strength)

Community Structurein Soil Porous System.

Soil is well recognized as a highly complex system. The interaction and coupled physical, chemical, and biological processes and phenomena occurring in the soil environment at different spatial and temporal scales are the main reasons for such complexity. There is a need for appropriate methodologies to characterize soil porous systems with an interdisciplinary character. Four different real soil samples, presenting different textures, have been modeled as heterogeneous complex networks, applying a model known as the heterogeneous preferential attachment. An analytical study of the degree distributions in the soil model shows a multiscaling behavior in the connectivity degrees, leaving an empirically testable signature of heterogeneity in the topology of soil pore networks. We also show that the power-law scaling in the degree distribution is a robust trait of the soil model. Last, the detection of spatial pore communities, as densely connected groups with only sparser connections between them, has been studied for the first time in these soil networks. Our results show that the presence of these communities depends on the parameter values used to construct the network. These findings could contribute to understanding the mechanisms of the diffusion phenomena in soils, such as gas and water diffusion, development and dynamics of microorganisms, among others.

Nitrogen uptake dynamics, yield and quality as influenced by nitrogen fertilization in "Piel de sapo" melon.

The need to reduce nitrogen (N) fertilizer pollution strengthens the importance of improving the utilization efficiency of applied N to crops. This requires knowledge of crop N uptake characteristics and how fertilization management affects it. A three-year field experiment was conducted from May to September in central Spain to investigate the influence of different N rates, which ranged from 11 to 393 kg ha-1, applied through drip irrigation, on the dynamics of N uptake, nitrogen use efficiency (NUE), fruit yield and quality of a ?Piel de sapo? melon crop (Cucumis melo L. cv. Sancho). Both N concentration and N content increased in different plant parts with the N rate. Leaves had the highest N concentration, which declined by 40-50% from 34-41 days after transplanting (DAT), while the highest N uptake rate was observed from 30-35 to 70-80 DAT, coinciding with fruit development. In each year, NUE declined with increasing N rate. With N fertilizer applications close to the optimum N rate of 90-100 kg ha-1, the fruits removed approximately 60 kg N ha-1, and the amount of N in the crop residue was about 80 kg N ha-1; this serves to replenish the organic nutrient pool in the soil and may be used by subsequent crops following mineralization.

Modeling GHG emissions, N and C dynamics in Spanish agricultural soils.

To date, only few initiatives have been carried out in Spain in order to use mathematical models (e.g. DNDC, DayCent, FASSET y SIMSNIC) to estimate nitrogen (N) and carbon (C) dynamics as well as greenhouse gases (GHG) in Spanish agrosystems. Modeling at this level may allow to gain insight on both the complex relationships between biological and physicochemical processes, controlling the processes leading to GHG production and consumption in soils (e.g. nitrification, denitrification, decomposing, etc.), and the interactions between C and N cycles within the different components of the continuum plant-soil-environment. Additionally, these models can simulate the processes behind production, consumition and transport of GHG (e.g. nitrous oxide, N2O, and carbon dioxide, CO2) in the short and medium term and at different scales. Other sources of potential pollution from soils can be identified and quantified using these process-based models (e.g. NO3 y NH3).

Initializing the DSSAT-Century model: inverse calibration of carbon pools from apparent soil N mineralization

The CENTURY soil organic matter model was adapted for the DSSAT (Decision Support System for Agrotechnology Transfer), modular format in order to better simulate the dynamics of soil organic nutrient processes (Gijsman et al., 2002). The CENTURY model divides the soil organic carbon (SOC) into three hypothetical pools: microbial or active material (SOC1), intermediate (SOC2) and the largely inert and stable material (SOC3) (Jones et al., 2003). At the beginning of the simulation, CENTURY model needs a value of SOC3 per soil layer which can be estimated by the model (based on soil texture and management history) or given as an input. Then, the model assigns about 5% and 95% of the remaining SOC to SOC1 and SOC2, respectively. The model performance when simulating SOC and nitrogen (N) dynamics strongly depends on the initialization process. The common methods (e.g. Basso et al., 2011) to initialize SOC pools deal mostly with carbon (C) mineralization processes and less with N. Dynamics of SOM, SOC, and soil organic N are linked in the CENTURY-DSSAT model through the C/N ratio of decomposing material that determines either mineralization or immobilization of N (Gijsman et al., 2002). The aim of this study was to evaluate an alternative method to initialize the SOC pools in the DSSAT-CENTURY model from apparent soil N mineralization (Napmin) field measurements by using automatic inverse calibration (simulated annealing). The results were compared with the ones obtained by the iterative initialization procedure developed by Basso et al., 2011.

Soil aggregate geometry: Measurements and morphology

Aggregates provide physical microenvironments for microorganisms, the vital actors of soil systems, and thus play a major role as both, an arena and a product of soil carbon stabilization and dynamics. The surface of an aggregate is what enables exchange of the materials and air and water fluxes between aggregate exterior and interior regions. We made use of 3D images from X-ray CT of aggregates and mathematical morphology to provide an exhaustive quantitative description of soil aggregate morphology that includes both intra-aggregate pore space structure and aggregate surface features. First, the evolution of Minkowski functionals (i.e. volume, boundary surface, curvature and connectivity) for successive dilations of the solid part of aggregates was investigated to quantify its 3D geometrical features. Second, the inner pore space was considered as the object of interest. We devised procedures (a) to define the ends of the accessible pores that are connected to the aggregate surface and (b) to separate accessible and inaccessible porosity. Geometrical Minkowski functionals of the intra-aggregate pore space provide the exhaustive characterization of the inner structure of the aggregates. Aggregates collected from two different soil treatments were analyzed to explore the utility of these morphological tools in capturing the impact on their morphology of two different soil managements, i.e. conventional tillage management, and native succession vegetation treatment. The quantitative tools of mathematical morphology distinguished differences in patterns of aggregate structure associated to the different soil managements.