Frequency Domain Response at Pacific Coast Harbors to Major Tsunamis of 2005–2011

Tsunamis waves caused by submarine earthquake or landslide might contain large wave energy, which could cause significant human loss and property damage locally as well as in distant region. The response of three harbors located at the Pacific coast (i.e. Crescent City Harbor, Los Angeles/Long Beach Port, and San Diego Harbor) to six well-known tsunamis events generated (both near-field and far-field) between 2005 and 2011 are examined and simulated using a hybrid finite element numerical model in frequency domain. The model incorporated the effects of wave refraction, wave diffraction, partial wave reflection from boundaries, entrance and bottom energy dissipation. It can be applied to harbor regions with arbitrary shapes and variable water depth. The computed resonant periods or modes of oscillation for three harbors are in good agreement with the energy spectral analysis of the time series of water surface elevations recorded at tide gauge stations inside three harbors during the six tsunamis events. The computed wave induced currents based on the present model are also in qualitative agreement with some of the reported eye-witness accounts absence of reliable current data. The simulated results show that each harbor responded differently and significantly amplified certain wave period(s) of incident wave trains according to the shape, topography, characteristic dimensions and water depth of the harbor basins.     

Numerical study on transient harbor oscillations induced by successive solitary waves

Tsunamis are traveling waves which are characterized by long wavelengths and large amplitudes close to the shore. Due to the transformation of tsunamis, undular bores have been frequently observed in the coastal zone and can be viewed as a sequence of solitary waves with different wave heights and different separation distances among them. In this article, transient harbor oscillations induced by incident successive solitary waves are first investigated. The transient oscillations are simulated by a fully nonlinear Boussinesq model, FUNWAVE-TVD. The incident successive solitary waves include double solitary waves and triple solitary waves. This paper mainly focuses on the effects of different waveform parameters of the incident successive solitary waves on the relative wave energy distribution inside the harbor. These wave parameters include the incident wave height, the relative separation distance between adjacent crests, and the number of elementary solitary waves in the incident wave train. The relative separation distance between adjacent crests is defined as the ratio of the distance between adjacent crests in the incident wave train to the effective wavelength of the single solitary wave. Maximum oscillations inside the harbor excited by various incident waves are also discussed. For comparison, the transient oscillation excited by the single solitary wave is also considered. The harbor used in this paper is assumed to be long and narrow and has constant depth; the free surface movement inside the harbor is essentially one-dimensional. This study reveals that, for the given harbor and for the variation ranges of all the waveform parameters of the incident successive solitary waves studied in this paper, the larger incident wave heights and the smaller number of elementary solitary waves in the incident tsunami lead to a more uniform relative wave energy distribution inside the harbor. For the successive solitary waves, the larger relative separation distance between adjacent crests can cause more obvious fluctuations of the relative wave energy distribution over different resonant modes. When the wave height of the elementary solitary wave in the successive solitary waves equals to that of the single solitary wave and the relative separation distance between adjacent crests is equal to or greater than 0.6, the maximum oscillation inside the harbor induced by the successive solitary waves is almost identical to that excited by the single solitary wave.     

Effects of Harbor Modification on Crescent City, California��s Tsunami Vulnerability

