Combined Effects of Tectonic and Landslide-Generated Tsunami Runup at Seward, Alaska During the M W 9.2 1964 Earthquake

We apply a recently developed and validated numerical model of tsunami propagation and runup to study the inundation of Resurrection Bay and the town of Seward by the 1964 Alaska tsunami. Seward was hit by both tectonic and landslide-generated tsunami waves during the $M_{\rm W}$ 9.2 1964 megathrust earthquake. The earthquake triggered a series of submarine mass failures around the fjord, which resulted in landsliding of part of the coastline into the water, along with the loss of the port facilities. These submarine mass failures generated local waves in the bay within 5?min of the beginning of strong ground motion. Recent studies estimate the total volume of underwater slide material that moved in Resurrection Bay to be about 211?million m³ (Haeussler et?al. in Submarine mass movements and their consequences, pp 269?C278, 2007). The first tectonic tsunami wave arrived in Resurrection Bay about 30?min after the main shock and was about the same height as the local landslide-generated waves. Our previous numerical study, which focused only on the local landslide-generated waves in Resurrection Bay, demonstrated that they were produced by a number of different slope failures, and estimated relative contributions of different submarine slide complexes into tsunami amplitudes (Suleimani et?al. in Pure Appl Geophys 166:131?C152, 2009). This work extends the previous study by calculating tsunami inundation in Resurrection Bay caused by the combined impact of landslide-generated waves and the tectonic tsunami, and comparing the composite inundation area with observations. To simulate landslide tsunami runup in Seward, we use a viscous slide model of Jiang and LeBlond (J Phys Oceanogr 24(3):559?C572, 1994) coupled with nonlinear shallow water equations. The input data set includes a high resolution multibeam bathymetry and LIDAR topography grid of Resurrection Bay, and an initial thickness of slide material based on pre- and post-earthquake bathymetry difference maps. For simulation of tectonic tsunami runup, we derive the 1964 coseismic deformations from detailed slip distribution in the rupture area, and use them as an initial condition for propagation of the tectonic tsunami. The numerical model employs nonlinear shallow water equations formulated for depth-averaged water fluxes, and calculates a temporal position of the shoreline using a free-surface moving boundary algorithm. We find that the calculated tsunami runup in Seward caused first by local submarine landslide-generated waves, and later by a tectonic tsunami, is in good agreement with observations of the inundation zone. The analysis of inundation caused by two different tsunami sources improves our understanding of their relative contributions, and supports tsunami risk mitigation in south-central Alaska. The record of the 1964 earthquake, tsunami, and submarine landslides, combined with the high-resolution topography and bathymetry of Resurrection Bay make it an ideal location for studying tectonic tsunamis in coastal regions susceptible to underwater landslides.     

Ries crater and suevite revisited—Observations and modeling Part I: Observations

We report results of an interdisciplinary project devoted to the 26 km‐diameter Ries crater and to the genesis of suevite. Recent laboratory analyses of “crater suevite” occurring within the central crater basin and of “outer suevite” on top of the continuous ejecta blanket, as well as data accumulated during the past 50 years, are interpreted within the boundary conditions imposed by a comprehensive new effort to model the crater formation and its ejecta deposits by computer code calculations (Artemieva et al. 2013). The properties of suevite are considered on all scales from megascopic to submicroscopic in the context of its geological setting. In a new approach, we reconstruct the minimum/maximum volumes of all allochthonous impact formations (108/116 km³), of suevite (14/22 km³), and the total volume of impact melt (4.9/8.0 km³) produced by the Ries impact event prior to erosion. These volumes are reasonably compatible with corresponding values obtained by numerical modeling. Taking all data on modal composition, texture, chemistry, and shock metamorphism of suevite, and the results of modeling into account, we arrive at a new empirical model implying five main consecutive phases of crater formation and ejecta emplacement. Numerical modeling indicates that only a very small fraction of suevite can be derived from the “primary ejecta plume,” which is possibly represented by the fine‐grained basal layer of outer suevite. The main mass of suevite was deposited from a “secondary plume” induced by an explosive reaction (“fuel‐coolant interaction”) of impact melt with water and volatile‐rich sedimentary rocks within a clast‐laden temporary melt pool. Both melt pool and plume appear to be heterogeneous in space and time. Outer suevite appears to be derived from an early formed, melt‐rich and clast‐poor plume region rich in strongly shocked components (melt ? clasts) and originating from an upper, more marginal zone of the melt pool. Crater suevite is obviously deposited from later formed, clast‐rich and melt‐poor plumes dominated by unshocked and weakly shocked clasts and derived from a deeper, central zone of the melt pool. Genetically, we distinguish between “primary suevite” which includes dike suevite, the lower sublayer of crater suevite, and possibly a basal layer of outer suevite, and “secondary suevite” represented by the massive upper sublayer of crater suevite and the main mass of outer suevite.     

Geochemical studies of the SUBO 18 (Enkingen) drill core and other impact breccias from the Ries crater,Germany

