High resolution mapping of Earth tide response based on GPS data in Japan

We observe the Earth tidal fields at diurnal and semi-diurnal periods using Kinematic Precise Point Positioning (KPPP) GPS analysis. Our KPPP GPS solutions compare well with super-conducting gravimeter (SG) observations and a theoretical Earth tidal model, that includes both ocean tide loading model and body tides. We make a high resolution map of the observed Earth tidal response fields using the Japanese GEONET GPS network which consists of 1200 sites. We find that: (1) the average phase of GPS data lags 0.11 ± 0.04° from our theoretical Earth tidal model, (2) the average amplitude ratio between GPS and the theoretical Earth tidal model is 1.007 ± 0.003, (3) the amplitude in the Kyushu district is about 1.0–1.5 ± 0.3% larger than in the Hokkaido district, and (4) the amplitude at the Japan Sea side is about 0.5 ± 0.2% larger than that at the Pacific Ocean side. These results suggest that we may be able to place constraints on Earth structure using GPS-derived tidal information.     

Applying local Green's functions to study the influence of the crustal structure on hydrological loading displacements

The influence of the elastic Earth properties on seasonal or shorter periodic surface deformations due to atmospheric surface pressure and terrestrial water storage variations is usually modeled by applying a local half-space model or an one dimensional spherical Earth model like PREM from which a unique set of elastic load Love numbers, or alternatively, elastic Green's functions are derived. The first model is valid only if load and observer almost coincide, the second model considers only the response of an average Earth structure. However, for surface loads with horizontal scales less than 2500 km², as for instance, for strong localized hydrological signals associated with heavy precipitation events and river floods, the Earth elastic response becomes very sensitive to inhomogeneities in the Earth crustal structure.We derive a set of local Green's functions defined globally on a 1° × 1° grid for the 3-layer crustal structure TEA12. Local Green's functions show standard deviations of ±12% in the vertical and ±21% in the horizontal directions for distances in the range from 0.1° to 0.5°. By means of Green's function scatter plots, we analyze the dependence of the load response to various crustal rocks and layer thicknesses. The application of local Green's functions instead of a mean global Green's function introduces a variability of 0.5–1.0 mm into the hydrological loading displacements, both in vertical and in horizontal directions. Maximum changes due to the local crustal structures are from −25% to +26% in the vertical and −91% to +55% in the horizontal displacements. In addition, the horizontal displacement can change its direction significantly. The lateral deviations in surface deformation due to local crustal elastic properties are found to be much larger than the differences between various commonly used one-dimensional Earth models.     

A collaborative study on temperature diurnal tide in the midlatitude mesopause region (41°N, 105°W) with Na lidar and TIMED/SABER observations

The na lidar-observed temperature diurnal tidal perturbations, based on full-diurnal-cycle observations from 2002 to 2008, are compared with tidal wave measurements by the TIMED/SABER instrument to elucidate the nature of diurnal tidal-period perturbations observed locally. The diurnal amplitude and phase profiles deduced by the two instruments are in very good agreement most of the year. However, the lidar-observed diurnal amplitudes during winter months and early spring are considerably larger than SABER observations, leading to the existence of a significant amplitude maximum of 12 K near 90 km in February and a different seasonal structure of temperature diurnal amplitude from the two instruments. The lidar-observed diurnal phase shows propagating wave characteristics during equinoctial months, but exhibit “evanescent wave” behavior in winter months, whereas SABER diurnal tidal phase exhibits propagating diurnal tidal character all year long with small seasonal variation. This anomalous tidal characteristic from the lidar observations repeats almost every winter. The exact mechanism behind this tidal feature is not fully understood, therefore further investigation and more experimental observations are necessary.     

Elastic anisotropy of ferromagnesian post-perovskite in Earth's D” layer

We report experimental observation of a sizable elastic anisotropy in a polycrystalline sample of ferromagnesian silicate in post-perovskite (ppv) structure. Using a novel composite X-ray transparent gasket to contain and synthesize ppv in a panoramic diamond-anvil cell along with oblique X-ray diffraction geometry, we observed the anisotropic lattice strain and {1 0 0} or {1 1 0} slip-plane texture in the sample at 140 GPa. We deduced the elasticity tensor (c_ij), orientation-dependent compressional wave velocities, polarization-dependent shear-wave velocities, and the velocity anisotropy of the silicate ppv. Our results are consistent with calculations and indicate that with sufficient preferred orientation, the elastic anisotropy of this phase can produce large shear-wave splitting.     

