A survey of Earth strain tides in Great Britain

Summary . We report 35 measurements of Earth strain tides at 16 near-surface sites in Great Britain. This is the first widespread survey undertaken specifically to examine the problem of the inhomogeneity of elastic strain fields near the Earth's surface. Some sites were instrumented intensively in order to examine variations of tidal admittance over distances of several hundred metres, while measurements from single instruments at other sites were compared with theoretically predicted strain tides. After allowing for cavity and topographic effects, our data show variations of up to 50 per cent in tidal admittance. We interpret such large anomalies as being due to variations in the regional elastic parameters of wavelength about a 100 m or less. The data indicate that strain measurements from single instrument sites must be interpreted with caution.     

Coseismic rotation changes from the 2004 Sumatra earthquake: the effects of Earth's compressibility versus earthquake induced topography

Polar motion is modelled for the large 2004 Sumatra earthquake via dislocation theory for an incompressible elastic earth model, where inertia perturbations are due to earthquake-triggered topography of density–contrast interfaces, and for a compressible model, where inertia perturbation due to compression-dilatation of Earth's material is included; density and elastic parameters are based on a multilayered reference Earth. Both models are based on analytical Green's functions, propagated from the centre to the Earth's surface. Preliminary and updated seismological solutions are considered in elucidating the effects of improving earthquake parameters on polar motion. The large Sumatra thrust earthquake was particularly efficient in driving polar motion since it was responsible for large material displacements occurring orthogonally to the strike of the earthquake and to the Earth's surface, as imaged by GRACE gravity anomalies over the earthquake area. The effects of earthquake-induced topography are four times larger than the effects of Earth's compressibility, for l = 2 geopotential components. For varying compressional Earth properties and seismic solution, modelled polar motion ranges from 8.6 to 9.4 cm in amplitude and between 117° and 130° east longitude in direction. The close relationship between polar motion direction, earthquake longitude and thrust nature of the event, are established in terms of basic physical concepts.     

Method for UAV-based 3D topography reconstruction of tidal creeks

It is common to obtain the topography of tidal flats by the Unmanned Aerial Vehicle(UAV) photogrammetry,but this method is not applicable in tidal creeks.The residual water will lead to inaccurate depth inversion results,and the topography of tidal creeks mainly de-pends on manual survey.The present study took the tidal creek of Chuandong port in Jiangsu Province,China,as the research area and used UAV oblique photogrammetry to reconstruct the topography of the exposed part above the water after the ebb tide.It also proposed a Trend Prediction Fitting (TPF) method for the topography of the unexposed part below the water to obtain a complete 3D topography.The topography above the water measured by UAV has the vertical precision of 12 cm.When the TPF method is used,the cross-section should be perpendicular the central axis of the tidal creek.A polynomial function can be adapted to most shape of sections,while a Fourier function obtains better results in asym-metrical sections.Compared with the two-order function,the three-order function lends itself to more complex sections.Generally,the TPF method is more suitable for small,straight tidal creeks with clear texture and no vegetation cover.     

On the excitation of resonant double Kelvin waves in the Barents Sea Opening

In the northern Barents Sea Opening (BSO) the K₁ tidal energy is predominant in the diurnal tidal frequency band, suggesting the generation of a topographic wave with the K₁ tidal frequency. Tidal energy of the K₁ component becomes strong where bottom topography undulates in the BSO and the scale of the undulation is close to the wavelength of the K₁ wave. An analytical model is developed to investigate the energy enhancement mechanism of the tidally induced topographic wave due to a resonance between tidal current, a topographic wave and periodic topography. The wave excited by the resonance is identified as a resonant double Kelvin wave (DKW) and the significant K₁ energy in the BSO could be due to the excitation of the resonant DKW.     

An application-driven approach to terrain model construction