Suevite and melt breccia compositions in the boreholes Enkingen and Polsingen are compared with compositions of suevites from other Ries boreholes and surface locations and discussed in terms of implications for impact breccia genesis. No significant differences in average chemical compositions for the various drill cores or surface samples are noted. Compositions of suevite and melt breccia from southern and northeastern sectors of the Ries crater do not significantly differ. This is in stark contrast to the published variations between within‐crater and out‐of‐crater suevites from northern and southern sectors of the Bosumtwi impact structure, Ghana. Locally occurring alteration overprint on drill cores—especially strong on the carbonate‐impregnated suevite specimens of the Enkingen borehole—does affect the average compositions. Overall, the composition of the analyzed impact breccias from Ries are characterized by very little macroscopically or microscopically recognized sediment‐clast component; the clast populations of suevite and impact melt breccia are dominated consistently by granitic and intermediate granitoid components. The Polsingen breccia is significantly enriched in a dioritic clast component. Overall, chemical compositions are of intermediate composition as well, with dioritic‐granodioritic silica contents, and relatively small contributions from mafic target components. Selected suevite samples from the Enkingen core have elevated Ni, Co, Cr, and Ir contents compared with previously analyzed suevites from the Ries crater, which suggest a small meteoritic component. Platinum‐group element (PGE) concentrations for some of the enriched samples indicate somewhat elevated concentrations and near‐chondritic ratios of the most immobile PGE, consistent with an extraterrestrial contribution of 0.1–0.2% chondrite‐equivalent.     

Geomagnetic control of polar mesosphere summer echoes

Using observations with the ALOMAR SOUSY radar near Andenes (69.3°N, 16.0°E) from 1994 until 1997 polar mesosphere summer echoes (PMSE) have been investigated in dependence on geomagnetic K indices derived at the Auroral Observatory Tromsø (69.66°N, 18.94°E). During night-time and morning hours a significant correlation between the signal-to-noise ratio (SNR) of the radar results and the geomagnetic K indices could be detected with a maximum correlation near midnight. The correlation becomes markedly smaller in the afternoon and early evening hours with a minimum near 17 UT. This diurnal variation is in reasonable agreement with riometer absorption at Ivalo (68.55°N, 27.28°E) and can be explained by the diurnal variation of ionization due to precipitating high energetic particles. Therefore, a part of the diurnal PMSE variation is caused by this particle precipitation. The variability of the solar EUV variation, however, has no significant influence on the PMSE during the observation period.     

安徽巢湖地区下三叠统综合层序

印度阶一奥伦尼克阶界线的全球层型候选剖面之一位于安徽巢湖地区;巢湖地区的下三叠统也是国际上同期地层中生物地层序列最有代表性、多重地层学手段应用齐备、研究效果最好的地层序列之一。根据巢湖地区3条代表性下三叠统剖面的岩石地层、生物地层和碳同位素地层的最新研究成果归纳出本区早三叠世综合地层序列,作为区域地层对比研究的标准．巢湖的下三叠统明确包含8个牙形石带和6个菊石带。它们具有区域甚至全球对比意义;巢湖地区早三叠世碳同位素δ^13Ccarb的演变呈现2个显著的漂移周期,这种有特色的碳同位素漂移,不仅具有地层学价值,而且可能对于三叠纪初的生物复苏和生态系演变具有指导意义．在综合地层序列基础上,将3条剖面的古地磁学研究成果链接,形成了巢湖地区完整的早三叠世磁性地层序列,包括5个主要的正向极性带和5个反向极性带,这也是目前在本区乃至华南获得的最完整的早三叠世磁极性序列之一．这些为该区域乃至全球相关地层研究树立了一个基本格架。     

Predicting Soil-Water and Soil-Air Transport Properties and Their Effects on Soil-Vapor Extraction Efficiency

Accurate prediction of water and air Iran sport parameters in variably saturated soil is necessary for modeling of soil-vapor extraction (SVE) at soil sites contaminated with volatile organic chemicals (VOCs). An expression for predicting saturated water permeability (k_l,s) in undisturbed soils from the soil total porosity and the field capacity soil-water content was developed by fitting a tortuous-tube fluid flow model to measured water permeability and gas diffusivity data. The new k_l,s expression gave accurate predictions when tested against independent k_l,s data. The k_l,s expression was implemented in the Campbell relative water permeability model to yield a predictive model for water permeability in variably saturated, undisturbed soil. The water permeability model, together with recently developed predictive equations for gas permeability and gas diffusivity, was used in a two-dimensional numerical SVE model that also included non-equilibrium mass transfer of VOC from a separate phase (nonaqueous phase liquid [NAPL]) to the air phase. SVE: calculations showed that gas permeability is likely the most important factor controlling VOC migration and vapor extraction efficiency. Water permeability and gas diffusivity effects became significant at water contents near and above field capacity. The NAPL-air mass transfer coefficient also had large impacts on simulated vapor extraction efficiency. The calculations suggest that realistic SVE models need to include predictive expressions for both conveciive, diffusive. and phase-partitioning processes in natural, undisturbed soils.     

西准噶尔库什库都克石英闪长岩地球化学、年代学及地质意义

