The seismological characteristics of the 15 February 2013 Chelyabinsk bolide explosion are investigated based on seismograms recorded at 50 stations with epicentral distances ranging from 229 to 4324 km. By using 8–25 s vertical-component Rayleigh waveforms, we obtain a surface-wave magnitude of 4.17±0.31 for this event. According to the relationship among the Rayleigh-wave magnitude, burst height and explosive yield, the explosion yield is estimated to be 686 kt. Using a single-force source to fit the observed Rayleigh waveforms, we obtain a single force of 1.03×1012 N, which is equivalent to the impact from the shock wave generated by the bolide explosion.
The medium-small earthquakes that occurred in the middle part of Tibetan Plateau (32°N–36°N, 90°E–93°E) from August 2016 to June 2017 were relocated using the absolute earthquake location method Hypo2000. Compared to the reports of Chinese Seismological Networks, our relocation results are more clustered on the whole, the horizontal location differences exceed 10 km, and the focal depths are concentrated in 0–8 km, which indicates that the upper crust inside the Tibetan Plateau is tectonically active. In June 2017 altogether eight earthquakes above magnitude 3.0 took place; their relocated epicenters are concentrated around Gêladaindong. The relocation results of M<3.0 small earthquakes also showed obvious differences. Therefore, we used the CAP method to invert for the focal mechanisms of theM ≥3.0 earthquakes; results generally tally with the surface geological structures, indicating that the Tibetan Plateau is still under the strong compressional force from the India Plate. Among them the eight earthquakes that occurred near Gêladaindong in June 2017 are all of normal fault type or with some strike-slip at the same time; based on previous research results we conjecture that these events are intense shallow crust responses to deep crust-mantle activities.
Considering the uncertainty of the electrical axis for two-dimensional audo-magnetotelluric (AMT) data processing, an AMT inversion method with the Central impedance tensor was presented. First, we present a calculation expression of the Central impedance tensor in AMT, which can be considered as the arithmetic mean of TE-polarization mode and TM-polarization mode in the two-dimensional geo-electrical model. Second, a least-squares iterative inversion algorithm is established, based on a smoothness-constrained model, and an improved L-curve method is adopted to determine the best regularization parameters. We then test the above inversion method with synthetic data and field data. The test results show that this two-dimensional AMT inversion scheme for the responses of Central impedance is effective and can reconstruct reasonable two-dimensional subsurface resistivity structures. We conclude that the Central impedance tensor is a useful tool for two-dimensional inversion of AMT data.
Exohiss is a low-frequency structureless whistler-mode emission potentially contributing to the precipitation loss of radiation belt electrons outside the plasmasphere. Exohiss is usually considered the plasmaspheric hiss leaked out of the dayside plasmapause. However, the evolution of exohiss after the leakage has not been fully understood. Here we report the prompt enhancements of exohiss waves following substorm injections observed by Van Allen Probes. Within several minutes, the energetic electron fluxes around 100 keV were enhanced by up to 5 times, accompanied by an up to 10-time increase of the exohiss wave power. These substorm-injected electrons are shown to produce a new peak of linear growth rate in the exohiss band (< 0.1fce). The corresponding path-integrated growth rate of wave power within 10° latitude of the magnetic equatorial plane can reach 13.4, approximately explaining the observed enhancement of exohiss waves. These observations and simulations suggest that the substorm-injected energetic electrons could amplify the preexisting exohiss waves.
Climate system models are useful tools for understanding the interactions among the components of the climate system and predicting/projecting future climate change. The development of climate models has been a central focus of the State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences (LASG/IAP) since the establishment of the laboratory in 1985. In China, many pioneering component models and fully coupled models of the climate system have been developed by LASG/IAP. The fully coupled climate system developed in the recent decade is named FGOALS (Flexible Global Ocean-Atmosphere-Land System Model). In this paper, an application-oriented review of the LASG/IAP FGOALS model is presented. The improved model performances are demonstrated in the context of cloud-radiation processes, Asian monsoon, ENSO phenomena, Atlantic Meridional Overturning Circulation (AMOC) and sea ice. The FGOALS model has contributed to both CMIP5 (Coupled Model Intercomparison Project-phase 5) and IPCC (Intergovernmental Panel on Climate Change) AR5 (the Fifth Assessment Report). The release of FGOALS data has supported the publication of nearly 500 papers around the world. The results of FGOALS are cited ~106 times in the IPCC WG1 (Working Group 1) AR5. In addition to the traditional long-term simulations and projections, near-term decadal climate prediction is a new set of CMIP experiment, progress of LAGS/IAP in the development of near-term decadal prediction system is reviewed. The FGOALS model has supported many Chinese national-level research projects and contributed to the national climate change assessment report. The crucial role of FGOALS as a modeling tool for supporting climate sciences is highlighted by demonstrating the model’s performances in the simulation of the evolution of Earth’s climate from the past to the future.
We present preliminary results of a new global Magnetohydrodynamics (MHD) simulation model of the Jovian magnetosphere. The model incorporates mass loading from Jupiter's satellite Io, the planet's fast corotation, and electrostatic coupling between its magnetosphere and ionosphere (M-I coupling). The basic configuration of the Jovian magnetosphere including the equatorial plasma flow pattern, the corotation enforcement current system, and the field aligned currents (FACs) in the ionosphere are presented under an antiparallel interplanetary magnetic field (IMF) condition. The simulation model results for equatorial density and pressure profiles are consistent with results from data-based empirical models. It is also found that there are similarities between the FACs distribution in the ionosphere and the observed aurora features, showing the potential application of the simple ionospheric model to the complicated M-I coupling. This model will help deepen our understanding of the global dynamics of the Jovian magnetosphere.
On November 13, 2016, an MW7.8 earthquake struck Kaikoura in South Island of New Zealand. By means of back-projection of array recordings, ASTFs-analysis of global seismic recordings, and joint inversion of global seismic data and co-seismic InSAR data, we investigated complexity of the earthquake source. The result shows that the 2016 MW7.8 Kaikoura earthquake ruptured about 100 s unilaterally from south to northeast (~N28°–33°E), producing a rupture area about 160 km long and about 50 km wide and releasing scalar moment 1.01×1021 Nm. In particular, the rupture area consisted of two slip asperities, with one close to the initial rupture point having a maximal slip value ~6.9 m while the other far away in the northeast having a maximal slip value ~9.3 m. The first asperity slipped for about 65 s and the second one started 40 s after the first one had initiated. The two slipped simultaneously for about 25 s. Furthermore, the first had a nearly thrust slip while the second had both thrust and strike slip. It is interesting that the rupture velocity was not constant, and the whole process may be divided into 5 stages in which the velocities were estimated to be 1.4 km/s, 0 km/s, 2.1 km/s, 0 km/s and 1.1 km/s, respectively. The high-frequency sources distributed nearly along the lower edge of the rupture area, the high-frequency radiating mainly occurred at launching of the asperities, and it seemed that no high-frequency energy was radiated when the rupturing was going to stop.
The potential for devastating earthquakes in the Himalayan orogeny has long been recognized. The 2015 MW7.8 Gorkha, Nepal earthquake has heightened the likelihood that major earthquakes will occur along this orogenic belt in the future. Reliable seismic hazard assessment is a critical element in development of policy for seismic hazard mitigation and risk reduction. In this study, we conduct probabilistic seismic hazard assessment using three different seismogenic source models (smoothed gridded, linear, and areal sources) based on the complicated tectonics of the study area. Two sets of ground motion prediction equations are combined in a standard logic tree by taking into account the epistemic uncertainties in hazard estimation. Long-term slip rates and paleoseismic records are also incorporated in the linear source model. Peak ground acceleration and spectral acceleration at 0.2 s and 1.0 s for 2% and 10% probabilities of exceedance in 50 years are estimated. The resulting maps show significant spatial variation in seismic hazard levels. The region of the Lesser Himalaya is found to have high seismic hazard potential. Along the Main Himalayan Thrust from east to west beneath the Main Central Thrust, large earthquakes have occurred regularly in history; hazard values in this region are found to be higher than those shown on existing hazard maps. In essence, the combination of long span earthquake catalogs and multiple seismogenic source models gives improved seismic hazard constraints in Nepal.
Eolian dust preserved in deep-sea sediments of the North Pacific Ocean (NPO) is an important recorder of paleoclimatic and paleoenvironmental changes in the Asian inland. To better understand changes in the dust provenances, in this study diffuse reflectance spectroscopy (DRS) was used to extract the eolian signal recorded in sediments of ODP Hole 885A recovered from the NPO. First, we systematically investigated sieving effects on the DRS data; then band positions of hematite (obtained from the second order derivative curves of the K-M remission function spectrum derived from the DRS) were used to distinguish different provenances of the eolian dust preserved in the pelagic sediments of this hole. Our results show that the sieving (38 μm) process can suppress effectively the experimental errors. Eolian signatures from Chinese Loess Plateau (CLP) sources and non-CLP-sources have been identified in the pelagic sediments of ODP Hole 885A from the late Pliocene to the early Pleistocene. The provenance differences account for the discrepancies in the eolian records recovered from the pelagic sediments in the NPO and profiles in the CLP. Temporal changes in dust provenances are caused by the latitudinal movement of the westerly jet mainstream. The hematite DRS band position is a useful tool to distinguish the provenance of eolian components preserved in pelagic sediments.