A Unified Variational Principle of Fluid Mechanics and Application on Solitary Subdomain or Point

-According to basic equations of fluid mechanics, this paper presents a unified variational principle of fluid mechanics (UVPFM) by using the optimization method of weighted residuals (OMWR). The advantages are as follows, the establishment of the functional and the variational principle is easy, it can change various problems of fluid mechanics derived by basic equations into a unified optimization problem, and the solution is the optimum one in some sense. According to the OMWR for the solitary subdomain, this paper uses UVPFM onto any solitary subdomain and gives the solution of the hydrodynamics equation which is suitable only for that solitary subdomain. According to the OMWR for solitary point, this paper uses UVPFM to any solitary point and gives the solution of the hydrodynamics equation (point solution) which is suitable only for that solitary point. As the solution for the solitary subdomain or solitary point is developed independently, the compatibility with other subdomain or other points, do     

Variational analysis of dynamics of fissured poroelastic rocks

Rocks can be modeled in a continuum framework as fissured, poroelastic materials, i.e., materials with two degrees of porosity, one due to the fissures and another one due to the pores. The governing equations of motion of fissured poroelastic rocks established by Beskos are rederived here by establishing a variational statement which also provides the boundary conditions of the problem. This is accomplished by considering strain, dissipation and kinetic energies as well as the work of external forces. The above statement is also derived here by employing the method of weighted residuals.     

Calibration of paired watersheds: Utility of moving sums in presence of externalities

Historically, paired watershed studies have been used to quantify the hydrological effects of land use and management practices by concurrently monitoring 2 similar watersheds during calibration (pretreatment) and post‐treatment periods. This study characterizes seasonal water table and flow response to rainfall during the calibration period and tests a change detection technique of moving sums of recursive residuals (MOSUM) to select calibration periods for each control–treatment watershed pair when the regression coefficients for daily water table elevation were most stable to minimize regression model uncertainty. The control and treatment watersheds were 1 watershed of 3–4‐year‐old intensely managed loblolly pine (Pinus taeda L.) with natural understory, 1 watershed of 3–4‐year‐old loblolly pine intercropped with switchgrass (Panicum virgatum), 1 watershed of 14–15‐year‐old thinned loblolly pine with natural understory (control), and 1 watershed of switchgrass only. The study period spanned from 2009 to 2012. Silvicultural operational practices during this period acted as external factors, potentially shifting hydrologic calibration relationships between control and treatment watersheds. MOSUM results indicated significant changes in regression parameters due to silvicultural operations and were used to identify stable relationships for water table elevation. None of the calibration relationships developed using this method were significantly different from the classical calibration relationship based on published historical data. We attribute that to the similarity of historical and 2010–2012 leaf area index on control and treatment watersheds as moderated by the emergent vegetation. Although the MOSUM approach does not eliminate the need for true calibration data or replace the classic paired watershed approach, our results show that it may be an effective alternative approach when true data are unavailable, as it minimizes the impacts of external disturbances other than the treatment of interest.     

Travel time residuals in southeastern Europe

The Pn travel time relative residuals, in respect to a crustal model of the Aegean area, have been determined for 103 permanent seismological stations in southeastern Europe, western Turkey and the Middle East. The values of these residuals are considered to depend mainly on the crustal thickness beneath the seismological stations. Based on these values seven regions with different crustal thickness, varying between 31 Km and 42 Km, have been defined. The crust in these regions is continental. A region with very high negative residuals has been defined in the Middle East (Egypt, Israel, Lebanon). These negative residuals are attributed to different crustal structure of the eastern Mediterranean (oceanic crust with an extra thick sedimentary layer) and not the crustal thickness at the station sites.Independently from the interpretation, these Pn residuals can be used successfully to considerably improve (up to 2 Km) the determination of the earthquake foci locations.     

Travel-time residuals and three-dimensional velocity structure of Italy

Jeffreys-Bullen P and PKP travel-time residuals observed at more than 50 seismic stations distributed along Italy and surrounding areas in the time interval 1962–1979, indicate the complex velocity pattern of this region. Strong lateral velocity inhomogeneities and low velocity zones are required to explain the observed pattern of residuals. In particular, late arrivals of about 1 sec are observed in the Apenninic mountain range, requiring both greater crustal thickness and low velocity layers, coherent with seismic refraction data and surface wave dispersion measurements. The seismic stations located in the Western and Eastern Alps indicate the presence of high velocities. In the Western Alps the strong azimuthal variation of residuals and the high values of early arrivals have a close relationship to the Ivrea body, an intrusive crustal complex characterized by a velocity as high as 7–7.2 km/sec.A travel-time inversion performed with theAki et al. (1977) block model, confirms the peculiar characteristics and the sharp variations in the lithosphere of the whole Italian region, with values of velocity perturbations between many adjacent blocks, ranging in size from 50 to 100 km, and independent from the earth parametrization chosen, reaching values up to 10% in the lithospheric part and 5% in the asthenosphere. 3-D inversion requires also high velocity along the Tyrrhenian coastal margin, equivalent to an uprise of major crustal and lithospheric discontinuities along this part of the Italian peninsula. Moreover low velocity material must be present in the northern part of the Adriatic foreland, in the lithosphere-asthenosphere system, closely related to the stress and seismicity pattern, and the lateral bending of the lithosphere in the same region.     

Crustal transit time monitoring using PKP

The arrival time difference for the AB branch of PKP from deep Tonga earthquakes is constant over years with a standard deviation of ±0.05 seconds at seismographs located 10 to 50 km from each other. If published travel time curves are used to calculate the relative residuals of PKP the standard deviation from the constant mean is improved by approximately 0.01 seconds for AB branch data. For the BC branch, standard deviations of relative travel times of ±0.06 seconds are reduced to less than ±0.05 seconds by calculating relative residuals. We conclude that changes of crustal transit time forP-waves could be resolved, based on careful PKP arrival time measurement at two or more neighboring stations if the changes exceed 0.05 sec and last for more than one year. The conditions for achieving this result are that PKP from Tonga earthquakes is clearly recorded, and that time-keeping is accurate. The data on which these conclusions are based were obtained from the Graefenberg seismograph array, which is located in West Gemany and consists of 13 stations separated by distances of 10 km to 100 km. We propose that relative arrival times of PKP from Tonga could be used in the Mediterranean - Middle East area to search for precursory travel time changes before large earthquakes.     

联邦滤波的改进算法

分析并综述了联邦滤波算法,以及该算法存在的问题.针对观测信息不足的情况下,结合抗差估计和自适应估计理论对联邦滤波算法进行了改进,给出了联邦滤波在局部传感器观测信息不足情况下的改进算法.由模拟计算结果可知,在局部传感器观测信息不足的情况下,改进的联邦滤波器能较好地抑制载体观测异常和状态扰动异常对动态系统参数估值的影响,具有很强的的抗差性和自适应性.