Horizontal rotation signals detected by “G-Pisa” ring laser for the M w = 9.0, March 2011, Japan earthquake

We report the observation of the ground rotation induced by the M _w = 9.0, 11th of March 2011, Japan earthquake. The rotation measurements have been conducted with a ring laser gyroscope operating in a vertical plane, thus detecting rotations around the horizontal axis. Comparison of ground rotations with vertical accelerations from a co-located force balance accelerometer shows excellent ring laser coupling at periods longer than 100?s. Under the plane wave assumption, we derive a theoretical relationship between horizontal rotation and vertical acceleration for Rayleigh waves. Due to the oblique mounting of the gyroscope with respect to the wave direction of arrival, apparent velocities derived from the acceleration/rotation rate ratio are expected to be always larger than or equal to the true wave propagation velocity. This hypothesis is confirmed through comparison with fundamental mode, Rayleigh-wave phase velocities predicted for a standard Earth model.     

Performance of “G-Pisa” ring laser gyro at the Virgo site

The ring laser gyroscope ??G-Pisa?? has been taking data inside the Virgo interferometer central area with the aim of performing high sensitivity measurements of rotations in the vertical as well as in the horizontal orientation. We discuss the main characteristics of the instrument, describing its mechanical design and presenting the measured sensitivity limit. By applying a simple effective model for the laser gyroscope, we show that the stability of the sensor above 10?s of integration time is mainly limited by backscattering effects. The horizontal rotation rate signal is also compared with the signals recorded by the Virgo environmental monitoring system and by a biaxial mechanical tiltmeter rigidly fixed on top of the gyrolaser mounting frame.     

Examining ambient noise using colocated measurements of rotational and translational motion

In the past decade, a number of studies have reported the observation of rotational motion associated with seismic events. We report a first observation of rotational motion in the microseismic ambient noise band. A striking feature of rotational motion measurements is that the information about the seismic phase velocity and source back azimuth is contained in the amplitude ratio of a point measurement of rotation rate and transverse acceleration. We investigate the possibility of applying this method to ambient noise measured with a ring laser and a broadband seismometer at the Wettzell Geodetic Observatory in Germany. Using data in the secondary microseismic band, we recover local phase velocities as well as the back azimuth of the strongest noise source for two different time periods. In order to confirm these findings, we additionally compare the results with classical array processing techniques of the Gr?fenberg array located nearby.     

Investigation of array-derived rotation in TAIPEI 101

We investigate rotational motions derived from measurements by arrays of translational seismometers??hereafter called array-derived rotations (ADRs)??and compare these to measurements made by a commercially available point rotation sensor (eentec? R-1?). We focus on two aspects of the array problem: (1) the requisite conditions for calculating an ADR well and (2) the effect of array configuration on the result. Our data set consists of translational accelerations and rotation rates recorded by an array of Kinemetrics? EpiSensor? accelerometers and two R-1? rotational sensors in the TAIPEI 101 building in Taipei, Taiwan. Our results indicate that (1) array configuration affects the accuracy of ADRs about orthogonal components in horizontal plane, (2) coherence between two point rotation measurements (two R-1?) can determine the maximum frequency of translations viably used for the calculating ADRs, and (3) the performance of the R-1? is adequate, at least above a frequency of 0.12?Hz (periods shorter than 8?s). We also discuss deriving strain from the same array.     

Effect of support rotation on triple friction pendulum bearing behavior

This paper presents numerical modeling techniques to capture the behavior of the triple friction pendulum (TFP) isolation bearing when rotation is permitted about the horizontal axes of the top and bottom components. This paper builds on a previous model for the TFP bearing presented by the authors that is based on the kinematic and constitutive relationships of the individual components of the TFP bearing. The effect of rotation on cyclic bearing behavior and seismic system behavior are investigated numerically for two cases: constant support rotation and variable support rotation. It is found that constant support rotations should be limited in amplitude to ensure standard TFP bearing behavior. Results suggest that flexible supports may not have a large effect on global structure performance as long as typical deformation limits for the supporting members are met. In cases of both constant support rotations and flexible supports, the hardening stages of TFP bearing behavior are diminished. Copyright © 2013 John Wiley & Sons, Ltd.     

A robust channel network extraction method combining discrete curve evolution and the skeleton construction technique

