@article{ZornLeCorvecVarleyetal.2019, author = {Zorn, Edgar Ulrich and Le Corvec, Nicolas and Varley, Nick R. and Salzer, Jacqueline T. and Walter, Thomas R. and Navarro-Ochoa, Carlos and Vargas-Bracamontes, Dulce M. and Thiele, Samuel T. and Ar{\´a}mbula Mendoza, Ra{\´u}l}, title = {Load stress controls on directional lava dome growth at Volcan de Colima, Mexico}, series = {Frontiers in Earth Science}, volume = {7}, journal = {Frontiers in Earth Science}, publisher = {Frontiers Media}, address = {Lausanne}, issn = {2296-6463}, doi = {10.3389/feart.2019.00084}, pages = {18}, year = {2019}, abstract = {During eruptive activity of andesitic stratovolcanoes, the extrusion of lava domes, their collapse and intermittent explosions are common volcanic hazards. Many lava domes grow in a preferred direction, in turn affecting the direction of lava flows and pyroclastic density currents. Access to active lava domes is difficult and hazardous, so detailed data characterizing lava dome growth are typically limited, keeping the processes controlling the directionality of extrusions unclear. Here we combine TerraSAR-X satellite radar observations with high-resolution airborne photogrammetry to assess morphological changes, and perform finite element modeling to investigate the impact of loading stress on shallow magma ascent directions associated with lava dome extrusion and crater formation at Volcan de Colima, Mexico. The TerraSAR-X data, acquired in similar to 1-m resolution spotlight mode, enable us to derive a chronology of the eruptive processes from intensity-based time-lapse observations of the general crater and dome evolution. The satellite images are complemented by close-range airborne photos, processed by the Structure-from-Motion workflow. This allows the derivation of high-resolution digital elevation models, providing insight into detailed loading and unloading features. During the observation period from Jan-2013 to Feb-2016, we identify a dominantly W-directed dome growth and lava flow production until Jan-2015. In Feb-2015, following the removal of the active summit dome, the surface crater widened and elongated along a NE-SW axis. Later in May-2015, a new dome grew toward the SW of the crater while a separate vent developed in the NE of the crater, reflecting a change in the direction of magma ascent and possible conduit bifurcation. Finite element models show a significant stress change in agreement with the observed magma ascent direction changes in response to the changing surface loads, both for loading (dome growth) and unloading (crater forming excavation) cases. These results allow insight into shallow dome growth dynamics and the migration of magma ascent in response to changing volcano summit morphology. They further highlight the importance of detailed volcano summit morphology surveillance, as changes in direction or location of dome extrusion may have major implications regarding the directions of potential volcanic hazards, such as pyroclastic density currents generated by dome collapse.}, language = {en} } @article{ZhuChenLiuetal.2019, author = {Zhu, Hui and Chen, Lunjin and Liu, Xu and Shprits, Yuri}, title = {Modulation of locally generated equatorial noise by ULF wave}, series = {Journal of geophysical research : Space physics}, volume = {124}, journal = {Journal of geophysical research : Space physics}, number = {4}, publisher = {American Geophysical Union}, address = {Washington}, issn = {2169-9380}, doi = {10.1029/2018JA026199}, pages = {2779 -- 2787}, year = {2019}, abstract = {In this paper we report a rare and fortunate event of fast magnetosonic (MS, also called equatorial noise) waves modulated by compressional ultralow frequency (ULF) waves measured by Van Allen Probes. The characteristics of MS waves, ULF waves, proton distribution, and their potential correlations are analyzed. The results show that ULF waves can modulate the energetic ring proton distribution and in turn modulate the MS generation. Furthermore, the variation of MS intensities is attributed to not only ULF wave activities but also the variation of background parameters, for example, number density. The results confirm the opinion that MS waves are generated by proton ring distribution and propose a new modulation phenomenon.}, language = {en} } @article{ZhuPilzCotton2019, author = {Zhu, Chuanbin and Pilz, Marco and Cotton, Fabrice}, title = {Which is a better proxy, site period or depth to bedrock, in modelling linear site response in addition to the average shear-wave velocity?}, series = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, volume = {18}, journal = {Bulletin of earthquake engineering : official publication of the European Association for Earthquake Engineering}, number = {3}, publisher = {Springer}, address = {Dordrecht}, issn = {1570-761X}, doi = {10.1007/s10518-019-00738-6}, pages = {797 -- 820}, year = {2019}, abstract = {This study aims to identify the best-performing site characterization proxy alternative and complementary to the conventional 30 m average shear-wave velocity V-S30, as well as the optimal combination of proxies in characterizing linear site response. Investigated proxies include T-0 (site fundamental period obtained from earthquake horizontal-to-vertical spectral ratios), V-Sz (measured average shear-wave velocities to depth z, z = 5, 10, 20 and 30 m), Z(0.8) and Z(1.0) (measured site depths to layers having shear-wave velocity 0.8 and 1.0 km/s, respectively), as well as Z(x-infer) (inferred site depths from a regional velocity model, x = 0.8 and 1.0, 1.5 and 2.5 km/s). To evaluate the performance of a site proxy or a combination, a total of 1840 surface-borehole recordings is selected from KiK-net database. Site amplifications are derived using surface-to-borehole response-, Fourier- and cross-spectral ratio techniques and then are compared across approaches. Next, the efficacies of 7 single-proxies and 11 proxy-pairs are quantified based on the site-to-site standard deviation of amplification residuals of observation about prediction using the proxy or the pair. Our results show that T-0 is the best-performing single-proxy among T-0, Z(0.8), Z(1.0) and V-Sz. Meanwhile, T-0 is also the best-performing proxy among T-0, Z(0.8), Z(1.0) and Z(x-infer) complementary to V-S30 in accounting for the residual amplification after V-S30-correction. Besides, T-0 alone can capture most of the site effects and should be utilized as the primary site indicator. Though (T-0, V-S30) is the best-performing proxy pair among (V-S30, T-0), (V-S30, Z(0.8)), (V-S30, Z(1.0)), (V-S30, Z(x-infer)) and (T-0, V-Sz), it is only slightly better than (T-0, V-S20). Considering both efficacy and engineering utility, the combination of T-0 (primary) and V-S20 (secondary) is recommended. Further study is needed to test the performances of various proxies on sites in deep sedimentary basins.}, language = {en} } @article{ZhuCottonPilz2019, author = {Zhu, Chuanbin and Cotton, Fabrice and Pilz, Marco}, title = {Testing the Depths to 1.0 and 2.5 km/s Velocity Isosurfaces in a Velocity Model for Japan and Implications for Ground-Motion Modeling}, series = {Bulletin of the Seismological Society of America}, volume = {109}, journal = {Bulletin of the Seismological Society of America}, number = {6}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0037-1106}, doi = {10.1785/0120190016}, pages = {2710 -- 2721}, year = {2019}, abstract = {In the Next Generation Attenuation West2 (NGA-West2) project, a 3D subsurface structure model (Japan Seismic Hazard Information Station [J-SHIS]) was queried to establish depths to 1.0 and 2.5 km/s velocity isosurfaces for sites without depth measurement in Japan. In this article, we evaluate the depth parameters in the J-SHIS velocity model by comparing them with their corresponding site-specific depth measurements derived from selected KiK-net velocity profiles. The comparison indicates that the J-SHIS model underestimates site depths at shallow sites and overestimates depths at deep sites. Similar issues were also identified in the southern California basin model. Our results also show that these underestimations and over-estimations have a potentially significant impact on ground-motion prediction using NGA-West2 ground-motion models (GMMs). Site resonant period may be considered as an alternative to depth parameter in the site term of a GMM.}, language = {en} } @article{ZhelayskayaVasileShpritsetal.2019, author = {Zhelayskaya, Irina S. and Vasile, Ruggero and Shprits, Yuri and Stolle, Claudia and Matzka, J{\"u}rgen}, title = {Systematic Analysis of Machine Learning and Feature Selection Techniques for Prediction of the Kp Index}, series = {Space Weather: The International Journal of Research and Applications}, volume = {17}, journal = {Space Weather: The International Journal of Research and Applications}, number = {10}, publisher = {American Geophysical Union}, address = {Washington}, issn = {1542-7390}, doi = {10.1029/2019SW002271}, pages = {1461 -- 1486}, year = {2019}, abstract = {The Kp index is a measure of the midlatitude global geomagnetic activity and represents short-term magnetic variations driven by solar wind plasma and interplanetary magnetic field. The Kp index is one of the most widely used indicators for space weather alerts and serves as input to various models, such as for the thermosphere and the radiation belts. It is therefore crucial to predict the Kp index accurately. Previous work in this area has mostly employed artificial neural networks to nowcast Kp, based their inferences on the recent history of Kp and on solar wind measurements at L1. In this study, we systematically test how different machine learning techniques perform on the task of nowcasting and forecasting Kp for prediction horizons of up to 12 hr. Additionally, we investigate different methods of machine learning and information theory for selecting the optimal inputs to a predictive model. We illustrate how these methods can be applied to select the most important inputs to a predictive model of Kp and to significantly reduce input dimensionality. We compare our best performing models based on a reduced set of optimal inputs with the existing models of Kp, using different test intervals, and show how this selection can affect model performance.}, language = {en} } @article{ZhangHuangZhangetal.2019, author = {Zhang, Yang and Huang, Wentao and Zhang, Yuanyuan and Poujol, Marc and Guillot, Stephane and Roperch, Pierrick and Dupont-Nivet, Guillaume and Guo, Zhaojie}, title = {Detrital zircon provenance comparison between the Paleocene-Eocene Nangqian-Xialaxiu and Gongjue basins: New insights for Cenozoic paleogeographic evolution of the eastern Tibetan Plateau}, series = {Palaeogeography, palaeoclimatology, palaeoecology : an international journal for the geo-sciences}, volume = {533}, journal = {Palaeogeography, palaeoclimatology, palaeoecology : an international journal for the geo-sciences}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0031-0182}, doi = {10.1016/j.palaeo.2019.109241}, pages = {19}, year = {2019}, abstract = {Paleogeographic reconstructions of terranes can greatly benefit from the provenance analysis of sediments. A series of Cenozoic basins provide key sedimentary archives for investigating the growth of the Tibetan Plateau, yet the provenance of the sediments in these basins has never been constrained robustly. Here we report sedimentary petrological and detrital zircon geochronological data from the Paleocene-Eocene Nangqian-Xialaxiu and Gongjue basins. Sandstone detrital modes and zircon morphology suggest that the samples collected in these two basins were sourced from recycled orogen. Detrital zircon geochronology indicates that sediments in the Nangqian-Xialaxiu Basin are characterized by two distinct age populations at 220-280 Ma and 405-445 Ma. In contrast, three predominant age populations of 207-256 Ma, 423-445 Ma, and 1851-1868 Ma, and two subordinate age populations of similar to 50 Ma and similar to 2500 Ma, are recognized in the Gongjue Basin. Comparison with detrital zircon ages from the surrounding terranes suggests that sediments in the Nangqian-Xialaxiu Basin come from the neighboring thrust belts, whereas sediments from the Gongjue Basin are predominantly derived from the distant Songpan-Ganzi Terrane with minor contribution from the surrounding areas. A three-stage Cenozoic evolution of the eastern Tibetan Plateau is proposed. During the Paleocene, the Nangqian-Xialaxiu Basin appeared as a set of small intermontane sub-basins and received plentiful sediments from the neighboring mountain belts; during the Eocene, the Gongjue Basin kept a relatively low altitude and was a depression at the edge of a proto-Plateau; since the Oligocene, the Tibetan Plateau further uplifted and the marginal Gongjue Basin was involved in the Tibetan interior orogeny, indicating the eastward propagation of the Tibetan Plateau.}, language = {en} } @article{ZhangNajmanMeietal.2019, author = {Zhang, Peng and Najman, Yani and Mei, Lianfu and Millar, Ian and Sobel, Edward and Carter, Andrew and Barfod, Dan and Dhuime, Bruno and Garzanti, Eduardo and Govin, Gwladys and Vezzoli, Giovanni and Hu, Xiaolin}, title = {Palaeodrainage evolution of the large rivers of East Asia, and Himalayan-Tibet tectonics}, series = {Earth science reviews}, volume = {192}, journal = {Earth science reviews}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-8252}, doi = {10.1016/j.earscirev.2019.02.003}, pages = {601 -- 630}, year = {2019}, abstract = {A number of sedimentary provenance studies have been undertaken in order to determine whether the palaeo-Red River was once a river of continental proportions into which the upper reaches of the Yangtze, Salween, Mekong, Irrawaddy, and Yarlung drained. We have assessed the evidence that the Yarlung originally flowed into the palaeo-Red river, and then sequentially into the Irrawaddy and Brahmaputra, connecting to the latter first via the Lohit and then the Siang. For this river system, we have integrated our new data from the Paleogene-Recent Irrawaddy drainage basin (detrital zircon U-Pb with Hf and fission track, rutile U-Pb, mica Ar-Ar, bulk rock Sr-Nd, and petrography) with previously published data, to produce a palaeodrainage model that is consistent with all datasets. In our model, the Yarlung never flowed into the Irrawaddy drainage: during the Paleogene, the Yarlung suture zone was an internally drained basin, and from Neogene times onwards the Yarlung drained into the Brahmaputra in the Bengal Basin. The Central Myanmar Basin, through which the Irrawaddy River flows today, received predominantly locally-derived detritus until the Middle Eocene, the Irrawaddy initiated as a through-going river draining the Mogok Metamorphic Belt and Bomi-Chayu granites to the north sometime in the Late Eocene to Early Oligocene, and the river was dominated by a stable MMB-dominated drainage throughout the Neogene to present day. Existing evidence does not support any connection between the Yarlung and the Red River in the past, but there is a paucity of suitable palaeo-Red River deposits with which to make a robust comparison. We argue that this limitation also precludes a robust assessment of a palaeo-connection between the Yangtze/ Salween/Mekong and the Red River; it is difficult to unequivocally interpret the recorded provenance changes as the result of specific drainage reorganisations. We highlight the palaeo-Red River deposits of the Hanoi Basin as a potential location for future research focus in view of the near-complete Cenozoic record of palaeo-Red River deposits at this location. A majority of previous studies consider that if a major continental-scale drainage ever existed at all, it fragmented early in the Cenozoic. Such a viewpoint would agree with the growing body of evidence from palaeoaltitude studies that large parts of SE Tibet were uplifted by this period. This then leads towards the intriguing question as to the mechanisms which caused the major period of river incision in the Miocene in this region.}, language = {en} } @article{ZapataSobelPapaetal.2019, author = {Zapata, Sebastian Henao and Sobel, Edward and Papa, Cecilia Eugenia del and Muruaga, Claudia and Zhou, R.}, title = {Miocene fragmentation of the Central Andean foreland basins between 26 and 28 degrees S}, series = {Journal of South American earth sciences}, volume = {94}, journal = {Journal of South American earth sciences}, publisher = {Elsevier}, address = {Oxford}, issn = {0895-9811}, doi = {10.1016/j.jsames.2019.102238}, pages = {18}, year = {2019}, abstract = {We present new U-Pb LA-ICP-MS data from the Central Andean foreland basins combined with new and published stratigraphic information in order to reconstruct the Miocene fragmentation of the Andean foreland between 26 and 28 degrees S. The disruption of this foreland basin and the subsequent development of elevated intermountain basins have been the focus of several studies. However, the absence of temporal constraints in the Miocene to Pliocene sedimentary record of the low elevation Choromoro and Tucuman foreland basins has presented an obstacle for precise paleogeographic reconstructions. We describe 11 discontinuous stratigraphic sections and use the U-Pb LA-ICP-MS method to date 10 pyroclastic-bearing sediments in order to reconstruct the stratigraphic evolution of the Choromoro and Tucuman basins. We combine our results with published strati graphic and thermochronologic data from adjacent basins to present a refined Miocene paleogeographic model. In a first stage, a continuous Early Miocene foreland lacustrine basin developed, filling up the preexisting Paleogene topography. The second stage is characterized by basin unroofing around similar to 12 Ma; the easily eroded sedimentary cover was removed, leading to the uplift of the underlying basement rocks and the segmentation of the lacustrine system. In the third stage, relief increase took place after similar to 6 Ma due to the low erodibility of the basement blocks; as a result, stable fluvial systems developed. Progressive relief development caused pronounced unconformities in the basins and the development of proximal fluvial-gravitational depositional systems after 3 Ma. This model emphasizes on the relations between tectonics, climate, and erodibility, and their control on the evolution of the depositional systems and relief.}, language = {en} } @article{ZaccarelliBindiStrolloetal.2019, author = {Zaccarelli, Riccardo and Bindi, Dino and Strollo, Angelo and Quinteros, Javier and Cotton, Fabrice}, title = {Stream2segment: An Open-Source Tool for Downloading, Processing, and Visualizing Massive Event-Based Seismic Waveform Datasets}, series = {Seismological research letters}, volume = {90}, journal = {Seismological research letters}, number = {5}, publisher = {Seismological Society of America}, address = {Albany}, issn = {0895-0695}, doi = {10.1785/0220180314}, pages = {2028 -- 2038}, year = {2019}, abstract = {The task of downloading comprehensive datasets of event-based seismic waveforms has been made easier through the development of standardized webservices but is still highly nontrivial because the likelihood of temporary network failures or subtle data errors naturally increases when the amount of requested data is in the order of millions of relatively short segments. This is even more challenging because the typical workflow is not restricted to a single massive download but consists of fetching all possible available input data (e.g., with several repeated download executions) for a processing stage producing any desired user-defined output. Here, we present stream2segment, a highly customizable Python 2+3 package helping the user in the entire workflow of downloading, inspecting, and processing event-based seismic data by means of a relational database management system as archiving storage, which has clear performance and usability advantages, and an integrated processing subroutine requiring a configuration file and a single Python function to produce user-defined output. Stream2segment can also produce diagnostic maps or user-defined plots, which, unlike existing tools, do not require external software dependencies and are not static images but instead are interactive browser-based applications ideally suited for data inspection or annotation tasks and subsequent training of classifiers in foreseen supervised machine-learning applications. Stream2segment has already been used as a data quality tool for datasets within the European Integrated Data Archive and to create a weak-motion database (in the form of a so-called flat file) for the stable continental region of Europe in the context of the European Ground Shaking Intensity Model service, in turn an important building block for seismic hazard studies.}, language = {en} } @article{YuanBraunGueritetal.2019, author = {Yuan, Xiaoping and Braun, Jean and Guerit, Laure and Simon, Brendan and Bovy, Beno{\^i}t and Rouby, Delphine and Robin, C{\´e}cile and Jiao, R.}, title = {Linking continental erosion to marine sediment transport and deposition: A new implicit and O(N) method for inverse analysis}, series = {Earth \& planetary science letters}, volume = {524}, journal = {Earth \& planetary science letters}, publisher = {Elsevier}, address = {Amsterdam}, issn = {0012-821X}, doi = {10.1016/j.epsl.2019.115728}, pages = {15}, year = {2019}, abstract = {The marine sedimentary record contains unique information about the history of erosion, uplift and climate of the adjacent continent. Inverting this record has been the purpose of many numerical studies. However, limited attention has been given to linking continental erosion to marine sediment transport and deposition in large-scale surface process evolution models. Here we present a new numerical method for marine sediment transport and deposition that is directly coupled to a landscape evolution algorithm solving for the continental fluvial and hillslope erosion equations using implicit and O(N) algorithms. The new method takes into account the sorting of grain sizes (e.g., silt and sand) in the marine domain using a non-linear multiple grain-size diffusion equation and assumes that the sediment flux exported from the continental domain is proportional to the bathymetric slope. Specific transport coefficients and compaction factors are assumed for the two different grain sizes to simulate the stratigraphic architecture. The resulting set of equations is solved using an efficient (O(N) and implicit) algorithm. It can thus be used to invert stratigraphic geometries using a Bayesian approach that requires a large number of simulations. This new method is used to invert the sedimentary geometry of a natural example, the Ogooue Delta (Gabon), over the last similar to 5 Myr. The objective is to unravel the set of erosional histories of the adjacent continental domain compatible with the observed geometry of the offshore delta. For this, we use a Bayesian inversion scheme in which the misfit function is constructed by comparing four geometrical parameters between the natural and the simulated delta: the volume of sediments stored in the delta, the surface slope, the initial and the final shelf lengths. We find that the best-fit values of the transport coefficients for silt in the marine domain are in the range of 300 - 500 m(2)/yr, in agreement with previous studies on offshore diffusion. We also show that, in order to fit the sedimentary geometry, erosion rate on the continental domain must have increased by a factor of 6 to 8 since 5.3 Ma. (C) 2019 Elsevier B.V. All rights reserved.}, language = {en} }