Entropy Analysis of Multiple Scale Causality and Qualitative Causal Shifts in Spatial Systems*

Spatial systems are typically characterized by multiple controlling factors and processes operating at different spatial and temporal scales (multiple scale causality [MSC]). An entropy decomposition‐based approach to MSC is presented here in two contexts. First, given maps or distributions of an observed phenomenon at two or more scales, the contribution at more local or global (relative to the primary scale of observation) controls to the observed entropy can be estimated. Second, a theoretical treatment of the entropy decomposition equations shows that as the range of scale is increased by broadening or narrowing resolutions or by incorporating more controls, the influence of larger or smaller‐scale influences not only changes, but may change qualitatively, e.g., in terms of having positive (entropy‐increasing) or negative (information‐increasing) effects. Such qualitative causal shifts have implications for efforts to use any single causal explanation across the molecular to planetary spatial and instantaneous to geological range of scales relevant to physical geography. The entropy decomposition method is illustrated with an application to soil landscapes in the Ouachita Mountains, Arkansas.     

Effect of Fe on garnet–melt trace element partitioning: experiments in FCMAS and quantification of crystal-chemical controls in natural systems

Garnet–melt trace element partitioning experiments were performed in the system FeO–CaO–MgO–Al₂O₃–SiO₂ (FCMAS) at 3 GPa and 1540°C, aimed specifically at studying the effect of garnet Fe²⁺ content on partition coefficients (D^Grt/Melt). D^Grt/Melt, measured by SIMS, for trivalent elements entering the garnet X-site show a small but significant dependence on garnet almandine content. This dependence is rationalised using the lattice strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal–melt partition coefficients from elastic moduli. Nature 372, 452–454], which describes partitioning of an element i with radius r_i and valency Z in terms of three parameters: the effective radius of the site r₀(Z), the strain-free partition coefficient D₀(Z) for a cation with radius r₀(Z), and the apparent compressibility of the garnet X-site given by its Young's modulus E_X(Z). Combination of these results with data in Fe-free systems [Van Westrenen, W., Blundy, J.D., Wood, B.J., 1999. Crystal-chemical controls on trace element partitioning between garnet and anhydrous silicate melt. Am. Mineral. 84, 838–847] and crystal structure data for spessartine, andradite, and uvarovite, leads to the following equations for r₀(3+) and E_X(3+) as a function of garnet composition (X) and pressure (P):

r₀(3+) [Å]=0.930X_Py+0.993X_Gr+0.916X_Alm+0.946X_Spes+1.05(X_And+X_Uv)−0.005(P [GPa]−3.0)(±0.005 Å)

E_X(3+) [GPa]=3.5×10¹²(1.38+r₀(3+) [Å])^−26.7(±30 GPa)

Accuracy of these equations is shown by application to the existing garnet–melt partitioning database, covering a wide range of P and T conditions (1.8 GPa<P<5.0 GPa; 975°C<T<1640°C). D^Grt/Melt for all 3+ elements entering the X-site (REE, Sc and Y) are predicted to within 10–40% at given P, T, and X, when D^Grt/Melt for just one of these elements is known. In the absence of such knowledge, relative element fractionation (e.g. D_Sm^Grt/Melt/D_Nd^Grt/Melt) can be predicted. As an example, we predict that during partial melting of garnet peridotite, group A eclogite, and garnet pyroxenite, r₀(3+) for garnets ranges from 0.939±0.005 to 0.953±0.009 Å. These values are consistently smaller than the ionic radius of the heaviest REE, Lu. The above equations quantify the crystal-chemical controls on garnet–melt partitioning for the REE, Y and Sc. As such, they represent a major advance en route to predicting D^Grt/Melt for these elements as a function of P, T and X.     

The effect of Ca-Tschermaks component on trace element partitioning between clinopyroxene and silicate melt

