Collision of the Ganges–Brahmaputra Delta with the Burma Arc: Implications for earthquake hazard

We take a fresh look at the topography, structure and seismicity of the Ganges–Brahmaputra Delta (GBD)–Burma Arc collision zone in order to reevaluate the nature of the accretionary prism and its seismic potential. The GBD, the world's largest delta, has been built from sediments eroded from the Himalayan collision. These sediments prograded the continental margin of the Indian subcontinent by  400 km, forming a huge sediment pile that is now entering the Burma Arc subduction zone. Subduction of oceanic lithosphere with > 20 km sediment thickness is fueling the growth of an active accretionary prism exposed on land. The prism starts at an apex south of the GBD shelf edge at  18°N and widens northwards to form a broad triangle that may be up to 300 km wide at its northern limit. The front of the prism is blind, buried by the GBD sediments. Thus, the deformation front extends 100 km west of the surface fold belt beneath the Comilla Tract, which is uplifted by 3–4 m relative to the delta. This accretionary prism has the lowest surface slope of any active subduction zone. The gradient of the prism is only  0.1°, rising to  0.5° in the forearc region to the east. This low slope is consistent with the high level of overpressure found in the subsurface, and indicates a very weak detachment. Since its onset, the collision of the GBD and Burma Arc has expanded westward at  2 cm/yr, and propagated southwards at  5 cm/yr. Seismic hazard in the GBD is largely unknown. Intermediate-size earthquakes are associated with surface ruptures and fold growth in the external part of the prism. However, the possibility of large subduction ruptures has not been accounted for, and may be higher than generally believed. Although sediment-clogged systems are thought to not be able to sustain the stresses and strain-weakening behavior required for great earthquakes, some of the largest known earthquakes have occurred in heavily-sedimented subduction zones. A large earthquake in 1762 ruptured  250 km of the southern part of the GBD, suggesting large earthquakes are possible there. A large, but poorly documented earthquake in 1548 damaged population centers at the northern and southern ends of the onshore prism, and is the only known candidate for a rupture of the plate boundary along the subaerial part of the GBD–Burma Arc collision zone.     

Differential developmental toxicity of naphthoic acid isomers in medaka (Oryzias latipes) embryos

Polycyclic aromatic hydrocarbons (PAHs) are widespread persistent pollutants that readily undergo biotic and abiotic conversion to numerous transformation products in rivers, lakes and estuarine sediments. Here we characterize the developmental toxicity of four PAH transformation products each structural isomers of hydroxynaphthoic acid: 1H2NA, 2H1NA, 2H3NA, and 6H2NA. Medaka fish (Oryzias latipes) embryos and eleutheroembryos were used to determine toxicity. A 96-well micro-plate format was used to establish a robust, statistically significant platform for assessment of early life stages. Individual naphthoic acid isomers demonstrated a rank order of toxicity with 1H2NA > 2H1NA > 2H3NA > 6H2NA being more toxic. Abnormalities of circulatory system were most pronounced including pericardial edema and tube heart. To determine if HNA isomers were AhR ligands, spatial-temporal expression and activity of CYP1A was measured via in vivo EROD assessments. qPCR measurement of CYP1A induction proved different between isomers dosed at respective concentrations affecting 50% of exposed individuals (EC50s). In vitro, all ANH isomers transactivated mouse AhR using a medaka CYP1A promoter specific reporter assay. Circulatory abnormalities followed P450 induction and response was consistent with PAH toxicity. A 96-well micro-plates proved suitable as exposure chambers and provided statistically sound evaluations as well as efficient toxicity screens. Our results demonstrate the use of medaka embryos for toxicity analysis thereby achieving REACH objectives for the reduction of adult animal testing in toxicity evaluations.     

Use of GRACE determined secular gravity rates for glacial isostatic adjustment studies in North-America

Monthly geopotential spherical harmonic coefficients from the GRACE satellite mission are used to determine their usefulness and limitations for studying glacial isostatic adjustment (GIA) in North-America. Secular gravity rates are estimated by unweighted least-squares estimation using release 4 coefficients from August 2002 to August 2007 provided by the Center for Space Research (CSR), University of Texas. Smoothing is required to suppress short wavelength noise, in addition to filtering to diminish geographically correlated errors, as shown in previous studies. Optimal cut-off degrees and orders are determined for the destriping filter to maximize the signal to noise ratio. The halfwidth of the Gaussian filter is shown to significantly affect the sensitivity of the GRACE data (with respect to upper mantle viscosity and ice loading history). Therefore, the halfwidth should be selected based on the desired sensitivity.It is shown that increase in water storage in an area south west of Hudson Bay, from the summer of 2003 to the summer of 2006, contributes up to half of the maximum estimated gravity rate. Hydrology models differ in the predictions of the secular change in water storage, therefore even 4-year trend estimates are influenced by the uncertainty in water storage changes. Land ice melting in Greenland and Alaska has a non-negligible contribution, up to one-fourth of the maximum gravity rate.The estimated secular gravity rate shows two distinct peaks that can possibly be due to two domes in the former Pleistocene ice cover: west and south east of Hudson Bay. With a limited number of models, a better fit is obtained with models that use the ICE-3G model compared to the ICE-5G model. However, the uncertainty in interannual variations in hydrology models is too large to constrain the ice loading history with the current data span. For future work in which GRACE will be used to constrain ice loading history and the Earth's radial viscosity profile, it is important to include realistic uncertainty estimates for hydrology models and land ice melting in addition to the effects of lateral heterogeneity.     