Terrain is a surface phenomenon that is measured, modelled, and mapped. However, it is continuously variable and must be simulated by points or mathematical equations that are inherently approximations. The error induced by digitally represented terrain can propagate to surface derivatives and geographical information science (GIS) applications where topography is considered. This can lead to uncertainty in model predictions and the use of data that are unfit for the application to which they are intended. This article outlines the problem of uncertainty in terrain representation and demonstrates the consequences for volcanic mudflow modelling. The response of a simple least-cost single flow algorithm to input parameters was investigated in order to assess output variation from the different sources of input variation. Elevation error was modelled with a probability density function (PDF) and propagated through stochastic simulation (Monte Carlo). Such combined uncertainty and sensitivity analyses enabled a qualitative judgement of the relative significance of elevation error on the flow model prediction. Different methods for terrain model construction were considered and show that supplementing global positioning system (GPS) measurements with information from field notes and reconnaissance photographs greatly improved the model performance and reduced the uncertainty. It is concluded that in terms of validity of model results, there is no substitute for constructing an elevation model that is informed by the terrain.     

Acoustic-Gravity Wave Propagation in a Measured Atmosphere

Summary An atmosphere is modelled by averaging measurements taken at eight locations for 1968 September 9. The effectiveness of the modelling is evaluated by comparing the atmospheric dispersion properties determined from an experimentally recorded signal resulting from a nuclear explosion on that date and those computed from the model atmosphere. The effects of atmospheric winds are shown to change the dispersion properties.
A theoretical barogram is synthesized for a receiver located 7930 km from the source and is compared to the recorded signal for the event. Barograms are synthesized for two different source functions and the source functions are discussed.     

Tidal gravity and ocean tide loading in Europe

b
The results are presented from tidal gravity measurements at five sites in Europe using LaCoste and Romberg ET gravimeters. Improvements that we have made to the accuracies of these gravimeters are discussed. It is shown that the 'standard' calibration of the International Center for Earth Tides, used for worldwide tidal gravity profiles, is 1.2 per cent too high. The M₂ and O₁ observations are compared with model calculations of the Earth's body tide and ocean tide loading and it is shown that there is a very significant improvement in the agreement between observations and models compared to that obtained with previous tidal gravity measurements. For O₁, where the ocean tide loading and attraction in central Europe is only 0.4 per cent of the body tide, our measurements verify that the Dehant-Wahr anelastic body tide model gravimetric factor is accurate to 0.2 per cent. It is also shown that the effects of lateral heterogeneities in Earth structure on tidal gravity are too small to explain the large anomalies in previously published tidal gravity amplitudes. The observations clearly show the importance of conserving tidal mass in the Schwiderski ocean tide model. For sites in central Europe, the M₂ and O₁ observations and the models are in agreement at the 0.1 μgal (10⁻⁹ m s⁻²) level and tidal corrections to this accuracy can now be made to absolute gravity measurements.     

流动沙丘涡度相关通量的坐标旋转订正与对比研究

复杂下垫面导致近地面大气运动具有特殊的动力学和热力学特征,进而影响物质和能量的交换过程。以科尔沁流动沙丘7月份涡度数据为例,分析由于地形引起的风速、旋转角度β的变化,采用二次旋转法（DR）和平面拟合法（PF）对水热通量、数据质量对比分析,结果表明：（1）地形平缓区（180°~330°）垂直风速>0,β在1.677 3°左右;地形陡峭区（330°~180°）垂直风速<0,β在-1.648 6°左右。（2）DR法和PF法对潜热通量影响很大,分别增加17%和减小15%,DR法订正前后水热通量相关性均高于PF法。（3）DR法和PF法订正后的数据分别有89.6%和89.3%通过湍流发展的充分性检验;DR法订正前后摩擦风速相关性高于PF法;能量闭合率分别提高了3.47%和2.78%。（4）整个生长季中,6、8月订正后涡度数据质量有较好的提高。因此,DR法更适用于起伏较大的流动沙丘。研究结果充分考虑地形对涡度数据影响,为科尔沁沙地及具有相似地形的地区涡度数据的合理使用提供参考。     

Tidal deformation of a viscoelastic body

Summary. The tidal deformation of a homogeneous viscoelastic sphere due to the gravitational attraction of an external body is calculated. The sphere is modelled as an incompressible Kelvin-Voigt solid. An equation for the displacement field is obtained assuming that strains are small and inertia is negligible. This equation has a series solution in terms of Legendre polynomials. The resulting expression for the displacement field reduces to that for an elastic solid and a viscous fluid in the appropriate limits of the material constants. The first term in the viscoelastic solution is used to calculate the moments induced by tidal deformation assuming a circular orbit. In the absence of obliquity and precession, these moments reduced to a torque about the spin axis. This torque is compared to that predicted by a phase lag analysis. These two approaches are formally equivalent if the tidal dissipation function Q ⁻¹ depends in a specific way on the difference of the spin and orbital angular velocities.     

