Geochemical Characterization of Bolivian Formative Earthen Architecture by Wavelength‐Dispersive X‐ray Fluorescence

In this study, we employ wavelength‐dispersive X‐ray fluorescence (WDXRF) to characterize construction materials from Formative civic architecture (1000 B.C.E.–C.E. 400), ethnographic mudbricks, and off‐site controls from the Taraco Peninsula, Bolivia. The preparation of earthen construction materials for civic buildings can shed light on aspects of community development such as labor organization, resource management, and architectural technologies. We apply geochemical results to reconstructing how public buildings were made as communities moved toward socio‐political complexity in this region. However, there are few geochemical studies in the Andes, and little prior scientific analysis of earthen architecture. We therefore tested the efficacy of WDXRF for this region, and developed control materials. Our archaeological samples were selected from two Formative villages, Chiripa and Kala Uyuni. In this study, we performed WDXRF analyses on 63 archaeological and control samples including archaeological floors, walling, plasters, and mortars, as well as contemporary ethnographic walling and topsoils. Elemental signatures for 28 elements clearly distinguished the archaeological flooring, walling, plaster, and mortars from ethnographic and off‐site controls. More subtle variations were detected that distinguish study sites and the different material types. Laboratory‐calibrated multi‐element XRF effectively supports efforts to reconstruct the pathways to social complexity in the Titicaca Basin.     

Geoarchaeological Investigations in Mesoamerica Move into the 21st Century: A Review

Over the past thirty years, geoarchaeology has moved from the fringe to mainstream status within Mesoamerican archaeological investigations. This review focuses on works published since the year 2000. Five themes are identified as central to recent studies: (1) the correlation of environmental change and cultural history; (2) anthropogenic environmental impacts; (3) ancient land cover, land use, and diet; (4) archaeological prospection; and (5) provenance studies. These themes are often interwoven in the application of complex systems approaches that allow scientists to more accurately model the intricacies of ancient human–environment interactions.     

Variability in bedform morphology and kinematics with transport stage

Bedform geometry is widely recognized to be a function of transport stage. Bedform aspect ratio (height/length) increases with transport stage, reaches a maximum, then decreases as bedforms washout to a plane bed. Bedform migration rates are also linked to bedform geometry, in so far as smaller bedforms in coarser sediment tend to migrate faster than larger bedforms in finer sediment. However, how bedform morphology (height, length and shape) and kinematics (translation and deformation) change with transport stage and suspension have not been examined. A series of experiments is presented where initial flow depth and grain size were held constant and the transport stage was varied to produce bedload dominated, mixed‐load dominated and suspended‐load dominated conditions. The results show that the commonly observed pattern in bedform aspect ratio occurs because bedform height increases then decreases with transport stage, against a continuously increasing bedform length. Bedform size variability increased with transport stage, leading to less uniform bedform fields at higher transport stage. Total translation‐related and deformation‐related sediment fluxes all increased with transport stage. However, the relative contribution to the total flux changed. At the bedload dominated stage, translation‐related and deformation‐related flux contributed equally to the total flux. As the transport stage increased, the fraction of the total load contributed by translation increased and the fraction contributed by deformation declined because the bedforms got bigger and moved faster. At the suspended‐load dominated transport stage, the deformation flux increased and the translation flux decreased as a fraction of the total load, approaching one and zero, respectively, as bedforms washed out to a plane bed.     

Significance of shallow core transects for reef models and sea‐level curves,Heron Reef,Great Barrier Reef

