Rice crop phenology mapping at high spatial and temporal resolution using downscaled MODIS time-series

Satellite data holds considerable potential as a source of information on rice crop growth which can be used to inform agronomy. However, given the typical field sizes in many rice-growing countries such as China, data from coarse spatial resolution satellite systems such as the Moderate Resolution Imaging Spectroradiometer (MODIS) are inadequate for resolving crop growth variability at the field scale. Nevertheless, systems such as MODIS do provide images with sufficient frequency to be able to capture the detail of rice crop growth trajectories throughout a growing season. In order to generate high spatial and temporal resolution data suitable for mapping rice crop phenology, this study fused MODIS data with lower frequency, higher spatial resolution Landsat data. An overall workflow was developed which began with image preprocessing, calculation of multi-temporal normalized difference vegetation index (NDVI) images, and spatiotemporal fusion of data from the two sensors. The Spatial and Temporal Adaptive Reflectance Fusion Model was used to effectively downscale the MODIS data to deliver a time-series of 30 m spatial resolution NDVI data at 8-day intervals throughout the rice-growing season. Zonal statistical analysis was used to extract NDVI time-series for individual fields and signal filtering was applied to the time-series to generate rice phenology curves. The downscaled MODIS NDVI products were able to characterize the development of paddy rice at fine spatial and temporal resolutions, across wide spatial extents over multiple growing seasons. These data permitted the extraction of key crop seasonality parameters that quantified inter-annual growth variability for a whole agricultural region and enabled mapping of the variability in crop performance between and within fields. Hence, this approach can provide rice crop growth data that is suitable for informing agronomic policy and practice across a wide range of scales.     

Unveiling buried aeolian landscapes: reconstructing a late Holocene dune environment using 3D ground-penetrating radar

Across the UK, sandy beaches and dunes protect coastal infrastructure from waves and extreme water levels during large-scale storms, while providing important habitats and recreational opportunities. Understanding their long-term evolution is vital in managing their condition in a changing climate. Recently, ground-penetrating radar (GPR) methods have grown in popularity in geomorphological applications, yielding centimetre-scale resolution images of near-surface stratigraphy and structure, thus allowing landscape evolution to be reconstructed. Additionally, abrupt changes in palaeo-environments can be visualized in three dimensions. Although often complemented by core data, GPR allows interpretations to be extended into areas with minimal ground-truth control. Nonetheless, GPR data interpretation can be non-intuitive and ambiguous, and radargrams may not initially resemble the expected subsurface geometry. Interpretation can be made yet more onerous when handling the large 3D data volumes that are facilitated with modern GPR technology. Here we describe the development of novel semi-automated GPR feature-extraction tools, based on ‘edge detection’ and ‘thresholding’ methods, which detect regions of increased GPR reflectivity which can be applied to aid in the reconstruction of a range Quaternary landscapes. Since reflectivity can be related to lithological and/or pore fluid changes, the 3D architecture of the palaeo-landscape can be reconstructed from the features extracted from a geophysical dataset. We present 500 MHz GPR data collected over a buried Holocene coastal dune system in North Wales, UK, now reclaimed for use as an airfield. Core data from the site, reaching a maximum depth 2 m, suggest rapid vertical changes from sand to silty-organic units, and GPR profiles suggest similar lateral complexity. By applying thresholding methods to GPR depth slices, these lateral complexities are effectively and automatically mapped. Furthermore, automatic extraction of the local reflection power yields a strong correlation with the depth variation of organic content, suggesting it is a cause of reflectivity contrast. GPR-interpolated analyses away from core control thus offer a powerful proxy for parameters derived from invasive core logging. The GPR data collected at Llanbedr airfield highlight a complex dune system to a depth of 2.8 m, probably deposited in several phases over ~700 years, similar to elsewhere in North Wales.     

Rapid U‐Pb Geochronology by Laser Ablation Multi‐Collector ICP‐MS

