A Note On Surface Humidity Measurements In The Cold Canadian Environment

In the cold Canadian environment, humidity measurements can be verydifficult to conduct. In this brief communication, humidity observationstaken by two different sensors at six remote Canadian Arctic locations arecompared. The observations collected by Vaisala HMP35CF sensors displaya strong tendency toward the ice saturation point whereas dew cell instrumentsexhibit significantly lower values of relative humidity with respect to ice (RH_i).Humidity data collected by HMP35CF hygrometers are therefore unreliable sincethey are subject to persistent icing that lead them to record values of RH_i near 100%, irrespective of the air temperature. The high humidity bias emerges atmeteorological stations mounted with the HMP35CF probes since these instrumentsare usually neither sheltered nor heated, and are not attended to at regular intervals. Thus, great care must be taken in utilizing humidity data recorded by HMP35CF sensors across the network of climate autostations in Canada.     

离散裂隙渗流方法与裂隙化渗透介质建模

流体渗流模拟的连续介质方法通常适用于多孔地质体，并不一定适用于裂隙岩体，由于裂隙分布及其特征与孔隙差异较大。若流体渗流主要受裂隙的控制，对于一定尺寸的裂隙岩体，多孔介质假设则较难刻划裂隙岩体的渗流特征。离散裂隙渗流方法不但可直接用于模拟裂隙岩体非均质性和各向异性等渗流特征，而且可用其确定所研究的裂隙岩体典型单元体及其水力传导（渗透）张量大小。主要讨论了以下问题：（1）饱和裂隙介质中一般的离散流体渗流模拟；（2）裂隙岩体中的REV（典型单元体）及其水力传导（渗透）张量的确定；（3）利用离散裂隙网络流体渗流模型研究裂隙方向几何参数对水力传导系数和REV的影响；（4）在二维和三维离散裂隙流体渗流模型中对区域大裂隙和局部小裂隙的处理方法。调查结果显示离散裂隙流体渗流数学模型可用来评价不同尺度上的裂隙岩体的水力特征，以及裂隙方向对裂隙化岩体的水力特征有着不可忽视的影响。同时，局部小裂隙、区域大裂隙应当区别对待，以便据其所起的作用及水力特征，建立裂隙化岩体相应的流体渗流模型。     

Stability of earthquake moment tensor inversions: effect of the double-couple constraint

We show that spurious large non-double-couple components can be obtained in inversions for the full deviatoric moment tensor for shallow crustal earthquakes due to inaccurate Earth models. The traditional “best double-couple” solution does not in general provide an optimal estimate of a double-couple mechanism, and is only reliable when the non-double-couple component of the full deviatoric solution is small. The inverse problem for the moment tensors of the 1998 Antarctic Plate and 2000 Wharton Basin strike-slip earthquakes is shown in each case to have two well-fitting minima in the misfit function of pure double-couple solutions. Such pairs of solutions are most likely to exist for earthquakes which are close either to vertical strike-slip or to dip-slip on a fault plane dipping at 45°. It is shown theoretically that these pairs of solutions arise from the combination of the pure double-couple constraint and the instability of two elements of the moment tensor. No significant non-double-couple component is found for the shallow thrusting 1996 Biak, Indonesia earthquake.     

Mechanism of the June 4, 2000 Southern Sumatra, Indonesia, earthquake

Introduction The June 4, 2000 southern Sumatra, Indonesia, earthquake (16h28min26.2s UTC, 4.72, 102.09, 33 km, MS=8.0 [IRIS]) occurred under the Indian Ocean, near the Mentawai fault, along the well-known Sumatran subduction zone and the great Sumatran fault, all of which trend northwest-southeast. 1 800 houses were totally destroyed, 10 196 were heavily damaged, and 18 378 were slightly damaged by the earthquake. At least 97 people were killed, 1 900 were injured, and 122 000 were left…     

Moment-tensor statistics of the 1997 Long Valley microearthquake swarm

