Local dynamical equivalence of certain food webs

An important challenge in theoretical ecology is to find good coarse-grained representations of complex food webs. Here, we use the approach of generalized modeling to show that it may be possible to formulate a coarse-graining algorithm that conserves the local dynamics of the model exactly. We show examples of food webs with a different number of species that have exactly identical local bifurcation diagrams. Based on these observations, we formulate a conjecture governing which populations of complex food webs can be grouped together into a single variable without changing the local dynamics. As an illustration, we use this conjecture to show that chaotic regions generically exist in the parameter space of a class of food webs with more than three trophic levels. While our conjecture is at present only applicable to relatively special cases, we believe that its applicability could be greatly extended if a more sophisticated mapping of parameters were used in the model reduction.     

基于SD和CA的城镇土地扩展模拟模型——以江苏省南通地区为例

充分利用系统动力学模型(System Dynamics,SD)在情景模拟和宏观因素反映上的优势和元胞自动机模型(Cellular Automata,CA)在微观土地利用空间格局反映上的优势,构建一个耦合SD和CA的城镇土地扩展模拟模型,并以江苏省南通地区为例,对模型的实证应用做了进一步的验证。结果表明,这种耦合模型不仅能够对研究区域未来城镇土地扩展数量给予一个比较好的预测,而且还对其空间分布效果做了一定精度上的模拟,这使得城市规划在土地利用预测方面有一个相对科学的依据。     

美国罗德河流域地下水中铝形成动力学

In the Rhode River estuarine/watershed area, increased acidity of precipitation from atmospheric acidic deposition has deleterious effects on the freshwater ecosystems. One of the characteristics of an acidic watershed is the mobilization of aluminum from soils to aquatic environment. Increased concentrations of aluminum in surface waters are toxic to living organisms. Detection of long-term changes of acidity and elevated Al concentrations in surface and soil waters is a central issue. In this present paper, the dynamics ofAl speciation in the shallow ground waters from Rhode River watershed 110 (la) were investigated. This research provides a unique, regional analysis of important controls on the transport of Al speciation through the rural watershed of the Chesapeake Bay.     

A multidisciplinary study of a slow-slipping fault for seismic hazard assessment: the example of the Middle Durance Fault (SE France)

Assessing seismic hazard in continental interiors is difficult because these regions are characterized by low strain rates and may be struck by infrequent destructive earthquakes. In this paper, we provide an example showing that interpretations of seismic cross sections combined with other kinds of studies such as analysis of microseismicity allow the whole seismogenic source area to be imaged in this type of region. The Middle Durance Fault (MDF) is an 80-km-long fault system located southeastern France that has a moderate but regular seismicity and some palaeoseismic evidence for larger events. It behaves as an oblique ramp with a left-lateral-reverse fault slip and has a low strain rate. MDF is one of the rare slow active fault system monitored by a dedicated dense velocimetric short period network. This study showed a fault system segmented in map and cross section views which consists of staircase basement faults topped by listric faults ramping off Triassic evaporitic beds. Seismic sections allowed the construction of a 3-D structural model used for accurate location of microseismicity. Southern part of MDF is mainly active in the sedimentary cover. In its northern part and in Alpine foreland, seismicity deeper than 8 km was also recorded meaning active faults within the crust cannot be excluded. Seismogenic potential of MDF was roughly assessed. Resulting source sizes and estimated slip rates imply that the magnitude upper limit ranges from 6.0 to 6.5 with a return period of a few thousand years. The present study shows that the coupling between 3-D fault geometry imaging and accurate location of microseismicity provides a robust approach to analyse active fault sources and consequently a more refined seismic hazard assessment.     

易门式大型铜矿床构造成矿动力学模型

针对金属矿床隐伏矿定位预测的关键问题,在提出构造成矿动力学的分支学术方向的基础上,以易门式大型铜矿床为例,从成矿地质背景、矿田构造、构造地球化学和成矿构造应力场等方面进行构造成矿动力学研究,阐述在构造应力场控制下成矿流体运移和聚集的规律。构造地球化学研究认为,构造地球化学异常反映矿体的原生晕;构造地球化学异常受构造控制,其分带特征可指示成矿流体的流向,并提供矿床成因和隐伏矿(化)体相对埋深的信息。构造应力场的控矿特征主要表现在:构造应力场导致控矿构造的形成,并驱动成矿流体的运移;构造应力场控制了成矿能量场;应力和能量的高值集中区分布特点可反映某些构造型式。在此基础上建立构造成矿动力学模型,概括隐伏矿成矿预测准则,据此提出若干重点找矿靶区和靶位,其中部分靶区得到工程验证,表明构造成矿动力学方法对隐伏矿定位预测和评价具有重要的指导意义。     