基于一种新智能优化算法与谱相关方法对重力固体潮的分析

根据月球与太阳相对于地球位置的关系,建立一个三维正交分解模型来揭示重力固体潮中包含的谐波成分。在此基础上,提出一种基于单形邻域与多角色进化策略智能优化算法的独立分量分析方法,用以提高算法效率和全局寻优能力;并结合谱相关方法,对重力固体潮观测数据进行处理和分析,观察其频域和时域的变化。最后引入理论信号和理论计算值作为参考背景,以凸显其异常变化特征。研究表明,利用该方法可以有效地反映出重力固体潮潮汐的日波、半日波、长周期波等特征。同时,发现重力固体潮潮汐的异常变化中隐含有地震前兆信息,通过对云南地区观测数据的分析,在排除测量误差和环境影响外,发现在地震前2~5个月内,重力固体潮的长周期波均呈现出规律性的异常变化。这反映出地震的前兆异常信息,可能对地震时间和地点的预测提供依据。     

Sea level rise in the north-western part of the Arabian Gulf

Relative sea level variations in the north-western part of the Arabian Gulf have been estimated in the past using no more than 10 to 15 years of observations. In this study, we have almost doubled the period to 28.7 years by examining all available tide gauge data in the area and constructing a mean gauge time-series from seven coastal tide gauges. We found for the period 1979–2007 a relative sea level rise of 2.2 ± 0.5 mm/year. Using the subsidence observed at 6 GPS stations within a radius of 100 km of the tide gauges as an indication of the vertical land motion, the corresponding absolute sea level rise is 1.5 ± 0.8 mm/year that is in agreement with the global estimate of 1.9 ± 0.1 mm/year (Church and White, 2011) for the same studied period. By taking into account the temporal correlations we conclude that previous published results underestimate the true sea level rate uncertainty in this area by a factor of 5–10.     

Reducing atmospheric noise in RST analysis of TIR satellite radiances for earthquakes prone areas satellite monitoring

Space–time fluctuations of the Earth’s emitted Thermal Infrared (TIR) radiation observed from satellite from months to weeks before an earthquake are reported in several studies. Among the others, a Robust Satellite data analysis Technique (RST) was proposed (and applied to different satellite sensors in various geo-tectonic contexts) to discriminate anomalous signal transients possibly associated with earthquake occurrence from normal TIR signal fluctuations due to other possible causes (e.g. solar diurnal–annual cycle, meteorological conditions, changes in observational conditions, etc.). Variations in satellite view angle depending on satellite’s passages (for polar satellites) and atmospheric water vapour fluctuations were recognized in the past as the main factors affecting the residual signal variability reducing the overall Signal-to-Noise (S/N) ratio and the potential of the RST-based approach in identifying seismically related thermal anomalies. In this paper we focus on both factors for the first time, applying the RST approach to geostationary satellites (which guarantees stable view angles) and using Land Surface Temperature (LST) data products (which are less affected by atmospheric water vapour variability) instead of just TIR radiances at the sensor.The first results, obtained in the case of the Abruzzo earthquake (6 April 2009, M_W ∼ 6.3) by analyzing 6 years of SEVIRI (Spinning Enhanced Visible and Infrared Imager on board the geostationary Meteosat Second Generation satellite) LST products provided by EUMETSAT, seem to confirm the major sensitivity of the proposed approach in detecting perturbations of the Earth’s thermal emission a few days before the main shock. The results achieved in terms of increased S/N ratio (in validation) and reduced “false alarms” rate (in confutation) are discussed comparing results obtained by applying RST to LST products with those achieved by applying an identical RST analysis (using the same MSG-SEVIRI 2005–2010 data-set) to the simple TIR radiances at the sensor.     

Inferred fault geometry and slip distribution of the 2010 Jiashian,Taiwan, earthquake is consistent with a thick-skinned deformation model