Identification of geologic contrasts from landscape dissection pattern: An application to the Cascade Range, Oregon, USA

This paper demonstrates the plausibility of inferring the spatial variability of geology from topographically derived landscape dissection patterns. This enables surveying large regions for spatial variability in geology, for which direct remote sensing is not feasible, by studying variability in dissection pattern, a feature extracted straight off from digital elevation model data. Dissection pattern is obtained automatically by a novel algorithm, especially designed to delineate the valleys with high accuracy in order to reflect spatial variability in dissection density. The dissection pattern is encapsulated by a continuous map of drainage density, a raster variable best suited for showing spatial variability of dissection. Such a map, constructed for the study area in the Cascade Range, Oregon, USA, shows a sharp discontinuity in the dissection pattern, indicating change in underlying geology. Possible factors controlling the dissection pattern such as climate, local and regional slopes, vegetation, and geology are examined, and geology has been found to be the dominant controlling factor. The dissection contrast coincides with the boundary between the Western and High Cascades, two geologic provinces with different rock ages and types. The older and less permeable Western Cascades are associated with denser dissection pattern, whereas the younger and more permeable High Cascades correspond to less dissected pattern. This new mapping method can be applied to locations where topography is the only readily available data, and the generated map could be used to extract previously unknown geologic or environmental information.     

Land use and climate control the spatial distribution of soil types in the grasslands of Inner Mongolia

The spatial distribution of soil types is controlled by a set of environmental factors such as climate, organisms, parent material and topography as well as time and space. A change of these factors will lead to a change in the spatial distribution of soil types. In this study, we use a digital soil mapping approach to improve our knowledge about major soil type distributing factors in the steppe regions of Inner Mongolia (China) which currently undergo tremendous environmental change, e.g. climate and land use change. We use Random Forests in an effort to map Reference Soil Groups according to the World Reference Base for Soil Resources (WRB) in the Xilin River catchment. We benefit from the superior prediction capabilities of RF and additional interpretive results in order to identify the major environmental factors that control spatial patterns of soil types. The nine WRB soil groups that were identified and spatially predicted for the study area are Arenosol, Calcisol, Cambisol, Chernozem, Cryosol, Gleysol, Kastanozem, Phaeozem and Regosol.Model and prediction performances of the RF model are high with an Out-of-Bag error of 51.6% for the model and a misclassification error for the predicted map of 28.9%. The main controlling factors of soil type distribution are land use, a set of topographic variables, geology and climate. However, land use and climate are of major importance and topography and geology are of minor importance. The visualizations of the predictions, the variable importance measures as result of RF and the comparisons of these with the spatial distribution of the environmental factors delivered additional, quantitative information of these controlling factors and revealed that intensively grazed areas are subjected to soil degradation. However, most of the area is still governed by natural soil forming processes which are driven by climate, topography and geology. Most importantly though, our study revealed that a shift towards warmer temperatures and lower precipitation regimes will lead to a change of the spatial distribution of RSGs towards steppe soils that store less carbon, i.e. a decrease of spatial extent of Phaeozems and an increase of spatial extent of Chernozems and Kastanozems.     

A finite element analysis of tidal deformation of the entire earth with a discontinuous outer layer

Tidal deformation of the Earth is normally calculated using the analytical solution with some simplified assumptions, such as the Earth is a perfect sphere of continuous media. This paper proposes an alternative way, in which the Earth crust is discontinuous along its boundaries, to calculate the tidal deformation using a finite element method. An in-house finite element code is firstly introduced in brief and then extended here to calculate the tidal deformation. The tidal deformation of the Earth due to the Moon was calculated for an geophysical earth model with the discontinuous outer layer and compared with the continuous case. The preliminary results indicate that the discontinuity could have different effects on the tidal deformation in the local zone around the fault, but almost no effects on both the locations far from the fault and the global deformation amplitude of the Earth. The localized deformation amplitude seems to depend much on the relative orientation between the fault strike direction and the loading direction (i.e. the location of the Moon) and the physical property of the fault.     

