北斗GEO卫星载波相位单差观测值的周期性分析

高精度的卫星定位测量中,由于卫星在测站上空的重复周期性运动,周期性误差是影响测量精度的重要误差源。恒星日滤波利用周期性误差的特点,提取前一天的周期性误差来改正后一天的定位结果,可滤波处理诸如多路径误差。本文从北斗GEO卫星载波相位观测值站间单差的角度来探讨其周期性特点,通过前后两天站间相位单差观测值进行对比,采用相关性分析方法,验证北斗GEO卫星载波相位单差观测值的周期性特征。     

Using micro-catchment experiments for multi-local scale modelling of nature-based solutions

The Q-natural flood management project has co-developed with the Environment Agency 18 monitored micro-catchments (~1 km²) in Cumbria, UK installing calibrated flumes aimed at quantifying the potential shift in observed flows resulting from a range of nature-based-solutions installed by local organizations. The small-scale reduces the influence of variability characterizing larger catchments that would otherwise mask any such shifts, which we attempt to relate to a shift in model parameters. This paper demonstrates an approach to applying donor-parameter-shifts obtained from modelling two of the paired micro-catchments to a much larger scale, in order to understand the potential for improved distributed modelling of nature-based solutions in the form of additional tree-planting. The models include a rainfall-runoff model, Dynamic Topmodel, and a 2D hydrodynamic model, JFlow, permitting analysis of changes in hillslope processes and channel hydrodynamics resulting from a range of distributed measures designed to emulate natural hydrological processes that evaporate, store or infiltrate flows. We report on attempts to detect shift in hydrological response using one of the paired-micro-catchment moorland versus forestry sites in Lorton using Dynamic Topmodel. A donor-parameter-shift approach is used in a hypothetical experiment to represent new woodland in a much larger catchment, although testing all combinations of spatial planting strategies, responses to multiple-extremes, failure-modes and changes to synchronization becomes intractable to support good decision making. We argue that the problem can be re-framed to use donor-parameter-shifts at multi-local-scale catchments above communities known to be at risk, commensurate with most of the evidence of NbS impacts being effective at the small scale (ca. 10 km²). This might lead to more effective modelling to help catchment managers prioritize those communities-at-risk where there is more evidence that NbS might be effective.     

Source Parameters of the Deadly M<Subscript>w</Subscript> 7.6 Kashmir Earthquake of 8 October, 2005

During the last six years, National Geophysical Research Institute, Hyderabad has established a semi-permanent seismological network of 5–8 broadband seismographs and 10–20 accelerographs in the Kachchh seismic zone, Gujarat with a prime objective to monitor the continued aftershock activity of the 2001 M_w 7.7 Bhuj mainshock. The reliable and accurate broadband data for the 8 October M_w 7.6 2005 Kashmir earthquake and its aftershocks from this network as well as Hyderabad Geoscope station enabled us to estimate the group velocity dispersion characteristics and one-dimensional regional shear velocity structure of the Peninsular India. Firstly, we measure Rayleigh-and Love-wave group velocity dispersion curves in the period range of 8 to 35 sec and invert these curves to estimate the crustal and upper mantle structure below the western part of Peninsular India. Our best model suggests a two-layered crust: The upper crust is 13.8 km thick with a shear velocity (Vs) of 3.2 km/s; the corresponding values for the lower crust are 24.9 km and 3.7 km/sec. The shear velocity for the upper mantle is found to be 4.65 km/sec. Based on this structure, we perform a moment tensor (MT) inversion of the bandpass (0.05–0.02 Hz) filtered seismograms of the Kashmir earthquake. The best fit is obtained for a source located at a depth of 30 km, with a seismic moment, Mo, of 1.6 × 10²⁷ dyne-cm, and a focal mechanism with strike 19.5°, dip 42°, and rake 167°. The long-period magnitude (M_A ~ M_w) of this earthquake is estimated to be 7.31. An analysis of well-developed sPn and sSn regional crustal phases from the bandpassed (0.02–0.25 Hz) seismograms of this earthquake at four stations in Kachchh suggests a focal depth of 30.8 km.     

A Least-squares Window Curves Method for Interpretation of Magnetic Anomalies Caused by Dipping Dikes