A sequence of shallow reef cores from Heron Reef, Great Barrier Reef, provides new insights into Holocene reef growth models. Isochron analysis of a leeward core transect suggests that the north‐western end of Heron Reef reached current sea‐level by ca 6·5 kyr bp and then prograded leeward at a rate of ca 19·6 m/kyr between 5·1 kyr and 4·1 kyr bp (pre‐1950) to the present reef margin. A single short core on the opposing margin of the reef is consistent with greater and more recent progradation there. Further to the east, one windward core reached modern sea‐level by ca 6·3 kyr bp , suggesting near ‘keep‐up’ behaviour at that location, but the opposing leeward margin behind the lagoon reached sea‐level much more recently. Hence, Heron Reef exhibited significantly different reef growth behaviour on different parts of the same margin. Mean reef accretion rates calculated from within 20 m of one another in the leeward core transect varied between ca 2·9 m and 4·7 m/kyr depending on relative position in the prograding wedge. These cores serve as a warning regarding the use of isolated cores to inform reef growth rates because apparent aggradation at any given location on a reef varies depending on its location relative to a prograding margin. Only transects of closely spaced cores can document reef behaviour adequately so as to inform reef growth models and sea‐level curves. The cores also emphasize potential problems in U‐series dates for corals within a shallow (ca 1·5 m) zone beneath the reef flat. Apparent age inversions restricted to that active diagenetic zone may reflect remobilization and concentration of Th in irregularly distributed microbialites or biofilms that were missed during sample vetting. Importantly, the Th‐containing contaminant causes ages to appear too old, rather than too young, as would be expected from younger cement.     

Thresholds of intrabed flow and other interactions of turbidity currents with soft muddy substrates

Controlled laboratory experiments reveal that the lower part of turbidity currents has the ability to enter fluid mud substrates, if the bed shear stress is higher than the yield stress of the fluid mud and the density of the turbidity current is higher than the density of the substrate. Upon entering the substrate, the turbidity current either induces mixing between flow‐derived sediment and substrate sediment, or it forms a stable horizontal flow front inside the fluid mud. Such ‘intrabed’ flow is surrounded by plastically deformed mud; otherwise it resembles the front of a ‘bottom‐hugging’ turbidity current. The ‘suprabed’ portion of the turbidity current, i.e. the upper part of the flow that does not enter the substrate, is typically separated from the intrabed flow by a long horizontal layer of mud which originates from the mud that is swept over the top of the intrabed flow and then incorporated into the flow. The intrabed flow and the mixing mechanism are specific types of interaction between turbidity currents and muddy substrates that are part of a larger group of interactions, which also include bypass, deposition, erosion and soft sediment deformation. A classification scheme for these types of interactions is proposed, based on an excess bed shear stress parameter, which includes the difference in the bed shear stress imposed by the flow and the yield stress of the substrate and an excess density parameter, which relies on the density difference between the flow and the substrate. Based on this classification scheme, as well as on the sedimentological properties of the laboratory deposits, an existing facies model for intrabed turbidites is extended to the other types of interaction involving soft muddy substrates. The physical threshold of flow‐substrate mixing versus stable intrabed flow is defined using the gradient Richardson number, and this method is validated successfully with the laboratory data. The gradient Richardson number is also used to verify that stable intrabed flow is possible in natural turbidity currents, and to determine under which conditions intrabed flow is likely to be unstable. It appears that intrabed flow is likely only in natural turbidity currents with flow velocities well below ca 3·5 m s^?1, although a wider range of flows is capable of entering fluid muds. Below this threshold velocity, intrabed flow is stable only at high‐density gradients and low‐velocity gradients across the upper boundary of the turbidity current. Finally, the gradient Richardson number is used as a scaling parameter to set the flow velocity limits of a natural turbidity current that formed an inferred intrabed turbidite in the deep‐marine Aberystwyth Grits Group, West Wales, United Kingdom.     

Development of large‐scale seismites in Upper Cretaceous fluvial sandstones in a fault‐proximal setting