More damaging tsunamis have impacted Crescent City, California in historic times than any other location on the West Coast of the USA. Crescent City??s harbor has undergone significant modification since the early 20th century, including construction of several breakwaters, dredging, and a 200?×?300?m² small boat basin. In 2006, a M _w 8.3 earthquake in the Kuril Islands generated a moderate Pacific-wide tsunami. Crescent City recorded the highest amplitudes of any tide gauge in the Pacific and was the only location to experience structural damage. Strong currents damaged docks and boats within the small boat basin, causing more than US?$20 million in damage and replacement costs. We examine how modifications to Crescent City??s harbor may have affected its vulnerability to moderate tsunamis such as the 2006 event. A bathymetric grid of the basin was constructed based on US Army Corps of Engineers soundings in 1964 and 1965 before the construction of the small boat basin. The method of splitting tsunamis was used to estimate tsunami water heights and current velocities at several locations in the harbor using both the 1964?C1965 grid and the 2006 bathymetric grid for the 2006 Kuril event and a similar-sized source along the Sanriku coast of Japan. Model velocity outputs are compared for the two different bathymetries at the tide gauge location and at six additional computational sites in the harbor. The largest difference between the two grids is at the small boat basin entrance, where the 2006 bathymetry produces currents over three times the strength of the currents produced by the 1965 bathymetry. Peak currents from a Sanriku event are comparable to those produced by the 2006 event, and within the boat basin may have been higher. The modifications of the harbor, and in particular the addition of the small boat basin, appear to have contributed to the high current velocities and resulting damage in 2006 and help to explain why the 1933 M _w 8.4?C8.7 Sanriku tsunami caused no damage at Crescent City.     

日本东北M_W9.0地震海啸在港池及邻近区域诱发的涡流危险性计算与评估分析

海啸造成的灾害与损失并非都与淹没有关,特别是港口中海啸诱导的强流会对船只及海事设施产生重要的影响及损害.由于海啸流观测数据稀缺及海啸诱导涡流机制的不确定性,过去60年海啸科学主要集中于对海啸波特征及淹没过程的研究与分析,海啸流模拟及验证工作开展较少,导致对海啸流基本特征及其造成灾害现象的曲解.开展海啸诱导的涡流研究及预警服务显得尤为重要及紧迫.考虑快速海啸预警需要,综合对比海啸诱导涡流的物理框架及模型方法,探索兼顾效率与计算精度的海啸流模拟方法是本文的核心工作及出发点.通过分析浅层湍流相干结构(TCS)产生的主要物理耗散机制,确定了考虑2D水平耗散机制的非线性浅水方程可用于海啸涡流的模拟分析.基于高精、高分辨率有限体积模型Geoclaw建立了三个精细化的港口海啸流模型,模型分辨率为5m.利用基于海啸浮标反演的海啸源模型作为初始条件,模拟分析了日本东北地震海啸在远场的海啸波流特征.海啸波流特征模拟结果与观测吻合较好,结果可信.对比发现:波驱动的自由表面流,小的位相或波幅误差就会导致大的流速误差,流的模拟和预报相对波幅来说更具挑战性.研究了海啸波流能量在港池中的分布特征,得到:港池入口及防波堤两端常被强流控制,具有极高的危险性;相对于波幅的空间变化,海啸流具有更强的空间敏感性;所建立的高分辨率海啸模型模拟再现了日本海啸在近场的涡旋结构,给出了与观测基本一致的涡流特征.最后,引入海啸流危险等级标准,分析了港口海啸流危险性等级分布、船只疏散的安全深度及回港的时间周期.针对港口、海湾同时考虑海啸波流特征的海啸预警与评估对于港口应急管理者科学决策具有重要意义.     

Numerical investigation of infragravity wave amplifications during harbor oscillations influenced by variable offshore topography

The infragravity (IG) period oscillations inside an elongated rectangular harbor near the offshore fringing reef induced by normal-incident bichromatic short wave groups are simulated using a fully nonlinear Boussinesq model, FUNWAVE 2.0. Based on an IG wave separation procedure, this article presents a systematical investigation on how the maximum IG period component amplitude, the bound and free IG waves, and their relative components inside the harbor change with respect to the plane reef-face slope and the incident short wave amplitude under the condition of the 2nd to the 5th modes. For the given harbor and the ranges of the reef-face slope and the incident short wave amplitude studied in this paper, it is shown that both the maximum IG period component amplitude and the free IG wave component amplitude inside the harbor fluctuate widely with the reef-face slope, and their changing trends with the reef-face slope are almost identical with each other, while the bound IG waves inside the harbor seem insensitive to it. Both the maximum IG period component amplitude and those of the bound and free IG standing waves inside the harbor change cubically with the incident short wave amplitude.     