We invert measurements of coseismic displacements from 139 continuously recorded GPS sites from the 2010, Jiashian, Taiwan earthquake to solve for fault geometry and slip distribution using an elastic uniform stress drop inversion. The earthquake occurred at a depth of ~ 23 km in an area between the Western Foothills fold-and-thrust belt and the crystalline high mountains of the Central Range, providing an opportunity to examine the deep fault structure under Taiwan. The inferred rupture plane is oblique to the prominent orientation of thrust faults and parallel to several previously recognized NW-striking transfer zones that appear to connect stepping thrusts. We find that a fault striking 318°–344° with dip of 26°–41° fits the observations well with oblique reverse-sinistral slip under a low stress drop of about 0.5 MPa. The derived geodetic moment of 2.92 × 10¹⁸ N-m is equivalent to a M_w = 6.24 earthquake. Coseismic slip is largely concentrated within a circular patch with a 10-km radius at the depth between 10 and 24 km and maximum slip of 190 mm. We suggest this earthquake ruptured the NW-striking Chishan transfer fault zone, which we interpret as a listric NE-dipping lateral ramp with oblique slip connecting stepping thrust faults (ramps). The inferred slip on the lateral ramp is considerably deeper than the 7–15 km deep detachment identified in previous studies of western Taiwan. We infer an active basal detachment under western Taiwan at a depth of at least ~ 20–23 km based on these inversion results. The earthquake may have nucleated at the base of the lateral ramp near the intersection with the basal detachment. Coulomb stress change calculations suggest that this earthquake moved several NE-striking active thrust faults in western Taiwan nearer to failure.     

Three-dimensional numerical thermal and rheological modelling in the central Fennoscandian Shield

A three-dimensional model for the central Fennoscandian Shield was constructed for analysing the thermal, the rheological and the structural conditions in the lithosphere. The mesh covers a rectangular area in the southern Finland with horizontal dimensions of 500 km × 400 km and a depth extent of 100 km. Structural boundaries are derived from the several deep seismic soundings carried out in the area. Constructed model is first used in the calculation of the thermal and the rheological models and secondly in analysing the stress and the deformational conditions with the obtained rheology. Thermal and structural models are solved with the finite element method. The calculated surface HFD is between 40 and 48 mW m⁻² in the Proterozoic southern part and below 40 mW m⁻² in the older and northern Archaean part of the model. The calculated rheological strength shows a layered structure with two individual rheologically weak layers in the crust and strong layer in the upper part of the lower crust. The minimum brittle–ductile transition (BDT) depth is around 10 km in the southern part of the model while in the north and north-eastern parts the BDT depth is around 45–50 km. Comparison with the focal depth data shows that as most of the earthquakes occur no deeper than the depth of 10 km are they located in the brittle regime. Resulting stress conditions and possible regions of deformation after the model is subjected to pressure of 50 MPa reveals that the stress field is quite uniformly distributed in different crustal layers and that the elastic parameters control more the state of the stress than the applied rheological structure. In the upper crust, the stress intensity has values between 42 and 45 MPa whereas in the middle crust the values are around 50 MPa. Comparison of the 3-D model with earlier 2-D models shows that some differences in the results are to be expected.     

The GGOS as the backbone for global observing and local monitoring: A user driven perspective

A key geodetic contribution to both the three Global Observing Systems and initiatives like the European Global Monitoring for Environment and Security is an accurate, long-term stable, and easily accessible reference frame as the backbone. Many emerging scientific as well as non-scientific high-accuracy applications require access to an unique, technique-independent reference frame decontaminated for short-term fluctuations due to global Earth system processes. Such a reference frame can only be maintained and made available through an observing system such as the Global Geodetic Observing System (GGOS), which is currently implemented and expected to provide sufficient information on changes in the Earth figure, its rotation and its gravity field. Based on a number of examples from monitoring of infrastructure, point positioning, maintenance of national references frames to global changes studies, likely future accuracy requirements for a global terrestrial reference frame are set up as function of time scales. Expected accuracy requirements for a large range of high-accuracy applications are less than 5 mm for diurnal and sub-diurnal time scales, 2–3 mm on monthly to seasonal time scales, better than 1 mm/year on decadal to 50 years time scales. Based on these requirements, specifications for a geodetic observing system meeting the accuracy requirements can be derived.     

Complex demodulation in VLBI estimation of high frequency Earth rotation components

