Energy flow in shallow depth based on vertical array records during recent strong earthquakes (supplement)

Based on subsurface wave energy flows calculated in the previous paper for strong motion records at 30 vertical array sites by assuming vertical propagation of SH waves, upward energy ratios between layers have been correlated to corresponding impedance ratios. This indicates that the upward energy tends to decrease almost in proportion to the power of 0.70 of the impedance ratio from deep-seated stiff rock (equivalent to seismological bedrock) to soft soil near ground surface. The upward energies at the seismological bedrock have been estimated from the observed energy at the deepest level for all vertical array sites by using this relationship and plotted versus hypocentral distances.     

The 1957 great Aleutian earthquake

The 9 March 1957 Aleutian earthquake has been estimated as the third largest earthquake this century and has the longest aftershock zone of any earthquake ever recorded—1200 km. However, due to a lack of high-quality seismic data, the actual source parameters for this earthquake have been poorly determined. We have examined all the available waveform data to determine the seismic moment, rupture area, and slip distribution. These data include body, surface and tsunami waves. Using body waves, we have estimated the duration of significant moment release as 4 min. From surface wave analysis, we have determined that significant moment release occurred only in the western half of the aftershock zone and that the best estimate for the seismic moment is 50–100×10²⁰ Nm. Using the tsunami waveforms, we estimated the source area of the 1957 tsunami by backward propagation. The tsunami source area is smaller than the aftershock zone and is about 850 km long. This does not include the Unalaska Island area in the eastern end of the aftershock zone, making this area a possible seismic gap and a possible site of a future large or great earthquake. We also inverted the tsunami waveforms for the slip distribution. Slip on the 1957 rupture zone was highest in the western half near the epicenter. Little slip occurred in the eastern half. The moment is estimated as 88×10²⁰ Nm, orM _w =8.6, making it the seventh largest earthquake during the period 1900 to 1993. We also compare the 1957 earthquake to the 1986 Andreanof Islands earthquake, which occurred within a segment of the 1957 rupture area. The 1986 earthquake represents a rerupturing of the major 1957 asperity.     

多层球体对电磁场散射响应的求解方法

为建立多层球体对电磁场散射响应的求解方法,将球外入射电磁场分为TM和TE极化波两部分,通过求解Debye位方程可得到多层球体各层中两种极化场的通解形式,利用边界条件建立了TM波（或TE波）球面波阻抗（或导纳）的递推关系,借助于阻抗（或寻纳）导出了球体对TM波（或TE波）所产生散射场的一般表达式．文中具体求解了球体对平面波径向和横向偶极子场的散射响应．     

Influence of wave activity on the composition of the polar stratosphere

The planetary wave impact on the polar vortex stability, polar stratosphere temperature, and content of ozone and other gases was simulated with the global chemical–climatic model of the lower and middle atmosphere. It was found that the planetary waves propagating from the troposphere into the stratosphere differently affect the gas content of the Arctic and Antarctic stratosphere. In the Arctic region, the degree of wave activity critically affects the polar vortex formation, the appearance of polar stratospheric clouds, the halogen activation on their surface, and ozone anomaly formation. Ozone anomalies in the Arctic region as a rule are not formed at high wave activity and can be registered at low activity. In the Antarctic Regions, wave activity affects the stability of polar vortex and the depth of ozone holes, which are formed at almost any wave activity, and the minimal ozone values depend on the strong or weak wave activity that is registered in specific years.     

Gravity waves in the mesopause region observed by meteor radar, 2: Climatologies of gravity waves in the Antarctic and Arctic