Measurement of long waves at the harbor of Marina di Carrara,Italy

This paper presents the analysis of measurements of long waves at the harbor of Marina di Carrara, Italy. Nine pressure gauges (eight in the inner harbor and one outside the harbor) were used to record continuously the water surface elevation, for about 3 years; incoming short waves were measured using a directional buoy at about 500 m offshore the breakwaters. The analysis is carried out splitting the energy of the long waves into two bands, i.e., very long waves (VLW; f < 0.003 Hz) and long waves (LW; 0.003 < f < 0.030 Hz); it is found that LW energy is strongly correlated with the energy of the incoming short waves while no correlation (or very poor) is found for VLW. During energetic sea states, surface elevation spectra of VLW and LW, normalized against the wave energy of each frequency band, appear to be autosimilar; this suggests that the spectra of the incoming long waves have a mostly constant shape. Simple numerical computations of harbor resonance, carried out using a mild-slope equation linearized model, indicate that the amplification of the height of the incoming LW generated by short waves (calculated using empirical formulae from previous researches by Melito et al. 2007) inside the harbor can be reasonably predicted using such simple approach, assuming that the LW spectra are rectangular. The shape of the amplification diagram calculated by the numerical model shows reasonable (but not perfect) agreement with the shape of the normalized measured spectra inside the harbor.     

In situ 2-D and 3-D measurements of radiation patterns of half-wave dipole GPR antennas

Dipole antennas for ground-penetrating radar (GPR) radiate and receive electromagnetic waves with a strong directional dependence. Thus, experiments to measure in situ antenna radiation as functions of direction and polarization are of practical interest. Three field experiments were performed. One experiment was over a layered fluvial/eolian sequence; the other two used controlled targets (buried pipes and a metal ball). The radiation patterns were sampled by incrementally varying the antenna orientations and separations while recording reflections from the known targets. The results show qualitative, but systematic, correspondence with approximate theoretical far-field radiation patterns. Slow variations of amplitude with antenna azimuth and dip indicate that antenna orientations within 20° in the standard TE and TM acquisition geometries are adequate for most field applications, but not if detailed amplitude analysis is to be performed. Variations in antenna orientation or height (particularly for heights less than one-quarter wavelength) above the free surface introduce corresponding biases or uncertainties into recorded amplitudes. The variance within any suite of measurements is, in part, a consequence of differences in ground impedance at each antenna location. The theoretically predicted sensitivity to antenna height is mediated by surface roughness at high frequencies. It is necessary to include, or compensate for, the antenna radiation pattern in analysis of field data amplitudes, in experiment design, and in selecting appropriate antennas for specific applications.     

Improved determination of heights using a conversion surface by combining gravimetric quasi-geoid/geoid and GPS-levelling height differences

The quasi-geoid/geoid can be determined from the Global Positioning System (GPS) ellipsoidal height and the normal/orthometric heights derived from levelling (GPS-levelling). In this study a gravimetric quasigeoid and GPS-levelling height differences are combined to develop a new surface, suitable for “levelling” by GPS. This new surface provides better conversion of GPS ellipsoidal heights to the national normal heights. Different combining procedures, a four-parameter solution, linear and cubic splines interpolations, as well as the least-squares collocation method were investigated and compared over entire Norway. More than 1700 GPS-levelling stations were used in this study. The combined surface provides significant accuracy improvement for the normal height transformation of GPS height data, as demonstrated by the post-fitting residuals. The best solution, based on the least-squares collocation, provided a conversion surface for the transformation of GPS heights into normal height in Norway with an accuracy of about 5 cm.     

Effect of multilayered groundwater mounds on water dynamics beneath a recharge basin: Numerical simulation and assessment of surface injection