The automatic mapping of drainage networks from terrain representation has been an interesting topic in hydrological and geomorphological modeling. However, the existing methods often suffer from high sensitivity to terrain noise or lose significant stream branches and accurate channel paths. In this paper, we propose a contour-based framework in drainage network extraction. The proposed framework incorporates discrete curve evolution (DCE) to eliminate the noise influence by dynamically segmenting the contour lines (CLs) into valley bends, and to detect the valley feature points. The skeleton construction technique is then applied to distill more accurate channel paths in complex terrain. Finally, a linking step is undertaken to generate the channel network. The proposed method was tested on a series of elevation datasets, with varied resolution, region size, and local relief. The experiments verified that the proposed method can achieve highly accurate channel networks and is robust, even in regions with high-contrast relief, and/or in cases with significant terrain noise and irregularities.     

Comparison of point and array-computed rotations for the TAIGER explosions of 4 March 2008

Two large explosions were recorded by a dense array of strong-motion accelerometers and rotational seismometers in northeastern Taiwan associated with a Taiwan Integrated Geodynamics Research long-range refraction experiment. The objective of this experiment was to test the response of the experimental eentec rotational seismometers against calculated array rotations. Computed array rotation rates are seen to have little variation across the array, but point rotation rate measurements obtained from individual rotational seismometers show significant deviations with each other and with the array rotation rates in the ranges of 3?C5?Hz and, especially, 3?C50?Hz. A cross-correlation method was used to compare array-computed rotation rates and point rotation rate measurements in the frequency band of 3?C5?Hz with the result that the absolute value of the normalized maximum correlation coefficient for each station set varied from 0.45 to 0.97 with an average of 0.84. Amplitude differences of the point rotation rate measurements are seen to be factors of 0.2 to 1.8 times the array rotations as well. It is not likely that the differences seen in the point and array-computed rotation rates are due to nonlinear or heterogeneous site conditions under each array element since these effects should also be seen in the acceleration data used to determine rotation rate. A rigorous method for accurately calibrating rotation rate instruments is needed to understand their response in the field.     

Vertical and horizontal components of the lake Baikal electrotelluric field and their relation to the electric conductivity

Synchronous observations of the electrotelluric field were performed on the lake ice and the coast. The ice-based vertical/horizontal field ratio at periods of a few seconds to a few minutes ranges from a few hundredths to a few thousandths, which is comparable to the noise level. The results of these observations favor a plane wave model. At the coastal station, the vertical and horizontal fields are similar, which is related to geoelectric heterogeneity of the medium. An electric tipper is used to control the behavior of heterogeneity in the spectral-time domain. Anomalous changes associated with earthquakes of K ≥ 12 are revealed from data of the tipper monitoring at epicentral distances of up to 300 km. A longitudinal MTS curve free from the effect of local geoelectric heterogeneities is obtained from ice-based data. Deep conducting layers are recognized as minimums in this curve.     

The Large Ring Laser G for Continuous Earth Rotation Monitoring

Ring Laser gyroscopes exploit the Sagnac effect and measure rotations absolute. They do not require an external reference frame and therefore provide an independent method to monitor Earth rotation. Large-scale versions of these gyroscopes promise to eventually provide a similar high resolution for the measurement of the variations in the Earth rotation rate as the established methods based on VLBI and GNSS. This would open the door to a continuous monitoring of LOD (Length of Day) and polar motion, which is not yet available today. Another advantage is the access to the sub-daily frequency regime of Earth rotation. The ring laser “G” (Grossring), located at the Geodetic Observatory Wettzell (Germany) is the most advanced realization of such a large gyroscope. This paper outlines the current sensor design and properties.     

Measurement of particle rotation in a saltation layer