We have developed a least-squares method to determine simultaneously the depth and the width of a buried thick dipping dike from residualized magnetic data using filters of successive window lengths. The method involves using a relationship between the depth and the half-width of the source and a combination of windowed observations. The relationship represents a family of curves (window curves). For a fixed window length, the depth is determined for each half-width value by solving one nonlinear equation of the form f (z) = 0 using the least-squares method. The computed depths are plotted against the width values representing a continuous curve. The solution for the depth and the width of the buried dike is read at the common intersection of the window curves. The method involves using a dike model convolved with the same moving average filter as applied to the observed data. As a result, this method can be applied to residuals as well as to measured magnetic data. Procedures are also formulated to estimate the amplitude coefficient and the index parameter. The method is applied to theoretical data with and without random errors. The validity of the method is tested on airborne magnetic data from Canada and on a vertical component magnetic anomaly from Turkey. In all cases examined, the model parameters obtained are in good agreement with the actual ones and with those given in the published literature.     

大姚地热动态特征分析

通过大姚高精度水温观测井水地质、井孔条件的分析，震例观测资料、实验观测数据的处理，以及井水水质分析结果，表明大姚水温周期性变化的突出位置是井深75m左右，根本原因是该地层石膏脉稳定释热，含水层热水与井内冷水混合过程中，在通道内形成钙盐类沉淀物（如CaSO4CaCO3），堵塞对流通道，含水层内外压差突破被阻塞通道，热水反复侵入井内，便形成水温周期性变化。地震的孕育、发生和调整过程，由于应力场作用，产生附加地热场，水温基值发生变化，周期畸变或消失，这种附加地热动态是可以恢复的。所以认为大姚地热动态是井区特殊地层地热信息的反映。     

ANALYTICAL HAYAMI SOLUTION FOR THE DIFFUSIVE WAVE FLOOD ROUTING PROBLEM WITH LATERAL INFLOW

The diffusive wave equation is generally used in flood routing in rivers. The two parameters of the equation, celerity and diffusivity, are usually taken as functions of the discharge. If these two parameters can be assumed to be constant without lateral inflow, the diffusive wave equation may have an analytical solution: the Hayami model. A general analytical method, based on ‘Hayami’s hypothesis, is developed here which resolves the diffusive wave flood routing equation with lateral inflow or outflow uniformly distributed over a channel reach. Flood routing parameters are then identified using observed inflow and outflow and the Hayami model used to simulate outflow. Two examples are discussed. Firstly, the prediction of the hydrograph at a downstream section on the basis of a knowledge of the hydrograph at an upstream section and the lateral inflow. The second example concerns lateral inflow identification between an upstream and a downstream section on the basis of a knowledge of hydrographs at the upstream and downstream sections. The new general Hayami model was applied to flood routing simulation and for lateral inflow identification of the River Allier in France. The major advantages of the method relate to computer simulation, real-time forecasting and control applications in examples where numerical instabilities, in the solution of the partial differential equations must be avoided.     

Analytical model of well leakage pressure perturbations in a closed aquifer system

Deep saline aquifers are commonly used for disposal and storage of various surface fluids. The target injection zone must be hydraulically isolated from overlying zones in order to ensure containment of the injected fluids. Improperly plugged nonoperational abandoned wells that penetrate the injection zone are the main potential leakage pathways. Leakage through such wells may cause an observable pressure signal in a zone overlying the injection zone; such a signal can be used to detect the leakage. In this paper we develop an analytical model to evaluate the pressure change induced by leakage through a well in a multilayer system. Unlike previous analytical models on the topic, our model uses a closed system, which may significantly affect the strength and behavior of the pressure signal induced by leakage. The analytical model is first presented for a two-layer system centered at the leaky well location. We evaluate the leakage-induced pressure change using the Laplace transform of Duhamel’s superposition integral, yielding the solution in the Laplace domain. We then derive a late-time asymptotic solution using the final value theorem, which suggests that the leakage rate becomes constant after sufficient time. We then obtain the multilayer solution by extending the two-layer solution and presenting it in matrix form in the Laplace domain. We apply the solution to three examples. In the first example, we apply the analytical model to a two-layer system, investigating its behavior and comparing the results with a numerical solution. In order to demonstrate behavior and potential applications of the multilayer analytical model, we present two multilayer examples: one with identical layers and another, replicating a CO₂ storage site, with dissimilar layers. The leakage-induced pressure change does not necessarily decrease as the distance increases from the injection zone toward the surface.     