A simple technique has been used to investigate gravity waves in the Antarctic and Arctic mesosphere. Data from two meteor radars are used, one at Rothera in the Antarctic (68°S) and one at Esrange in the Arctic (68°N). Observations have been made from 2005 to 2008 in the Antarctic and from 2000 to 2008 in the Arctic. The results reveal a seasonal behaviour with a semi-annual cycle. Over both sites wave activity maximises at the solstices and minimises at the equinoxes. Monthly-mean gravity-wave activity increases with height in all seasons except in summer when gravity-wave variances show little or no increase with height at heights below about 90 km. Above about 90 km in summer there is a rapid increase in wave activity. We suggest that this summer-time behaviour is a consequence of rapid wave growth as gravity waves ascend into the more stable air of the lower thermosphere. Differences between the Arctic and Antarctic gravity-wave field are also evident. In particular, a higher level of gravity-wave activity is observed in the Antarctic spring compared to the Arctic spring. We suggest that these inter-hemispheric differences may be the result of differences in the generation of gravity waves in the polar troposphere and stratosphere.     

Detection of crustal inhomogeneity in the Nojima Fault zone using earth current noise

Abstract A resistivity survey method using artificial telluric noise was examined and applied to a field of a fault zone. The electric earth current was measured at 50 sites in the Nojima Fault zone, which is in the northwestern part of Awaji Island, southwestern Japan. The dominant component of the observed electric field is supposed to be leakage currents from DC electric railways running outside the island. Amplitude and polarization of the stray current were systematically investigated and were revealed to represent the subsurface electrical structure of the study area. Some features on the fault zone's electrical structure have been pointed out, including: (i) an electrical boundary that corresponds to a geological one between granite (resistive) and sediments (conductive); and (ii) a low resistivity spot on the surface rupture of the earthquake fault. The structure estimated in the present study is both qualitatively and quantitatively consistent with previous resistivity surveys done using other methods pursued in the same area. It shows the validity of the 'stray current method' as one that is easy and uses low-cost resistivity exploration tools in a region where the effect of artificial noise caused mainly by leakage currents from electrical railways cannot be ignored.     

Annual variations in the electromagnetic field of the earth and electrical conductivity of the geological environment

Synchronous annual variations in the geoelectric and geomagnetic field are studied on the basis of long-term electromagnetic monitoring. It is shown that the annual geoelectric variations have intraterrestrial origin and are not related to the annual geomagnetic variations. Temporal variations in the magnetotelluric impedance and magnetic tipper, which characterize the electrical conductivity of the geological environment, are analyzed. It is established that annual variations in the magnetotelluric impedance mainly describe the variations in the electrical conductivity of surface crustal layers and are less sensitive to the deep electrical conductivity of the Earth. The annual variations in the imaginary magnetic tipper at the periods of 1000–3000 s probably reflect the changes in conductivity of a deep transversal low-resistive zone (the fault). It is suggested that annual variations in the geoelectrical and geomagnetic fields, as well as in the electrical conductivity of the geological environment, arise as a response to the changes in the geodynamical processes caused by the revolution of the Earth around the Sun.

     

Generalized polarization relationships for acoustic gravity waves in the nonisothermal atmosphere with the wind

The nonlinear equation of the first order of the Riccati type has been obtained for the wave impedance of acoustic gravity waves in the nonisothermal atmosphere. The vertically nonuniform horizontal wind and the effect of viscosity on the horizontal components of the velocity field have been taken into account in the calculations. The boundary-value problem for the Riccati equation is defined by the boundary emission condition at high altitudes. Upon finding the wave impedance along with the generalized polarization relationship, all remaining disturbances of the atmospheric parameters related to acoustic gravity waves are found with the help of a simple integration. The results of using a developed formalism are illustrated by the numerical computation of acoustic gravity wave fields in the atmosphere with real vertical profiles of temperature and horizontal field velocity.     

'Green' enclosure of ocean space--déjà vu?