Two computational methods to measure particle rotations from shadow images of sand particles saltating in a wind tunnel are presented. One method calculates the maximum of the cross‐correlations through multiple angular rotations of an imaged particle. The second method polar transforms both images and then calculates the correlation coefficient for multiple pixel displacements in the θ axis, corresponding to particle rotations. The results from both methods were analysed as a function of height above sand bed (3.7–33.4 mm) and particle size (0.32–0.93 mm equivalent mean diameter). Our results indicate little evidence that particle rotation speeds depend on either their size or height above the sand bed. Though similar results were obtained from both methods, there existed different advantages and disadvantages between the methods. Erroneous results likely arose from particles that were inadequately described by a 2‐D rotation axis, or from poorly imaged particles. At a wind tunnel speed of about 12 m/s, most particles rotated at around 300–400 rev/s. Negative rotations were also found, and their proportion was approximately 15% within the total range of ?450 to 850 rev/s. The ratio of displacement kinetic energy to rotation energy was compared across the various groups and had values between 15 and 40. The quotient showed little dependence on height, though decreased with increasing particle size. Wider applicability of the measurement methodology to study snow particle rotation is also discussed. Copyright © 2014 John Wiley & Sons, Ltd.     

Nonsmooth dynamics prediction of measured bridge response involving deck‐abutment pounding

Earthquake‐induced deck‐abutment contact alters the boundary conditions at the deck level and might activate a different mechanical system than the one assumed during the design of the bridge. Occasionally this discrepancy between the assumed and the actual seismic behavior has detrimental consequences, for example, pier damage, deck unseating, or even collapse. Recently, an insightful shake‐table testing of a scaled deck‐abutment bridge model 1 , showed unexpected in‐plane rotations even though the deck was straight. These contact‐induced rotations produced significant residual displacements and damage to the piers and the bents. The present paper utilizes that experimental data to examine the validity and the limitations of a proposed nonsmooth dynamic analysis framework. The results show that the proposed approach satisfactorily captures the planar rigid‐body dynamics of the deck which is characterized by deck‐abutment contact. The analysis brings forward the role of friction on the physical mechanism behind the rotation of the deck, and underlines the importance of considering the frictional contact forces during deck‐abutment interaction even for straight bridges, which typically are neglected. Finally, the paper investigates the sensitivity of the rotation with respect to macroscopic contact parameters (i.e., the coefficient of friction and the coefficient of restitution). Copyright © 2017 John Wiley & Sons, Ltd.     

Range imaging: a new method for high‐resolution topographic measurements in small‐ and medium‐scale field sites

Topographic measurements are essential for the study of earth surface processes. Three‐dimensional data have been conventionally obtained through terrestrial laser scanning or photogrammetric methods. However, particularly in steep and rough terrain, high‐resolution field measurements remain challenging and often require new creative approaches. In this paper, range imaging is evaluated as an alternative method for obtaining surface data in such complex environments. Range imaging is an emerging time‐of‐flight technology, using phase shift measurements on a multi‐pixel sensor to generate a distance image of a surface. Its suitability for field measurements has yet not been tested. We found ambient light and surface reflectivity to be the main factors affecting error in distance measurements. Low‐reflectivity surfaces and strong illumination contrasts under direct exposure to sunlight lead to noisy distance measurements. However, regardless of lighting conditions, the accuracy of range imaging was markedly improved by averaging multiple images of the same scene. For medium ambient lighting (shade) and a light‐coloured surface the measurement uncertainty was approximately 9 mm. To further test the suitability of range imaging for field applications we measured a reach of a steep mountain stream with a horizontal resolution of approximately 1 cm (in the focal plane of the camera), allowing for the interpolation of a digital elevation model on a 2 cm grid. Comparison with an elevation model obtained from terrestrial laser scanning for the same site revealed that both models show similar degrees of topographic detail. Despite limitations in measurement range and accuracy, particularly at bright ambient lighting, range imaging offers three‐dimensional data in real time and video mode without the need of post‐processing. Therefore, range imaging is a useful complement or alternative to existing methods for high‐resolution measurements in small‐ to medium‐scale field sites. Copyright © 2012 John Wiley & Sons, Ltd.     

Structural Setting of the Northeastern Biga Peninsula (Turkey) from Tilt Derivatives of Gravity Gradient Tensors and Magnitude of Horizontal Gravity Components

The tilt derivatives (or angles) are determined using the gravity gradient tensors (GGT) from two horizontal components (Gx, Gy) and vertical component (Gz), and the magnitude of the horizontal components (MHC). We show that the tilt derivatives from GGT and MHC are highly suitable for mapping linear geological structures or edges of target geology. The results obtained from theoretical data, with and without random noise, have been analyzed in this study. The tilt derivatives from GGT and MHC allow imaging the horizontal boundaries of gravity sources with a high resolution and show an improvement performance as edge detectors since they are ultra sensitive for detecting source boundaries. The tilt derivatives from GGT and MHC over the northeast of the Biga Peninsula in northwestern, Turkey are interpreted to image the edges of geological structures. Tilt derivatives provide a quick method which works well to obtain a detailed structural image of the complex area, with the most precision. The results indicate that the most predominant structural trends are in a NE-SW direction. Therefore, a good correspondence is recognized between structures well known by surface geology and many other lineaments which are unknown or only partially known.     