Interactions between groundwater mounds caused by a geologic layer contrast affect the efficiency of managed aquifer recharge in arid areas. However, research has rarely examined the roles of groundwater mounding size variations on soil water dynamics in a stratified vadose zone in response to a sustained infiltration source. Numerical experiments were conducted on a two-dimensional vertical-section domain using HYDRUS software to simulate the behaviours of two adjacent (upper and lower) groundwater mounds underlying an infiltration basin subjected to clay loam and sandy alternately-layered soil profiles. The model successfully predicted the volume and extent of perched water and approximated vertical travel times during events generating downward fluxes from the surface injection. The response time of the mounding width (lateral extension) to the surface injection was delayed as compared to that of the mounding height (vertical extension), especially for the lower water mound. The mounding heights and widths show a strongly positive correlation with the infiltration rates of both high- and low-permeability layers where the injected water mounded, while the water storage amounts in the high- and low-permeability layers were governed by the mounding height and width, respectively. Exploratory simulations were then employed to assess the dependence of groundwater mounding behaviours and recharge performances on surface injection strategies. Results suggest that, by reducing injection rate or shortening injection duration, the near-term fraction of the surface injection converted to deep recharge is likely to be increased due to the narrowed groundwater mounding size, which would be limited by the water-retarding effect of layer contrasts. This study has important implications for predicting and understanding multilayered groundwater mounding behaviours and associated water mass balance under the geologic stratification, and is expected to aid in optimizing the infiltration basin operation for aquifer recharge.     

Response of the nighttime midlatitude ionosphere to the passage of an atmospheric gravity wave

Using the data of the ionospheric vertical sounding in Almaty, the response of various parameters of the nighttime F layer to the passage of an atmospheric gravity wave, generated during the large magnetic storm on July 24–25, 2004, is studied. The analysis of the phase relations between the variations in the electron density at the F layer maximum (NmF), the layer maximum height (hmF), and the layer half-thickness showed that they are determined by the slope of the wave phase front. It is shown that the half-thickness of the layer changes in anti-phase with the variations in NmF2. The known fact that the amplitudes of variations in the critical frequencies of the F ₂ layer are smaller than the amplitudes of electron density variations at fixed heights is explained.     

Reaction of the ionospheric <Emphasis Type="Italic">F</Emphasis>1 region to disturbances in October and November 2003 (analysis from Irkutsk digisonde measurements)

We consider the ionospheric response at heights of 120–220 km to geomagnetic disturbances in October and November 2003, which caused a strongly pronounced decrease in the electron content in the noted height range. For this disturbance period, using a technique of the authors, midday variations of the relative content of atomic and molecular oxygen at a height of 120 km were estimated. This estimation was performed with the help of ionospheric measurements (by a digital ionosonde) conducted at the Institute of Solar-Terrestrial Physics of the Russian Academy of Sciences from 2003 to 2006. Comparison of these values with similar values from the mass spectrometer-incoherent scatter (MSIS)-86 model showed that our estimates during disturbance days were two times less. This study continues research into the field of using ionospheric measurements to estimate the relative gas composition of the thermosphere at heights below the maximum of the F2 layer.     

港湾中海啸流的模拟及其时空演变特征对输入条件不确定性的响应分析