Large‐scale soft‐sediment deformation structures occur within fluvial sandstone bodies of the Upper Cretaceous Wahweap Formation in the Kaiparowits basin, southern Utah, USA. These structures represent an exceptional example of metre‐scale fault‐proximal, seismogenic load structures in nearly homogenous sandstones. The load structures consist of two types: large‐scale load casts and wedge‐shaped load structures. Large‐scale load casts penetrate up to 4·5 m into the underlying sandstone bed. Wedge‐shaped load structures include metre‐scale, parallel, sub‐vertical features and decimetre‐scale features along the periphery of the large‐scale load casts or other wedge‐shaped load structures. Wedge‐shaped load structures contain well‐developed, medial cataclastic shear deformation bands. All load structures contain pervasive well‐defined millimetre‐thick to centimetre‐thick internal laminae, oriented parallel to the outside form of the load structures and asymptotic to deformation bands. Both types of load structures formed because of an inverted density profile, earthquake‐triggered liquefaction and growth of irregularities (a Rayleigh–Taylor instability) on the sandstone–sandstone erosional contact. The internal laminae and deformation bands formed during deformation and clearly demonstrate polyphase deformation, recording a transition from liquefied to hydroplastic to brittle modes of deformation. Decimetre‐scale wedge‐shaped load structures on the edge of the large‐scale load casts probably formed towards the end of a seismic event after the sediment dewatered and increased the frictional contact of grains enough to impart strength to the sands. Metre‐scale wedge‐shaped load structures were created as the tips of downward foundering sediments were driven into fractures, which widened incrementally with seismic pulsation. With each widening of the fracture, gravity and a suction effect would draw additional sediment into the fracture. Superimposed laminae indicate a secondary syndeformational origin for internal laminae, probably by flow‐generated shearing and vibrofluidization mechanisms. Large‐scale and wedge‐shaped load structures, polyphase deformation and secondary laminae may characterize soft‐sediment deformation in certain fault‐proximal settings.     

Progradational and backstepping shoreface deposits in the Ladinian to Early Norian Snadd Formation of the Barents Sea

Shoreface sandstone deposits within the Early Carnian part of the Snadd Formation of the Norwegian Barents Sea can be traced for hundreds of kilometres in the depositional strike direction and for tens of kilometres in the depositional‐dip direction. This study uses three‐dimensional seismic attribute mapping and two‐dimensional regional seismic profiles to visualize the seismic facies of these shoreface deposits and to map their internal stratigraphic architecture at a regional scale. The shoreface deposits are generally elongate but show variable width from north‐east to south‐west, which corresponds to a sediment source in the northern part of the basin and a southward decrease in longshore sediment transport. The Snadd Formation presents an example of how large‐scale progradational shoreface deposits develop. The linear nature of its shoreface deposits contrasts with more irregular, cuspate wave‐dominated deltaic shorelines that contain river outlets, and instead implies longshore drift as the main sediment source. In map view, discrete sets of linear features bounded by truncation surfaces scale directly to beach ridge sets in modern counterparts. The shoreface deposits studied here are characteristic in terms of scale and basin‐wide continuity, and offer insight into the contrast between shallow marine deposition under stable Triassic Greenhouse and fluctuating Holocene Icehouse climates. Findings presented herein are also important for hydrocarbon exploration in the Barents Sea, because they describe a hitherto poorly understood reservoir play in the Triassic interval, wherein the most prominent reservoir plays have so far been considered to be found in channelized deposits in net‐progradational delta‐plain strata that form the topsets to shelf‐edge clinoforms. The documented presence of widespread wave‐dominated shoreface deposits also has implications for how the relative importance of different sedimentary processes is considered within the basin during this period.     

Iron isotopic systematics of UHP eclogites respond to oxidizing fluid during exhumation