Remote Sensing Data in Wind Velocity Field Modelling: a Case Study from the Sudetes (SW Poland)

The phenomena of wind-field deformation above complex (mountainous) terrain is a popular subject of research related to numerical modelling using GIS techniques. This type of modelling requires, as input data, information on terrain roughness and a digital terrain/elevation model. This information may be provided by remote sensing data. Consequently, its accuracy and spatial resolution may affect the results of modelling. This paper represents an attempt to conduct wind-field modelling in the area of the ?nie?nik Massif (Eastern Sudetes). The modelling process was conducted in WindStation 2.0.10 software (using the computable fluid dynamics solver Canyon). Two different elevation models were used: the Global Land Survey Digital Elevation Model (GLS DEM) and Digital Terrain Elevation Data (DTED) Level 2. The terrain roughness raster was generated on the basis of Corine Land Cover 2006 (CLC 2006) data. The output data were post-processed in ArcInfo 9.3.1 software to achieve a high-quality cartographic presentation. Experimental modelling was conducted for situations from 26 November 2011, 25 May 2012, and 26 May 2012, based on a limited number of field measurements and using parameters of the atmosphere boundary layer derived from the aerological surveys provided by the closest meteorological stations. The model was run in a 100-m and 250-m spatial resolution. In order to verify the model’s performance, leave-one-out cross-validation was used. The calculated indices allowed for a comparison with results of former studies pertaining to WindStation’s performance. The experiment demonstrated very subtle differences between results in using DTED or GLS DEM elevation data. Additionally, CLC 2006 roughness data provided more noticeable improvements in the model’s performance, but only in the resolution corresponding to the original roughness data. The best input data configuration resulted in the following mean values of error measure: root mean squared error of velocity = 1.0 m/s and mean absolute error of direction = 30°. The author concludes that, within specific meteorological conditions (relatively strong and constant synoptic forcing) and using the aforementioned input data, the Canyon model provides fairly acceptable results. Similarly, the quality of the presented remote sensing data is suitable for wind velocity modelling in the proposed resolution. However, CLC 2006 land use data should be first verified with a higher-resolution satellite or aerial imagery.     

日本琵琶湖流域大气边界层的三维数值模拟

本文在地形坐标下引入一种水平导数差分近似方法,建立一个适合日本琵琶湖复杂地形、满足静力平衡条件的三维非定常大气边界层数值模式。在弱的背景风场下,模式加入了琵琶湖流域的实际地形,计算结果表明:模式运行十分稳定,且占机内存小,节省计算时间,同时模式还做了流域地形、湖泊本身存在与否等因子对边界层风场影响的数值试验。在强的背景风场下,模式加入了经过适当平滑的琵琶湖流域地形,得出了一些有益的看法。     

Triangulation hand‐held laser‐scanning (TriHaLaS) for micro‐ and meso‐habitat surveys in streams

A novel hand‐held laser‐based stream bed survey system is presented. The system facilitates the capture of detailed 3D mapping of shallow (< 0.7 m) riverbed topography in sections approximately 4 m by 2 m. The system includes a stationary reference system, which projects three laser sheets (two at offset angles), within which a hand‐held monitoring pole is moved. The unique configuration of the light sheets intercepts with the pole as it moves within the survey area providing an exact horizontal location. Pole tilt is compensated for by an inertial measurement unit on the pole, and the height above the bed of the pole and submerged scanning laser are monitored relative to the horizontal laser sheet. Verification and application measurements demonstrate high resolution and accuracy in the horizontal (~5 mm) and vertical (~1 mm) direction. The system can be applied at sites where a free view is blocked and other optical through‐water methods fail. It is appropriate for studies on riverbed statistics and dynamics, which necessitate non‐invasive in‐situ surveys. Copyright © 2017 John Wiley & Sons, Ltd.     

