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.     

Travel times for snowmelt-dominated headwater catchments: Influences of wetlands and forest harvesting,and linkages to stream water quality

The time it takes water to travel through a catchment, from when it enters as rain and snow to when it leaves as streamflow, may influence stream water quality and catchment sensitivity to environmental change. Most studies that estimate travel times do so for only a few, often rain-dominated, catchments in a region and use relatively short data records (<10 years). A better understanding of how catchment travel times vary across a landscape may help diagnose inter-catchment differences in water quality and response to environmental change. We used comprehensive and long-term observations from the Turkey Lakes Watershed Study in central Ontario to estimate water travel times for 12 snowmelt-dominated headwater catchments, three of which were impacted by forest harvesting. Chloride, a commonly used water tracer, was measured in streams, rain, snowfall and as dry atmospheric deposition over a 31 year period. These data were used with a lumped convolution integral approach to estimate mean water travel times. We explored relationships between travel times and catchment characteristics such as catchment area, slope angle, flowpath length, runoff ratio and wetland coverage, as well as the impact of harvesting. Travel time estimates were then used to compare differences in stream water quality between catchments. Our results show that mean travel times can be variable for small geographic areas and are related to catchment characteristics, in particular flowpath length and wetland cover. In addition, forest harvesting appeared to decrease mean travel times. Estimated mean travel times had complex relationships with water quality patterns. Results suggest that biogeochemical processes, particularly those present in wetlands, may have a greater influence on water quality than catchment travel times.     

Construction of ray synthetic seismograms using interpolation of travel times and ray amplitudes

The seismic wave field, in its high-frequency asymptotic approximation, can be interpolated from a low- to a high-resolution spatial grid of receivers and, possibly, point sources by interpolating the eikonal (travel time) and the amplitude. These quantities can be considered as functions of position only. The travel time and the amplitude are assumed to vary in space only slowly, otherwise the validity conditions of the theory behind would be violated. Relatively coarse spatial sampling is then usually sufficient to obtain their reasonable interpolation. The interpolation is performed in 2-D models of different complexity. The interpolation geometry is either 1-D, 2-D, or 3-D according to the source-receiver distribution. Several interpolation methods are applied: the Fourier interpolation based on the sampling theorem, the linear interpolation, and the interpolation by means of the paraxial approximation. These techniques, based on completely different concepts, are tested by comparing their results with a reference ray-theory solution computed for gathers and grids with fine sampling. The paraxial method holds up as the most efficient and accurate in evaluating travel times from all investigated techniques. However, it is not suitable for approximation of amplitudes, for which the linear interpolation has proved to be universal and accurate enough to provide results acceptable for many seismological applications.     

Grid travel-time tracing: Second-order method for the first arrivals in smooth media

A new computational scheme for calculating the first-arrival travel times on a rectangular grid of points is proposed. The new proposed method is of second-order accuracy. This means that the error of the calculated travel time is proportional to the second power of the grid spacing. The method should be sufficiently accurate for all applications in smooth seismic models. On the other hand, the method is not, in its present form, proposed for models with structural interfaces which make the method unstable and generate travel-time errors of the first order. Equations are also presented for the appropriate evaluation of the errors of calculated travel times to check their accuracy, and the proposed method is compared with other numerical methods. The method is developed, described and demonstrated in 2-D, but may also be extended to 3-D models and to general models with structural interfaces.     

Application of Snell's law in reflection raytracing using the multistage fast marching method

Accurate calculations of travel times and raypaths of reflection waves are important for reflection travel time tomography. The multistage shortest path method(MSPM) and multistage fast marching method(MFMM) have been widely used in reflection wave raytracing, and both of them are characterized by high efficiency and accuracy. However, the MSPM does not strictly follow Snell's law at the interface because it treats the interface point as a sub-source, resulting in a decrease in accuracy. The MFMM achieves high accuracy by solving the Eikonal equation in local triangular mesh. However, the implementation process is complex. Here we propose a new method which uses linear interpolation to compute the incident travel time of interface points and then using Snell's law to compute the reflection travel time of grid points just above the interface. Our new method is much simpler than the MFMM; furthermore, numerical simulations show that the accuracy of the MFMM and our new method are basically the same, thus the reflection tomography algorithms which use our new method are easier to implement without decreasing accuracy. Besides, our new method can be extended easily to other grid-based raytracing methods.     

A hybrid approach to compute seismic travel times in three-dimensional tetrahedral meshes

