Are peatlands in different states with respect to their thermodynamic behaviour? A simple test of peatland energy and entropy budgets

Whilst all ecosystems must obey the second law of thermodynamics, these physical bounds and controls on ecosystem evolution and development are largely ignored across the ecohydrological literature. To unravel the importance of these underlying restraints on ecosystem form and function, and their power to inform our scientific understanding, we have calculated the entropy budget of a range of peat ecosystems. We hypothesize that less disturbed peatlands are ‘near equilibrium’ with respect to the second law of thermodynamics and thus respond to change by minimizing entropy production. This ‘near equilibrium’ state is best achieved by limiting evaporative losses. Alternatively, peatlands ‘far-from-equilibrium’ respond to a change in energy inputs by maximizing entropy production which is best achieved by increasing evapotranspiration. To test these alternatives this study examined the energy balance time series from seven peatlands across a disturbance gradient. We estimate the entropy budgets for each and determine how a change in net radiation (ΔR_n) was transferred to a change in latent heat flux (ΔλE). The study showed that: (i) The transfer of net radiation to latent heat differed significantly between peatlands. One group transferred up to 64% of the change in net radiation to a change in latent heat flux, while the second transferred as little as 27%. (ii) Sites that transferred the most energy to latent heat flux were those that produced the greatest entropy. The study shows that an ecosystem could be ‘near equilibrium’ rather than ‘far from equilibrium’.     

Transport of water in frozen soil IV. Analysis of experimental results on the effects of ice content

Effects of ice content on the transport of water in frozen soil are studied experimentally and theoretically under isothermal conditions. A physical law, that the flux of water in unsaturated frozen soil is proportional to the gradient of total water content is proposed. Theoretical justification is made by the use of the two-phase flow theory. The experimental results are shown to support the proposed physical law. The results of this study are presented in two parts and this is the second paper describing the theoretical aspects of the study.     

A parallel evolutionary strategy based simulation–optimization approach for solving groundwater source identification problems

Groundwater characterization involves the resolution of unknown system characteristics from observation data, and is often classified as an inverse problem. Inverse problems are difficult to solve due to natural ill-posedness and computational intractability. Here we adopt the use of a simulation–optimization approach that couples a numerical pollutant-transport simulation model with evolutionary search algorithms for solution of the inverse problem. In this approach, the numerical transport model is solved iteratively during the evolutionary search. This process can be computationally intensive since several hundreds to thousands of forward model evaluations are typically required for solution. Given the potential computational intractability of such a simulation–optimization approach, parallel computation is employed to ease and enable the solution of such problems. In this paper, several variations of a groundwater source identification problem is examined in terms of solution quality and computational performance. The computational experiments were performed on the TeraGrid cluster available at the National Center for Supercomputing Applications. The results demonstrate the performance of the parallel simulation–optimization approach in terms of solution quality and computational performance.     

A unifying framework for watershed thermodynamics: constitutive relationships

The balance equations for mass and momentum, averaged over the scale of a watershed entity, need to be supplemented with constitutive equations relating flow velocities, pressure potential differences, as well as mass and force exchanges within and across the boundaries of a watershed. In this paper, the procedure for the derivation of such constitutive relationships is described in detail. This procedure is based on the method pioneered by Coleman and Noll through exploitation of the second law of thermodynamics acting as a constraint-type relationship. The method is illustrated by its application to some common situations occurring in real world watersheds. Thermodynamically admissible and physically consistent constitutive relationships for mass exchange terms among the subregions constituting the watershed (subsurface zones, overland flow regions, channel) are proposed. These constitutive equations are subsequently combined with equations of mass balance for the subregions. In addition, constitutive relationships for forces exchanged amongst the subregions are also derived within the same thermodynamic framework. It is shown that, after linearisation of the latter constitutive relations in terms of the velocity, a watershed-scale Darcy's law governing flow in the unsaturated and saturated zones can be obtained. For the overland flow, a second order constitutive relationship with respect to velocity is proposed for the momentum exchange terms, leading to a watershed-scale Chezy formula. For the channel network REW-scale Saint–Venant equations are derived. Thus, within the framework of this approach new relationships governing exchange terms for mass and momentum are obtained and, moreover, some well-known experimental results are derived in a rigorous manner.     

A Self-tuning Robust Control System for nonlinear real-time hybrid simulation

In a real-time hybrid simulation, a transfer system is used to enforce the interface interaction between computational and physical substructures. A model-based, multilayer nonlinear control system is developed to accommodate extensive performance variations and uncertainties in a physical substructure. The aim of this work is to extend the application of real-time hybrid simulation to investigating failure, nonlinearity, and nonstationary behavior. This Self-tuning Robust Control System (SRCSys) consists of two layers: robustness and adaptation. The robustness layer synthesizes a nonlinear control law such that the closed-loop dynamics perform as intended under a broad range of parametric and nonparametric uncertainties. Sliding mode control is employed as the control scheme in this layer. Then, the adaptation layer reduces uncertainties at run time through slow and controlled learning of the control plant. The tracking performance of the SRCSys is evaluated in two experiments that have highly uncertain physical specimens.     

