An extreme wave event in eastern Yucatán,Mexico: Evidence of a palaeotsunami event during the Mayan times

The Yucatán Peninsula, Mexico, has typically been considered a tectonically stable region with little significant seismic activity. The region though, is one that is regularly affected by hurricanes. A detailed survey of ca 100 km of the eastern Yucatán and Cozumel coast identified the presence of ridges containing individual boulders measuring >1 m in length. The boulder ridges reach 5 m in height and their origin is associated with extreme wave event activity. Previously modelled tsunami waves from known seismically active zones in the region (Muertos Thrust Belt and South Caribbean Deformed Belt) are not of sufficient scale in the area of the Yucatán Peninsula to have produced the boulder ridges recorded in this study. The occurrence of hurricanes in this region is more common, but two of the most destructive (Hurricane Gilbert 1988 and Hurricane Wilma 2005) produced coastal waves too small to have created the ridges recorded here. In this paper, a new tsunami model with a source area located on the Motagua/Swan Island Fault System has been generated that indicates a tsunami event may have caused the extreme wave events that resulted in the deposition of the boulder ridges.     

Recent progress on combining geomorphological and geochronological data with ice sheet modelling,demonstrated using the last British–Irish Ice Sheet

Palaeo-ice sheets are important analogues for understanding contemporary ice sheets, offering a record of ice sheet behaviour that spans millennia. There are two main approaches to reconstructing palaeo-ice sheets. Empirical reconstructions use the available glacial geological and chronological evidence to estimate ice sheet extent and dynamics but lack direct consideration of ice physics. In contrast, numerically modelled simulations implement ice physics, but often lack direct quantitative comparison with empirical evidence. Despite being long identified as a fruitful scientific endeavour, few ice sheet reconstructions attempt to reconcile the empirical and model-based approaches. To achieve this goal, model-data comparison procedures are required. Here, we compare three numerically modelled simulations of the former British–Irish Ice Sheet with the following lines of evidence: (a) position and shape of former margin positions, recorded by moraines; (b) former ice-flow direction and flow-switching, recorded by flowsets of subglacial bedforms; and (c) the timing of ice-free conditions, recorded by geochronological data. These model–data comparisons provide a useful framework for quantifying the degree of fit between numerical model simulations and empirical constraints. Such tools are vital for reconciling numerical modelling and empirical evidence, the combination of which will lead to more robust palaeo-ice sheet reconstructions with greater explicative and ultimately predictive power.     

北方泥炭沼泽中酚类物质释放的热敏感性研究:以英国威尔士地区雨养和矿养泥炭沼泽为例

利用高精度多梯度热培养系统,对英国威尔士地区两类典型的北方泥炭沼泽中酚类物质和溶解有机碳释放的热敏感性进行了为期1年的调查研究.研究结果表明,雨养泥炭沼泽(Bog)酚类物质释放的热敏感性Q10(总酚)月变化为0.92～1.57,而矿养泥炭沼泽(Fen)的月变化范围为0.93～1.30.酚类释放的热敏感性与土壤温度大致呈正相关关系.此外,土壤温度与溶解有机碳释放的热敏感性Q10(DOC)也呈正相关关系.总体上,雨养泥炭沼泽的热敏感性比矿养泥炭沼泽略高.除温度外,水文条件和植被等多种环境因子对酚类物质和溶解有机碳释放的热敏感性也有较大影响.研究结果还表明酚碳热敏感性比值Q10o(总酚)/Q10(DOC)相对稳定,土壤温度的变化对其影响有限.     

Stable isotope study of the Archaean rocks of the Vredefort impact structure,central Kaapvaal Craton,South Africa

The Vredefort dome in the Kaapvaal Craton was formed as a result of the impact of a large meteorite at 2.02 Ga. The central core of Archaean granitic basement rocks is surrounded by a collar of uplifted and overturned strata of the Witwatersrand Supergroup, exposing a substantial depth section of the Archaean crust. Orthogneisses of the core show little variation in whole-rock δ ¹⁸O value, with the majority being between 8 and 10‰, with a mean of 9.2‰ (n = 35). Quartz and feldspar have per mil differences that are consistent with O-isotope equilibrium at high temperatures, suggesting minimal interaction with fluids during subsequent cooling. These data refute previous suggestions that the Outer Granite Gneiss (OGG) and Inlandsee Leucogranofels (ILG) of the core represent middle and lower crust, respectively. Granulite-facies greenstone remnants from the ILG have δ ¹⁸O values that are on average 1.5‰ higher than the ILG host rocks and are unlikely, therefore, to represent the residuum from the partial melting event that formed the host rock. Witwatersrand Supergroup sedimentary rocks of the collar, which were metamorphosed at greenschist-to amphibolite-facies conditions, generally have lower δ ¹⁸O values than the core rocks with a mean value for metapelites of 7.7‰ (n = 45). Overall, through an ∼20 km thick section of crust, there is a general increase in whole-rock δ ¹⁸O value with increasing depth. This is the reverse of what is normal in the crust, largely because the collar rocks have δ ¹⁸O values that are unusually low in comparison with metamorphosed sedimentary rocks worldwide. The collar rocks have δD values ranging from −35 to −115‰ (average −62‰, n = 29), which are consistent with interaction with water of meteoric origin, having a δD of about −25 to −45‰. We suggest that fluid movement through the collar rocks was enhanced by impact-induced secondary permeability in the dome structure. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Gaseous mercury emissions from urban surfaces: Controls and spatiotemporal trends

