Parameters of the normal gravity field deduced from satellite observations

Summary The parameters of the normal gravity field were deduced from the harmonic coefficients[3, 4] upto n=6 and compared with the parameters used hitherto. The symbols used are the same as in papers[5, 6, 8] with which this paper connects up.     

On the wave propagation in inhomogeneous non-isotropic media

Summary FollowingEason, we have discussed here the propagation of elastic waves in non-homogeneous spheres and cylinders when the curved surface is given a uniform normal loading. The material is assumed to be transversely isotropic with respect to a direction of symmetry, the stress and displacement components within the body may be assumed to depend on one space co-ordinate and time alone. The particular case in which the elastic parameters are proportional to (radius) ⁿ has been considered.     

引黄水源水库平-丰水期有色可溶性有机物(CDOM)的来源、组成及差异分析

为探究引黄水源水库——门楼水库平水期和丰水期有色可溶性有机物(CDOM)的组成特征、来源及差异,运用紫外—可见光谱技术(UV-vis)和三维荧光光谱(EEMs)技术,结合平行因子分析法(PARAFAC)分析2022年5月(平水期)和2022年7月(丰水期)有色可溶性有机物含量及组分变化。研究结果表明:PARAFAC识别出2类荧光组分,分别是C1(Ex=355 nm,Em=476 nm,类腐殖质组分)和组分C2(Ex=225 nm,Em=320 nm,类蛋白组分);丰水期CDOM组分荧光强度显著高于平水期。CDOM光谱参数表明,门楼水库水体处于中营养状态,水体CDOM受新生内源和外源输入共同影响,以自生源为主;水库CDOM具有相对分子量小、腐殖化程度较弱的特点。丰水期水库富营养化水平和CDOM相对浓度低于平水期;丰水期CDOM疏水性组分比例和芳香化程度高于平水期。水质理化指标、CDOM组分和光谱参数相关性分析结果表明SUVA₂₆₀和SUVA₂₈₀与DOC呈显著负相关,说明紫外—可见光谱参数在一定程度上可以用来估算DOC的浓度;Chl.a浓度作为...     

NON-NORMAL INCIDENCE INVERSION: EXISTENCE OF SOLUTION*

An inverse problem is one in which the parameters of a model are determined from measured seismic data. Important to the solution of inverse problems is the issue of whether or not a solution exists. In this paper we show, in a constructive manner, that a solution does exist to the specific inverse problem of determining the parameters of a horizontally stratified, lossless, isotropic and homogeneous layered system that is excited by a non-normal incidence (NNI) plane wave. Mode conversion between P- and S-waves is included. We develop a seven-step layer-recursive procedure for determining all of the parameters for layer j. These parameters are P-wave and S-wave velocities and angles of incidence, density, thickness, traveltimes, and reflection- and transmission-coefficient matrices. Downward continuation of data from the top of one layer to the top of the next lower layer is an important step in our procedure, just as it is in normal incidence (NI) inversion. We show that, in order to compute all parameters of layer j, we need to (and can) compute some parameters for layer j+ 1. This is a non-causal phenomenon that seems to be necessary in NNI inversion but is not present in NI inversion.     

Effects of stress reduction on geomechanical and acoustic relationship of overconsolidated sands

Relationship between different geomechanical and acoustic properties measured from seven laboratory-tested unconsolidated natural sands with different mineralogical compositions and textures were presented. The samples were compacted in the uniaxial strain configuration from 0.5 to 30 MPa effective stress. Each sand sample was subjected to three loading–unloading cycles to study the influence of stress reduction. Geomechanical, elastic and acoustic parameters are different between normal compaction and overconsolidation (unloaded and reloaded). Stress path (K₀) data differ between normal consolidated and overconsolidated sediments. The K₀ value of approximately 0.5 is founded for most of the normal consolidated sands, but varies during unloading depending on mineral compositions and textural differences. The K₀ and overconsolidation ratio relation can be further simplified and can be influenced by the material compositions. K₀ can be used to estimate horizontal stress for drilling applications. The relationship between acoustic velocity and geomechanical is also found to be different between loading and unloading conditions. The static moduli of the overconsolidated sands are much higher than normal consolidated sands as the deformation is small (small strain). The correlation between dynamic and static elastic moduli is stronger for an overconsolidation stage than for a normal consolidation stage. The results of this study can contribute to geomechanical and acoustic dataset which can be applied for many seismic-geomechanics applications in shallow sands where mechanical compaction is the dominant mechanism.     