On the relationship between MODIS Land Surface Temperature and topography in Iran

Land Surface Temperature (LST) is considered important in monitoring the energy flux between the land surface and atmosphere. Due to the diversity of topography in Iran and its effect on the climate diversity, we decided to study the effect of topography on the LST variations. To this end, the LST digital data derived from the observations of the MODIS Terra and Aqua were used. The results indicated that, during the daytime, from sea level up to a height of 400 meters, the LST increased, and then the temperature decreased with increasing altitude, and up to a height of 3000 meters, there was a strong correlation between the two. LST lapse rate was more during the daytime compared with that of the night time and it was more during the winter compared with the summer. LST lapse rate showed larger variability in diurnal cycle, but its monthly patterns were similar in different aspects. The aspect had substantial effect on LST inversion elevation. Furthermore, the inverse relationship between LST and slope was strong in slopes up to 20°.     

洋山港海域潮动力特征及其对工程的响应

距上海陆域外约32 km的洋山港是以小洋山岛链为基础建造的中国首个离岸式深水港口, 其港区陆域是通过封堵北岛链三个主要汊道并抛填成陆而成。大规模的封堵陆连工程改变了该海域的地形边界条件, 海域潮动力特征也随之响应。通过对洋山港海域堵汊陆连工程前后潮位及潮流实测资料的分析研究, 发现洋山港海域的潮动力特征主要表现在:海域潮波倍潮及复合潮增幅随工程的进行而增大, 潮动力表现出明显的非线性特征;海域中部发生壅水;海域潮波向驻波性质变异并于憩流时段引发环状流场。     

Reconstructing ancient topography through erosion modelling

One of the main aims of geomorphology is to understand how geomorphic processes change topography over long time scales. Over the last decades several landscape evolution models have been developed in order to study this question. However, evaluation of such models has often been very limited due to the lack of necessary field data. In this study we present a topography based hillslope erosion and deposition model that is based on the WATEM/SEDEM model structure and works on a millennial time scale. Soil erosion, transport and deposition are calculated using slope and unit contributing area. The topography is iteratively rejuvenated by taking into account modelled erosion and deposition rates, thereby simulating topographic development backwards in time. A first attempt has been made to spatially evaluate the model, using detailed estimates for historical soil erosion and sediment deposition volumes, obtained from an augering campaign in a small catchment in the Belgian Loess Belt. The results show that the model can simulate realistic soil redistribution patterns. However, further research is necessary in order to deal with artificial flaws that cause routing problems and significantly influence results. Common problems and issues related to this type of backward modelling are also discussed.     

影响深埋长隧道岩爆的主要因素分析

岩爆是深埋长隧道的主要工程地质问题之一。影响岩石发生岩爆的主要因素有围岩性质、埋深、区域构造应力、地形地貌、活断层和地震等。文中对上述因素对岩爆的影响进行了分析,并指出精细的地面地质调查是合理评价深埋长隧道岩爆的基础。     

Landslide susceptibility mapping using GIS and the weight-of-evidence model

The weights-of-evidence model (a Bayesian probability model) was applied to the task of evaluating landslide susceptibility using GIS. Using landslide location and a spatial database containing information such as topography, soil, forest, geology, land cover and lineament, the weights-of-evidence model was applied to calculate each relevant factor's rating for the Boun area in Korea, which had suffered substantial landslide damage following heavy rain in 1998. In the topographic database, the factors were slope, aspect and curvature; in the soil database, they were soil texture, soil material, soil drainage, soil effective thickness and topographic type; in the forest map, they were forest type, timber diameter, timber age and forest density; lithology was derived from the geological database; land-use information came from Landsat TM satellite imagery; and lineament data from IRS satellite imagery. Tests of conditional independence were performed for the selection of factors, allowing 43 combinations of factors to be analysed. For the analysis of mapping landslide susceptibility, the contrast values, W ⁺ and W ^-, of each factor's rating were overlaid spatially. The results of the analysis were validated using the previous landslide locations. The combination of slope, curvature, topography, timber diameter, geology and lineament showed the best results. The results can be used for hazard prevention and land-use planning.