Toxicities of antifouling biocide Irgarol 1051 and its major degraded product to marine primary producers

Irgarol 1051 (2-methythiol-4-tert-butylamino-6-cyclopropylamino-s-triazine) is an algaecide commonly used in antifouling paints. It undergoes photodegradation which yields M1 (2-methylthio-4-tert-butylamino-6-amino-s-triazine) as its major and most stable degradant. Elevated levels of both Irgarol and M1 have been detected in coastal waters worldwide; however, ecotoxicity effects of M1 to various marine autotrophs such as cyanobacteria are still largely unknown. This study firstly examined and compared the 96 h toxicities of Irgarol and M1 to the cyanobacterium Chroococcus minor and two marine diatom species, Skeletonema costatum and Thalassiosira pseudonana. Our results suggested that Irgarol was consistently more toxic to all of the three species than M1 (96 h EC50 values: C. minor, 7.71 microug L(-1) Irgarol vs. > 200 microg L(-1) M1; S. costatum, 0.29 microg L(-1) Irgarol vs. 11.32 microg L(-1)M1; and T. pseudonana, 0.41 microg L(-1) Irgarol vs. 16.50 microg L(-1)M1). Secondly, we conducted a meta-analysis of currently available data on toxicities of Irgarol and M1 to both freshwater and marine primary producers based on species sensitivity distributions (SSDs). Interestingly, freshwater autotrophs are more sensitive to Irgarol than their marine counterparts. For marine autotrophs, microalgae are generally more sensitive to Irgarol than macroalgae and cyanobacteria. With very limited available data on M1 (i.e. five species), M1 might be less toxic than Irgarol; nonetheless this finding warrants further confirmation with additional data on other autotrophic species.     

Evolution of the variability of surface temperature and vegetation density in the great plains

This study focuses on how the variability of land surface temperature and vegetation density at the SGP ARM-CART site changes over episodic (day to day) and seasonal time scales using AVHRR satellite data. Four drying periods throughout the year are analyzed. Land surface temperature had an erratic relationship with time exhibiting no deterministic pattern from day-to-day or season-to-season. Furthermore, it did not exhibit spatial pattern persistence. On the other hand, vegetation density had a consistent spatial pattern and temporal decay during average length drying periods (less than 7 days) as well as within a season. However, there were distinct differences in the seasonal pattern of variation between winter and growing seasons. In addition, the paper highlights a methodology to quantify the relationships that exist at the land surface between the primary parameter of interest and the controlling variables.     

Flow and deposition of pyroclastic granular flows: A type example from the 1975 Ngauruhoe eruption,New Zealand

Small-volume pyroclastic density currents (PDCs) are generated frequently during explosive eruptions with little warning. Assessing their hazard requires a physical understanding of their transport and sedimentation processes which is best achieved by the testing of experimental and numerical models of geophysical mass flows against natural flows and/or deposits. To this end we report on one of the most detailed sedimentological studies ever carried out on a series of pristine small-volume PDC deposits from the 1975 eruption of Ngauruhoe volcano, whose emplacement were also witnessed during eruption. Using high-resolution GPS surveys, a series of lateral excavations across the deposits, and bulk sedimentological analysis we constrained the geomorphology, internal structure and texture of the deposits with respect to laterally varying modes of deposition.     

Spatial and temporal analysis of hillslope–channel coupling and implications for the longitudinal profile in a dryland basin

The long‐term evolution of channel longitudinal profiles within drainage basins is partly determined by the relative balance of hillslope sediment supply to channels and the evacuation of channel sediment. However, the lack of theoretical understanding of the physical processes of hillslope–channel coupling makes it challenging to determine whether hillslope sediment supply or channel sediment evacuation dominates over different timescales and how this balance affects bed elevation locally along the longitudinal profile. In this paper, we develop a framework for inferring the relative dominance of hillslope sediment supply to the channel versus channel sediment evacuation, over a range of temporal and spatial scales. The framework combines distinct local flow distributions on hillslopes and in the channel with surface grain‐size distributions. We use these to compute local hydraulic stresses at various hillslope‐channel coupling locations within the Walnut Gulch Experimental Watershed (WGEW) in southeast Arizona, USA. These stresses are then assessed as a local net balance of geomorphic work between hillslopes and channel for a range of flow conditions generalizing decadal historical records. Our analysis reveals that, although the magnitude of hydraulic stress in the channel is consistently higher than that on hillslopes, the product of stress magnitude and frequency results in a close balance between hillslope supply and channel evacuation for the most frequent flows. Only at less frequent, high‐magnitude flows do channel hydraulic stresses exceed those on hillslopes, and channel evacuation dominates the net balance. This result suggests that WGEW exists mostly (~50% of the time) in an equilibrium condition of sediment balance between hillslopes and channels, which helps to explain the observed straight longitudinal profile. We illustrate how this balance can be upset by climate changes that differentially affect relative flow regimes on slopes and in channels. Such changes can push the long profile into a convex or concave condition. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.