The iron stable isotope compositions (δ⁵⁶Fe) and iron valence states of ultrahigh‐pressure eclogites from Bixiling in the Dabie orogen belt, China, were measured to trace the changes of geochemical conditions during vertical transportation of earth materials, for example, oxygen fugacity. The bulk Fe³⁺/ΣFe ratios of retrograde eclogites, determined by Mössbauer spectroscopy, are consistently higher than those of fresh eclogites, suggesting oxidation during retrograde metamorphism and fluid infiltration. The studied eclogites (five samples) display limited mid‐ocean ridge basalts (MORB)‐like (~0.10‰) δ⁵⁶Fe values, which are indistinguishable from their protoliths, that is, gabbro cumulates formed through differentiation of mantle‐derived basaltic magma. This suggests that Fe isotope fractionation during continental subduction is limited. Garnet separates display limited δ⁵⁶Fe variation ranging from ?0.08 ± 0.07 ‰ to 0.02 ± 0.07‰, whereas coexisting omphacite displays a large variation of δ⁵⁶Fe values from 0.15 ± 0.07‰ to 0.47 ± 0.07‰. Omphacite also has highly variable Fe³⁺/ΣFe ratios from 0.367 ± 0.025 to 0.598 ± 0.024, indicating modification after peak metamorphism. Omphacite from retrograde eclogites has elevated Fe³⁺/ΣFe ratios (0.54–0.60) compared to that from fresh eclogites (~0.37), whereas garnet displays a narrow range of ferric iron content with Fe³⁺/ΣFe ratios from 0.039 ± 0.013 to 0.065 ± 0.022. The homogenous δ⁵⁶Fe values and Fe³⁺/ΣFe ratios of garnet suggest that it survived the retrograde metamorphism and preserved its Fe‐isotopic features and ferric contents of peak metamorphism. Because of similar diffusion rates of Fe and Mg in garnet and omphacite, and constant Δ²⁶Mg_{omphacite‐garnet} values (1.14 ± 0.04‰), equilibrium iron isotope fractionation between garnet and omphacite was probably achieved during peak metamorphism. Elevated Fe³⁺/ΣFe ratios of omphacite from retrograde eclogites and variant Δ⁵⁶Fe_{omphacite‐garnet} values of the studied eclogites (0.13 ± 0.10‰ to 0.48 ± 0.10‰) indicate that oxidized geofluid infiltration resulted in the elevation of δ⁵⁶Fe values of omphacite during retrograde metamorphism.     

Modeling Urban Land Use Changes Using Support Vector Machines

Support Vector Machines (SVM) is a machine learning (ML) algorithm commonly applied to the classification of remotely sensing data and more recently for modeling land use changes. However, in most geospatial applications the current literature does not elaborate on specifications of the SVM method with respect to data sampling, attribute selection and optimal parameters choices. Therefore the main objective of this study is to present and investigate the SVM technique for modeling urban land use change. The SVM model building procedure is presented together with the detailed evaluation of the output results with respect to the choice of datasets, attributes and the change of SVM parameters. Geospatial datasets containing nine land use classes and spatial attributes for the Municipality of Zemun, Republic of Serbia were used for years 2001, 2003, 2007 and 2011. The Correlation‐based Feature Subset method, kappa coefficient, Area Under Receiver Operating Characteristic Curve (AUC) and kappa simulation were used to perform the model evaluation and compare the model outputs with the real land use datasets. The obtained results indicate that the SVM‐based models perform better when implementing balanced data sampling, reduced data sets to informative subsets of attributes and properly identify the optimal learning parameters.     

A hybrid approach to modeling territorial control in violent armed conflicts

Territorial control is central to the understanding of violent armed conflicts, yet reliable and valid measures of this concept do not exist. We argue that geospatial analysis provides an important perspective to measure the concept. In particular, measuring territorial control can be seen as an application of calculating service areas around points of control. The modeling challenge is acute for areas with limited road infrastructure, where no complete network is available to perform the analysis, and movements largely occur off road. We present a new geospatial approach that applies network analysis on a hybrid transportation network with both actual road data and hexagon‐fishnet‐based artificial road data representing on‐road and off‐road movements, respectively. Movement speed or restriction can be readily adjusted using various input data. Simulating off‐road movement with hexagon‐fishnet‐based artificial road data has a number of advantages including scalability to small or large study areas and flexibility to allow all‐directional travel. We apply this method to measuring territorial control of armed groups in Sub‐Saharan Africa where inferior transport infrastructure is the norm. Based on the Uppsala Conflict Data Program's (UCDP) Georeferenced Event Data (GED) as well as spatial data on terrain, population locations, and limited transportation networks, we enhance the delineation of the specific areas directly controlled by each warring party during civil wars within a given travel time.