We propose an optimized method to compute travel times for seismic inversion problems. It is a hybrid method combining several approaches to deal with travel time computation accuracy in unstructured meshes based on tetrahedral elementary cells. As in the linear travel time interpolation method, the proposed approach computes travel times using seismic ray paths. The method operates in two sequential steps: At a first stage, travel times are computed for all nodes of the mesh using a modified version of the shortest path method. The difference with the standard version is that additional secondary nodes (called tertiary nodes) are added temporarily around seismic sources in order to improve accuracy with a reasonable increase in computational cost. During the second step, the steepest travel time gradient method is used to trace back ray paths for each source–receiver pair. Travel times at each receiver are then recomputed using slowness values at the intersection points between the ray path and the traversed cells. A number of numerical tests with an array of different velocity models, mesh resolutions and mesh topologies have been carried out. These tests showed that an average relative error in the order of 0.1% can be achieved at a computational cost that is suitable for travel time inversion.     

Stochastic ground water flow simulation with a fracture zone continuum model

A method is presented for incorporating the hydraulic effects of vertical fracture zones into two-dimensional cell-based continuum models of ground water flow and particle tracking. High hydraulic conductivity features are used in the model to represent fracture zones. For fracture zones that are not coincident with model rows or columns, an adjustment is required for the hydraulic conductivity value entered into the model cells to compensate for the longer flowpath through the model grid. A similar adjustment is also required for simulated travel times through model cells. A travel time error of less than 8% can occur for particles moving through fractures with certain orientations. The fracture zone continuum model uses stochastically generated fracture zone networks and Monte Carlo analysis to quantify uncertainties with simulated advective travel times. An approach is also presented for converting an equivalent continuum model into a fracture zone continuum model by establishing the contribution of matrix block transmissivity to the bulk transmissivity of the aquifer. The methods are used for a case study in west-central Florida to quantify advective travel times from a potential wetland rehydration site to a municipal supply wellfield. Uncertainties in advective travel times are assumed to result from the presence of vertical fracture zones, commonly observed on aerial photographs as photolineaments.     

Scattering objects in the lower mantle beneath northeastern China observed with a short-period seismic array

Precursor and coda portions of short-period PcP waves (reflected P wave from the core-mantle boundary, CMB) recorded at J-array stations in Japan were analyzed in order to extract weak scattered signals originating from small-scale heterogeneities in the lowermost mantle beneath northeastern China. Two nuclear explosions at Lop Nor in China detonated on 21 May 1992 (Mb=6.5) and 8 June 1996 (Mb=5.9) were used for our analysis.Three-dimensional grids above the CMB were defined in the area around the PcP bounce points beneath northeastern China to calculate theoretical travel times of scattered waves which propagate from the sources to each grid point and arrive at each station based on the IASP91 model. Subsequently the waveforms were aligned with respect to the theoretical travel times and the semblance (an amplitude dependent measure of coherency) was calculated for each grid point. In order to obtain a more accurate travel time correction, we applied a cross correlation method to PcP waveforms in order to reduce picking error of the PcP onset time. A cross convolution method was also applied so that the two events could be analyzed simultaneously without using unstable deconvolutions.We could identify regions with relative high semblance values in semblance contour maps at about 200 and 375 km above the CMB. Stacking waveforms with respect to the theoretical travel times for the grid points with relative high semblance values indicate coherent wavelets originating at those grid points, that is, they correspond to scattered waves originating from small-scale heterogeneities in the lowermost mantle. Our results indicate the existence of small-scale scattering objects in the D″ layer, especially in the depth range of 200 and 375 km above the CMB beneath northeastern China. Considering recent tomographic images of high velocity anomalies in this area, these scattering objects could be fragments of old oceanic crusts which have subducted through the lower mantle and have accumulated in the D″ layer beneath northeastern China.     

SKKS震相记录特征及其对极远震深度的测定

就穿透地球内核的SKKS震相的传播路径、走时以及随深度变化的特性进行了分析讨论。对银川台记录SKKS震相的特征和震相识别做了研究总结，同时给出了SKKS震相求极远震深度列线图。经过震例检验，用SKKS震相测得的深度在单台定位中相对比较准确，为应用SKKS震相测准极远震深度，提高银川单台月报质量，提供了有利工具。     

Determination of upward/downward groundwater fluxes using transient variations of soil profile temperature: test of the method with Voyons (Aube,France) experimental data