一直以来,海啸波特征作为表征海啸潜在破坏性的参数指标得到了广泛应用,特别是针对近场极端海啸事件造成的灾害来说,这种表征具有较好的适用性.然而总结分析历史海啸事件造成的损失发现：在远场近岸及港湾系统中,海啸诱导的强流却是造成损失的主要原因.陆架或港湾振荡导致海啸波幅快速升降诱发强流,可能促使港工设施受到威胁及损害,进而对海啸预警服务及海事应急管理提出了新的挑战.因此,全面理解与评估海啸在港湾中诱发的灾害特征,探索港湾中海啸流的数值模拟方法,发展针对港湾尺度的海啸预警服务指导产品尤为迫切.受限于海啸流验证数据的缺乏及准确模拟海啸流技术方法的诸多不确定性,大部分海啸数值模拟研究工作主要是针对水位特征的研究及验证,可能导致对港湾中海啸灾害危险性认识的曲解与低估.本研究基于非线性浅水方程,针对夏威夷群岛三个典型港湾建立了精细化海啸数值模型（空间分辨率达到10 m）,并联合有限断层破裂模型计算分析了日本东北地震海啸在三个港湾及其邻近区域的海啸特征,波、流计算结果与实测结果吻合较好,精细化的海啸港湾模型模拟结果可信.模拟发现港湾中较小的波幅,同样可以产生强流.综合分析日本东北地震海啸波、流特征对输入条件不确定性的响应结果发现：港湾中海啸波-流能量的空间分布特征差异较大,这与港湾系统中海啸波的驻波特性相关;相比海啸波幅空间特征,海啸流特征具有更强的空间敏感性;海啸流时空分布特征对输入条件的不确定性响应比海啸波幅对这些不确定性的响应更强,海啸流的模拟与预报更有挑战性;不确定性对海啸流计算精度的影响会进一步传导放大港湾海啸流危险性的评估及对港工设施产生的应力作用的误差,合理的输入条件对海啸流的精确模拟至关重要.最后,希望通过本文的研究可以从海啸波-流特征角度更加全面认识近岸海啸灾害特征,拓展海啸预警服务的广度与深度,从而为灾害应急管理部门提供更加科学合理的辅助决策产品.     

无破碎带断层场地对SH波的散射

在层状半空间精确动力刚度矩阵和斜线荷载动力格林函数的基础上建立间接边界元方法,在频域内求解无破碎带断层场地对入射平面SH波的散射。为方便求解,将总波场分解为自由波场和散射波场,自由波场由直接刚度法求得,断层两侧的散射波场通过在断层面上分别对两侧施加均布斜线荷载产生的动力响应来模拟,虚拟荷载的密度可通过引入断层表面的边界条件确定,最后叠加自由波场和散射波场求得总波场。以有落差断层和无落差断层模型为例进行数值计算,分析断层落差、断层倾角以及断层两侧介质的刚度比对散射效应的影响。研究表明,断层落差与波长相当时,断层对SH波的放大作用最大;地表位移幅值随着断层倾角的增大逐渐增大;若断层无落差且其两侧刚度不同时,一般刚度较小一侧地表位移幅值较大且振荡更为剧烈,波从刚度较小一侧入射时位移幅值放大尤为显著。     

中国不同纬度背景Na层夜间和季节变化特征的激光雷达研究

基于子午工程的北京、合肥和海南三个Na荧光激光雷达对中国不同纬度上空夜间背景Na层的长期观测,分析了我国沿东经120°N上空Na层的夜间变化和季节变化特征.对照2010年12月17日夜间三个雷达站的观测结果,发现三个地方Na层的夜间变化并不具有相关性.Na层长期变化的年加半年变化拟合结果显示,北京和合肥上空Na层柱密度具有明显的年变化特点,而海南上空Na层柱密度的半年变化特征更明显; Na层的质心高度和RMS宽度具有明显的半年变化特点,但海南地区Na层的RMS宽度的长期变化不具有半年变化特征.Na层参数的统计和对照显示,Na层柱密度的季节变化与大气温度季节变化相关,在冬季最大,在夏季最小.Na层柱密度随纬度升高而增大,同时年变化性逐渐增强;质心高度随纬度变化趋势不明显,但三个地方Na层质心高度的长期变化都具有较明显的半年变化特征;在各个月份中,北京地区Na层RMS宽度最大,合肥地区Na层RMS宽度最小,海南地区居中.     

基于电离层反射的袖珍云闪(CID)三维定位研究