The spectrum of high frequency Earth rotation variations contains strong harmonic signal components mainly excited by ocean tides along with much weaker non-harmonic fluctuations driven by irregular processes like the diurnal thermal tides in the atmosphere and oceans. In order to properly investigate non-harmonic phenomena a representation in time domain is inevitable. We present a method, operating in time domain, which is easily applicable within Earth rotation estimation from Very Long Baseline Interferometry (VLBI). It enables the determination of diurnal and subdiurnal variations, and is still effective with merely diurnal parameter sampling. The features of complex demodulation are used in an extended parameterization of polar motion and universal time which was implemented into a dedicated version of the Vienna VLBI Software VieVS. The functionality of the approach was evaluated by comparing amplitudes and phases of harmonic variations at tidal periods (diurnal/semidiurnal), derived from demodulated Earth rotation parameters (ERP), estimated from hourly resolved VLBI ERP time series and taken from a recently published VLBI ERP model to the terms of the conventional model for ocean tidal effects in Earth rotation recommended by the International Earth Rotation and Reference System Service (IERS). The three sets of tidal terms derived from VLBI observations extensively agree among each other within the three-sigma level of the demodulation approach, which is below 6 μas for polar motion and universal time. They also coincide in terms of differences to the IERS model, where significant deviations primarily for several major tidal terms are apparent. An additional spectral analysis of the as well estimated demodulated ERP series of the ter- and quarterdiurnal frequency bands did not reveal any significant signal structure. The complex demodulation applied in VLBI parameter estimation could be demonstrated a suitable procedure for the reliable reproduction of high frequency Earth rotation components and thus represents a qualified tool for future studies of irregular geophysical signals in ERP measured by space geodetic techniques.     

Taiwan strait current in winter

Four bottom-mounted current profilers were deployed across the Taiwan Strait from September 28 to December 14 of 1999 to monitor the current velocity when the northeast monsoon was strong. Results indicate both diurnal and semidiurnal tidal currents were primarily barotropic. The barotropic diurnal tide might be explained by a single Kelvin wave propagating along the Mainland China coast from north to south. However, the barotropic semidiurnal tide manifested as a more complicated form in the Taiwan Strait.The subtidal current generally fluctuated with the northeast winds. When the northeast wind was weak, the along- and cross-strait subtidal current flowed primarily against the wind and toward Taiwan, respectively. As the northeast wind intensified, the along-strait current flowed downwind, brought the cold China coastal water southward, and formed a baroclinic velocity front in the western portion of the Taiwan Strait. The Ekman effect forced the cross-strait current toward Mainland China in the upper water column and toward Taiwan in the lower water column, respectively. The along-strait volume transport, estimated from interpolated current velocity, varied from −5 to 2 Sv with a mean value of 0.12±0.33 Sv. Similar transport was also estimated from the sea level difference across the Taiwan Strait.Although the local wind played a dominant role for the fluctuations of current velocity and transport in the Taiwan Strait, it could be not the only important factor. The current or transport directed frequently against the wind could be related to the northward current, which was consistently observed in the Penghu Channel.     

On the tidal heating of Enceladus

Enceladus, one of Saturn's moons, shows significant volcanic activity identified by the Cassini spacecraft. The aim of the present study is to investigate – with the adaptation of mathematical tools used in geodynamics – the extent of tidal heating due to the mean motion resonance with Dione. For the purpose of calculations a two-layer model of Enceladus was used. The inner part of the model is a “rocky core” with a relative radius 0.55, while the outer part is composed of water ice. The results of model calculations show that the effective tidal heating is not uniformly distributed within Enceladus. It was found for the selected model of Enceladus, that the tidal heating is maximum within the depth interval (25–75) km. Due to the inhomogeneity within Enceladus, 85% of the tidal energy is generated in a volume that contains just 39% of its mass. In time intervals of 3.0 × 10⁸ and 5.3 × 10⁸ years the temperature increase in the relative depth range 0.70 ≤ r/a_E ≤ 0.90 is approximately 270 and 370 K, respectively.     

Large-scale sedimentary bedforms and sediment dynamics on a glaciated tectonic continental shelf: Examples from the Pacific margin of Canada