The zero-line method of interpreting total field magnetic anomalies of spherical ore-bodies

Summary A method called the zero-line method for interpreting the total field magnetic anomalies of spherical bodies is described. The contour of zero-anomaly is controlled by a linear equation of the typeAx ² +By ² +Cxy +Dx +Ey +F = 0, where the coefficientsA toF are related to the depth of the sphere and the parameters relating to its position and magnetisation. From the coordinates (x, y) of the various points lying on the zero-contour, five normal equations for the above equation are set and solved for the various coefficients. The various parameters relating to the position and magnetisation of the sphere are then determined from these coefficients. It is expected that similar methods can be developed for interpreting magnetic anomalies of other three dimensional bodies also.     

Variability of parameters of the <Emphasis Type="Italic">F</Emphasis>2-layer maximum in the quiet midlatitude ionosphere under low solar activity: 1. Statistical properties

Results of statistical analysis of the properties of variability of F2-layer maximum parameters (critical frequency foF2 and the height hmF2) in quiet midlatitude ionosphere under low solar activity in the daytime (1000–1500 LT) and nighttime (2200–0300 LT) hours are presented on the basis of Irkutsk station data for 2007–2008. It is found that the distribution density of δfoF2 could be presented as consisting of two distinctly different normal laws of this distribution, one of which corresponds to weak (|δfoF2| < 10%) fluctuations in foF2 and the other corresponds to strong (30% > |δfoF2| > 10%) fluctuations. Weak fluctuations in foF2 to a substantial degree are related to ionospheric variability at times less of than 1–3 h and determine the δfoF2 variability in the daytime hours. Strong fluctuations in foF2 are mainly related to day-to-day variability of the ionosphere at a fixed local time, the variability increasing by approximately a factor of 3 during the transition from day to night and determining the δfoF2 variability in the nighttime hours. The distribution density of ΔhmF2 is close to the normal distribution law. An interpretation of the different character of the distribution densities of δfoF2 and ΔhmF2 is given.     

Transect study on solute transport in a macroporous soil

Solute transport experiments using a non-reactive tracer were conducted on short, undisturbed, saturated columns of a sandy loam soil. All columns, 20 cm in diameter and 20 cm long, were collected along a transect of 35 m. Most of the soil columns had pre-existing macropores. The columns were leached at a steady flow-rate under ponding conditions. The resulting breakthrough curves (BTCs) showed a large heterogeneity. Several of the BTCs displayed early breakthrough and long tailing. All the data were interpreted in terms of dimensional time moments, the classical convection-dispersion equation (CDE) and the mobile-immobile transport model (MIM). Experimental time moments were found to vary significantly among the different BTCs. Analysis of the time moments also revealed that the variance of the field-scale BTC was several times larger than the average of the local-scale variance. The pore water velocity v and dispersion coefficient D were obtained by fitting the CDE to the local-scale BTCs, resulting in an average dispersivity of 7·4 cm. Frequency distributions for the CDE parameters v and D were equally well described by a normal or log-normal probability density function (pdf). When a log-normal pdf for D is considered, the variance of the log_e transformed D values (σ_{ln D}²) was found to be 2·1. For the MIM model, two additional parameters were fitted: the fraction of mobile water, θ_m/θ, and the first-order mass transfer coefficient, α. The MIM was more successful in describing the data than the CDE transport model. For the MIM model, the average dispersivity was about 2 cm. The MIM parameters v, D and θ_m/θ were best described by a log-normal pdf rather than a normal pdf. Only the parameter α was better described by a normal pdf. Mobile water fractions, θ_m/θ ranged from 0·01 to 0·98, with a mean of 0·43 (based on a log-normal pdf). When the CDE and MIM were applied to the data, the fitted pore water velocities, v, compared favourably with the effective pore water velocities, v_eff, obtained from moment analysis.     

