Observations and modelling of the travel time delay and leading negative phase of the 16 September 2015 Illapel,Chile tsunami

The systematic discrepancies in both tsunami arrival time and leading negative phase (LNP) were identified for the recent transoceanic tsunami on 16 September 2015 in Illapel, Chile by examining the wave characteristics from the tsunami records at 21 Deep-ocean Assessment and Reporting of Tsunami (DART) sites and 29 coastal tide gauge stations. The results revealed systematic travel time delay of as much as 22 min (approximately 1.7％ of the total travel time) relative to the simulated long waves from the 2015 Chilean tsunami. The delay discrepancy was found to increase with travel time. It was difficult to identify the LNP from the near-shore observation system due to the strong background noise, but the initial negative phase feature became more obvious as the tsunami propagated away from the source area in the deep ocean. We determined that the LNP for the Chilean tsunami had an average duration of 33 min, which was close to the dominant period of the tsunami source. Most of the amplitude ratios to the first elevation phase were approximately 40％, with the largest equivalent to the first positive phase amplitude. We performed numerical analyses by applying the corrected long wave model, which accounted for the effects of seawater density stratification due to compressibility, self-attraction and loading (SAL) of the earth, and wave dispersion compared with observed tsunami waveforms. We attempted to accurately calculate the arrival time and LNP, and to understand how much of a role the physical mechanism played in the discrepancies for the moderate transoceanic tsunami event. The mainly focus of the study is to quantitatively evaluate the contribution of each secondary physical effect to the systematic discrepancies using the corrected shallow water model. Taking all of these effects into consideration, our results demonstrated good agreement between the observed and simulated waveforms. We can conclude that the corrected shallow water model can reduce the tsunami propagation speed and reproduce the LNP, which is observed for tsunamis that have propagated over long distances frequently. The travel time delay between the observed and corrected simulated waveforms is reduced to <8 min and the amplitude discrepancy between them was also markedly diminished. The incorporated effects amounted to approximately 78％ of the travel time delay correction, with seawater density stratification, SAL, and Boussinesq dispersion contributing approximately 39％, 21％, and 18％, respectively. The simulated results showed that the elastic loading and Boussinesq dispersion not only affected travel time but also changed the simulated waveforms for this event. In contrast, the seawater stratification only reduced the tsunami speed, whereas the earth's elasticity loading was responsible for LNP due to the depression of the seafloor surrounding additional tsunami loading at far-field stations. This study revealed that the traditional shallow water model has inherent defects in estimating tsunami arrival, and the leading negative phase of a tsunami is a typical recognizable feature of a moderately strong transoceanic tsunami. These results also support previous theory and can help to explain the observed discrepancies.     

Duality properties of Gorringe–Leach equations

In the category of motions preserving the angular momentum direction, Gorringe and Leach exhibited two classes of differential equations having elliptical orbits. After enlarging slightly these classes, we show that they are related by a duality correspondence of the Arnold–Vassiliev type. The specific associated conserved quantities (Laplace–Runge–Lenz vector and Fradkin–Jauch–Hill tensor) are then dual reflections of each other.     

Analytical solutions for the equations of motion of a space vehicle during the atmospheric re-entry phase on a 2-D trajectory

A practical and important problem encountered during the atmospheric re-entry phase is to determine analytical solutions for the space vehicle dynamical equations of motion. The author proposes new solutions for the equations of trajectory and flight-path angle of the space vehicle during the re-entry phase in Earth’s atmosphere. Explicit analytical solutions for the aerodynamic equations of motion can be effectively applied to investigate and control the rocket flight characteristics. Setting the initial conditions for the speed, re-entering flight-path angle, altitude, atmosphere density, lift and drag coefficients, the nonlinear differential equations of motion are linearized by a proper choice of the re-entry range angles. After integration, the solutions are expressed with the Exponential Integral, and Generalized Exponential Integral functions. Theoretical frameworks for proposed solutions as well as, several numerical examples, are presented.     

Solution of Lane-Emden type equations using Bernstein operational matrix of differentiation

The purpose of this paper is to propose an efficient numerical method for solving Lane-Emden type equations arising in astrophysics using Bernstein polynomials. First Bernstein operational matrix of differentiation is derived using Bernstein polynomials and then applied to solve the linear and nonlinear differential equations of Lane-Emden type. Some illustrative examples are given to demonstrate the efficiency and validity of the proposed algorithm.     

二维圆柱坐标下FDTD法对多频电磁波测井的数值模拟

在介质围绕Z轴旋转对称的情况下,三维直角坐标下的麦克斯韦方程可变成两组分别对应TM波和TE波的二维柱坐标下的偏微分方程.采用线圈作为电磁波测井的发射器和接收器的时候,可采用TE波对应的方程,这时相当于磁极子的情况.而采用电偶极子作为发射和接收器的时候,则对应TM波的方程.采用时间和空间均为二阶精确的有限差分方法,将偏微分方程进行差分化.这样,空间的电磁场可由时间域有限差分法(FDTD)来求解.     