库什库都克石英闪长岩位于新疆西准噶尔的达尔布特晚古生代中酸性岩浆岩带中。通过LA-ICP-MS锆石U-Pb定 年,结合锆石阴极发光图像和U、Th元素特征,获得石英闪长岩的年龄为305.3±2.8 Ma（n=19,MSWD=0.17）,时代属晚石 炭世晚期（C2）。石英闪长岩呈半自形不等粒结构,块状构造,矿物成分主要为斜长石、石英、角闪石和黑云母等,局部发 育含柱状辉石的石英辉长闪长岩。岩石的SiO2含量为59.38%~63.37%,全碱变化于6.55%~7.51%;里特曼指数( σ )变化于 2.62~2.88,铝饱和指数(A/CNK值)介于0.83~0.90,为准铝质;岩石轻稀土富集且轻重稀土元素分馏明显,（La/Yb）N 为2.99~ 3.66,负Eu异常（δEu为0.64~0.79）;岩石相对富集Ba、Nd、Zr、K和LREE等元素,而相对亏损Rb、Nb、Th、Sr、P、Ti 和HREE等元素,属高钾钙碱性准铝质I型花岗岩类岩石。石英闪长岩地球化学特征表明其为上地壳变质玄武岩熔融的产 物,在岩浆演化过程中发生了以斜长石、辉石、金红石和磷灰石为主的分离结晶作用,而角闪石基本未发生分离结晶。区 域构造环境演化分析认为西准噶尔不仅在晚石炭世早期存在与俯冲作用相关的岛弧花岗岩,在晚石炭世晚期-早二叠世早 期也存在持续俯冲作用,在此环境下形成了与俯冲作用相关的岛弧型库什库都克石英闪长岩。     

Importance of geological information for assessing drain flow in a Danish till landscape

There has, in recent years, been an increasing interest in developing nutrient load mitigation measures focussing on tile drains. To plan the location of such tile drain measures, it is important to know where in the landscape drain flow is generated and to understand the key factors governing drain flow dynamics. In the present study, we test two approaches to assess spatial patterns in drain flow generation and thereby assess the importance of including geological information. The approaches are the widely used topographical wetness index (TWI), based solely on elevation data, and hydrological models that include the subsurface geology. We set‐up an ensemble of 20 hydrological models based on 20 stochastically generated geological models to predict drain flow dynamics in the clay till Norsminde catchment in Denmark and test the results against TWI. We find that the hydrological models predict observed daily drain flow reasonably well. High drain flow volumes were found in stream valleys and in wetlands and lower drain flow volumes in the more hilly parts of the catchment. In spite of the apparent connection to the landscape, there was no statistically significant correlation between TWI and drain flow at grid scale (100 × 100 m). TWI was therefore not found to be a sufficient index on its own to assess where drain flow is generated, especially in the highlands of the catchment. The geology below 3 m was found to have a large impact on the drain flow, and correlations between sand percentage in the subsurface geology and drain flow volume were found to be statistically significant. Geological uncertainty therefore give rise to uncertainty on simulated drain flow, and this uncertainty was found to be high at the model grid scale but decreasing with increasing scale.     

Asymptotic performance of a pulse-to-pulse incoherent Doppler sonar in an oceanic environment

Active Doppler sonar systems are gaining popularity as a method to profile remotely the oceanic velocity field. Estimation of the variance with which an ideal single.beam sonar system measures relative Doppler motion is the subject of this paper. Theoretical formulations, simulated reverberation data, and sea data are considered to determine this variance. The sea data considered in this paper were selected carefully to avoid performance degrading physical/instrumental phenomena and bias corrected to remove further instrumental effects. Asymptotic performance, determined in this research, is realized when instrumental and physical phenomena do not measurably affect experimental performance. Oceanic and simulated data sets are processed with several velocity-estimation algorithms that are implemented in existing sonar systems. These algorithms are compared against each other and the velocity-estimation standard-deviation results tabulated. The algorithms are shown to possess similar performance characteristics when simulated, and sea data are processed in spite of the fact that some of these algorithms are suboptimal from a theoretical viewpoint. The results of this study show that the asymptotic performance of the Doppler system examined is approximately 25 percent above a theoretical lower bound.     

Oceanic incoherent Doppler sonar spectral analysis by conventional and finite-parameter modeling methods

The spectral structure of oceanic incoherent Doppler sonar data is analyzed with Fourier and two finite-parameter modeling methods (Spectral Moment Estimation (SME) and AutoRegressive (AR)). Fourier spectral analysis, applied to long-range segments of echo return data, validates a theoretical point scattering model. It is demonstrated that Fourier spectral techniques, applied to data from incoherent sonar systems, cannot be used to probe oceanic spectral and velocity fields on the scales that frequently are important in oceanographic applications. The finite-parameter methods model the spectral structure on closely spaced range intervals. The SME method calculates the first moment of the Doppler spectrum at each range interval and the AR method provides spectral estimates at these same intervals. Trade-offs in estimating scatterer radial velocity using the various algorithms are considered. It is shown that a frequency-integrated AR technique has velocity estimation performance comparable to the SME method. In addition, the AR technique reveals that several regions of data possess asymmetric spectral structure. The implications of this spectral structure for oceanic velocity estimates are discussed.