Calibration of the Tibetan Plateau Using Regional Seismic Waveforms

We use the recordings from 51 earthquakes produced by a PASSCAL deployment in Tibet to develop a two-layer crustal model for the region. Starting with their ISC locations, we iteratively fit the P-arrival times to relocate the earthquakes and estimate mantle and crustal seismic parameters. An average crustal P velocity of 6.2–6.3 km/s is obtained for a crustal thickness of 65 km while the P velocity of the uppermost mantle is 8.1 km/s. The upper layer of the model is further fine-tuned by obtaining the best synthetic SH waveform match to an observed waveform for a well-located event. Green's functions from this model are then used to estimate the source parameters for those events using a grid search procedure. Average event relocation relative to the ISC locations, excluding two poorly located earthquakes, is 16 km. All but one earthquake are determined by the waveform inversion to be at depths between 5 and 15 km. This is 15 km shallower, on average, than depths reported by the ISC. The shallow seismicity cut-off depth and low crustal velocities suggest high temperatures in the lower crust. Thrust faulting source mechanisms dominate at the margins of the plateau. Within the plateau, at locations with surface elevations less than 5 km, source mechanisms are a mixture of strike-slip and thrust. Most events occurring in the high plateau where elevations are above 5 km show normal faulting. This indicates that a large portion of the plateau is under EW extension.     

A semi-empirical method for predicting hydrological drought magnitudes in the Canadian prairies

Abstract

A hydrological drought magnitude (M _T ) expressed in standardized terms is predicted on annual, monthly and weekly time scales for a sampling period of T years in streamflow data from the Canadian prairies. The drought episodes are considered to follow the Poisson law of probability and, when coupled with the gamma probability distribution function (pdf) of drought magnitude (M) in the extreme number theorem, culminate in a relationship capable of evaluating the expected value, E(M _T ). The parameters of the underlying pdf of M are determined based on the assumption that the drought intensity follows a truncated normal pdf. The E(M _T ) can be evaluated using only standard deviation (σ), lag-1 autocorrelation (ρ) of the standardized hydrological index (SHI) sequence, and a weighting parameter Φ (ranging from 0 to 1) to account for the extreme drought duration (L _T ), as well as the mean drought duration (L_m ), in a characteristic drought length (L_c ). The SHI is treated as standard normal variate, equivalent to the commonly-used standardized precipitation index. A closed-form relationship can be used for the estimation of first-order conditional probabilities, which can also be estimated from historical streamflow records. For all rivers, at the annual time scale, the value of Φ was found equal to 0.5, but it tends to vary (in the range 0 to 1) from river to river at monthly and weekly time scales. However, for a particular river, the Φ value was nearly constant at monthly and weekly time scales. The proposed method estimates E(M _T ) satisfactorily comparable to the observed counterpart. At the annual time scale, the assumption of a normal pdf for drought magnitude tends to yield results in close proximity to that of a gamma pdf. The M _T , when transformed into deficit-volume, can form a basis for designing water storage facilities and for planning water management strategies during drought periods.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Sharma, T.C. and Panu, U.S., 2013. A semi-empirical method for predicting hydrological drought magnitudes in the Canadian prairies. Hydrological Sciences Journal, 58 (3), 549–569.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

Anisotropic azimuthal surface-wave inversion for a vertically fractured and transverse isotropy medium: Theory and examples from a controlled field experiment