The spatial and temporal variability of Hg emissions from urban paved surfaces was assessed through repeated measurements under varying environmental conditions at six sample sites in Toronto, Ontario, Canada. The results show significant spatial variability of the Hg emissions with median values ranging from below detection limit to 5.2 ng/m²/h. Two of the sites consistently had higher Hg emissions (on several occasions >20 ng/m²/h) than the other 4, which were equivalently low (maximum emission: 2.1 ng/m²/h). A surrogate measure of the pavement Hg concentrations was obtained during each day of sampling through the collection of street dust. The median street dust concentration also showed significant spatial variability (ranging from 9.6 to 44.5 ng/g). Regression analysis showed that the spatial variability of the Hg emissions was significantly related to the street dust concentrations. Controlled experiments using Hg amended street dust confirmed the relationship between Hg surface concentration and emission magnitude. Within a given sample site, Hg emissions varied temporally and multiple regression analysis showed that within-site variability was significantly influenced by changes in solar radiation with only a minor effect from surface temperature. Controlled experiments using shade cloths confirmed that solar radiation can have a large influence on the magnitude of Hg emissions within a given site. The emissions measured in Toronto were contextualized through comparison sampling in Austin, Texas. The Hg emissions measured in Austin were within the range detected in Toronto and also showed significant correlation with Hg street dust concentrations between sites. To provide a holistic assessment of Hg emissions from urban environments, samples were also collected from other common urban surfaces (soil, roofs, and windows). Soils consistently had higher emissions than all the other surfaces (7.3 ng/m²/h, n = 39).     

Assimilation of carbonate country rock by the parent magma of the Panzhihua Fe-Ti-V deposit（SW China）：Evidence from stable isotopes

The Panzhihua intrusion in southwest China is part of the Emeishan Large Igneous Province and host of a large Fe-Ti-V ore deposit.During emplacement of the main intrusion,multiple generations of mafic dykes invaded carbonate wall rocks,producing a large contact aureole.We measured the oxygen-isotope composition of the intrusions,their constituent minerals,and samples of the country rock.Magnetite and plagioclase from Panzhihua intrusion haveδ¹⁸O values that are consistent with magmatic equilibrium, and formed from magmas withδ¹⁸O values that were 1-2‰higher than expected in a mantle-derived magma.The unmetamorphosed country rock has highδ¹⁸O values,ranging from 13.2‰（sandstone） to 24.6-28.6‰（dolomite）.The skarns and marbles from the aureole have lowerδ¹⁸O andδ¹³C values than their protolith suggesting interaction with fluids that were in exchange equilibrium with the adjacent mafic magmas and especially the numerous mafic dykes that intruded the aureole.This would explain the alteration ofδ¹⁸O of the dykes which have significantly higher values than expected for a mantle-derived magma.Depending on the exactδ¹⁸O values assumed for the magma and contaminant, the amount of assimilation required to produce the elevatedδ¹⁸O value of the Panzhihua intrusion was between 8 and 13.7 wt.%,assuming simple mixing.The exact mechanism of contamination is unclear but may involve a combination of assimilation of bulk country rock,mixing with a melt of the country rock and exchange with CO₂-rich fluid derived from decarbonation of the marls and dolomites.These mechanisms,particularly the latter,were probably involved in the formation of the Fe-Ti-V ores.     

Wetland hydroperiod classification in the western prairies using multitemporal synthetic aperture radar