Abstract. A simple, fast, moment-tensor inversion method using bandpass-filtered P-amplitudes was used to study the moment-tensor statistics of Long Valley caldera microearthquakes. The events were recorded in the summer of 1997, during a swarm in the caldera. The swarm was associated with geodetic extension, uplift, and subsequent moderate earthquake activity. The moment tensor solutions for 1,993 events were calculated using the new method. The majority of the resulting focal mechanisms appear to be explained in terms of double couple mechanisms. Since some events did exhibit considerable deviation from double-couples, the moment data were studied for their statistical significance. The moments of the actual data were compared to the moments of synthetic data with varying degrees of random noise in their spectra. The results of this study suggested that unless data from more than 20 stations are used and the earthquake epicenter is located inside or very close to the network area, moment-tensor inversion does not correctly resolve the non-double-couple components of microearthquakes. Analysis of the inversion residuals shows that the average noise in the P-wave spectra was close to 20%. The fluctuations of the volumetric components of the moment-tensor are in good agreement with those of the synthetic pure double-couples with 20% of added noise. Thus the moment-tensor statistics suggests that little if any volume change is required to explain the observed seismic energy release in the swarm. However, the statistics do show that a significant compensated-linear-vector-dipole component maybe present in the bulk of the seismicity. Given the network used in this study, such a component could not be precisely resolved for individual earthquakes. This possibility deserves further investigation because of its bearing on the nature of fluid-fault-earthquake processes in swarms.     

Relation of the Wave-Propagation Metric Tensor to the Curvatures of the Slowness and Ray-Velocity Surfaces

The contravariant components of the wave-propagation metric tensor equal half the second-order partial derivatives of the selected eigenvalue of the Christoffel matrix with respect to the slowness-vector components. The relations of the wave-propagation metric tensor to the curvature matrix and Gaussian curvature of the slowness surface and to the curvature matrix and Gaussian curvature of the ray-velocity surface are demonstrated with the help of ray-centred coordinates.     

Inferring erosional resistance of bedrock units in the east Tennessee mountains from digital elevation data

Measures of local relief, regional relief, and slope were calculated from digital elevation models (DEMs) for 50 bedrock units in the Ridge and Valley and Blue Ridge provinces of Tennessee. Each of these measures was normalized and the three were then averaged to produce the erosional resistance index (ERI). Bedrock units with higher ERI values include coarse clastics, intermediate clastics, and metaplutonics. Units with lower values include shales, limestones, limestones plus dolostones, and carbonates plus fine clastics. Dolostones tend to have intermediate values. The calculated ERI values were compared with subjective ratings by a geologist with decades of field experience in east Tennessee. Generally, the agreement between the two ratings was good, the most glaring exception being several shales with improbably high ERI values. These turned out to be thin units cropping out beneath very hard sandstones, allowing them to stand higher and steeper than would otherwise be possible. A systematic method for detecting such erroneously high ERI values is suggested. Inspection of a drainage map superimposed on the geology map shows that in a given area, streams tend to flow on rock units with the lowest ERI values. In addition, statistical analysis shows that bedrock units with the lowest ERI values are, on average, almost three times closer to the nearest stream and six times as likely to have streams flowing on them than are units with highest values. Further, the effect of ERI on stream location is strongest for streams with drainage areas between 1 and 30 km². Thus, small streams appear to be subject to greater lithologic control than are larger streams.     

An Mw based earthquake Catalogue for central, northern and northwestern Europe using a hierarchy of magnitude conversions

Data from 25 local catalogues and 30special studies of earthquakes in central,northern and northwestern Europe have beenincorporated into a Databank. The dataprocessing includes discriminating eventtypes, eliminating fake events and dupletsand converting different magnitudes andintensities to M_w if this is not givenby the original source. The magnitudeconversion is a key task of the study andimplies establishment of regressionequations where no local relations exist.The Catalogue contains tectonic events fromthe Databank within the area44°N–72°N,25°W–32°E and the time period1300–1993. The lower magnitude level forthe Catalogue entries is setat M_w == 3.50. The area covered by thedifferent catalogues are associated withpolygons. Within each polygon only datafrom one or a small number of the localcatalogues, supplemented by data fromspecial studies, enter the Catalogue. Ifthere are two or more such catalogues orstudies providing a solution for an event,a priority algorithm selects one entry forthe Catalogue. Then M_w is calculatedfrom one of the magnitude types, or frommacroseismic data, given by the selectedentry according to another priority scheme.The origin time, location, M_w magnitude and reference are specified for eachentry of the Catalogue. So is theepicentral intensity, I₀, if providedby the original source. Following thesecriteria, a total of about 5,000earthquakes constitute the Catalogue.Although originally derived for the purposeof seismic hazard calculation within GSHAP,the Catalogue provides a data base for manytypes of seismicity and seismic hazardstudies.