Thermal history and large scale differentiation of the Saturn’s satellite Rhea

Thermal history of Rhea from the beginning of accretion is investigated. We developed a numerical model of convection combined with the parameterized theory. Large scale melting of the satellite’s matter and gravitational differentiation of silicates from ices are included. The results are confronted with observational data from Cassini spacecraft that indicate minor differentiation of the satellite’s interior. We suggest that partial differentiation of the satellite’s interior is accompanied (or followed) by the process of light fraction uprising to the surface. The calculation indicates that the partial differentiation of the matter of the satellite’s interior is possible only for narrow range of parameters. In particular, we found that the time from the formation of CAI (calciumaluminum rich inclusions in chondrites) to the end of accretion of Rhea is in the range of 3–4 My.     

Inclined lidar observations of boundary layer aerosol particles above the Kongsfjord,Svalbard

An inclined lidar with vertical resolution of 0.4 m was used for detailed boundary layer studies and to link observations at Zeppelin Mountain (474 m) and Ny-Ålesund, Svalbard. We report on the observation of aerosol layers directly above the Kongsfjord. On 29 April 2007, a layer of enhanced backscatter was observed in the lowest 25 m above the open water surface. The low depolarization ratio indicated spherical particles. In the afternoon, this layer disappeared. The ultrafine particle concentration at Zeppelin and Corbel station (close to the Kongsfjord) was low. On 1 May 2007, a drying process in the boundary layer was observed. In the morning, the atmosphere up to Zeppelin Mountain showed enhanced values of the backscatter coefficient. Around noon, the top of the highly reflecting boundary layer decreased from 350 to 250 m. The top of the boundary layer observed by lidar was confirmed by radiosonde data.     

Seismic vulnerability assessment of a mixed masonry–RC building aggregate by linear and nonlinear analyses

From the beginning of the twentieth century, and due to the rapid increase of reinforced concrete (RC) usage, mixed masonry–RC buildings have emerged. In Lisbon, Portugal, old mixed masonry–RC buildings appeared between 1920 and 1960, representing the transition period between masonry and proper RC. These buildings are often integrated in blocks, and frequently share the side-walls, implying, thus, the need to assess the seismic vulnerability of building aggregates. The present paper approaches the seismic vulnerability assessment of a specific type of old mixed masonry–RC buildings in Lisbon. The study comprises the analysis of a building, both as an isolated structure and inserted in its aggregate, using two approaches: (1) linear dynamic analysis with SAP2000 and (2) nonlinear static analysis by means of 3Muri/Tremuri software. A comparison of both approaches derives a good matching between the obtained results. However, a nonlinear analysis is required to identify, in an adequate manner, the critical areas of the structure requiring strengthening.     

Interaction of geometry and mechanical property of trapezoidal sedimentary basins with incident <Emphasis Type="Italic">SH</Emphasis> waves

This paper investigated the effects of basin geometry and material property on the response of 2D trapezoidal sediment-filled basin to incident plane SH waves. Ten basin configurations with different geometries were developed, and then their seismic responses to both Ricker wavelets and seismic records were simulated by using an explicit finite difference scheme. The definition of deep/shallow basin, the precondition for the observation of prominent surface waves and the influential area of edge effects of the shallow basin were discussed quantitatively in this study. The followings were concluded: in the common velocity contrast range (v _s1/v _s2 < 10), the fundamental frequency a basin with W/H > 3.0 can be estimated approximately by 1D theory. The complexity of peak ground acceleration distribution pattern, the width of the most affected section as well as the amplitude of ground motion in the Edge Region increase with incident frequency. Prominent surface waves can only be observed when the incident wavelength is shorter than the critical wavelength λ _c . The interaction between incident wave and basin dynamic property plays a dominant role on the peak ground acceleration amplitude while the interaction between incident wave and geometry plays a more significant role on the peak ground acceleration distribution. For very shallow basin, different areas along the basin width are affected to different extents. Only a limited area close to the basin edge is influenced significantly. It is more feasible to propose spectral aggravation factor for different surface zones respectively than a uniform constant as a tool to calibrate the 1D-based design spectrum so as to take the basin effects into account.     

Modelling Freshwater Resources at the Global Scale: Challenges and Prospects

Quantification of spatially and temporally resolved water flows and water storage variations for all land areas of the globe is required to assess water resources, water scarcity and flood hazards, and to understand the Earth system. This quantification is done with the help of global hydrological models (GHMs). What are the challenges and prospects in the development and application of GHMs? Seven important challenges are presented. (1) Data scarcity makes quantification of human water use difficult even though significant progress has been achieved in the last decade. (2) Uncertainty of meteorological input data strongly affects model outputs. (3) The reaction of vegetation to changing climate and CO₂ concentrations is uncertain and not taken into account in most GHMs that serve to estimate climate change impacts. (4) Reasons for discrepant responses of GHMs to changing climate have yet to be identified. (5) More accurate estimates of monthly time series of water availability and use are needed to provide good indicators of water scarcity. (6) Integration of gradient-based groundwater modelling into GHMs is necessary for a better simulation of groundwater–surface water interactions and capillary rise. (7) Detection and attribution of human interference with freshwater systems by using GHMs are constrained by data of insufficient quality but also GHM uncertainty itself. Regarding prospects for progress, we propose to decrease the uncertainty of GHM output by making better use of in situ and remotely sensed observations of output variables such as river discharge or total water storage variations by multi-criteria validation, calibration or data assimilation. Finally, we present an initiative that works towards the vision of hyperresolution global hydrological modelling where GHM outputs would be provided at a 1-km resolution with reasonable accuracy.     

Earth’s Electromagnetic Environment