Detrital zircon (DZ) U‐Pb laser ablation‐inductively coupled plasma‐mass spectrometry (LA‐ICP‐MS) has revolutionised the way geologists approach many Earth science questions. Although recent research has focused on rapid sample throughput, acquisition rates are limited to 100–300 analyses h^?1. We present a method to acquire zircon U‐Pb dates at rates of 120, 300, 600 and 1200 analyses h^?1 (30, 12, 6 and 3 s per analysis) by multi‐collector LA‐ICP‐MS. We demonstrate the efficacy of this method by analysing twelve zircon reference materials with dates from ~ 3465 to ~ 28 Ma. Mean offset from high‐precision dates increases with faster rates from 0.9% to 1.1%; mean random 1s uncertainty increases from 0.6% to 1.3%. We tested this new method on a sandstone sample previously characterised by large‐n DZ geochronology. Quantitative comparison shows increased correspondence among age distributions comprising > 300 dates. This new method holds promise for DZ geochronology because (a) it requires no major changes to hardware, but rather modifications to software; (b) it yields robust age distributions well‐suited for quantitative analysis and maximum depositional age calculations; (c) there is only a minor sacrifice of accuracy and measurement uncertainty; and (d) there is less burden to researchers in terms of time investment and analytical cost.     

When hazard avoidance is not an option: lessons learned from monitoring the postdisaster Oso landslide,USA

On 22 March 2014, a massive, catastrophic landslide occurred near Oso, Washington, USA, sweeping more than 1 km across the adjacent valley flats and killing 43 people. For the following 5 weeks, hundreds of workers engaged in an exhaustive search, rescue, and recovery effort directly in the landslide runout path. These workers could not avoid the risks posed by additional large-scale slope collapses. In an effort to ensure worker safety, multiple agencies cooperated to swiftly deploy a monitoring and alerting system consisting of sensors, automated data processing and web-based display, along with defined communication protocols and clear calls to action for emergency management and search personnel. Guided by the principle that an accelerating landslide poses a greater threat than a steadily moving or stationary mass, the system was designed to detect ground motion and vibration using complementary monitoring techniques. Near real-time information was provided by continuous GPS, seismometers/geophones, and extensometers. This information was augmented by repeat-assessment techniques such as terrestrial and aerial laser scanning and time-lapse photography. Fortunately, no major additional landsliding occurred. However, we did detect small headscarp failures as well as slow movement of the remaining landslide mass with the monitoring system. This was an exceptional response situation and the lessons learned are applicable to other landslide disaster crises. They underscore the need for cogent landslide expertise and ready-to-deploy monitoring equipment, the value of using redundant monitoring techniques with distinct goals, the benefit of clearly defined communication protocols, and the importance of continued research into forecasting landslide behavior to allow timely warning.

     

Estimation of peak relative velocity and peak absolute acceleration of linear SDOF systems

In the seismic analysis and design of structures, the true velocity and absolute acceleration are usually approximated by their corresponding pseudo-values. This approach is simple and works well for structures with small damping (say, less than 15%). When the damping of a structure is enhanced for the purpose of response reduction, it may result in large analysis and design errors. Based on theory of random vibration and the established mechanism of seismic response spectra analysis, a method is developed (1) to predict the relative velocity spectra with any damping ratio level directly from the 5% standard pseudo-acceleration spectrum; and (2) to estimate the peak absolute acceleration. The accuracy of both is validated by using two selected ensembles of ground motion records.     

Foreword

Destructive earthquakes have caused great damage in China and the United States and collapsing buildings have caused many deaths and injuries. The field of earthquake engineering studies     

Performance of pile-restrained flexible floating breakwaters

The overall performance of pile-restrained flexible floating breakwaters is investigated under the action of linear monochromatic incident waves in the frequency domain. The aforementioned floating breakwaters undergo only vertical structural deflections along their length and are held in place by means of vertical piles. The total number of degrees of freedom equals the six conventional body modes, when the breakwater moves as a rigid body, plus the extra bending modes. These bending modes are introduced to represent the structural deflections of the floating breakwater and are described by the Bernoulli–Euler flexible beam equation. The number of bending modes introduced is determined through an appropriate iterative procedure. The hydrostatic coefficients corresponding to the bending modes are also derived. The numerical analysis of the flexible floating breakwaters is based on a three-dimensional hydrodynamic formulation of the floating body. A parametric study is carried out for a wide range of structural stiffness parameters and wave headings, to investigate their effect on the performance of flexible floating breakwaters. Moreover, this performance is compared with that of the corresponding pile-restrained rigid floating breakwater. Results indicated that the degree of structural stiffness and the wave heading strongly affect the performance of flexible floating breakwaters. The existence of an “optimum” value of structural stiffness is demonstrated for the entire wave frequency range.     