井间地震波衰减成像的几种方法

在井间走时成像的基础上，进行衰减成像，可以克服单一物性反演的局限性，从而更全面地认识井间研究区的地质结构与物性分布，本文列了四种具有代表性的井间地震波衰减成像方法，并就各种方法的稳定性及其在井间成像中的适用情况进行了讨论。     

3D numerical simulation of seagrass movement under waves and currents with GPUSPH

The current study tries a new approach to simulating interactions between waves and seagrass through Smoothed Particle Hydrodynamics (SPH). In this model, the plants are defined as a solid that respects Hooke's law, and are assumed to have direct interaction with the fluid. Given the characteristics of the SPH method, especially in terms of computational time, the dimensions of the simulations were limited. The first goal of the current study was to optimize the approach to avoid reaching certain limits such as the rupture of the simulated plant. Plant movements under waves and/or currents have been studied by several authors in various in-situ, physical, and numerical experiments concerning various vegetation species, thus proving that plant movements can be successfully reproduced by SPH 2D/3D. Manning's roughness coefficient, n, was calculated to confirm that the results were in accordance with what had been measured in flume studies. Even though there is still room for improvement, it is shown that this method can be used to estimate Manning's coefficient for coastal vegetation (seagrass and saltmarsh vegetation) and to greatly improve the modeling and forecasting of coastal erosion and storm surge risks by including the effects of vegetation in integrated models.     

Wave-induced dynamic response of saturated multi-layer porous media: Analytical solutions and validity regions of various formulations in non-dimensional parametric space

In this paper, dynamic response of saturated-layered porous media under harmonic waves is evaluated through a semi-analytical solution. The coupled differential equations governing the dynamics of saturated or nearly saturated porous media such as soils containing all the inertial terms of solid and fluid phases are presented for a multi-layer system. Possible simplifications of the equations which are called formulations are introduced based upon the presence of inertial terms associated with the phases. The semi-analytical solutions to the response of multiple layers for all the formulations are presented in terms of pore water pressure and stress variations considering a set of non-dimensional parameters and their respective ratios. Validity of the formulations is presented in a non-dimensional parametric space. The maximum discrepancies in the pore pressure response of the formulations leading to validity regions are illustrated for typical dynamic problems. Subsequently, the effects of layering and drainage conditions on these regions are also presented. The proposed semi-analytical solution may be served as a benchmark one for validating the coupled numerical solutions, which can be used to deal with real scientific and geo-engineering problems in the emerging field of computational geomechanics.     

A simple strategy for dynamic substructuring: I. Concept and development

Dynamic substructuring refers to physical testing with computational models in the loop. This paper presents a new strategy for such testing. The key feature of this strategy is that it decouples the substructuring controller from the physical subsystem. Unlike conventional approaches, it does not explicitly include a tracking controller. Consequently, the design and implementation of the substructuring controls are greatly simplified. This paper motivates the strategy and discusses the main concept along with details of the substructuring control design. The focus is on configurations that use shake tables and active mass drivers. An extensive experimental assessment of the new strategy is presented in a companion paper, where the influence of various factors such as virtual subsystem dynamics, control gains, and nonlinearities is investigated, and it is shown that robustly stable and accurate substructuring is achieved.     

The Lagrangian governing equations of a rotating,compressible, viscous atmosphere

In this paper, the Lagrangian equations of motion, thermodynamics, continuity and diffusion of a rotating, compressible, viscous atmosphere are derived. It is shown that the pressure force and Laplacian terms involving velocity, temperature and concentration in the Eulerian system become nonlinear in the Lagrangian system. In the case that these Laplacian terms can be neglected, the governing equations in the Lagrangian system can be greatly simplified, and particle dynamics and dispersion can be investigated.     

Biofilm growth in porous media: Experiments,computational modeling at the porescale,and upscaling

Biofilm growth changes many physical properties of porous media such as porosity, permeability and mass transport parameters. The growth depends on various environmental conditions, and in particular, on flow rates. Modeling the evolution of such properties is difficult both at the porescale where the phase morphology can be distinguished, as well as during upscaling to the corescale effective properties. Experimental data on biofilm growth is also limited because its collection can interfere with the growth, while imaging itself presents challenges.In this paper we combine insight from imaging, experiments, and numerical simulations and visualization. The experimental dataset is based on glass beads domain inoculated by biomass which is subjected to various flow conditions promoting the growth of biomass and the appearance of a biofilm phase. The domain is imaged and the imaging data is used directly by a computational model for flow and transport. The results of the computational flow model are upscaled to produce conductivities which compare well with the experimentally obtained hydraulic properties of the medium. The flow model is also coupled to a newly developed biomass–nutrient growth model, and the model reproduces morphologies qualitatively similar to those observed in the experiment.     