The Pacific margin of Canada has been subjected to tectonism, dramatic sea level change and vigorous storm and tidal energy since glacial times resulting in a complex seafloor. Extensive multibeam mapping of this shelf has provided an opportunity to understand how these processes have impacted sedimentology and morphology. Bathymetric restriction of the tidally dominated flow between the inland seas and the open Pacific has resulted in the development of very large subaqueous dune fields and terrace moats. For example, in the southern Strait of Georgia nearly symmetrical dunes with wavelengths between 100 and 300 m, dune heights up to 28 m, cover the seafloor in 170–210 m water depth. In northern Hecate Strait a 72 km² area of large 2D dunes occurs at the transition with Dixon Entrance which opens to the Pacific Ocean and steep (>10°) wave-cut terraces and drowned spits, a result of sea level changes during the Holocene, are now being undercut to generate moats 7 m deep, in a narrowing shelf trough. Currents, with velocities ranging between 0.2 and 2.2 m s^?1, are dominated by semi-diurnal tidal streams that are continually modified by wind and estuarine circulation. There appears to be a clear association of grain size, water depth and flow velocity controlling the size of the subaqueous dunes.     

The thermal state and strength of the lithosphere in the Spanish Central System and Tajo Basin from crustal heat production and thermal isostasy

In this work we have modeled the thermal structure of the lithosphere of the Spanish Central System and the Tajo Basin, and their implications for lithospheric strength. For this, we have used refined heat-producing elements (HPE) values to obtain new estimates of heat production rates in the Spanish Central System and Tajo Basin areas, which have been used joined to the relation between topography and thermal structure of the lithosphere to calculate the best-fit surface heat flows in the study area. Moreover, we have implemented a temperature-dependent thermal conductivity (appropriate for olivine) for the lithospheric mantle to improve the calculations of temperature profiles in the mantle. The geotherms so obtained, together with the implementation of a new rheological law for the upper lithospheric mantle, have been used to calculate refined estimations of the strength and effective elastic thickness of the lithosphere. We have obtained surface heat flow values of 84 mW m⁻² and ∼82 mW m⁻² for the Spanish Central System and the Tajo Basin, respectively. The thermal state of the lithosphere affects mantle temperatures, and hence may be playing an important role in the uplift and maintenance of the Spanish Central System.     

New data on the Late Cenozoic basaltic volcanism in Syria,applied to its origin

New data on geology and 21 K–Ar dates of the Late Oligocene–Quaternary basalts in Syria, combined with analysis of the new and previous data are used to reconstruct the volcanic history and relations between it and tectonic events. Volcanism began at the end of Oligocene (26–24 Ma) and was concentrated in the Late Oligocene–Early Miocene along a N-trending band, which stretches from the Jebel Arab (Harrat Ash Shaam) up to Kurd Dagh and southern Turkey. Activity waned in the Middle Miocene (17–12 Ma), but was resumed in the same band in the Tortonian and increased in the Messinian and Early Pliocene (6.3–4 Ma), when volcanism spread to the Shin Plateau and its coastal extension. After a brief hiatus ～ 4–3.5 Ma, volcanism became still more intensive and spread from the N-trending band to the east into the northern margin of the Mesopotamian Foredeep and to the west into the Dead Sea Transform zone. Additional eruptions continued into the Holocene.Volcanism lasted > 25 million years in the Jebel Arab Highland and > 15 million years in the Aleppo Plateau. The long duration of volcanism in the same parts of the moving Arabian plate and absence of records of one-way migration of the activity mean that the magmatic sources moved together with the plate, i.e., they were situated within the lithosphere mantle. Coincidence of the tectonic and volcanic stages of the Arabian plate development proves that volcanic activity depended on the geodynamic situation, caused by the plate motion. Situated within the lithosphere, magmatic sources within this transverse band were possibly caused by thermal and deforming influences of the asthenospheric lateral flow, moved laterally from the Ethiopia–Afar deep superplume.     

Single-crystal elasticity of stishovite: New experimental data obtained using high-frequency resonant ultrasound spectroscopy and a Gingham check structure model