Clinoform architecture and along‐strike facies variability through an exhumed erosional to accretionary basin margin transition

Exhumed basin margin‐scale clinothems provide important archives for understanding process interactions and reconstructing the physiography of sedimentary basins. However, studies of coeval shelf through slope to basin‐floor deposits are rarely documented, mainly due to outcrop or subsurface dataset limitations. Unit G from the Laingsburg depocentre (Karoo Basin, South Africa) is a rare example of a complete basin margin scale clinothem (>60 km long, 200 m‐high), with >10 km of depositional strike control, which allows a quasi‐3D study of a preserved shelf‐slope‐basin floor transition over a ca. 1,200 km² area. Sand‐prone, wave‐influenced topset deposits close to the shelf‐edge rollover zone can be physically mapped down dip for ca. 10 km as they thicken and transition into heterolithic foreset/slope deposits. These deposits progressively fine and thin over tens of km farther down dip into sand‐starved bottomset/basin‐floor deposits. Only a few km along strike, the coeval foreset/slope deposits are bypass‐dominated with incisional features interpreted as minor slope conduits/gullies. The margin here is steeper, more channelized and records a stepped profile with evidence of sand‐filled intraslope topography, a preserved base‐of‐slope transition zone and sand‐rich bottomset/basin‐floor deposits. Unit G is interpreted as part of a composite depositional sequence that records a change in basin margin style from an underlying incised slope with large sand‐rich basin‐floor fans to an overlying accretion‐dominated shelf with limited sand supply to the slope and basin floor. The change in margin style is accompanied with decreased clinoform height/slope and increased shelf width. This is interpreted to reflect a transition in subsidence style from regional sag, driven by dynamic topography/inherited basement configuration, to early foreland basin flexural loading. Results of this study caution against reconstructing basin margin successions from partial datasets without accounting for temporal and spatial physiographic changes, with potential implications on predictive basin evolution models.     

Stable isotope constraints on the fluid source of hydrothermal breccia pipes in the Tankwa Karoo depocentre,South Africa: Breakdown of authigenic minerals during sill intrusion

The Karoo Basin covers much of South Africa and is an area of prospective shale gas exploration, with the Whitehill Formation the target shale unit. However, the sedimentary succession, including the Whitehill, has been intruded by a series of sills and dykes associated with the Karoo Large Igneous Province (~183 Ma), which are expected to have modified the thermal history of the basin dramatically. Here, we investigate a secondary effect of these intrusions: a series of hydrothermal vent complexes, or breccia pipes, focusing on using O, H, and C isotopes to constrain the origin and evolution of fluids produced during the intrusion of basaltic sills. A cluster of breccia pipes have been eroded down to the level of the Ecca Group at Luiperdskop on the western edge of the Karoo basin; a small isolated pipe of similar appearance crops out 13 km to the east. The Luiperdskop pipes are underlain by a Karoo dolerite sill that is assumed to provide the heat driving fluidization. The pipes consist of fine‐grained matrix and about 8% clasts, on average, of mostly sedimentary material; occasional large rafts of quartzite and dolerite are also present. The presence of clasts apparently from the Dwyka Group is consistent with the depth of formation of the pipes being at, or near, the base of the Karoo Supergroup, between 400 and 850 m below present surface. The presence of chlorite as the dominant hydrous mineral is consistent with an emplacement temperature between 300 and 350°C. The major and trace element, and O‐ and H‐isotope composition of the Tankwa breccias is homogenous, consistent with them being derived from the same source. The δ¹⁸O values (vsVSMOW) of the breccias are relatively uniform (7.1‰–8.7‰), and are similar to that of the country rock shale, and both are lower than expected for shale. The water content of the breccia is between 2.7 and 3.1 wt.% and the δD values range from ?109‰ to ?144‰. Calcite in vesicles has δ¹³C and δ¹⁸O (VSMOW) values of ?4.2‰ and 24.0‰, respectively. The low δD value of the breccia rocks does not appear to be due to the presence of methane in the fluid. Instead, it is proposed that low δD and δ¹⁸O values are the result of the fluid being derived from the breakdown of clay minerals that formed and were deposited at a time of cold climate at ~290 Ma.