袖珍云闪是一类区别于常规闪电放电过程的特殊放电现象,能够同时产生极强的高频和低频辐射信号,其低频辐射信号在电离层与地面之间反射后能够在电场变化波形上形成电离层反射脉冲对.电离层反射信号与原信号的时间差包含着放电三维位置和电离层高度的信息,而借助于多站闪电探测网络的同步观测就能够反演这些信息.基于这一规律,本文发展了一种对袖珍云闪实现三维定位的新方法.这种方法不仅能够对大范围内的袖珍云闪实现准确的三维定位,同时还能够反演电离层的高度,是一种潜在的研究电离层相关性质的有效手段.通过将定位结果与雷达回波比较,证明这种方法具有较高的精度.利用这种方法,计算了5489例正极性袖珍云闪和1400例负极性袖珍云闪的放电高度,发现正极性袖珍云闪主要集中在7~14 km,而负极性袖珍云闪达到了15~18 km.负极性袖珍云闪的放电高度总体上与对流层顶高度相当,其数量相比于正极性袖珍云闪明显偏少,因此很可能产生于较为罕见的极旺盛的雷暴过程中.     

Variations in the phase of the semidiurnal tide over Davis,Antarctica

The amplitude and time of maximum of the semidiurnal tide in the mesosphere and lower thermosphere above Davis varies through the summer months. In particular, the time of maximum can oscillate around a fixed value or, at some heights, shift by a whole cycle of local time. Recent studies of the semi-diurnal tide at the South Pole have suggested that an s=1 mode is common at high-southern latitudes. At mid-latitudes, the s=2 mode is thought to dominate suggesting that a region where these modes overlap is a possibility. Although it is not possible to discern which of these modes are present using single station data, the concept that more than one 12-h wave is present allows a new interpretation of the Davis semi-diurnal amplitudes and times of maximum. In this study, wind data obtained at Davis, Antarctica, during the summer of 1996/97 are used to show that the complex variations in time of maximum and amplitude can be explained by a sum of two waves. The simplest of models combines an invariant and a varying 12-h wave and yields characteristics for each component that are less complex than the original observation. There is also potential for changes to the interpretation of time of maximum vs. height profiles.     

Experimental investigations on the seiches of Kake Constance

The seiches of Lake Constance and their possible exciting mechanisms were studied. The fundamental oscillation with a period of 55.6 minutes and other modes were observed at several stations. The dependance of the oscillation period on the water level is proved. By a model calculation it is shown what amplitudes can be excited by a wavelike air pressure disturbance which travels over the lake.     

Effect of soil water-repellent layer depth on post-wildfire hydrological processes

Soil water repellency induced by wildfires can alter hydraulic properties and hydrologic processes; however, the persistence and vertical position (i.e., depth) of water-repellent layers can vary between systems and fires, with limited understanding of how those variations affect infiltration processes. This study occurred in two forested locations in the south-central Appalachian Mountains that experienced wildfires in late 2016: Mount Pleasant Wildfire Refuge, Virginia, and Chimney Rock State Park, North Carolina. In each location, sites were selected to represent unburned conditions and low to moderate burn intensities. At each site, we measured the soil water repellency at the surface (ash layer or O horizon) and ~2 cm below the surface (A horizon) using the water drop penetration time method (n = 10–14). Soil water content was also measured over the upper 10 cm of the soil (n = 10), and infiltration tests were conducted using a tension infiltrometer (n = 6–8). The results showed that soil repellency was highest in the surface layer at the Mount Pleasant location and was highest in the subsurface layer at the Chimney Rock location. Soil water content was lower in unburned soil than in burned soil, especially for measurements taken immediately postfire, with soil water content negatively correlated with water repellency. Water repellency in the surface layer significantly reduced relative infiltration rates (estimated as differences between initial and steady-state rates), whereas subsurface water repellency did not affect relative infiltration. As a result, water repellency persisted longer in sites with surface as opposed to subsurface water repellency. Finally, differences between burned and unburned sites showed that although the wildfires increased the occurrence of water repellency, they did not alter the underlying relationship between relative infiltration and water repellency of the surface soil.