Multi-wavelength simulations of atmospheric radiation from Io with a 3-D spherical-shell backward Monte Carlo radiative transfer model

Conflicting observations regarding the dominance of either sublimation or volcanism as the source of the atmosphere on Io and disparate reports on the extent of its spatial distribution and the absolute column abundance invite the development of detailed computational models capable of improving our understanding of Io’s unique atmospheric structure and origin. Improving upon previous models, Walker et al. (Walker, A.C., Gratiy, S.L., Levin, D.A., Goldstein, D.B., Varghese, P.L., Trafton, L.M., Moore, C.H., Stewart, B. [2009]. Icarus) developed a fully 3-D global rarefied gas dynamics model of Io’s atmosphere including both sublimation and volcanic sources of SO₂ gas. The fidelity of the model is tested by simulating remote observations at selected wavelength bands and comparing them to the corresponding astronomical observations of Io’s atmosphere. The simulations are performed with a new 3-D spherical-shell radiative transfer code utilizing a backward Monte Carlo method. We present: (1) simulations of the mid-infrared disk-integrated spectra of Io’s sunlit hemisphere at 19 μm, obtained with TEXES during 2001-2004; (2) simulations of disk-resolved images at Lyman-α obtained with the Hubble Space Telescope (HST), Space Telescope Imaging Spectrograph (STIS) during 1997-2001; and (3) disk-integrated simulations of emission line profiles in the millimeter wavelength range obtained with the IRAM-30 m telescope in October-November 1999. We found that the atmospheric model generally reproduces the longitudinal variation in band depth from the mid-infrared data; however, the best match is obtained when our simulation results are shifted ∼30° toward lower orbital longitudes. The simulations of Lyman-α images do not reproduce the mid-to-high latitude bright patches seen in the observations, suggesting that the model atmosphere sustains columns that are too high at those latitudes. The simulations of emission line profiles in the millimeter spectral region support the hypothesis that the atmospheric dynamics favorably explains the observed line widths, which are too wide to be formed by thermal Doppler broadening alone.     

How tides and river flows determine estuarine bathymetries

For strongly tidal, funnel-shaped estuaries, we examine how tides and river flows determine size and shape. We also consider how long it takes for bathymetric adjustment, both to determine whether present-day bathymetry reflects prevailing forcing and how rapidly changes might occur under future forcing scenarios.Starting with the assumption of a 'synchronous' estuary (i.e., where the sea surface slope resulting from the axial gradient in phase of tidal elevation significantly exceeds the gradient in tidal amplitude ), an expression is derived for the slope of the sea bed. Thence, by integration we derive expressions for the axial depth profile and estuarine length, L, as a function of and D, the prescribed depth at the mouth. Calculated values of L are broadly consistent with observations. The synchronous estuary approach enables a number of dynamical parameters to be directly calculated and conveniently illustrated as functions of and D, namely: current amplitude Û, ratio of friction to inertia terms, estuarine length, stratification, saline intrusion length, flushing time, mean suspended sediment concentration and sediment in-fill times.Four separate derivations for the length of saline intrusion, L_I, all indicate a dependency on (U_o is the residual river flow velocity and f is the bed friction coefficient). Likely bathymetries for `mixed' estuaries can be delineated by mapping, against and D, the conditions L<sub>I</sub>/L<1,E<sub>X</sub>/L<1 (E<sub>X</sub> is the tidal excursion) alongside the Simpson-Hunter criteria D/U<sup>3</sup><50 m<sup>−2</sup> s<sup>3</sup>. This zone encompasses 24 out of 25 `randomly' selected UK estuaries.However, the length of saline intrusion in a funnel-shaped estuary is also sensitive to axial location. Observations suggest that this location corresponds to a minimum in landward intrusion of salt. By combining the derived expressions for L and L_I with this latter criterion, an expression is derived relating D_i, the depth at the centre of the intrusion, to the corresponding value of U_o. This expression indicates U_o is always close to 1 cm s⁻¹, as commonly observed. Converting from U_o to river flow, Q, provides a morphological expression linking estuarine depth to Q (with a small dependence on side slope gradients).These dynamical solutions are coupled with further generalised theory related to depth and time-mean, suspended sediment concentrations (as functions of and D). Then, by assuming the transport of fine marine sediments approximates that of a dissolved tracer, the rate of estuarine supply can be determined by combining these derived mean concentrations with estimates of flushing time, F_T, based on L_I. By further assuming that all such sediments are deposited, minimum times for these deposition rates to in-fill estuaries are determined. These times range from a decade for the shortest, shallowest estuaries to upwards of millennia in longer, deeper estuaries with smaller tidal ranges.     