Fractures in elastic media add compliance to a rock in the direction normal to the fracture strike. Therefore, elastic wave velocities in a fractured rock will vary as a function of the energy propagation direction relative to the orientation of the aligned fracture set. Anisotropic Thomson–Haskell matrix Rayleigh-wave equations for a vertically transverse isotropic media can be used to model surface-wave dispersion along the principal axes of a vertically fractured and transversely isotropic medium. Furthermore, a workflow combining first-break analysis and azimuthal anisotropic Rayleigh-wave inversion can be used to estimate P-wave and S-wave velocities, Thomsen's ε, and Thomsen's δ along the principal axes of the orthorhombic symmetry. In this work, linear slip theory is used to map our inversion results to the equivalent vertically fractured and transversely isotropic medium coefficients. We carried out this inversion on a synthetic example and a field example. The synthetic data example results show that joint estimation of S-wave velocities with Thomsen's parameters ε and δ along normal and parallel to the vertical fracture set is reliable and, when mapped to the corresponding vertically fractured and transversely isotropic medium, provides insight into the fracture compliances. When the inversion was carried out on the field data, results indicated that the fractured rock is more compliant in the azimuth normal to the visible fracture set orientation and that the in situ normal fracture compliance to tangential fracture compliance ratio is less than half, which implies some cementation may have occurred along the fractures. Such an observation has significant implications when modelling the transport properties of the rock and its strength. Both synthetic and field examples show the potential of azimuthal anisotropic Rayleigh-wave inversion as the method can be further expanded to a more general case where the vertical fracture set orientation is not known a priori.     

Verification of the TSMP-assisted digisonde topside profiling technique

The purpose of this paper is to evaluate the performance of the TSMP-assisted Digisonde (TaD) topside profiling technique. We present systematic comparisons between electron density profiles and TEC parameters extracted from TaD model with (a) CHAMP-derived TEC parameters, (b) CHAMP reconstructed profiles, (c) ground based GPSderived TEC parameters, and (d) profiles reconstructed from RPI/IMAGE plasmagrams. In all cases, TaD follows the general trend of plasmaspheric observations derived from the above datasets. Especially during storm cases, TaD shows remarkable agreement with the variations of the ground based GPS-derived TEC parameters. Overall, the comparison results shows that TaD method can be adopted by EURIPOS to provide the electron density distribution up to plasmaspheric heights in real-time.     

Linear combinations of order statistics to estimate the position and scale parameters of the Generalized Extreme Values distribution

 Ad hoc techniques for estimating the position and the scale parameters of the Generalized Extreme Values distribution are introduced. The estimators proposed are simple linear combinations of the order statistics: they provide valuable estimates of the parameters of interest both when the shape parameter is known and when it is unknown (this latter case being of great relevance in practical applications), and show a good performance as well when the sample size is small. In addition, weakly-consistent estimators of the position parameter are introduced, whose calculation does not require the knowledge of any of the remaining parameters. The procedures are tested on simulated data, and comparisons with other techniques are shown.     

The constant part of the tidal field in the theory of heights

Summary The problem of the constant part of the tidal field is still topical in view of the recommendations of IAG[1, 2] to eliminate the tidal effect of external masses from all geodetic measurements under preservation of the effect of the time-constant tidal deformation of the Earth. The paper discusses the consequences of accepting this recommendation for normal heights, and suggests a solution based on the new definition of the normal gravity field[3].     

Vertically fractured transversely isotropic media: dimensionality and deconstruction

A vertically fractured transversely isotropic (VFTI) elastic medium is one in which any number of sets of vertical aligned fractures (each set has its normal lying in the horizontal x₁, x₂‐plane) pervade the medium and the sets of aligned fractures are the only features of the medium disturbing the axi‐symmetry about the x₃‐axis implying that in the absence of fractures, the background medium is transversely isotropic (TI). Under the assumptions of long wavelength equivalent medium theory, the compliance matrix of a fractured medium is the sum of the background medium's compliance matrix and a fracture compliance matrix. For sets of parallel rotationally symmetric fractures (on average), the fracture compliance matrix is dependent on 3 parameters − its normal and tangential compliance and its strike direction. When one fracture set is present, the medium is orthorhombic and the analysis is straightforward. When two (non‐orthogonal) or more sets are present, the overall medium is in general elastically monoclinic; its compliance tensor components are subject to two equalities yielding an 11 parameter monoclinic medium. Constructing a monoclinic VFTI medium with n embedded vertical fracture sets, requires 5 TI parameters plus 3×n fracture set parameters. A deconstruction of such an 11 parameter monoclinic medium involves using its compliance tensor to find a background transversely isotropic medium and several sets of vertical fractures which, in the long wavelength limit, will behave exactly as the original 11 parameter monoclinic medium. A minimal deconstruction, would be to determine, from the 11 independent components, the transversely isotropic background (5 parameters) and two fracture sets (specified by 2 × 3 = 6 parameters). Two of the background TI medium's compliance matrix components are known immediately by inspection, leaving nine monoclinic components to be used in the minimal deconstruction of the VFTI medium. The use of the properties of a TI medium, which are linear relations on its compliance components, allows the deconstruction to be reduced to solving a pair of non‐linear equations on the orientations of two fracture sets. A single root yielding a physically meaningful minimum deconstruction yields a unique minimal representation of the monoclinic medium as a VFTI medium. When no such root exists, deconstruction requires an additional fracture set and uniqueness is lost. The boundary between those monoclinic media that have a unique minimal representation and those that do not is yet to be determined.     