Time-domain wavefield reconstruction inversion

Wavefield reconstruction inversion (WRI) is an improved full waveform inversion theory that has been proposed in recent years. WRI method expands the searching space by introducing the wave equation into the objective function and reconstructing the wavefield to update model parameters, thereby improving the computing efficiency and mitigating the influence of the local minimum. However, frequency-domain WRI is difficult to apply to real seismic data because of the high computational memory demand and requirement of time-frequency transformation with additional computational costs. In this paper, wavefield reconstruction inversion theory is extended into the time domain, the augmented wave equation of WRI is derived in the time domain, and the model gradient is modified according to the numerical test with anomalies. The examples of synthetic data illustrate the accuracy of time-domain WRI and the low dependency of WRI on low-frequency information.     

Dimensional analysis of the earthquake‐induced pounding between adjacent structures

In this paper the dynamic response of two and three pounding oscillators subjected to pulse‐type excitations is revisited with dimensional analysis. Using Buckingham's Π‐theorem the number of variables that govern the response of the system is reduced by three. When the response is presented in the dimensionless Π‐terms remarkable order emerges. It is shown that regardless of the acceleration level and duration of the pulse all response spectra become self‐similar and follow a single master curve. This is true despite the realization of finite duration contacts with increasing durations as the excitation level increases. All physically realizable contacts (impacts, continuous contacts, and detachments) are captured via a linear complementarity approach. The study confirms the existence of three spectral regions. The response of the most flexible among the two oscillators amplifies in the low range of the frequency spectrum (flexible structures); whereas, the response of the most stiff among the two oscillators amplifies at the upper range of the frequency spectrum (stiff structures). Most importantly, the study shows that pounding structures such as colliding buildings or interacting bridge segments may be most vulnerable for excitations with frequencies very different from their natural eigenfrequencies. Finally, by applying the concept of intermediate asymptotics, the study unveils that the dimensionless response of two pounding oscillators follows a scaling law with respect to the mass ratio, or in mathematical terms, that the response exhibits an incomplete self‐similarity or self‐similarity of the second kind with respect to the mass ratio. Copyright © 2008 John Wiley & Sons, Ltd.     

熵的原理与地震活动研究

熵是多种学科的一个重要交叉点。熵的原理广泛适用于物理、化学、生物、医学、气象乃至社会经济等科学领域,并已取得了令人注目的成果。本文以熵的基本含意出发,简要地阐述了熵概念的扩展,包括热温熵、信息熵以及条件熵等。将地震孕育过程作为一个复杂的开放系统,基于耗散结构理论的基本思想,为定量描述地震活动时、空、强分布的有序度,引入地震活动熵作为序参量,并将其应用于大地震前后地震活动演变特征的研究。初步结果表明,随着临近大地震的发生,在一定范围内地震活动熵值减小,即地震活动分布状态从无序趋向有序。可以预期,把熵的原理引入地震研究领域必将加深人们对地震发生物理过程本质的认识。同时也将为探索地震预测提供新的思路和新的途径。     

A new 3-D ray tracing method based on LTI using successive partitioning of cell interfaces and traveltime gradients

We present a new method of three-dimensional (3-D) seismic ray tracing, based on an improvement to the linear traveltime interpolation (LTI) ray tracing algorithm. This new technique involves two separate steps. The first involves a forward calculation based on the LTI method and the dynamic successive partitioning scheme, which is applied to calculate traveltimes on cell boundaries and assumes a wavefront that expands from the source to all grid nodes in the computational domain. We locate several dynamic successive partition points on a cell's surface, the traveltimes of which can be calculated by linear interpolation between the vertices of the cell's boundary. The second is a backward step that uses Fermat's principle and the fact that the ray path is always perpendicular to the wavefront and follows the negative traveltime gradient. In this process, the first-arriving ray path can be traced from the receiver to the source along the negative traveltime gradient, which can be calculated by reconstructing the continuous traveltime field with cubic B-spline interpolation. This new 3-D ray tracing method is compared with the LTI method and the shortest path method (SPM) through a number of numerical experiments. These comparisons show obvious improvements to computed traveltimes and ray paths, both in precision and computational efficiency.     

复阻尼地震反应谱的计算方法及其它

本文阐述了为什么要研究复阻尼地震反应谱问题,并给出了复阻尼地震反应谱计算的精确方法．实例计算结果表明,复阻尼反应港与常用的粘性阻尼反应谱之间存在着相当大的差异．文中还简述了复阻尼系统运动方程的建立方法。