Single-crystal elasticity of stishovite was examined using a new experimental technique and an empirical macroscopic model. Employing high-frequency resonant ultrasound spectroscopy, single-crystal elastic constants of stishovite were determined: C₁₁ = 443(3), C₃₃ = 781(4), C₁₂ = 193(2), C₂₃ = 199(2), C₄₄ = 256(2), and C₃₃ = 316(2) GPa. The frequency range of the resonant ultrasound spectroscopy techniques was 6–20 MHz, which is much lower than the ～10 GHz range of the Brillouin scattering technique. Of the elastic constants, the shear elastic constants C₄₄ and C₆₆ are consistent with the average value of the previously mentioned Brillouin scattering. Conversely, the four elastic constants, C₁₁, C₃₃, C₁₂, and C₂₃, slightly deviate outside the range of previous Brillouin scattering results. The present results, except those for C₁₂, are consistent with recent lattice dynamic analysis of inelastic X-ray scattering data. The adiabatic bulk modulus was calculated as 298 GPa, which is smaller and more consistent with the result of compression experiments than any other Brillouin scattering results (301–312 GPa). The present result shows greater P-wave velocity anisotropy (24.7%) than any preceding work. To understand the unique elastic properties of stishovite, the Gingham check model was proposed and examined. The result shows that the octahedron of 6-coordinated Si in stishovite crystal has stiffness comparable to that of diamond.     

Least squares self-coherence for sub-nGal signal detection in the superconducting gravimeter records

In this paper, we use superconducting gravimeter (SG) data recorded at three stations of the global geodynamics project (GGP) network, with good geographical distribution, to search for possible significant peaks in the gravity spectrum that are in the assumed period range of the Slichter triplet. Seven-year long series from Cantley (Canada), and four-year long series from both Canberra (Australia) and Moxa (Germany) stations are used. First, a solid Earth and ocean tide model is subtracted from the data, followed by a local atmospheric pressure correction based on a frequency-, and location-dependent admittance estimated by the least squares response method. Subsequently, the residual series are filtered with a Parzen-based bandpass filter with a passband (12 h–78 s). A sub-nGal detection level is confirmed by injecting an artificial sine wave of different amplitudes. The Least Squares Self-Coherency spectrum shows the existence of many apparently statistically significant peaks at the 95% confidence level in the band (3–8 h). Although a few peaks are close to the claimed Slichter periods in previous research, the large number of candidate peaks may be related to other mechanisms such as global pressure variations, or hydrology.     

Oxygen isotope constraints on the origin and differentiation of the Moon

We report new high-precision laser fluorination three-isotope oxygen data for lunar materials. Terrestrial silicates with a range of δ¹⁸O values (− 0.5 to 22.9‰) were analyzed to independently determine the slope of the terrestrial fractionation line (TFL; λ = 0.5259 ± 0.0008; 95% confidence level). This new TFL determination allows direct comparison of lunar oxygen isotope systematics with those of Earth. Values of Δ¹⁷O for Apollo 12, 15, and 17 basalts and Luna 24 soil samples average 0.01‰ and are indistinguishable from the TFL. The δ¹⁸O values of high- and low-Ti lunar basalts are distinct. Average whole-rock δ¹⁸O values for low-Ti lunar basalts from the Apollo 12 (5.72 ± 0.06‰) and Apollo 15 landing sites (5.65 ± 0.12‰) are identical within error and are markedly higher than Apollo 17 high-Ti basalts (5.46 ± 0.11‰). Evolved low-Ti LaPaz mare-basalt meteorite δ¹⁸O values (5.67 ± 0.05‰) are in close agreement with more primitive low-Ti Apollo 12 and 15 mare basalts. Modeling of lunar mare-basalt source composition indicates that the high- and low-Ti mare-basalt mantle reservoirs were in oxygen isotope equilibrium and that variations in δ¹⁸O do not result from fractional crystallization. Instead, these differences are consistent with mineralogically heterogeneous mantle sources for mare basalts, and with lunar magma ocean differentiation models that result in a thick feldspathic crust, an olivine–pyroxene-rich mantle, and late-stage ilmenite-rich zones that were convectively mixed into deeper portions of the lunar mantle. Higher average δ¹⁸O (WR) values of low-Ti basalts compared to terrestrial mid ocean ridge basalts (Δ=0.18‰) suggest a possible oxygen isotopic difference between the terrestrial and lunar mantles. However, calculations of the δ¹⁸O of lunar mantle olivine in this study are only 0.05‰ higher than terrestrial mantle olivine. These observations may have important implications for understanding the formation of the Earth–Moon system.