Sensitivity analysis of a no-crossflow model for the transient flowmeter test

 Logarithmic sensitivities and plausible relative errors are studied in a simple no-crossflow model of a transient flowmeter test (TFMT). This model is identical to the model of a constant-rate pumping test conducted on a fully penetrating well with wellbore storage, surrounded by a thick skin zone, and situated in a homogeneous confined aquifer. The sensitivities of wellbore drawdown and wellface flowrate to aquifer and skin parameters are independent of the pumping rate. However, the plausible relative errors in the aquifer and skin parameters estimated from drawdown and wellface flowrate data can be proportionally decreased by increasing the pumping rate. The plausible relative errors vary by many orders of magnitude from the beginning of the TFMT. The practically important flowrate and drawdown measurements in this test, for which the plausible relative errors vary by less than one order of magnitude from the minimum plausible relative errors, can begin approximately when the dimensionless wellface flowrate exceeds q _D =q/Q≈0.4. During most of this stage of the test, the plausible relative errors in aquifer hydraulic conductivity (K _a ) are generally an order of magnitude smaller than those in aquifer specific storativity. The plausible relative errors in the skin hydraulic conductivity (K _s ) are generally larger than the plausible relative errors in the aquifer specific storativity when the thick skin is normal (K _s >K _a ) and smaller when the thick skin is damaged (K _s <K _a ). The specific storativity of the skin zone would be so biased that one should not even attempt to estimate it from the TFMT. We acknowledge Wiebe H. van der Molen for recommending the De Hoog algorithm and sharing his code. This research was partially supported by the US Geological Survey, USGS Agreement #1434-HQ-96-GR-02689 and North Carolina Water Resources Research Institute, WRRI Project #70165.     

Testing geomorphic signal of active normal faulting: The case of the Cittanova Fault (Calabria,southern Italy)

Extracting tectonic signals from the landscape is an important challenge for constraining the style and rate of deformation associated with active faults, especially where their displacement history cannot be independently determined. Based on previous paleoseismological data coupled with new geomorphological field work and ¹⁴C dating of geomorphic markers, we analysed the geomorphic signal of the along‐strike differential throw of the Cittanova Fault in southern Calabria (Italy), the recent activity of which is already well documented and constrained. Through DTM‐derived stream power law parameters (SL and χ), we provide evidence of drainage network disequilibrium and reorganization in response to fault growth and deformation style. Furthermore, a methodological test of the reliability of the χ metric as a proxy for the differential throw along the strike of active normal faults provided good preliminary results, consistent with a strong inverse linear correlation with fault throw. Copyright © 2018 John Wiley & Sons, Ltd.     

On the grain size dependence of the behavicur of fine magnetic particles in rocks

The grain size dependence of the ratio of saturation remanent magnetization to saturation magnetization (J _R :J _S ), weak field susceptibility (X ₀), thermoremanent magnetization (TRM) and its stability against AF demagnetization are interpreted in terms of nucleation theory. It is concluded that each of these parameters exhibits grain size dependence due to two effects. The first is the increasing difficulty with which domain walls are neucleated as grain size decreases. The second is an intrinsic grain size dependence of the parameters in multidomain particles.