We have studied the influence of Ca-Tschermaks (Calcium Tschermaks or CaTs) content of clinopyroxene on the partitioning of trace elements between this phase and silicate melt at fixed temperature and pressure. Ion probe analyses of experiments carried out in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂, at 0.1 MPa and 1218°C, produced crystal-melt partition coefficients (D) of 36 trace elements (Li, Cl, Sc, Ti, V, Cr, Fe, Co, Ge, Sr, Y, Zr, Nb, Mo, Ru, Rh, In, Sn, Sb, Ba, La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu, Hf, Ta and W), for clinopyroxene compositions between 10 and 32 mol% CaTs. Partition coefficients for 2+ to 5+ cations show, for each charge, a near parabolic dependence of log D on ionic radius of the substituting cation, for partitioning into both the M1 and M2 sites of clinopyroxene. Fitting the results to the elastic strain model of Blundy and Wood [Blundy, J.D., Wood, B.J., 1994. Prediction of crystal-melt partition coefficients from elastic moduli. Nature 372, 452–454] we obtain results for the strain-free partition coefficients of theoretical cations (D₀), with site radius r₀, and for the site's Young's Modulus (E).

In agreement with earlier data our results show that increasing ^ivAl concentration in cpx is matched by increasing D, E_M1, E_M2 and D₀ for tri-, tetra- and pentavalent cations. The degree of fractionation between chemically similar elements (i.e. Ta/Nb, Zr/Hf) also increases. In contrast, D values for mono-, di- and hexavalent cations decrease with increasing ^ivAl in the cpx. The large suite of trace elements used has allowed us to study the effects of cation charge on D₀, r₀ and E. We have found that D₀ and r₀ decrease with increasing cation charge, e.g. r₀=0.66 Å for 4+ cations and 0.59 Å for 5+ cations substituting into M1. Values of E_M1 and E_M2 increase with cation charge as well as with increasing ^ivAl content. The increase in E_M2 is linear and close to the trend set by Hazen and Finger [Hazen, R.M., Finger, L.W., 1979. Bulk modulus-volume relationship for cation–anion polyhedra. J. Geophys. Res. 84 (10) 6723–6728] for oxides. E_M1 values are much higher and do not fit the trend predicted by the Hazen and Finger relationship.     

Application of airborne LiDAR data and airborne multispectral imagery to structural mapping of the upper section of the Troodos ophiolite, Cyprus

Structural maps are traditionally produced by mapping features such as faults, folds, fabrics, fractures and joints in the field. However, large map areas and the spatially limited ground perspective of the field geologist can potentially increase the likelihood that not all structural features will be identified within a given area. The ability to recognise and map both local and regional structural features using high-resolution remote sensing data provides an opportunity to complement field-based mapping to help generate more comprehensive structural maps. Nonetheless, vegetation cover can adversely affect the extraction of structural information from remotely sensed data as it can mask the appearance of subtle spectral and geomorphological features that correspond to geological structures. This study investigates the utility of airborne Light Detection And Ranging (LiDAR) data and airborne multispectral imagery for detailed structural mapping in vegetated ophiolitic rocks and sedimentary cover of a section of the northern Troodos ophiolite, Cyprus. Visual enhancement techniques were applied to a 4-m airborne LiDAR digital terrain model and 4-m airborne multispectral imagery to assist the generation of structural lineament maps. Despite widespread vegetation cover, dykes and faults were recognisable as lineaments in both data sets, and the predominant strike trends of lineaments in all resulting maps were found to be in agreement with field-based structural data. Interestingly, prior to fieldwork, most lineaments were assumed to be faults, but were ground-verified as dykes instead, emphasising the importance of ground-truthing. Dyke and fault trends documented in this study define a pervasive structural fabric in the upper Troodos ophiolite that reflects the original sea-floor spreading history in the Larnaca graben. This structural fabric has not previously been observed in such detail and is likely to be continuous in adjacent regions under sedimentary cover. This information may be useful to future exploration efforts in the region focused on identification of structurally controlled mineral and groundwater resources. Overall, our case study highlights the efficacy of airborne LiDAR data and airborne multispectral imagery for extracting detailed and accurate structural information in hard-rock terrain to help complement field-based mapping.     