The temperature variations recorded at several points of a vertical shallow‐depth profile are governed both by conductive and convective heat transfers and can be used to calculate the vertical component of the Darcy velocity and thermal diffusivity in the soil. This paper describes such calculations when transient variations over tens of days are considered and tests them using data collected at Voyons (Aube, France). The temperature was recorded during a year and a half period with a 1 h sampling time step at three different depths: 0·2, 0·4 and 0·75 m. By processing the annual variation of temperature, we obtained a value of the Darcy velocity in good agreement with the value of actual/potential evapotranspiration ratio. By processing transient variations, despite the limitation of the calculations due to the lack of sensitivity of the sensors, results obtained at Voyons were in good correlation with tensiometric data. Copyright © 2005 John Wiley & Sons, Ltd.     

Interpretation of first arrival travel times in seismic refraction work

The necessary condition for the seismic refraction method to succeed is that the refracted first arrivals from each layer in a multilayered earth system should be detected on a seismogram as first arrivals, and this is possible only when velocities of all underlying layers are successively greater. The usual procedure to interpret the refraction travel times is to fit such a data set with several intersecting straight lines by employing a visual technique which may lead to errors of subjective judgment, as the velocity model depends on the selection of various line segments through the data. To remove the visual fit we propose here a layer stripping method based on minimum intercept time, apparent velocity, rms residual, and maximum data points by least-squares fitting to yield several intersecting straight lines. Once data are segmented out, the conventional equations can be used to determine the velocity structure.     

考虑发射波形的瞬变电磁三维模型降阶快速正演算法

瞬变电磁响应受到发射波形的影响,正演计算中需考虑发射波形.基于位移逆Krylov子空间投影的模型降阶算法能够精确模拟三维全波形瞬变电磁正演响应,但该算法计算精确的on-time响应非常耗时.目前实际工作中多数情况下都是对off-time时间段的数据进行处理解释,因此可以将正演问题简化为精确模拟考虑发射波形的off-ti...     

A general model for Mt. Ruapehu lahars

A mathematical model of the motion of lahars is presented. Lahar flows and travel speeds are calculated using a kinematic wave model which equates gravitational accelerations to frictional losses. A chezyor Manning-type law of friction is assumed, in which lahar flow rate is a simple power function of lahar depth, multiplied by another simple power of the chanel slope. Use of the model requires knowledge of essentially only one parameter which appears to be relatively insensitive for flows down a given channel. Variable channel slope effects are removed by a longitudinal scaling which applies to all flows down a given channel. For lahars generated by a single explosive event it is unnecessary to perform numerical calculations to predict lahar flow and travel time, but for lahnars produced by multiple sources in which different lahar flows are interacting, numerical calculations appear necessary. The model is applied to all recorded lahar flows from Mt. Ruapehu, and satisfactorily described all lahar flows generated by a single explosive mechanism. Such flows depend essentially only on total lahar volume. The 1968 Mt. Ruapehu lahar, generated by a series of smaller eruptive mechanisms, was modelled as the interaction of seven point sources of fluid originating from positions mathematically extrapolated up the mountain. Good agreement was obtained between the predicted times of formation of these 1968 lahars, and the times of greatest seismic amplitude.     

Three-dimensional Travel-time Calculation Based on Fermat's Principle

— We present a new travel-time calculation method based on Fermat's principle. In the method, travel times are recursively calculated on horizontal planes of increasing depth. For typical configurations of exploration geophysics, the distance between the planes is on the order of 100 m. The travel times on the first plane are calculated by connecting straight ray segments from the source point to the grid points on the plane and integrating the slowness along each segment. The times on the other planes are calculated by finding the minimum of the combination of the times on the plane above plus the additional time along segments connecting grid points on the two planes. The travel-time calculation method is designed for calculating either the first arrival times, or the time of the shortest travel path arrival. The method is extended to handle vertically transverse isotropic (VTI) media by an approach which increases the computing time only slightly. The algorithm is tested against synthetic examples for isotropic and VTI wave propagation.     

Parsimonious wave‐equation travel‐time inversion for refraction waves

We present a parsimonious wave‐equation travel‐time inversion technique for refraction waves. A dense virtual refraction dataset can be generated from just two reciprocal shot gathers for the sources at the endpoints of the survey line, with N geophones evenly deployed along the line. These two reciprocal shots contain approximately 2N refraction travel times, which can be spawned into refraction travel times by an interferometric transformation. Then, these virtual refraction travel times are used with a source wavelet to create N virtual refraction shot gathers, which are the input data for wave‐equation travel‐time inversion. Numerical results show that the parsimonious wave‐equation travel‐time tomogram has about the same accuracy as the tomogram computed by standard wave‐equation travel‐time inversion. The most significant benefit is that a reciprocal survey is far less time consuming than the standard refraction survey where a source is excited at each geophone location.     

Waveform Retrieval and Phase Identification for Seismic Data from the CASS Experiment

