Earthquake-induced gravity variation refers to changes in the earth’s gravity field associated with seismic activities.In recent years,development in the theories has greatly promoted seismic deformation research,lay...Earthquake-induced gravity variation refers to changes in the earth’s gravity field associated with seismic activities.In recent years,development in the theories has greatly promoted seismic deformation research,laying a solid theoretical foundation for the interpretation and application of seismological gravity monitoring.Traditional terrestrial gravity measurements continue to play a significant role in studies of interseismic,co-seismic,and post-seismic gravity field variations.For instance,superconducting gravimeter networks can detect co-seismic gravity change at the sub-micro Gal level.At the same time,the successful launch of satellite gravity missions(e.g.,the Gravity Recovery and Climate Experiment or GRACE)has also facilitated applied studies of the gravity variation associated with large earthquakes,and several remarkable breakthroughs have been achieved.The progress in gravity observation technologies(e.g.,GRACE and superconducting gravimetry)and advances in the theories have jointly promoted seismic deformation studies and raised many new research topics.For example,superconducting gravimetry has played an important role in analyses of episodic tremor,slow-slip events,and interseismic strain patterns;the monitoring of transient gravity signals and related theories have provided a new perspective on earthquake early warning systems;the mass transport detected by the GRACE satellites several months before an earthquake has brought new insights into earthquake prediction methods;the use of artificial intelligence to automatically identify tiny gravity change signals is a new approach to accurate and rapid determination of earthquake magnitude and location.Overall,many significant breakthroughs have been made in recent years,in terms of the theory,application,and observation measures.This article summarizes the progress,with the aim of providing a reference for seismologists and geodetic researchers studying the phenomenon of gravity variation,advances in related theories and applications,and future research directions in this discipline.展开更多
The Earth’s surface kinematics and deformation are fundamental to understanding crustal evolution.An effective research approach is to estimate regional motion field and deformation fields based on modern geodetic ne...The Earth’s surface kinematics and deformation are fundamental to understanding crustal evolution.An effective research approach is to estimate regional motion field and deformation fields based on modern geodetic networks.If the discrete observed velocity field is obtained,the velocity related fields,such as dilatation rate and maximum shear strain rate,can be estimated by applying varied mathematical approaches.This study applied Akaike's Bayesian Information Criterion(ABIC)method to calculate strain rate fields constrained by GPS observations in the southeast Tibetan Plateau.Comparison with results derived from other three methods revealed that our ABIC-derived strain rate fields were more precise.The maximum shear strain rate highlighted the Xianshuihe–Xiaojiang fault system as the main boundary for the outward migration of material in southeastern Tibet,indicating rotation of eastern Tibet material around the eastern Himalaya rather than whole extrusion along a fixed channel.Additionally,distinct dilatation rate patterns in the northeast and southwest regions of the fault system were observed.The northeast region,represented by the Longmenshan area,exhibited negative dilatational anomalies;while the southwest region,represented by the Jinsha River area north of 29°N,displayed positive dilatational anomalies.This indicates compression in the former and extension in the latter.Combined with deep geophysical observations,we believe that the upper and lower crusts of the Jinsha River area north of 29°N are in an entire expanding state,probably caused by the escape-drag effect of material.The presence of a large,low-viscosity region south of 29°N may not enable the entire escape of the crust,but instead result in a differential escape of the lower crust faster than the upper crust.展开更多
Using a new vortex detection and tracing method, a dataset of the Southwest Vortex(SWV) is established based on Japanese 25-year Reanalysis(JRA-25) reanalysis data during 1979–2008. The spatiotemporal features of the...Using a new vortex detection and tracing method, a dataset of the Southwest Vortex(SWV) is established based on Japanese 25-year Reanalysis(JRA-25) reanalysis data during 1979–2008. The spatiotemporal features of the SWV are derived from the dataset. In comparison to other seasons, summer yields the least SWVs, but with the highest probability that they will migrate from their region of origin. SWVs mostly emerge in the southwest of the Sichuan Basin and the southeast of the Tibetan Plateau. Migratory SWVs mainly move along either an eastward or southeastward path. Detailed composite analysis of warm-season SWVs shows that the subtropical high is a key factor in determining the direction of migratory SWVs. Furthermore, the steering wind at 700 hPa dominates the moving direction of migratory SWVs. Potential stability diagnosed by pseudo-equivalent potential temperature ? se is of certain significance for the evolution and movement of SWVs. On the other hand, migratory SWVs possess relatively greater strength than stationary SWVs, due to a stronger low-level jet with enhanced baroclinicity and moisture transport providing more energy to support the growth of SWVs along their paths of movement.展开更多
Data obtained by GRACE(Gravity Recovery and Climate Experiment) have been used to invert for the seismic source parameters of megathrust earthquakes under the assumption of either uniform slip over an entire fault or ...Data obtained by GRACE(Gravity Recovery and Climate Experiment) have been used to invert for the seismic source parameters of megathrust earthquakes under the assumption of either uniform slip over an entire fault or a point-like seismic source.Herein, we further extend the inversion of GRACE long-wavelength gravity changes to heterogeneous slip distributions during the 2011 Tohoku earthquake using three fault models:(Ⅰ) a constant-strike and constant-dip fault,(Ⅱ) a variable dip fault, and(Ⅲ) a realistically varying strike fault. By removing the post-seismic signal from the time series, and taking the effect of ocean water redistribution into account, we invert for slip models I, II, and III using co-seismic gravity changes measured by GRACE, de-striped by DDK3 decorrelation filter. The total seismic moments of our slip models, with respective values of 4.9×10^(22) Nm, 5.1×10^(22) Nm, and 5.0×10^(22) Nm, are smaller than those obtained by other studies relying on GRACE data. The resulting centroids are also located at greater depths(20 km, 19.8 km,and 17.4 km, respectively). By combining onshore GPS, GPS-Acoustic, and GRACE data, we obtain a jointly inverted slip model with a seismic moment of 4.8×10^(22) Nm, which is larger than the seismic moment obtained using only the GPS displacements. We show that the slip inverted from low degree space-borne gravimetric data, which contains information at the ocean region, is affected by the strike of the arcuate trench. The space-borne gravimetric data help us constrain the source parameters of a megathrust earthquake within the frame of heterogeneous slip models.展开更多
Modern geodetic techniques have developed rapidly in recent years, providing reliable observation data and new effective approaches, and greatly enhancing studies of the Tibetan geodynamics. For instance, the well-kno...Modern geodetic techniques have developed rapidly in recent years, providing reliable observation data and new effective approaches, and greatly enhancing studies of the Tibetan geodynamics. For instance, the well-known GPS technique has been employed to measure seismic slips for many faults in the Tibetan Plateau. GPS data agree well with the hypothesis of a thickening crust and eastward mass flow. Moreover, absolute gravimetric data have been applied to interpret geophysical phenomena such as crust movement, co-seismic gravity change, GIA, and ground water change. The satellite gravity mission GRACE launched in 2002 provided global gravity models with unprecedentedly high precision and high spatial resolution. It has been used in implementing temporal gravity changes and improving our knowledge of the Earth's interior, including lithosphere dynamics, mantle viscosity and rheology, plateau uplift, and subduction processing. It is noteworthy that gravity presents unique advantages for the study of Tibetan geodynamics because of its sensitivity to mass migration and dynamic redistribution. To date, great advances have been made in applying modern geodetic data in studying dynamic changes of Tibetan plateau. For instance, the horizontal displacement field from GPS data revealed dynamical characteristics of the present-day Tibetan plateau. The combination of gravity anomalies and topographic data describe the tectonic characteristics of Tibetan plateau. The combination of gravity data and GPS data show present properties of the Tibetan plateau such as crust thickening, Moho's subsidence, and plateau uplift. GRACE data were used to estimate the distribution of ice/snow melting. These results demonstrate that mere application of integrated geodetic data as well as geophysical methods and numerical simulations can enhance our knowledge of Tibetan plateau dynamics. It must be pointed out that GRACE data include various geophysical signals such as crust vertical movement, denudation, ice and snow melting, GIA, ground water change, and permafrost degradation. To separate the tectonic information from other impulses, each physical signal must be evaluated and corrected carefully from the GRACE data. The Tibetan geodynamic problem is a complicated and synthetic issue that must be addressed through collaboration of workers in many fields. Succinctly put, although great achievements have been made in studying Tibetan plateau dynamics from each field, the dynamical process remains unclear. Some fundamental problems remain unresolved. They should be solved with modern geodetic data, such as GRACE, GPS, and absolute gravity data, combined with meteorological and geological data, for quantitative analysis of Tibetan plateau dynamics affected by respective geophysical sources. This review article introduces and discusses the scientific importance, advances, problems, and prospects of modern geodesy applied to the study of geodynamic changes of the Tibetan plateau.展开更多
Using the Simple Biosphere Model(SiB2),soil thermal properties(STP) were examined in a Tibetan prairie during the monsoon period to investigate ground surface temperature prediction.We improved the SiB2 model by incor...Using the Simple Biosphere Model(SiB2),soil thermal properties(STP) were examined in a Tibetan prairie during the monsoon period to investigate ground surface temperature prediction.We improved the SiB2 model by incorporating a revised force-restore method(FRM) to take the vertical heterogeneity of soil thermal diffusivity(k) into account.The results indicate that(1) the revised FRM alleviates daytime overestimation and nighttime underestimation in modeled ground surface temperature(Tg),and(2) its role in little rainfall events is significant because the vertical gradient of k increases with increasing surface evaporation.Since the original formula of thermal conductivity(λ) in the SiB2 greatly underestimates soil thermal conductivity,we compared five algorithms of λ involving soil moisture to investigate the cause of overestimation during the day and underestimation at night on the basis of the revised FRM.The results show that(1) the five algorithms significantly improve Tg prediction,especially in daytime,and(2) taking one of these five algorithms as an example,the simulated Tg values in the daytime are closer to the field measurements than those in the nighttime.The differences between modeled Tg and field measurements are mostly within the margin of error of ±2 K during 3 August to 4 September 1998.展开更多
Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the...Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the passive source surface-wave method is introduced into the site engineering seismic survey and practically applied in an engineering site of Shijingshan District. By recording the ubiquitous weak vibration on the earth surface, extract the dispersion curve from the surface-wave components using the SPAC method and obtain the shear-wave velocity structure from inversion. Over the depth of 42 m underground, it totally consists of five layers with interface depth of 3.31, 4.50, 7.23, 17.41, and 42.00 m; and shear-wave velocity of 144.0, 198.3, 339.4, 744.2, and 903.7 m/s, respectively. The inversion result is used to evaluate site classification, determine the maximum shear modulus of soil, provide basis for further seismic hazard analysis and site assessment or site zoning, etc. The result shows that the passive source surface-wave method is feasible in the site engineering seismic survey and can replace boreholes,shorten survey period, and reduce engineering cost to some extent.展开更多
The fracturing technology for shale gas reservoir is the key to the development of shale gas industrialization.It makes much sense to study the mechanical properties and deformation characteristics of shale,due to its...The fracturing technology for shale gas reservoir is the key to the development of shale gas industrialization.It makes much sense to study the mechanical properties and deformation characteristics of shale,due to its close relationship with the fracability of shale gas reservoir.This paper took marine shale in the Changning area,southern Sichuan Basin of China as the research object.Based on field profile and hand specimen observation,we analyzed the development of natural fractures and collected samples from Wufeng Formation and Longmaxi Formation.Combining with the indoor experiment,we investigated the macroscopic and microscopic structural features and the remarkable heterogeneity of shale samples.Then we illustrated the mechanics and deformation characteristics of shale,through uniaxial compression test and direct shear test.The shale has two types of fracture modes,which depend on the angular relation between loading direction and the bedding plane.Besides,the Wufeng shale has a higher value of brittleness index than the Longmaxi shale,which was calculated using two methods,mechanical parameters and mineral composition.Given the above results,we proposed a fracability evaluation model for shale gas reservoir using the analytic hierarchy process.Four influence factors,brittleness index,fracture toughness,natural fractures and cohesive force,are considered.Finally,under the control of normalized value and weight coefficient of each influence factor,the calculations results indicate that the fracability index of the Wufeng Formation is higher than that of the Longmaxi Formation in Changning area,southern Sichuan Basin.展开更多
The boundary value problem of deflections of vertical with ellipsoid boundary is studied in the paper.Based on spherical harmonic series,the ellipsoidal corrections for the boundary value problem are derived so that i...The boundary value problem of deflections of vertical with ellipsoid boundary is studied in the paper.Based on spherical harmonic series,the ellipsoidal corrections for the boundary value problem are derived so that it can be well solved.In addition,an imitation arithmetic is given for examining the accuracies of solutions for the boundary value problem as well as its spherical approximation problem,and the computational results illustrate that the boundary value problem has higher accuracy than its spherical approximation problem if deflection of the vertical are measured on geoid.展开更多
The Yellowstone volcano is one of the largest active volcanoes in the world,and its potential hazards demand detailed seismological and geodetic studies.Previous studies with travel time tomography and receiver functi...The Yellowstone volcano is one of the largest active volcanoes in the world,and its potential hazards demand detailed seismological and geodetic studies.Previous studies with travel time tomography and receiver functions have revealed a low-velocity layer in the crust beneath the Yellowstone volcano,suggesting the presence of a magma chamber at depth.We use ambient seismic noise from regional seismic stations to retrieve short-period surface waves and then study the shallow shear velocity structure of the Yellowstone region by surface wave dispersion analysis.We first obtained a crustal model of the area outside of the Yellowstone volcano and then constructed an absolute shear wave velocity structure in combination with receiver function results for the crust beneath the Yellowstone volcano.The velocity model shows a low-velocity layer with shear velocity at around1.3 km/s,suggesting that a large-scale magma chamber exists at shallow levels within the crust of the Yellowstone volcanic region.展开更多
In this paper,we study how coseismic deformations calculated in 1066 Earth models are affected by how the models treat Earth discontinuities.From the results of applying models 1066A(continuous)and 1066B(discontinuous...In this paper,we study how coseismic deformations calculated in 1066 Earth models are affected by how the models treat Earth discontinuities.From the results of applying models 1066A(continuous)and 1066B(discontinuous),we find that the difference in Love numbers of strike-slip and horizontal tensile sources are bigger than dip-slip and vertical tensile sources.Taken collectively,discontinuities have major effects on Green’s functions of four independent sources.For the near-field coseismic deformations of the 2013 Okhotsk earthquake(Mw 8.3),the overall differences between theoretical calculations in vertical displacement,geoid,and gravity changes caused by discontinuities are 10.52 percent,9.07 percent and 6.19 percent,with RMS errors of 0.624 mm,0.029 mm,and 0.063μGal,respectively.The difference in far-field displacements is small,compared with GPS data,and we can neglect this effect.For the shallow earthquake,2011 Tohoku-Oki earthquake(Mw 9.0),the differences in near-field displacements are 0.030 m(N-S),0.093 m(E-W),and 0.025 m(up-down)in our study area with the ARIA slip model,which gives results closer to GPS data than those from the USGS model.The difference in vertical displacements and gravity changes on the Earth’s surface caused by discontinuities are larger than 10 percent.The difference in the theoretical gravity changes at spatially fixed points truncated to degrees 60,as required by GRACE data,is 0.0016μGal and the discrepancy is 11 percent,with the theoretical spatial gravity changes from 1066B closer to observations than from 1066A.The results show that an Earth model with discontinuities in the medium has a large effect on the calculated coseismic deformations.展开更多
The complex geographical environment in China makes its gravity signals miscellaneous.This work gives a comprehensive representation and explanation in secular trend of gravity change in different regions,the key feat...The complex geographical environment in China makes its gravity signals miscellaneous.This work gives a comprehensive representation and explanation in secular trend of gravity change in different regions,the key features of which include positive trend in inner Tibet Plateau and South China and negative trend in North China plain and high mountain Asia(HMA).We also present the patterns of amplitudes and phases of annual and semiannual change.The mechanism underlying the semiannual period is explicitly discussed.The displacement in three directions expressed in terms of geo-potential spherical coefficients and load Love numbers are given.A case study applied with these equations is presented.The results show that Global Positioning System(GPS) observations can be used to compare with Gravity Recovery and Climate Experiment(GRACE) derived displacement and the vertical direction has a signal-noise-ratio of about one order of magnitude larger than the horizontal directions.展开更多
The distribution of Coulomb failure stress (CFS) change in the steep excavation slope is calculated by finite element method in this paper,and the failure mechanics under different conditions have been investigated.Co...The distribution of Coulomb failure stress (CFS) change in the steep excavation slope is calculated by finite element method in this paper,and the failure mechanics under different conditions have been investigated.Comparing the CFSs before and after the slope excavation (stress loading and unloading processes),the dangerous internal zone and the most likely failure external area are attained.Given the shear cracks on the top surface while tensile stress or cracks along the toe of the slope,we analyze the high cutting-angle steep slope in Kaixian county of the Three Gorges Reservoir region.We bring forward that the peak value of CFS after excavation can reach to the order of 0.1MPa,which is greatly higher than that of before.Our preliminary results are useful for optimizing the reinforcement structure during the steep slope stabilization engineering.展开更多
Nonlinear force-free magnetic field(NLFFF) extrapolation based on the observed photospheric magnetic field is the most important method to obtain the coronal magnetic field nowadays.However, raw photospheric magnetogr...Nonlinear force-free magnetic field(NLFFF) extrapolation based on the observed photospheric magnetic field is the most important method to obtain the coronal magnetic field nowadays.However, raw photospheric magnetograms contain magnetic forces and small-scale noises, and fail to be consistent with the force-free assumption of NLFFF models. The procedure for removing the forces and noises in observed data is called preprocessing. In this paper, we extend the preprocessing code of Jiang& Feng to spherical coordinates for a full sphere. We first smooth the observed data with Gaussian smoothing, and then split the smoothed magnetic field into a potential field and a non-potential field.The potential part is computed by a numerical potential field model, and the non-potential part is preprocessed using an optimization method to minimize the magnetic forces and magnetic torques. Applying the code to synoptic charts of the vector magnetic field from SDO/HMI, we find it can effectively reduce the noises and forces, and improve the quality of data for a better input which will be used for NLFFF extrapolations applied to the global corona.展开更多
Global mean sea level(GMSL) change is one of the important indicators of global climate change and is a crucial scientific issue of continuing interest. As satellite altimeter data, the Gravity Recovery and Climate Ex...Global mean sea level(GMSL) change is one of the important indicators of global climate change and is a crucial scientific issue of continuing interest. As satellite altimeter data, the Gravity Recovery and Climate Experiment(GRACE) and Argo continue to be updated, especially with the release of GRACE Follow-On(GRACE-FO) data, making it necessary to combine these latest data to estimate sea level change. Determinations on whether the GRACE and GRACE-FO observation systems provide unbiased global observation data have not been effectively evaluated. Therefore, this research mainly investigated the consistency of GRACE and GRACE-FO observation data in studying GMSL change. By comparing the sum of the GMSL calculated by the two gravity satellites and Argo data with the GMSL calculated by satellite altimeters, the discrepancy between GRACE-FO + Argo and satellite altimeter data is about 7.9 ± 2.3 mm, which is likely derived from the inconsistency between GRACE and GRACE-FO data.展开更多
Under the limited cultivated land area and the pursuit of sustainable agricultural development,it is essential for the safety of grain production to study agricultural management approaches on narrowing the winter whe...Under the limited cultivated land area and the pursuit of sustainable agricultural development,it is essential for the safety of grain production to study agricultural management approaches on narrowing the winter wheat yield gap and improving nitrogen use efficiency (NUE) in China.In this study,DSSAT-CERES-Wheat Model is used to simulate winter wheat yield under different agricultural treatments,and we analyze yield gaps and NUE with different management scenarios at regional scales and evaluate the suitable approaches for reducing yield gap and increasing NUE.The results show that,the potential of narrowing yield gap ranges 300–900 kg ha^(–1) with soil nutrients increase,400–1 200 kg ha^(–1) with sowing date adjustment and 0–400 kg ha^(–1) with planting density increase as well as 700–2 200 kg ha^(–1) with adding nitrogen fertilizer.Contribution rates of management measures of soil nutrients,sowing date adjusting,planting density,and nitrogen fertilizers are 5–15%,5–15%,0–4%,and 10–20%,respectively.Difference in nitrogen partial productivity ranges 3–10 kg kg^(–1) for soil nutrients,1–10 kg kg^(–1) for sowing date adjusting,1–5 kg kg^(–1) for planting density increase,and–12–0 kg kg^(–1) for adding nitrogen fertilizers,respectively.It indicates that four treatments can narrow yield gap and improve the NUE in varying degrees,but increasing nitrogen fertilizer leads to the decrease of NUE.展开更多
Global seismicity catalogs are sufficient for characterizing double seismic zones(DSZs)in subducting slab and facilitate to estimate layer separation without inconsistent uncertainties as local catalogs.Previous studi...Global seismicity catalogs are sufficient for characterizing double seismic zones(DSZs)in subducting slab and facilitate to estimate layer separation without inconsistent uncertainties as local catalogs.Previous studies have shown the correlation between DSZs layer separation and plate age while correlation for those younger than ～60 Ma is suspicious.The lacking of DSZs with layer separation less than 10 km further makes it difficult to precisely estimate such correlation.Thus,we incorporate eight DSZs data determined through local seismicity into globally-determined dataset and reexamine such correlation.The best fitting results show that both a linear model and a square root of plate age can mathematically fit the layer separation well.However,it is difficult to distinguish these two models when plate age is greater than～20 Ma since their difference is less than 2 km.However,if extrapolation is possible,both models should provide physical information that DSZs will not form if there is no subducting lithosphere.As a result,the DSZs cannot be produced until the oceanic lithospheric age becomes greater than 0.9 Ma in the square root model while the linear model gives a misleading result.As such the square root model demonstrates the relationship physically better than the linear one,it still needs further test in the future with more available data,nevertheless,our study might also provide evidence for the suggestion that the plate age is a primary control factor of the DSZs geometry as well as the subducting process which disregards any local tectonic stresses.展开更多
The relationships among an ocean wave spectrum,a fully polarimetric coherence matrix,and radar parameters are deduced with an electromagnetic wave theory.Furthermore,the relationship between the polarimetric entropy a...The relationships among an ocean wave spectrum,a fully polarimetric coherence matrix,and radar parameters are deduced with an electromagnetic wave theory.Furthermore,the relationship between the polarimetric entropy and ocean wave spectrum is established based on the definition of entropy and a twoscale scattering model of the ocean surface.It is the first time that the polarimetric entropy of the ocean surface is presented in theory.Meanwhile,the relationships among the fully polarimetric entropy and the parameters related to radar and ocean are discussed.The study is the basis of further monitoring targets on the ocean surface and deriving oceanic information with the entropy from the ocean surface.The contrast enhancement between human-made targets and the ocean surface with the entropy is presented with quad-pol airborne synthetic aperture radar(AIRSAR)data.展开更多
文摘Earthquake-induced gravity variation refers to changes in the earth’s gravity field associated with seismic activities.In recent years,development in the theories has greatly promoted seismic deformation research,laying a solid theoretical foundation for the interpretation and application of seismological gravity monitoring.Traditional terrestrial gravity measurements continue to play a significant role in studies of interseismic,co-seismic,and post-seismic gravity field variations.For instance,superconducting gravimeter networks can detect co-seismic gravity change at the sub-micro Gal level.At the same time,the successful launch of satellite gravity missions(e.g.,the Gravity Recovery and Climate Experiment or GRACE)has also facilitated applied studies of the gravity variation associated with large earthquakes,and several remarkable breakthroughs have been achieved.The progress in gravity observation technologies(e.g.,GRACE and superconducting gravimetry)and advances in the theories have jointly promoted seismic deformation studies and raised many new research topics.For example,superconducting gravimetry has played an important role in analyses of episodic tremor,slow-slip events,and interseismic strain patterns;the monitoring of transient gravity signals and related theories have provided a new perspective on earthquake early warning systems;the mass transport detected by the GRACE satellites several months before an earthquake has brought new insights into earthquake prediction methods;the use of artificial intelligence to automatically identify tiny gravity change signals is a new approach to accurate and rapid determination of earthquake magnitude and location.Overall,many significant breakthroughs have been made in recent years,in terms of the theory,application,and observation measures.This article summarizes the progress,with the aim of providing a reference for seismologists and geodetic researchers studying the phenomenon of gravity variation,advances in related theories and applications,and future research directions in this discipline.
基金supported by grants from the Ministry of Science and Technology(Grant Nos.2021FY100101,2019QZKK0901)the National Natural Science Foundation of China(Grant Nos.41941016,42230312,42020104007)China Geological Survey(Grant No.DD20221630).
文摘The Earth’s surface kinematics and deformation are fundamental to understanding crustal evolution.An effective research approach is to estimate regional motion field and deformation fields based on modern geodetic networks.If the discrete observed velocity field is obtained,the velocity related fields,such as dilatation rate and maximum shear strain rate,can be estimated by applying varied mathematical approaches.This study applied Akaike's Bayesian Information Criterion(ABIC)method to calculate strain rate fields constrained by GPS observations in the southeast Tibetan Plateau.Comparison with results derived from other three methods revealed that our ABIC-derived strain rate fields were more precise.The maximum shear strain rate highlighted the Xianshuihe–Xiaojiang fault system as the main boundary for the outward migration of material in southeastern Tibet,indicating rotation of eastern Tibet material around the eastern Himalaya rather than whole extrusion along a fixed channel.Additionally,distinct dilatation rate patterns in the northeast and southwest regions of the fault system were observed.The northeast region,represented by the Longmenshan area,exhibited negative dilatational anomalies;while the southwest region,represented by the Jinsha River area north of 29°N,displayed positive dilatational anomalies.This indicates compression in the former and extension in the latter.Combined with deep geophysical observations,we believe that the upper and lower crusts of the Jinsha River area north of 29°N are in an entire expanding state,probably caused by the escape-drag effect of material.The presence of a large,low-viscosity region south of 29°N may not enable the entire escape of the crust,but instead result in a differential escape of the lower crust faster than the upper crust.
基金supported by the National Natural Science Foundation of China (Grant Nos. 41275064 and 41175057)the Public Science and Technology Research Fund Projects of the Ocean (Grant No. 201005019)
文摘Using a new vortex detection and tracing method, a dataset of the Southwest Vortex(SWV) is established based on Japanese 25-year Reanalysis(JRA-25) reanalysis data during 1979–2008. The spatiotemporal features of the SWV are derived from the dataset. In comparison to other seasons, summer yields the least SWVs, but with the highest probability that they will migrate from their region of origin. SWVs mostly emerge in the southwest of the Sichuan Basin and the southeast of the Tibetan Plateau. Migratory SWVs mainly move along either an eastward or southeastward path. Detailed composite analysis of warm-season SWVs shows that the subtropical high is a key factor in determining the direction of migratory SWVs. Furthermore, the steering wind at 700 hPa dominates the moving direction of migratory SWVs. Potential stability diagnosed by pseudo-equivalent potential temperature ? se is of certain significance for the evolution and movement of SWVs. On the other hand, migratory SWVs possess relatively greater strength than stationary SWVs, due to a stronger low-level jet with enhanced baroclinicity and moisture transport providing more energy to support the growth of SWVs along their paths of movement.
基金supported by the National Natural Science Foundation of China(Grant Nos.41030422, 40972131,40772135 and 41202120)the National Basic Research Program of China(Grant Nos.2009CB219601 and 2006CB202201)the China Postdoctoral Science Foundation Funded Project(2012M510590)
基金supported financially by the National Natural Science Foundation of China (No.41574021,41474059,41331066,41774088,and 41174063)a research grant from the Institute of Crustal Dynamics,China Earthquake Administration (No.ZDJ2017-23)+4 种基金the CAS/CAFEA International Partnership Program for Creative Research Teams (No.KZZD-EW-TZ-19)the Key Research Program of Frontier Sciences CAS (Chinese Academy of Sciences) (QYZDY-SSW-SYS003)the SKLGED foundation (SKLGED2014-1-1-E)the GOCE Italy Project (the Italian Space Agency and the European Space Agency Endorsement)the China Postdoctoral Science Foundation (No.133014)
文摘Data obtained by GRACE(Gravity Recovery and Climate Experiment) have been used to invert for the seismic source parameters of megathrust earthquakes under the assumption of either uniform slip over an entire fault or a point-like seismic source.Herein, we further extend the inversion of GRACE long-wavelength gravity changes to heterogeneous slip distributions during the 2011 Tohoku earthquake using three fault models:(Ⅰ) a constant-strike and constant-dip fault,(Ⅱ) a variable dip fault, and(Ⅲ) a realistically varying strike fault. By removing the post-seismic signal from the time series, and taking the effect of ocean water redistribution into account, we invert for slip models I, II, and III using co-seismic gravity changes measured by GRACE, de-striped by DDK3 decorrelation filter. The total seismic moments of our slip models, with respective values of 4.9×10^(22) Nm, 5.1×10^(22) Nm, and 5.0×10^(22) Nm, are smaller than those obtained by other studies relying on GRACE data. The resulting centroids are also located at greater depths(20 km, 19.8 km,and 17.4 km, respectively). By combining onshore GPS, GPS-Acoustic, and GRACE data, we obtain a jointly inverted slip model with a seismic moment of 4.8×10^(22) Nm, which is larger than the seismic moment obtained using only the GPS displacements. We show that the slip inverted from low degree space-borne gravimetric data, which contains information at the ocean region, is affected by the strike of the arcuate trench. The space-borne gravimetric data help us constrain the source parameters of a megathrust earthquake within the frame of heterogeneous slip models.
基金supported by National Science and Technology Major Project(No2009 ZX 05039-003,2009 ZX 05039-004,2011ZX05060-005)the National Natural Science Foundation of China(No 2010CB428801-1)state-owned land resources investigation(1212010430351)
基金financially supported by the National Natural Science Foundation of China (Grant No.41174063)
文摘Modern geodetic techniques have developed rapidly in recent years, providing reliable observation data and new effective approaches, and greatly enhancing studies of the Tibetan geodynamics. For instance, the well-known GPS technique has been employed to measure seismic slips for many faults in the Tibetan Plateau. GPS data agree well with the hypothesis of a thickening crust and eastward mass flow. Moreover, absolute gravimetric data have been applied to interpret geophysical phenomena such as crust movement, co-seismic gravity change, GIA, and ground water change. The satellite gravity mission GRACE launched in 2002 provided global gravity models with unprecedentedly high precision and high spatial resolution. It has been used in implementing temporal gravity changes and improving our knowledge of the Earth's interior, including lithosphere dynamics, mantle viscosity and rheology, plateau uplift, and subduction processing. It is noteworthy that gravity presents unique advantages for the study of Tibetan geodynamics because of its sensitivity to mass migration and dynamic redistribution. To date, great advances have been made in applying modern geodetic data in studying dynamic changes of Tibetan plateau. For instance, the horizontal displacement field from GPS data revealed dynamical characteristics of the present-day Tibetan plateau. The combination of gravity anomalies and topographic data describe the tectonic characteristics of Tibetan plateau. The combination of gravity data and GPS data show present properties of the Tibetan plateau such as crust thickening, Moho's subsidence, and plateau uplift. GRACE data were used to estimate the distribution of ice/snow melting. These results demonstrate that mere application of integrated geodetic data as well as geophysical methods and numerical simulations can enhance our knowledge of Tibetan plateau dynamics. It must be pointed out that GRACE data include various geophysical signals such as crust vertical movement, denudation, ice and snow melting, GIA, ground water change, and permafrost degradation. To separate the tectonic information from other impulses, each physical signal must be evaluated and corrected carefully from the GRACE data. The Tibetan geodynamic problem is a complicated and synthetic issue that must be addressed through collaboration of workers in many fields. Succinctly put, although great achievements have been made in studying Tibetan plateau dynamics from each field, the dynamical process remains unclear. Some fundamental problems remain unresolved. They should be solved with modern geodetic data, such as GRACE, GPS, and absolute gravity data, combined with meteorological and geological data, for quantitative analysis of Tibetan plateau dynamics affected by respective geophysical sources. This review article introduces and discusses the scientific importance, advances, problems, and prospects of modern geodesy applied to the study of geodynamic changes of the Tibetan plateau.
基金supported by National Natural Science Foundation of China (Grant No.40874047)supported by National Natural Science Foundation of China (Grant No.40975009)supported by the National Key Basic Research Program (Grant No. 2012CB417203)
文摘Using the Simple Biosphere Model(SiB2),soil thermal properties(STP) were examined in a Tibetan prairie during the monsoon period to investigate ground surface temperature prediction.We improved the SiB2 model by incorporating a revised force-restore method(FRM) to take the vertical heterogeneity of soil thermal diffusivity(k) into account.The results indicate that(1) the revised FRM alleviates daytime overestimation and nighttime underestimation in modeled ground surface temperature(Tg),and(2) its role in little rainfall events is significant because the vertical gradient of k increases with increasing surface evaporation.Since the original formula of thermal conductivity(λ) in the SiB2 greatly underestimates soil thermal conductivity,we compared five algorithms of λ involving soil moisture to investigate the cause of overestimation during the day and underestimation at night on the basis of the revised FRM.The results show that(1) the five algorithms significantly improve Tg prediction,especially in daytime,and(2) taking one of these five algorithms as an example,the simulated Tg values in the daytime are closer to the field measurements than those in the nighttime.The differences between modeled Tg and field measurements are mostly within the margin of error of ±2 K during 3 August to 4 September 1998.
基金supported by National Natural Science Foundation of China (No. 41174085)Chinese Academy of Sciences (KZZD-EW-TZ-19)China Geological Survey (12120113101400)
文摘Site engineering seismic survey provides basic data for seismic effect analysis. As an important parameter of soil, shear-wave velocity is usually obtained through wave velocity testing in borehole. In this paper, the passive source surface-wave method is introduced into the site engineering seismic survey and practically applied in an engineering site of Shijingshan District. By recording the ubiquitous weak vibration on the earth surface, extract the dispersion curve from the surface-wave components using the SPAC method and obtain the shear-wave velocity structure from inversion. Over the depth of 42 m underground, it totally consists of five layers with interface depth of 3.31, 4.50, 7.23, 17.41, and 42.00 m; and shear-wave velocity of 144.0, 198.3, 339.4, 744.2, and 903.7 m/s, respectively. The inversion result is used to evaluate site classification, determine the maximum shear modulus of soil, provide basis for further seismic hazard analysis and site assessment or site zoning, etc. The result shows that the passive source surface-wave method is feasible in the site engineering seismic survey and can replace boreholes,shorten survey period, and reduce engineering cost to some extent.
基金financially supported by the National Natural Science Foundation of China(Grant No.41530315,41372213)the National Science and Technology Major Project of China(Grant No.2016ZX05066003,2016ZX05066006)the“Climate Change:Carbon Budget and Related Issues”Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDA05030100)
文摘The fracturing technology for shale gas reservoir is the key to the development of shale gas industrialization.It makes much sense to study the mechanical properties and deformation characteristics of shale,due to its close relationship with the fracability of shale gas reservoir.This paper took marine shale in the Changning area,southern Sichuan Basin of China as the research object.Based on field profile and hand specimen observation,we analyzed the development of natural fractures and collected samples from Wufeng Formation and Longmaxi Formation.Combining with the indoor experiment,we investigated the macroscopic and microscopic structural features and the remarkable heterogeneity of shale samples.Then we illustrated the mechanics and deformation characteristics of shale,through uniaxial compression test and direct shear test.The shale has two types of fracture modes,which depend on the angular relation between loading direction and the bedding plane.Besides,the Wufeng shale has a higher value of brittleness index than the Longmaxi shale,which was calculated using two methods,mechanical parameters and mineral composition.Given the above results,we proposed a fracability evaluation model for shale gas reservoir using the analytic hierarchy process.Four influence factors,brittleness index,fracture toughness,natural fractures and cohesive force,are considered.Finally,under the control of normalized value and weight coefficient of each influence factor,the calculations results indicate that the fracability index of the Wufeng Formation is higher than that of the Longmaxi Formation in Changning area,southern Sichuan Basin.
基金funded jointly by State's Key Project of Research and Development Plan(2016YFB0501702)National natural science fund of China(41274034)+1 种基金CAS/CAFEA international partnership for creative research teams(KZZD-EW-TZ-19)Beijing key laboratory of urban spatial information engineering (2016205)
文摘The boundary value problem of deflections of vertical with ellipsoid boundary is studied in the paper.Based on spherical harmonic series,the ellipsoidal corrections for the boundary value problem are derived so that it can be well solved.In addition,an imitation arithmetic is given for examining the accuracies of solutions for the boundary value problem as well as its spherical approximation problem,and the computational results illustrate that the boundary value problem has higher accuracy than its spherical approximation problem if deflection of the vertical are measured on geoid.
基金supported by China Postdoctoral Science Foundation 2012M510043,2013T60166,NSFC 41074032 and CAS KZCX2-EW-121
文摘The Yellowstone volcano is one of the largest active volcanoes in the world,and its potential hazards demand detailed seismological and geodetic studies.Previous studies with travel time tomography and receiver functions have revealed a low-velocity layer in the crust beneath the Yellowstone volcano,suggesting the presence of a magma chamber at depth.We use ambient seismic noise from regional seismic stations to retrieve short-period surface waves and then study the shallow shear velocity structure of the Yellowstone region by surface wave dispersion analysis.We first obtained a crustal model of the area outside of the Yellowstone volcano and then constructed an absolute shear wave velocity structure in combination with receiver function results for the crust beneath the Yellowstone volcano.The velocity model shows a low-velocity layer with shear velocity at around1.3 km/s,suggesting that a large-scale magma chamber exists at shallow levels within the crust of the Yellowstone volcanic region.
基金the National Natural Science Foundation of China(No.41604067,41974093,41331066,and 41774088)the Basic Research Fund of Chinese Academy of Surveying and Mapping(No.AR 1906)+1 种基金the special project of high-resolution Earth observation system(42-Y20A09-9001-17/18)the Key Research Program of Frontier Sciences Chinese Academy of Sciences(QYZDY-SSWSYS003).
文摘In this paper,we study how coseismic deformations calculated in 1066 Earth models are affected by how the models treat Earth discontinuities.From the results of applying models 1066A(continuous)and 1066B(discontinuous),we find that the difference in Love numbers of strike-slip and horizontal tensile sources are bigger than dip-slip and vertical tensile sources.Taken collectively,discontinuities have major effects on Green’s functions of four independent sources.For the near-field coseismic deformations of the 2013 Okhotsk earthquake(Mw 8.3),the overall differences between theoretical calculations in vertical displacement,geoid,and gravity changes caused by discontinuities are 10.52 percent,9.07 percent and 6.19 percent,with RMS errors of 0.624 mm,0.029 mm,and 0.063μGal,respectively.The difference in far-field displacements is small,compared with GPS data,and we can neglect this effect.For the shallow earthquake,2011 Tohoku-Oki earthquake(Mw 9.0),the differences in near-field displacements are 0.030 m(N-S),0.093 m(E-W),and 0.025 m(up-down)in our study area with the ARIA slip model,which gives results closer to GPS data than those from the USGS model.The difference in vertical displacements and gravity changes on the Earth’s surface caused by discontinuities are larger than 10 percent.The difference in the theoretical gravity changes at spatially fixed points truncated to degrees 60,as required by GRACE data,is 0.0016μGal and the discrepancy is 11 percent,with the theoretical spatial gravity changes from 1066B closer to observations than from 1066A.The results show that an Earth model with discontinuities in the medium has a large effect on the calculated coseismic deformations.
基金supported financially by the National Natural Science Foundation of China(41174063,41331066 and41474059)the CAS/CAFEA International Partnership Program for Creative Research Teams(KZZD-EW-TZ-19)the SKLGED Foundation(2014-1-1-E)
文摘The complex geographical environment in China makes its gravity signals miscellaneous.This work gives a comprehensive representation and explanation in secular trend of gravity change in different regions,the key features of which include positive trend in inner Tibet Plateau and South China and negative trend in North China plain and high mountain Asia(HMA).We also present the patterns of amplitudes and phases of annual and semiannual change.The mechanism underlying the semiannual period is explicitly discussed.The displacement in three directions expressed in terms of geo-potential spherical coefficients and load Love numbers are given.A case study applied with these equations is presented.The results show that Global Positioning System(GPS) observations can be used to compare with Gravity Recovery and Climate Experiment(GRACE) derived displacement and the vertical direction has a signal-noise-ratio of about one order of magnitude larger than the horizontal directions.
基金supported by Sinoprob-Deep Exploration Program in China sponsored by Ministry of Land and Resources of the People’s Republic of China (No.0819011A90)National Natural Science Foundation of China (No.D0408/4097409)
文摘The distribution of Coulomb failure stress (CFS) change in the steep excavation slope is calculated by finite element method in this paper,and the failure mechanics under different conditions have been investigated.Comparing the CFSs before and after the slope excavation (stress loading and unloading processes),the dangerous internal zone and the most likely failure external area are attained.Given the shear cracks on the top surface while tensile stress or cracks along the toe of the slope,we analyze the high cutting-angle steep slope in Kaixian county of the Three Gorges Reservoir region.We bring forward that the peak value of CFS after excavation can reach to the order of 0.1MPa,which is greatly higher than that of before.Our preliminary results are useful for optimizing the reinforcement structure during the steep slope stabilization engineering.
基金supported by the National Natural Science Foundation of China (41604140)the Specialized Research Fund for State Key Laboratories and China Postdoctoral Science Foundation Fund
文摘Nonlinear force-free magnetic field(NLFFF) extrapolation based on the observed photospheric magnetic field is the most important method to obtain the coronal magnetic field nowadays.However, raw photospheric magnetograms contain magnetic forces and small-scale noises, and fail to be consistent with the force-free assumption of NLFFF models. The procedure for removing the forces and noises in observed data is called preprocessing. In this paper, we extend the preprocessing code of Jiang& Feng to spherical coordinates for a full sphere. We first smooth the observed data with Gaussian smoothing, and then split the smoothed magnetic field into a potential field and a non-potential field.The potential part is computed by a numerical potential field model, and the non-potential part is preprocessed using an optimization method to minimize the magnetic forces and magnetic torques. Applying the code to synoptic charts of the vector magnetic field from SDO/HMI, we find it can effectively reduce the noises and forces, and improve the quality of data for a better input which will be used for NLFFF extrapolations applied to the global corona.
基金This research was supported financially by the NNSFC(41774088,41331066,42104084,and 41474059)the Key Research Program of Frontier Sciences CAS(Chinese Academy of Sciences)(QYZDY-SSW-SYS003)。
文摘Global mean sea level(GMSL) change is one of the important indicators of global climate change and is a crucial scientific issue of continuing interest. As satellite altimeter data, the Gravity Recovery and Climate Experiment(GRACE) and Argo continue to be updated, especially with the release of GRACE Follow-On(GRACE-FO) data, making it necessary to combine these latest data to estimate sea level change. Determinations on whether the GRACE and GRACE-FO observation systems provide unbiased global observation data have not been effectively evaluated. Therefore, this research mainly investigated the consistency of GRACE and GRACE-FO observation data in studying GMSL change. By comparing the sum of the GMSL calculated by the two gravity satellites and Argo data with the GMSL calculated by satellite altimeters, the discrepancy between GRACE-FO + Argo and satellite altimeter data is about 7.9 ± 2.3 mm, which is likely derived from the inconsistency between GRACE and GRACE-FO data.
基金financially supported by the National Key Research and Development Program of China(2016YFD0300110,2016YFD0300101)the National Natural Science Foundation of China(31671585,41871253)the Natural Science Foundation of Shandong Province,China(ZR2017ZB0422)。
文摘Under the limited cultivated land area and the pursuit of sustainable agricultural development,it is essential for the safety of grain production to study agricultural management approaches on narrowing the winter wheat yield gap and improving nitrogen use efficiency (NUE) in China.In this study,DSSAT-CERES-Wheat Model is used to simulate winter wheat yield under different agricultural treatments,and we analyze yield gaps and NUE with different management scenarios at regional scales and evaluate the suitable approaches for reducing yield gap and increasing NUE.The results show that,the potential of narrowing yield gap ranges 300–900 kg ha^(–1) with soil nutrients increase,400–1 200 kg ha^(–1) with sowing date adjustment and 0–400 kg ha^(–1) with planting density increase as well as 700–2 200 kg ha^(–1) with adding nitrogen fertilizer.Contribution rates of management measures of soil nutrients,sowing date adjusting,planting density,and nitrogen fertilizers are 5–15%,5–15%,0–4%,and 10–20%,respectively.Difference in nitrogen partial productivity ranges 3–10 kg kg^(–1) for soil nutrients,1–10 kg kg^(–1) for sowing date adjusting,1–5 kg kg^(–1) for planting density increase,and–12–0 kg kg^(–1) for adding nitrogen fertilizers,respectively.It indicates that four treatments can narrow yield gap and improve the NUE in varying degrees,but increasing nitrogen fertilizer leads to the decrease of NUE.
基金supported by the National Natural Science Foundation of China(grant Nos.40874047 and 40574047)
文摘Global seismicity catalogs are sufficient for characterizing double seismic zones(DSZs)in subducting slab and facilitate to estimate layer separation without inconsistent uncertainties as local catalogs.Previous studies have shown the correlation between DSZs layer separation and plate age while correlation for those younger than ～60 Ma is suspicious.The lacking of DSZs with layer separation less than 10 km further makes it difficult to precisely estimate such correlation.Thus,we incorporate eight DSZs data determined through local seismicity into globally-determined dataset and reexamine such correlation.The best fitting results show that both a linear model and a square root of plate age can mathematically fit the layer separation well.However,it is difficult to distinguish these two models when plate age is greater than～20 Ma since their difference is less than 2 km.However,if extrapolation is possible,both models should provide physical information that DSZs will not form if there is no subducting lithosphere.As a result,the DSZs cannot be produced until the oceanic lithospheric age becomes greater than 0.9 Ma in the square root model while the linear model gives a misleading result.As such the square root model demonstrates the relationship physically better than the linear one,it still needs further test in the future with more available data,nevertheless,our study might also provide evidence for the suggestion that the plate age is a primary control factor of the DSZs geometry as well as the subducting process which disregards any local tectonic stresses.
基金The National Natural Science Foundation of China under contract No.61001137the Project of Knowledge Innovative Program of the Chinese Academy of Sciences and other projects under contract Nos Y1530151A81530151G4 and Y15102EN00
文摘The relationships among an ocean wave spectrum,a fully polarimetric coherence matrix,and radar parameters are deduced with an electromagnetic wave theory.Furthermore,the relationship between the polarimetric entropy and ocean wave spectrum is established based on the definition of entropy and a twoscale scattering model of the ocean surface.It is the first time that the polarimetric entropy of the ocean surface is presented in theory.Meanwhile,the relationships among the fully polarimetric entropy and the parameters related to radar and ocean are discussed.The study is the basis of further monitoring targets on the ocean surface and deriving oceanic information with the entropy from the ocean surface.The contrast enhancement between human-made targets and the ocean surface with the entropy is presented with quad-pol airborne synthetic aperture radar(AIRSAR)data.