The effect of El Niño Southern Oscillation cycles on the decadal scale suspended sediment behavior of a coastal dry‐summer subtropical catchment

Rivers display temporal dependence in suspended sediment–water discharge relationships. Although most work has focused on multi‐decadal trends, river sediment behavior often displays sub‐decadal scale fluctuations that have received little attention. The objectives of this study were to identify inter‐annual to decadal scale fluctuations in the suspended sediment–discharge relationship of a dry‐summer subtropical river, infer the mechanisms behind these fluctuations, and examine the role of El Niño Southern Oscillation climate cycles. The Salinas River (California) is a moderate sized (11 000 km²), coastal dry‐summer subtropical catchment with a mean discharge (Q_mean) of 11.6 m³ s^?1. This watershed is located at the northern most extent of the Pacific coastal North America region that experiences increased storm frequency during El Niño years. Event to inter‐annual scale suspended sediment behavior in this system was known to be influenced by antecedent hydrologic conditions, whereby previous hydrologic activity regulates the suspended sediment concentration–water discharge relationship. Fine and sand suspended sediment in the lower Salinas River exhibited persistent, decadal scale periods of positive and negative discharge corrected concentrations. The decadal scale variability in suspended sediment behavior was influenced by inter‐annual to decadal scale fluctuations in hydrologic characteristics, including: elapsed time since small (~0.1 × Q_mean), and moderate (~10 × Q_mean) threshold discharge values, the number of preceding days that low/no flow occurred, and annual water yield. El Niño climatic activity was found to have little effect on decadal‐scale fluctuations in the fine suspended sediment–discharge relationship due to low or no effect on the frequency of moderate to low discharge magnitudes, annual precipitation, and water yield. However, sand concentrations generally increased in El Niño years due to the increased frequency of moderate to high magnitude discharge events, which generally increase sand supply. Copyright © 2014 John Wiley & Sons, Ltd.     

臭氧和平流层动力学的相互作用

讨论了给定的南极春季臭氧洞(取自1979—1985年臭氧减少的观测结果)对二维平流层—对流层模式中温度和环流的影响。11月份,南极上空约17km处,温度最多可降低6℃。这种温度变化引起的平均经向环流对臭氧洞起填塞作用,不过,这种影响很小,每年仅产生14DU的变化。观测事实表明,近年来10月份,南半球波活动减弱。为此,我们作了南半球波作用全年都减少一半,并考虑了臭氧洞的数值试验。结果表明,臭氧柱在11月份76°S减少了44DU,在赤道却增加了12DU。     

空间尺度转换与跨尺度信息链接:区域生态水文模拟研究空间尺度转换方法综述

空间尺度转换是近年来区域生态水文研究领域的一个基本研究问题。其需要主要是源于模型的输入数据与所能提供的数据空间尺度不一致以及模型所代表的地表过程空间尺度与所观测的地表过程空间尺度不吻合。综述了目前区域生态水文模拟研究中常用的空间尺度转换研究方法,包括向上尺度转换和向下尺度转换。详细论述了2种向下尺度转换方法: 统计学经验模型和动态模型。前者是通过将GCM大尺度数据与长期的历史观测数据比较从而建立统计学相关模型, 然后利用这个统计学经验模型进行向下的空间尺度转换. 然而动态模型并不直接对GCM数据进行向下尺度的转换,而是对与GCM进行动态耦合的区域气候模型(RCM) 的输出数据进行空间尺度转换. 通常后者所获得的数据精度要比前者高,但是一个主要缺点就是并不是全球所有的研究区域都有对应的RCM。还详细论述了2种向上尺度转换方法: 统计学经验模型和斑块模型。前者是建立一个能代表小尺度信息在大尺度上分布的密度分布概率函数, 然后利用这个函数在所需的大尺度上进行积分而求得大尺度所需的信息。而后者是根据相似性最大化原则将大尺度划分为若干个可操作的小尺度斑块,然后将计算的每个小尺度斑块的信息平均化得到大尺度所需的信息。通常在计算这种斑块化的小尺度信息的时候,对每个小尺度也会采用统计学经验模型来计算代表整个斑块小尺度的信息。建议用斑块模型与统计学经验模型相集合的方法来实现向上的空间尺度转换     