The natural spectrum of electromagnetic variations surrounding Earth extends across an enormous frequency range and is controlled by diverse physical processes. Electromagnetic (EM) induction studies make use of external field variations with frequencies ranging from the solar cycle which has been used for geomagnetic depth sounding through the 10\(^{-4}\)–10\(^4\) Hz frequency band widely used for magnetotelluric and audio-magnetotelluric studies. Above 10\(^4\) Hz, the EM spectrum is dominated by man-made signals. This review emphasizes electromagnetic sources at \(\sim\)1 Hz and higher, describing major differences in physical origin and structure of short- and long-period signals. The essential role of Earth’s internal magnetic field in defining the magnetosphere through its interactions with the solar wind and interplanetary magnetic field is briefly outlined. At its lower boundary, the magnetosphere is engaged in two-way interactions with the underlying ionosphere and neutral atmosphere. Extremely low-frequency (3 Hz–3 kHz) electromagnetic signals are generated in the form of sferics, lightning, and whistlers which can extend to frequencies as high as the VLF range (3–30 kHz).The roughly spherical dielectric cavity bounded by the ground and the ionosphere produces the Schumann resonance at around 8 Hz and its harmonics. A transverse resonance also occurs at 1.7–2.0 kHz arising from reflection off the variable height lower boundary of the ionosphere and exhibiting line splitting due to three-dimensional structure. Ground and satellite observations are discussed in the light of their contributions to understanding the global electric circuit and for EM induction studies.     

Present-Day Crustal Vertical Motion Around the Ordos Block Constrained by Precise Leveling and GPS Data

Precise leveling data observed in the period of 1970–2014 around the Ordos block were collected and processed to estimate present-day crustal vertical movement. Vertical rates of 6 GPS sites were employed as a priori constraints to define the reference frame. The velocity field shows that the interior of the Ordos block moves upward at a rate of 3 mm/a as a stable block. With respect to the central Ordos, the grabens and rifts around the Ordos block are undergoing subsidence, while the northeastern and southwestern Ordos uplift at the average rates of 2 and 1 mm/a, respectively. To the southeastern margin of the Ordos block, the Weihe and southern Shanxi grabens are subsiding at the rates of 4–6 mm/a. The subsidence of the Shanxi graben indicates that the graben is experiencing extensional movement on a long timescale. To the northwestern margin of the Ordos block, the Hetao and Yinchuan rifts are subsiding at the rates of 2–3 mm/a. A 2-D buried faulting model is used to infer the normal or reverse dip-slip rates. Our solution shows that most of the normal slip rates along the faults in the grabens and rifts are ~2 mm/a.     

Source parameters for the 2013–2015 earthquake sequence in Nógrád county,Hungary

Between 2013 June and 2015 January, 35 earthquakes with local magnitude M _L ranging from 1.1 to 4.2 occurred in Nógrád county, Hungary. This earthquake sequence represents above average seismic activity in the region and is the first one that was recorded by a significant number of three-component digital seismographs in the county. Using a Bayesian multiple-event location algorithm, we have estimated the hypocenters of 30 earthquakes with M _L ≥1.5. The events occurred in two small regions of a few squared kilometers: one to the east of Érsekvadkert and the other at Iliny. The uncertainty of the epicenters is about 1.5–1.7 km in the E-W direction and 1.8–2.1 km in the N-S direction at the 95 % confidence level. The estimated event depths are confined to the upper 3 km of the crust. We have successfully estimated the full moment tensors of 4 M _w ≥3.6 earthquakes using a probabilistic waveform inversion procedure. The non-double-couple components of the retrieved moment tensor solutions are statistically insignificant. The negligible amount of the isotropic component implies the tectonic nature of the investigated events. All of the analyzed earthquakes have strike-slip mechanism with either right-lateral slip on an approximately N-S striking or left-lateral movement on a roughly E-W striking nodal plane. The orientations of the obtained focal mechanisms are in good agreement with the main stress pattern published for the epicentral region. Both the P and T principal axes are horizontal, and the P axis is oriented along a NE-SW direction.     

Earthquake dynamics on circular faults: a review 1970–2015

We review the development of earthquake dynamics taken from the point of view of the origin of seismic radiation instead of the detailed study of rupture propagation on complex surfaces. Many features of seismic radiation can be explained by simple models that serve as elementary canonical problems. Some of these properties are very well known like the fact that at low frequencies, the seismic spectrum is proportional to the seismic moment. At high frequencies, on the other hand, radiation is generated by the motion of the rupture front, in particular stopping phases and geometrical obstacles (barriers). A rupture front moving at constant speed does not produce far-field radiation. For many practical applications, for determining source size and stress drop, for example, it is not necessary to determine geometrical details of the rupture. For such cases, a simple circular crack model is quite sufficient. An improvement on this method is to do dynamic inversion on simple, elliptical-shaped sources and letting the rupture start arbitrarily from a point on the fault. This problem can be solved nowadays with finite differences and a variety of inversion techniques.     

Three-dimensional admittance analysis of lithospheric elastic thickness over the Louisville Ridge

Using bathymetry and altimetric gravity anomalies, a 1° 9 1° lithospheric effective elastic thickness(Te) model over the Louisville Ridge and its adjacent regions is calculated using the moving window admittance technique. For comparison, three bathymetry models are used: general bathymetric charts of the oceans, SIO V15.1,and BAT_VGG. The results show that BAT_VGG is more suitable for calculating T e than the other two models. T e along the Louisville Ridge was re-evaluated. The southeast of the ridge has a medium Te of 10–20 km, while Te increases dramatically seaward of the Tonga-Kermadec trench as a result of the collision of the Pacific and IndoAustralian plates.