Wetlands represent one of the world's most biodiverse and threatened ecosystem types and were diminished globally by about two‐thirds in the 20th century. There is continuing decline in wetland quantity and function due to infilling and other human activities. In addition, with climate change, warmer temperatures and changes in precipitation and evapotranspiration are reducing wetland surface and groundwater supplies, further altering wetland hydrology and vegetation. There is a need to automate inventory and monitoring of wetlands, and as a study system, we investigated the Shepard Slough wetlands complex, which includes numerous wetlands in urban, suburban, and agricultural zones in the prairie pothole region of southern Alberta, Canada. Here, wetlands are generally confined to depressions in the undulating terrain, challenging wetlands inventory and monitoring. This study applied threshold and frequency analysis routines for high‐resolution, single‐polarization (HH) RADARSAT‐2, synthetic aperture radar mapping. This enabled a growing season surface water extent hyroperiod‐based wetland classification, which can support water and wetland resource monitoring. This 3‐year study demonstrated synthetic aperture radar‐derived multitemporal open‐water masks provided an effective index of wetland permanence class, with overall accuracies of 89% to 95% compared with optical validation data, and RMSE between 0.2 and 0.7 m between model and field validation data. This allowed for characterizing the distribution and dynamics of 4 marsh wetlands hydroperiod classes, temporary, seasonal, semipermanent, and permanent, and mapping of the sequential vegetation bands that included emergent, obligate wetland, facultative wetland, and upland plant communities. Hydroperiod variation and surface water extent were found to be influenced by short‐term rainfall events in both wet and dry years. Seasonal hydroperiods in wetlands were particularly variable if there was a decrease in the temporary or semipermanent hydroperiod classes. In years with extreme rain events, the temporary wetlands especially increased relative to longer lasting wetlands (84% in 2015 with significant rainfall events, compared with 42% otherwise).     

Sedimentology of the continental end‐Permian extinction event in the Sydney Basin,eastern Australia

Upper Permian to Lower Triassic coastal plain successions of the Sydney Basin in eastern Australia have been investigated in outcrop and continuous drillcores. The purpose of the investigation is to provide an assessment of palaeoenvironmental change at high southern palaeolatitudes in a continental margin context for the late Permian (Lopingian), across the end‐Permian Extinction interval, and into the Early Triassic. These basins were affected by explosive volcanic eruptions during the late Permian and, to a much lesser extent, during the Early Triassic, allowing high‐resolution age determination on the numerous tuff horizons. Palaeobotanical and radiogenic isotope data indicate that the end‐Permian Extinction occurs at the top of the uppermost coal bed, and the Permo‐Triassic boundary either within an immediately overlying mudrock succession or within a succeeding channel sandstone body, depending on locality due to lateral variation. Late Permian depositional environments were initially (during the Wuchiapingian) shallow marine and deltaic, but coastal plain fluvial environments with extensive coal‐forming mires became progressively established during the early late Permian, reflected in numerous preserved coal seams. The fluvial style of coastal plain channel deposits varies geographically. However, apart from the loss of peat‐forming mires, no significant long‐term change in depositional style (grain size, sediment‐body architecture, or sediment dispersal direction) was noted across the end‐Permian Extinction (pinpointed by turnover of the palaeoflora). There is no evidence for immediate aridification across the boundary despite a loss of coal from these successions. Rather, the end‐Permian Extinction marks the base of a long‐term, progressive trend towards better‐drained alluvial conditions into the Early Triassic. Indeed, the floral turnover was immediately followed by a flooding event in basinal depocentres, following which fluvial systems similar to those active prior to the end‐Permian Extinction were re‐established. The age of the floral extinction is constrained to 252.54 ± 0.08 to 252.10 ± 0.06 Ma by a suite of new Chemical Abrasion Isotope Dilution Thermal Ionization Mass Spectrometry U‐Pb ages on zircon grains. Another new age indicates that the return to fluvial sedimentation similar to that before the end‐Permian Extinction occurred in the basal Triassic (prior to 251.51 ± 0.14 Ma). The character of the surface separating coal‐bearing pre‐end‐Permian Extinction from coal‐barren post‐end‐Permian Extinction strata varies across the basins. In basin‐central locations, the contact varies from disconformable, where a fluvial channel body has cut down to the level of the top coal, to conformable where the top coal is overlain by mudrocks and interbedded sandstone–siltstone facies. In basin‐marginal locations, however, the contact is a pronounced erosional disconformity with coarse‐grained alluvial facies overlying older Permian rocks. There is no evidence that the contact is everywhere a disconformity or unconformity.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

A Truncated Pareto Model to Estimate the Under Recovery of Large Diamonds

The metallurgical recovery processes in diamond mining may, under certain circumstances, cause an under-recovery of large diamonds. In order to predict high quantiles or tail probabilities we use a Bayesian approach to fit a truncated Generalized Pareto Type distribution to the tail of the data consisting of the weights of individual diamonds. Based on the estimated tail probability, the expected number of diamonds larger than a specified weight can be estimated. The difference between the expected and observed frequencies of diamond weights above an upper threshold provides an estimate of the number of diamonds lost during the recovery process.