Landslide Research at the British Geological Survey:Capture,Storage and Interpretation on a National and Site-Specific Scale

Abstract: Landslide research at the British Geological Survey (BGS) is carried out through a number of activities, including surveying, database development and real-time monitoring of landslides. Landslide mapping across the UK has been carried out since BGS started geological mapping in 1835. Today, BGS geologists use a combination of remote sensing and ground-based investigations to survey landslides. The development of waterproof tablet computers (BGS·SIGMAmobile), with inbuilt GPS and GIS for field data capture provides an accurate and rapid mapping methodology for field surveys. Regional and national mapping of landslides is carried out in conjunction with site-specific monitoring, using terrestrial LiDAR and differential GPS technologies, which BGS has successfully developed for this application. In addition to surface monitoring, BGS is currently developing geophysical ground-imaging systems for landslide monitoring, which provide real-time information on subsurface changes prior to failure events. BGS’s mapping and monitoring activities directly feed into the BGS National Landslide Database, the most extensive source of information on landslides in Great Britain. It currently holds over 14?000 records of landslide events. By combining BGS’s corporate datasets with expert knowledge, BGS has developed a landslide hazard assessment tool, GeoSure, which provides information on the relative landslide hazard susceptibility at national scale.     

Devon's desert ‘worms’

Worm‐like trace fossils, sometimes of large size, have regularly been reported from the otherwise generally poorly‐fossiliferous Permo‐Triassic continental red beds of the East Devon coast, southwest England. Selected examples are discussed here to outline the difficulties involved in elucidating the true producers of these burrows and interpreting their significance in the local palaeoenvironment.     

Characterizing and classifying urban watersheds with compositional and structural attributes

Current land-use classifications used to assess urbanization effects on stream water quality date back to the 1980s when limited information was available to characterize watershed attributes that mediate non-point source pollution. With high resolution remote sensing and widely used GIS tools, there has been a vast increase in the availability and precision of geospatial data of built environments. In this study, we leverage geospatial data to expand the characterization of developed landscapes and create a typology that allows us to better understand the impact of complex developed landscapes across the rural to urban gradient. We assess the ability of the developed landscape typology to reveal patterns in stream water chemistry previously undetected by traditional land-cover based classification. We examine the distribution of land-cover, infrastructure, topography and geology across 3876 National Hydrography Dataset Plus catchments in the Piedmont region of North Carolina, USA. From this dataset, we generate metrics to evaluate the abundance, density and position of landscape features relative to streams, catchment outlets and topographic wetness metrics. While impervious surfaces are a key distinguishing feature of the urban landscape, sanitary infrastructure, population density and geology are better predictors of baseflow stream water chemistry. Unsupervised clustering was used to generate a distinct developed landscape typology based on the expanded, high-resolution landscape feature information. Using stream chemistry data from 37 developed headwater catchments, we compared the baseflow water chemistry grouped by traditional land-cover based classes of urbanization (rural, low, medium and high density) to our composition and structure-based classification (a nine-class typology). The typology based on 22 metrics of developed landscape composition and structure explained over 50% of the variation in NO₃⁻-N, TDN, DOC, Cl⁻, and Br⁻ concentration, while the ISC-based classification only significantly explained 23% of the variation in TDN. These results demonstrate the importance of infrastructure, population and geology in defining developed landscapes and improving discrete classes for water management.