Kinematic bending moments in pile foundations

In this paper the kinematic seismic interaction of single piles embedded in soil deposits is evaluated by focusing the attention on the bending moments induced by the transient motion. The analysis is performed by modeling the pile like an Euler–Bernoulli beam embedded in a layered Winkler-type medium. The excitation motion is obtained by means of a one-D propagation analysis. A comprehensive parametric analysis is carried out by varying the main parameters governing the dynamic response of piles like the soil properties, the bedrock location, the diameter and embedment in the bedrock of piles. On the basis of the parametric analysis, a new design formula for predicting the kinematic bending moments for both the cross-sections at the deposit–bedrock interface and at the pile head is proposed.     

http://www.sciencedirect.com/science/article/pii/S1674987114000905

Although there has been significant progress in the seismic imaging of mantle heterogeneity,the outstanding issue that remains to be resolved is the unknown distribution of mantle temperature anomalies in the distant geological past that give rise to the present-day anomalies inferred by global tomography models.To address this question,we present 3-D convection models in compressible and self-gravitating mantle initialised by different hypothetical temperature patterns.A notable feature of our forward convection modelling is the use of self-consistent coupling of the motion of surface tectonic plates to the underlying mantle flow,without imposing prescribed surface velocities(i.e.,plate-like boundary condition).As an approximation for the surface mechanical conditions before plate tectonics began to operate we employ the no-slip(rigid) boundary condition.A rigid boundary condition demonstrates that the initial thermally-dominated structure is preserved,and its geographical location is fixed during the evolution of mantle flow.Considering the impact of different assumed surface boundary conditions(rigid and plate-like) on the evolution of thermal heterogeneity in the mantle we suggest that the intrinsic buoyancy of seven superplumes is most-likely resolved in the tomographic images of present-day mantle thermal structure.Our convection simulations with a plate-like boundary condition reveal that the evolution of an initial cold anomaly beneath the Java-Indonesian trench system yields a long-term,stable pattern of thermal heterogeneity in the lowermost mantle that resembles the presentday Large Low Shear Velocity Provinces(LLSVPs),especially below the Pacific.The evolution of subduction zones may be,however,influenced by the mantle-wide flow driven by deeply-rooted and longlived superplumes since Archean times.These convection models also detect the intrinsic buoyancy of the Perm Anomaly that has been identified as a unique slow feature distinct from the two principal LLSVPs.We find there is no need for dense chemical ’piles’ in the lower mantle to generate a stable distribution of temperature anomalies that are correlated to the LLSVPs and the Perm Anomaly.Our tomography-based convection simulations also demonstrate that intraplate volcanism in the south-east Pacific may be interpreted in terms of shallow small-scale convection triggered by a superplume beneath the East Pacific Rise.     

莱州湾海湾扇贝养殖区海水中悬浮颗粒的动态变化

1997年10月和1998年6月，在莱州湾两个主要海湾扇贝养殖区－芙蓉岛海区和金城海区，各定点进行了24h连续观测和取样，并在实验室内分析了总颗粒物（TPM）、颗粒有机物（POM）和叶绿素α的含量。结果表明，6月份金城海区的颗粒有机物和叶绿素α的含量均高于荚蓉岛海区；在10月份，芙蓉岛海区颗粒有机物的含量高于金城海区，但其叶绿素α的含量却低于金城海区，表明金城海区颗粒有机物的质量高于芙蓉岛海区。