The little destruction to the deployment site and high repeatability of the Controlled Accurate Seismic Source (CASS) shows its potential for investigating seismic wave velocities in the Earth's crust. However, the difficulty in retrieving impulsive seismic waveforms from the CASS data and identifying the seismic phases substantially prevents its wide applications. For example, identification of the seismic phases and accurate measurement of travel times are essential for resolving the spatial distribution of seismic velocities in the crust. Until now, it still remains a challenging task to estimate the accurate travel times of different seismic phases from the CASS data which features extended wave trains, unlike processing of the waveforms from impulsive events such as earthquakes or explosive sources. In this study, we introduce a time-frequency analysis method to process the CASS data, and try to retrieve the seismic waveforms and identify the major seismic phases traveling through the crust. We adopt the Wigner-Ville Distribution (WVD) approach which has been used in signal detection and parameter estimation for linear frequency modulation (LFM) signals, and proves to feature the best time-frequency convergence capability. The Wigner-Hough transform (WHT) is applied to retrieve the impulsive waveforms from multi-component LFM signals, which comprise seismic phases with different arrival times. We processed the seismic data of the 40-ton CASS in the field experiment around the Xinfengjiang reservoir with the WVD and WHT methods. The results demonstrate that these methods are effective in waveform retrieval and phase identification, especially for high frequency seismic phases such as PmP and SmS with strong amplitudes in large epicenter distance of 80–120 km. Further studies are still needed to improve the accuracy on travel time estimation, so as to further promote applicability of the CASS for and imaging the seismic velocity structure.     

Fractal travel time estimates for dispersive contaminants

Alternative fractional models of contaminant transport lead to a new travel time formula for arbitrary concentration levels. For an evolving contaminant plume in a highly heterogeneous aquifer, the new formula predicts much earlier arrival at low concentrations. Travel times of contaminant fronts and plumes are often obtained from Darcy's law calculations using estimates of average pore velocities. These estimates only provide information about the travel time of the average concentration (or peak, for contaminant pulses). Recently, it has been shown that finding the travel times of arbitrary concentration levels is a straightforward process, and equations were developed for other portions of the breakthrough curve for a nonreactive contaminant. In this paper, we generalize those equations to include alternative fractional models of contaminant transport.     

核幔界面衍射波Pdiff震相最大 记录震中距探讨

智利M_S8.8特大地震在甘肃数字测震台网记录的最大震中距达179.7°, 通过读取该地震初至震相的观测走时, 计算该震相的平均观测慢度为4.104 s/°, 与P_diff震相的理论慢度4.439 s/°基本一致. 研究结果表明: 由于甘肃数字测震台网和甘东南野外观测流动台阵独特的地理位置, 记录到了至今全球最大震中距(179.7°)的核幔界面衍射波P_diff震相; P_diff震相起始观测走时明显滞后于IASP91理论走时; 核幔界面复杂的形态对不同方位传播的P_diff震相观测走时产生了一定的影响.      

Viruses as Groundwater Tracers: Using Ecohydrology to Characterize Short Travel Times in Aquifers

Viruses are attractive tracers of short (<3 year) travel times in aquifers because they have unique genetic signatures, are detectable in trace quantities, and are mobile in groundwater. Virus “snaphots” result from infection and disappearance in a population over time; therefore, the virus snapshot shed in the fecal wastes of an infected population at a specific point in time can serve as a marker for tracking virus and groundwater movement. The virus tracing approach and an example application are described to illustrate their ability to characterize travel times in high‐groundwater velocity settings, and provide insight unavailable from standard hydrogeologic approaches. Although characterization of preferential flowpaths does not usually characterize the majority of other travel times occurring in the groundwater system (e.g., center of plume mass; tail of the breakthrough curve), virus approaches can trace very short times of transport, and thus can fill an important gap in our current hydrogeology toolbox.     

甘肃岷县漳县6.6级地震发生日期的触发因素分析及对短临预测方法的思考

短临预报是目前地震预测的难点,通过对2013年甘肃岷县漳县MS6.6地震发生日期的研究,发现发生日期7月22日是节气、低点位移、异年倍七律和朔望共同作用的时间点。据此,本文用一些基于可能触发地震发生的非传统方法的时间预测方法,回顾性地讨论了甘肃岷县漳县MS6.6地震发生日期的临震日期预测问题。通过对岷县漳县MS6.6地震研究,认为在地震短临日期的预测上。一些非传统方法虽然在机理上还有待研究,但在统计上多次证明该方法具有预测准确性的效能。要基于传统方法与非传统方法结合,才能对临震预测的时间逼近,同时提出今后预测临震时间的思路:多因素不断拦截法。