The first wave of unilaterally declared environmental protection zones (starting with Canada's Arctic Waters Pollution Prevention Act of 1970) may have been halted by the 1982 UN Convention on the Law of the Sea. Yet, a second wave now seems to be gaining momentum: the French 'ecological protection zone in the Mediterranean' of April 2003 (promptly emulated by Croatia, Slovenia, and Italy), the British 'environment (protection and preservation) zone' around the Chagos Archipelago of September 2003, and the US 'marine national monument' around the Northwestern Hawaiian Islands of June 2006, all claiming new environmental regulatory powers in sea areas way beyond territorial waters. This viewpoint article questions the unilateral legitimation of those claims, and postulates a trusteeship ethos and fiduciary accountability for ocean governance in these areas.     

Time-dependence of electromagnetic transfer functions and their association with tectonic activity

Laboratory experiments and theoretical considerations have suggested that anomalous dilatant regions can develop in the earth's crust during the period of strain accumulation prior to an earthquake. For moderate and major earthquakes such anomalous regions could be tens or even hundreds of kilometers in extent and should be detectable at the surface with appropriate survey or sounding techniques. Since electrical resistivity is one of the rock properties likely to be strongly modified in a dilatant zone, magnetotelluric impedance and geomagnetic transfer functions might be expected to show time-dependent precursory effects if monitored over a period of time above the focal region of an impending earthquake. Such experiments have been conducted in Japan and in other parts of the world and several examples of resitivity changes in the crust prior to earthquake occurrence have been reported. These results and their association with local seismicity are reviewed in this paper. The available evidence indicates that transfer functions and impedance can display significant time-dependent response to changing crustal conditions in some regions. However, the correspondence between these effects and earthquake occurrence is usually not very clear.Contribution from the Earth Physics Branch No. 858.     

Effects of weathering on site response

Nearly all studies of site response under seismic excitation employ idealized site models in which uniform layers of homogeneous soils are separated by perfect interfaces. The impedance mismatch at each interface gives rise to wave reflections which strongly affect the calculated free surface response. In reality, these perfect interfaces may often not exist. Instead, a zone in which materials from both above and below are mingled may form a blurred interface between layers. This article examines the effects of such weathered zones. The zone is modelled by a linearly varying shear modulus so that no abrupt impedance change occurs. Significant differences in response of weathered and unweathered profiles are found.     

MT方法中Basokur视电阻率定义的研究

电法勘探中不同的视电阻率定义可以带来不同的反演效果,好的视电阻率定义可以提高反演分辨率。传统的Cagniar视电阻率(1953)只是利用了波阻抗的模值,而Basokur 1994年提出的视电阻率利用波阻抗的实部和虚部。本文做了4个不同模型的Basokur视电阻率响应及反演,结果显示新的视电阻率对薄层有很好的反应,对常规模型也有不错的表现。     

Groundwater modelling to quantify the effect of water harvesting structures in Wagarwadi watershed,Parbhani district,Maharashtra, India

Marathwada Agricultural University, Pharbani, has developed about 560 hectares of Wagarwadi watershed in Pharbani district since 1987. Groundwater monitoring on 16 observations wells at weekly intervals commenced in January 1992, and rainfall and pan evaporation has been measured daily at a hydrometeorological station situated in the nearby university campus. Aquifer parameters, namely, transmissivity and specific yield, have been estimated by carrying out a pumping test on a large diameter well. Groundwater recharge resulting from rainfall has been estimated using a water balance model of the soil moisture zone considering soil zone thickness and crops grown. The SCS (Soil Conservation Service) curve number method was used for surface runoff estimation. The groundwater flow model has been constructed using the nested squares, finite difference method. Nested square meshes of sizes 160 m×160 m and 80 m×80 m have been used and the steady-state condition of aquifer system was simulated in the model assuming the June 1992 water level configuration under equilibrium conditions. The model has been calibrated for transient conditions incorporating additional seepage from the water harvesting structures and their contribution to the groundwater regime has been assessed. © 1998 John Wiley & Sons, Ltd.     

Sea ice thickness measurement and its underside morphology analysis using radar penetration in the Arctic Ocean

Sea ice, as an important component of the Arctic climate system, has drawn significant sci-entific interest. Sea ice thickness and its morphology have dramatic impacts on ocean-atmos- phere-ice interactions[1—4], which directly affect the exchange proces…     

Experimental estimation of velocities and anisotropy of a series of Swedish crystalline rocks and ores

To provide a guide for future deep (<1.5 km) seismic mineral exploration and to better understand the nature of reflections imaged by surface reflection seismic data in two mining camps and a carbonatite complex of Sweden, more than 50 rock and ore samples were collected and measured for their seismic velocities. The samples are geographically from the northern and central parts of Sweden, ranging from metallic ore deposits, meta‐volcanic and meta‐intrusive rocks to deformed and metamorphosed rocks. First, ultrasonic measurements of P‐ and S‐wave velocities at both atmospheric and elevated pressures, using 0.5 MHz P‐ and S‐wave transducers were conducted. The ultrasonic measurements suggest that most of the measured velocities show positive correlation with the density of the samples with an exception of a massive sulphide ore sample that shows significant low P‐ and S‐wave velocities. The low P‐ and S‐wave velocities are attributed to the mineral texture of the sample and partly lower pyrite content in comparison with a similar type sample obtained from Norway, which shows significantly higher P‐ and S‐wave velocities. Later, an iron ore sample from the central part of Sweden was measured using a low‐frequency (0.1–50 Hz) apparatus to provide comparison with the ultrasonic velocity measurements. The low‐frequency measurements indicate that the iron ore sample has minimal dispersion and attenuation. The iron ore sample shows the highest acoustic impedance among our samples suggesting that these deposits are favourable targets for seismic methods. This is further demonstrated by a real seismic section acquired over an iron ore mine in the central part of Sweden. Finally, a laser‐interferometer device was used to analyse elastic anisotropy of five rock samples taken from a major deformation zone in order to provide insights into the nature of reflections observed from the deformation zone. Up to 10% velocity‐anisotropy is estimated and demonstrated to be present for the samples taken from the deformation zone using the laser‐interferometery measurements. However, the origin of the reflections from the major deformation zone is attributed to a combination of anisotropy and amphibolite lenses within the deformation zone.     

低速带的波场数值模拟及影响因素初步分析

低速带具有“自由表面、低速度、高吸收”的特征,能引起反射波能量减弱,产生很强的频散,造成近地表波场畸变,严重影响地震勘探的精度。而地震采集过程又不可避免遇到以下低速带问题:①低速带存在与否会对地震波传播造成什么影响;②近地表变化大,怎么才能得到精确的低速带参数;③存在低速带时,激发井深、最佳的药量(与子波频率有关)、最佳的耦合(与低速带速度相关)等采集参数会对采集资料的品质产生什么影响。本文通过数值模拟,集中对低速带与地震波激发的井深、速度、子波主频等核心采集参数的关系进行分析研究,同时在低速带背景下对数值模拟的网格、边界等问题也进行了讨论。模拟得到的波场特征及低速带参数影响的机理分析,对做好野外地震波激发,获取高品质的地震资料,具有理论意义和指导作用。      

Study of Regional Surface Waves and Frequency-dependent Ms:mb Discrimination in the European Arctic

—?Accurate discrimination of seismic events with a regional network requires detailed knowledge of the propagation characteristics of seismic waves in the region. At present, such propagation characteristics are reasonably well known for P and S waves in the European Arctic, however much work remains to be done regarding surface wave propagation and magnitude estimation.¶Regional long-period or broadband seismic data in digital form has been available in the European Arctic for only a few years. In order to assess regional surface wave propagation, and in particular to evaluate the M _s :m _b discriminant at regional distances, it is therefore necessary to take advantage of the historic analog recordings. The station APA in Apatity forms a unique source of such data, with high-quality long-period seismic recordings of regional earthquakes and nuclear explosions dating back about 30 years.¶This paper presents initial results from a project to digitize APA surface waves of selected regional events. The recordings for recent years have been compared to a colocated broadband Guralp three-component seismometer in order to verify the response characteristics and the quality of the digitization process. It turns out that the quality of the digitized records is excellent, and can be used over a spectral band ranging from 5?seconds to at least 30?seconds period.¶We demonstrate the capabilities of the APA surface wave recordings to provide a promising separation of earthquakes and explosions in the European Arctic over a range of frequencies using the M _s :m _b discriminant, although we note that additional work is required in regionalization of the propagation paths to take into account the major tectonic features in the region. We also note that the body-wave magnitudes provided by international agencies are not always reliable for events in this region, and must be reassessed in order to make full use of the earthquake-explosion discrimination potential.     

西吉县西南山区典型黄土地震滑坡高密度电法物探解译分析

本文运用高密度电法对宁夏西吉县西南山区典型的黄土地震滑坡进行了探测,并结合钻探资料进行验证分析,目的是查明滑坡区域的地层结构、黄土厚度、基岩埋深、富水地段以及空间展布等特征.结果表明:电阻率参数能够很好的反应研究区的地层分布特征,表层疏松干燥的黄土为高阻反映,随着埋深增加,含水量较大的黄土呈现出低阻反映,在黄土与泥岩接触带饱水区呈现低阻反映;滥泥河流域典型的黄土梁斜坡具有阳坡黄土沉积薄、富水层薄、基岩埋深浅的特征,而在阴坡则表现出黄土沉积厚、富水层厚、基岩埋深大的特征.探测结果可有效的应用于黄土地震滑坡的勘察,为进一步开展黄土地震滑坡的成因机理和防治研究提供参考数据.     

Predictability of storm wave heights in the ice-free Beaufort Sea

A predictability study on wave forecast of the Arctic Ocean is necessary to help identify hazardous areas and ensure sustainable shipping along the trans-Arctic routes. To assist with validation of the Arctic Ocean wave model, two drifting wave buoys were deployed off Point Barrow, Alaska for two months in September 2016. Both buoys measured significant wave heights exceeding 4 m during two different storm events on 19 September and 22 October. The NOAA-WAVEWATCH III^? model with 16-km resolution was forced using wind and sea ice reanalysis data and obtained general agreement with the observation. The September storm was reproduced well; however, model accuracy deteriorated in October with a negative wave height bias of around 1 m during the October storm. Utilising reanalysis data, including the most up-to-date ERA5, this study investigated the cause: grid resolution, wind and ice forcing, and in situ sea level pressure observations assimilated for reanalysis. The analysis has found that there is a 20% reduction of in situ SLP observations in the area of interest, presumably due to fewer ships and deployment options during the sea ice advance period. The 63-member atmospheric ensemble reanalysis, ALERA2, has shown that this led to a larger ensemble spread in the October monthly mean wind field compared to September. Since atmospheric physics is complex during sea ice advance, it is speculated that the elevated uncertainty of synoptic-scale wind caused the negative wave model bias. This has implications for wave hindcasts and forecasts in the Arctic Ocean.     

Ground surface response induced by shallow buried explosions

Based on the spherical cavity expansion theory in the elastic half space,the ground surface movement characteristics of shallowly buried explosions are analyzed.The results show that the induced seismic wave is a longitudinal wave in the near zone and a Rayleigh wave in the far zone.The maximum displacement(velocity) of the longitudinal wave and the Rayleigh wave are inversely proportional to the scaled distance,and can be described by exponential function with exponents equal to 1.4 and 0.5,respectively.The vibration frequencies of the waves have almost no change.The vibration frequency of the longitudinal wave approximates the natural vibration frequency of the cavity in the broken area,and the vibration frequency of the Rayleigh wave is about half that of the longitudinal wave.On the same reduced buried depth and reduced distance,the particle displacement is directly proportional to the product of the boundary loading and cavity radius,and is inversely proportional to the transversal wave velocity.Meanwhile,the particle velocity is directly proportional to the boundary loading and inversely proportional to the wave velocity ratio.In the far zone,the buried depth of the explosive only has a slight effect on the longitudinal wave,but has a larger effect on the Rayleigh wave.