Free‐field rotations during earthquakes: Relevance on buildings

Soil rotations around horizontal axes, during an earthquake, are studied through records collected by closely spaced arrays of strong motion accelerometers. The cross power spectrum of accelerations at nearby stations has been generally utilized to describe the spatial distribution of the motion. A number of cross spectra have been obtained during the training of these arrays. To take profit of these elaborations, a mathematical relation is established between the cross power spectrum and the power spectrum of the rotation. Rotation data presented by Liu et al, concerning 52 earthquake records collected at a single station in Taiwan, are compared with rotation data computed according to our procedure. The two series of data are suitably normalized to the peak horizontal acceleration. The data are shown in function of the distance from epicentre. The same ratio, computed according to our procedure, is in good agreement with the average value of these data. Direct measurements and the present approach have lead to evaluations of rotation higher than those predicted by mathematical investigations on the basis of the wave propagation theory, for comparable circumstances. The relevance of this input motion for relatively tall structures is examined, with reference to the structural effects that the horizontal motion concurrently provides. Meaningful will be ranked those effects of the order of magnitude of 20% or higher than those implied by the horizontal excitation. For understanding the relevance on building structures, the procedure has two areas of concern: 1) the coherence implicit in the cross power spectra, which depends on the interpolation process of the original records, collected in the arrays of instruments, and 2) the relative importance of the vertical to the horizontal input motion. As to the second item, the relevance of the rotation component on structures largely depends on the relative importance of the vertical to the horizontal input motion. When the records in an area a few km from the epicenter are considered, the response spectrum of vertical motion is comparable and in some records even higher, than that of horizontal motion, over the entire range of frequencies. This has been observed as well for the 2009 earthquake event of L'Aquila, Italy, and that at the Christchurch (New Zealand) 2011. When the response spectrum of vertical motion is comparable to that of horizontal motion, the effects of rotational motions on most engineering structures can be meaningful. Copyright © 2011 John Wiley & Sons, Ltd.     

Use of fused airborne scanning laser altimetry and digital map data for urban flood modelling

Flood modelling of urban areas is still at an early stage, partly because until recently topographic data of sufficiently high resolution and accuracy have been lacking in urban areas. However, digital surface models (DSMs) generated from airborne scanning laser altimetry (LiDAR) having sub‐metre spatial resolution have now become available, and these are able to represent the complexities of urban topography. This paper describes the development of a LiDAR post‐processor for urban flood modelling based on the fusion of LiDAR and digital map data. The map data are used in conjunction with LiDAR data to identify different object types in urban areas, though pattern recognition techniques are also employed. Post‐processing produces a digital terrain model (DTM) for use as model bathymetry, and also a friction parameter map for use in estimating spatially distributed friction coefficients. In vegetated areas, friction is estimated from LiDAR‐derived vegetation height, and (unlike most vegetation removal software) the method copes with short vegetation less than ～1 m high, which may occupy a substantial fraction of even an urban floodplain. The DTM and friction parameter map may also be used to help to generate an unstructured mesh of a vegetated urban floodplain for use by a two‐dimensional finite element model. The mesh is decomposed to reflect floodplain features having different frictional properties to their surroundings, including urban features (such as buildings and roads) and taller vegetation features (such as trees and hedges). This allows a more accurate estimation of local friction. The method produces a substantial node density due to the small dimensions of many urban features. Copyright © 2007 John Wiley & Sons, Ltd.     

起伏地形下黄河流域太阳散射辐射分布式模拟研究

针对天文、大气、宏观地理与局地地形等因子对起伏地形下太阳散射辐射的复杂作用,将影响起伏地形下太阳散射辐射的天空因素与地面因素分开处理.通过基于数字高程模型（DEM）数据的起伏地形下天文辐射模型和地形开阔度模型,综合考虑地面因素对散射辐射的影响;基于常规地面气象站观测资料建立的水平面散射辐射模型考虑天空因素对散射辐射的影响;依据各向异性散射机理,建立了起伏地形下太阳散射辐射分布式计算模型,探索出一条利用DEM数据和常规气象观测资料实现山区太阳散射辐射定量模拟的技术路线.以1 km×1 km分辨率的DEM数据作为地形的综合反映,实现了起伏地形下黄河流域1 km×1 km分辨率的太阳散射辐射分布式模拟.