CALCULATING THE DYNAMIC STIFFNESS MATRIX OF 2-D FOUNDATIONS BY DISCRETE WAVE NUMBER INDIRECT BOUNDARY ELEMENT METHODS

Apart from some special cases, calculating the dynamic stiffness matrix of foundations on a layered half-space, especially for embedded foundations, is computationally expensive. An efficient method for two-dimensional foundations in a horizontally layered soil media is presented in this paper. This method is based on indirect boundary element methods and uses discrete wave number solution methods for calculating Green's functions for displacements and analytical methods for the integrations over the boundary. For surface foundations, the present method applies at all frequencies. For embedded foundations or for constructing energy transmitting boundaries, because the free-field part is modelled by boundary elements and the excavated part is modelled by finite elements, the present method applies only at low frequencies for the spring coefficients (the real parts of the dynamic stiffness matrix) but applies at all frequencies for the damping coefficients (the imaginary part of the dynamic stiffness matrix) for undamped sites. The novelty of the method can be used for three-dimensional foundations. © 1997 by John Wiley & Sons, Ltd.     

On the frontal condition for finite amplitude α2Ω-dynamo wave trains in stellar shells

Abstract

In a previous paper, Bassom et al. (Proc. R. Soc. Lond. A, 455, 1443–1481, 1999) (BKS) investigated finite amplitude αΩ-dynamo wave trains in a thin turbulent, differentially rotating convective stellar shell; nonlinearity arose from α-quenching. There asymptotic solutions were developed based upon the small aspect ratio ε of the shell. Specifically, as a consequence of a prescribed latitudinally dependent α-effect and zonal shear flow, the wave trains have smooth amplitude modulation but are terminated abruptly across a front at some high latitude θ_F. Generally, the linear WKB-solution ahead of the front is characterised by the vanishing of the complex group velocity at a nearby point θ_f; this is essentially the Dee–Langer criterion, which determines both the wave frequency and front location.

Recently, Griffiths et al. (Geophys. Astrophys. Fluid Dynam. 94, 85–133, 2001) (GBSK) obtained solutions to the α²Ω-extension of the model by application of the Dee—Langer criterion. Its justification depends on the linear solution in a narrow layer ahead of the front on the short O(θ_f—θ_F) length scale; here conventional WKB-theory, used to describe the solution elsewhere, is inadequate because of mode coalescence. This becomes a highly sensitive issue, when considering the transition from the linear solution, which occurs when the dynamo number D takes its critical value D _c corresponding to the onset of kinematic dynamo action, to the fully nonlinear solutions, for which the Dee—Langer criterion pertains.

In this paper we investigate the nature of the narrow layer for α²Ω-dynamos in the limit of relatively small but finite α-effect Reynolds numbers R _α, explicitly ε^½ ? R ² _α ? 1. Though there is a multiplicity of solutions, our results show that the space occupied by the corresponding wave train is generally maximised by a solution with θ_f—θ_F small; such solutions are preferred as evinced by numerical simulations. This feature justifies the application by GBSK of the Dee—Langer criterion for all D down to the minimum D _min that the condition admits. Significantly, the frontal solutions are subcritical in the sense that |D _min| ≤ |D _c|; equality occurs as the α-effect Reynolds number tends to zero. We demonstrate that the critical linear solution is not connected by any parameter track to the preferred nonlinear solution associated with D _min. By implication, a complicated bifurcation sequence is required to make the connection between the linear and nonlinear states. This feature is in stark contrast to the corresponding results for αΩ-dynamos obtained by BKS valid in the limit R _{2 α} ? ε^½, which, though exhibiting a weak subcriticality, showed that the connection follows a clearly identifiable nonbifurcating track.     

2013年洱源M_S5.5地震序列M_L≥3.4地震的震源机制解研究

利用云南省地震台网的宽频带数据,采用CAP方法计算了2013年3月3日洱源MS5.5地震前后共10次ML≥3.4地震的震源机制解,分析了洱源MS5.5地震前后震源机制解的一致性变化特征,结合精定位结果、震后现场灾害评估报告及震区地质构造情况,判定维西—巍山断裂以西的苍山西麓炼铁盆地东缘主边界断裂带为洱源MS5.5地震的发震构造,震源机制解节面Ⅱ(走向160°、倾角48°、滑动解-99°)为其发震断层面。