Numerical study for vegetation effects on coastal wave propagation by using nonlinear Boussinesq model

The vegetation has important impacts on coastal wave propagation. In the paper, the sensitivities of coastal wave attenuation due to vegetation to incident wave height, wave period and water depth, as well as vegetation configurations are numerically studied by using the fully nonlinear Boussinesq model. The model is based on the implementation of drag resistances due to vegetation in the fully nonlinear Boussinesq equation where the drag resistance is provided by the Morison’s formulation for rigid structure induced drag stresses. The model is firstly validated by comparing with the experimental results for wave propagation in vegetation zones. Subsequently, the model is used to simulate waves with different height, period propagating on vegetation zones with different water depth and vegetation configurations. The sensitivities of wave attenuation to incident wave height, wave period, water depth, as well as vegetation configurations are investigated based on the numerical results. The numerical results indicate that wave height attenuation due to vegetation is sensitive to incident wave height, wave period, water depth, as well as vegetation configurations, and attenuation ratio of wave height is increased monotonically with increases of incident wave height and decreases of water depth, while it is complex for wave period. Moreover, more vegetation segments can strengthen the interaction of vegetation and wave in a certain range.     

A geomorphology‐based integrated stream–aquifer interaction model for semi‐gauged catchments

Many of the existing stream–aquifer interaction models available in the literature are very complex with limited applicability in semi‐gauged and ungauged catchments. In this study, to estimate the influent and effluent subsurface water fluxes under limited geo‐hydrometeorological data availability conditions, a simple stream–aquifer interaction model, namely, the variable parameter McCarthy–Muskingum (VPMM) hillslope‐storage Boussinesq (hsB) model, has been developed. This novel model couples the VPMM streamflow transport with the hsB groundwater flow transport modules in online mode. In this integrated model, the surface water–groundwater flux exchange process is modelled by the Darcian approach with the variable hydraulic heads between the river stage and groundwater table accounting for the rainfall forcing. Considering the exchange fluxes in the hyporheic zone and lateral overland flow contribution, this approach is field tested in a typical 48‐km stretch of the Brahmani River in eastern India to simulate the streamflow and its depth with the minimum Nash–Sutcliffe efficiency of 94% and 88%; the maximum root mean square error of 134 m³/s and 0.35 m; and the minimum index of agreement of 98% and 97%, respectively. This modelling approach could be very well utilized in data‐scarce world‐river basins to estimate the stream–aquifer exchange flux due to rainfall forcings.     

最小二乘定权方法对流动重力计算结果的影响

在处理流动重力观测数据时,采用最小二乘方法解算由权系数和观测值组成的超定方程组,得到最优权值,并利用最优权计算重力值。结果表明,等权方法计算得到的结果点值精度不高,但包含了较全的重力异常地震分析预报信息;不等权方法计算得到的结果点值精度较高,但削弱了有用的重力异常地震分析预报信息;最小二乘方法计算得到的结果不但有较高的点值精度,对M_S4.0及以上地震的分析预报效果也最优。     

扩展共线方程结合车载组合全景相机影像区域网平差

利用扩展共线方程实现组合相机单片影像的空中三角测量,通过该方法来提高全景影像的定位精度。由于组合全景相机中各子相机刚性固定,相机间的相对位姿参数可以精确标定,因此利用扩展共线方程模型,可以将组合全景相机的构像单元从单个子相机影像扩展为整个组合相机获取的多张影像。与已有的球模型方法,包括球理想模型和球严密模型方法相比,本文所用的扩展共线方程方法使用各子相机影像像点作为观测值,避免了球模型方法中由各子相机摄影中心不重合导致的球模型投影误差和摄站内重叠区域融合误差。试验结果表明:与球模型方法相比,本文方法无需模型投影和摄站内重叠区域融合等繁琐的处理过程,避免了由其引起的精度损失,保证了空三网稳定性的同时提高了定位精度。     

电力系统元件非线性微分-代数子系统模型的逆系统控制:一种新算法

针对电力系统元件非线性微分-代数子系统模型,本文提出一种新算法研究其逆系统控制问题.所提出的新算法不需要对控制输出及其高阶导数做复杂的变换,具有更好的应用性.本文的逆系统控制方法主要分为两步：第1步,利用所提出的新算法来判断被控元件的可逆性,若可逆,则基于状态反馈与动态补偿,构造出元件的α阶积分右逆系统,实现复合系统的线性化和解耦;第2步,利用线性控制的理论和方法设计闭环控制器,使得元件被控对象满足期望的性能指标.最后按照本文所提出的方法,研究了多机电力系统的分散非线性汽门控制问题.仿真结果验证了本文所提方法的有效性.