An estimation of groundwater storage variations from GRACE gravity satellites in the Heihe River Basin, northwestern China

There are only limited surface water resources available in the Heihe River Basin （HRB）, a typical inland river basin in the arid region of northwestern China, where groundwater overexploitation is a serious problem. Groundwater has become one of main resources of fresh water in the HRB. In this paper, temporal and spatial variations of groundwater in the HRB are estimated by the Gravity Recovery and Climate Experiment （GRACE） satellites. Our analysis shows that groundwater storage in the HRB reaches its highest in the summer of 2005, and then begins to decline in the following years and reaches steady status in 2008. Spatially, groundwater shows a decline in the upper HRB in the first two years and a slight increase in the following years, while this phenomenon is reversed in the middle HRB where groundwater slightly increases in 2005 and then declines in the following three years. In the lower HRB, GRACE detects a continual increase in the full six-year period. This approach is proven successful when employed in the HRB and thus offers a new insight into monitoring groundwater variations in a river basin with limited or even without any observed data.     

A probabilistic approach to the inversion of data from a seismic array and its application to volcanic signals

Array techniques are particularly well‐suited for detecting and quantifying the complex seismic wavefields associated with volcanic activity such as volcanic tremor and long‐period events. The methods based on the analysis of the signal in the frequency domain, or spectral methods, have the main advantages of both resolving closely spaced sources and reducing the necessary computer time, but may severely fail in the analysis of monochromatic, non‐stationary signals. Conversely, the time‐domain methods, based on the maximization of a multichannel coherence estimate, can be applied even for short‐duration pulses. However, for both the time and the frequency domain approaches, an exhaustive definition of the errors associated with the slowness vector estimate is not yet available. Such a definition become crucial once the slowness vector estimates are used to infer source location and extent. In this work we develop a method based on a probabilistic formalism, which allows for a complete definition of the uncertainties associated with the estimate of frequency–slowness power spectra from measurement of the zero‐lag cross‐correlation. The method is based on the estimate of the theoretical frequency–slowness power spectrum, which is expressed as the convolution of the true signal slowness with the array response pattern. Using a Bayesian formalism, the a posteriori probability density function for signal slowness is expressed as the difference, in the least‐squares sense, between the model spectrum and that derived from application of the zero‐lag cross‐correlation technique. The method is tested using synthetic waveforms resembling the quasi‐monochromatic signals often associated with the volcanic activity. Examples of application to data from Stromboli volcano, Italy, allow for the estimate of source location and extent of the explosive activity.     

Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity

The GRACE satellite mission has been measuring the Earth's gravity field and its temporal variations since 2002 April. Although these variations are mainly due to mass transfer within the geofluid envelops, they also result from mass displacements associated with phenomena including glacial isostatic adjustment and earthquakes. However, these last contributions are difficult to isolate because of the presence of noise and of geofluid signals, and because of GRACE's coarse spatial resolution (>400 km half-wavelength). In this paper, we show that a wavelet analysis on the sphere helps to retrieve earthquake signatures from GRACE geoid products. Using a wavelet analysis of GRACE geoids products, we show that the geoid variations caused by the 2004 December ( M _w= 9.2) and 2005 March ( M _w= 8.7) Sumatra earthquakes can be detected. At GRACE resolution, the 2004 December earthquake produced a strong coseismic decrease of the gravity field in the Andaman Sea, followed by relaxation in the area affected by both the Andaman 2004 and the Nias 2005 earthquakes. We find two characteristic timescales for the relaxation, with a fast variation occurring in the vicinity of the Central Andaman ridge. We discuss our coseismic observations in terms of density changes of crustal and upper-mantle rocks, and of the vertical displacements in the Andaman Sea. We interpret the post-seismic signal in terms of the viscoelastic response of the Earth's mantle. The transient component of the relaxation may indicate the presence of hot, viscous material beneath the active Central Andaman Basin.     

Iceocean mass balance during the Late Pleistocene glacial cycles in view of CHAMP and GRACE satellite missions

During the last glacial cycles, global sea level dropped several times by about 120 m and large ice sheets covered North America, northern Europe and Antarctica during the glacial stages. The changes in the iceocean mass balance have displaced mantle material mainly via viscous flow, and the perturbation of the equilibrium figure of the Earth by glacial isostatic adjustment is still observable today in timedependent changes of gravitational and rotational observations. Contemporary iceocean mass balance from volume changes of polar ice caps also contributes to secular variations of the Earth's gravitational field.
In the near future, several satellite gravity missions will significantly improve the accuracy of the observed timedependent gravitational field. In view of the expected improvements in the observations, we predict glacially induced perturbations of the gravitational field, induced by Late Pleistocene and contemporary ice volume changes, for a variety of radial mantle viscosity profiles. We assess the degree of uncertainty for the glacially induced contributions to gravitational and rotational parameters, both in the spectral and the spatial domain.
Predictions of power spectra for the glacially induced freeair gravity and geoid anomalies are about one order of magnitude lower than the observed values, and uncertainties arising from different plausible viscosity profiles are around 0.150.4 mGal and 0.21.5 m, respectively. Uncertainties from different ice models are of secondary importance for the predicted power spectra. Predicted secular changes in geoid anomalies in formerly glaciated areas are mainly controlled by the viscosity profile and contemporary ice volume changes. We also show that the simple threelayer viscosity profiles currently employed for the majority of postglacial rebound studies represent a limited subset for model predictions of the timedependent gravitational field.     

2002-2010年长江流域GRACE水储量时空变化特征

利用高斯平滑滤波对2002年4月-2010年12月逐月GRACE卫星的时变重力场数据反演得到长江流域大尺度陆地水储量变化,对其时空变化进行研究,并将结果与全球陆面同化数据(GLDAS)模拟结果进行比较。其结论为:根据GRACE数据反演与GGLDAS模拟得到的水储量结果在大多数区域变化趋势相同,两者具有一致性,相关性达到0.89(P<0.05)。GRACE水储量研究结果表明:①2002-2010年长江流域水储量呈增加趋势,平均增长速率为0.43mm/月,相当于约95.04亿m³/年。长江上游增长速率为0.53mm/月,相当于约67.13亿m³/年;中游增长速率为0.51mm/月,相当于25.73亿m³/年;下游增长速率为0.36mm/月,相当于9.14亿m³/年。近9年长江流域水储量共增加约855.33亿m³。②从多年平均水储量空间分布来看,长江流域冬季月份(12、1、2、3月)水储量处于亏损状态,7-9月水储量处于盈余状态,4-6月下游至上游地区由亏损向盈余状态过渡,而10-11月则从上游至下游地区由盈余向亏损状态过渡。③全流域、上游及中游水储量逐月增长速率最大值出现在9月,分别为1.01cm/a、1.37cm/a、1.05cm/a;而下游地区则出现在7月,增长速率为1.62cm/a。     

联合卫星重力和卫星测高确定南极绕极流

联合基于GRACE重力卫星观测资料恢复的重力场模型(EIGEN-GL04S1)和卫星测高推求的平均海面高模型(KMSS04)来构造南极绕极流区域的平均海面动力地形,并利用小波滤波方法去掉短波及噪声信号,进而推算大、中尺度的绕极流。与非卫星重力场模型、同化资料及海洋水文资料确定相应结果的验证分析表明:基于新的卫星重力场模型推算的南极绕极流区域的海面动力地形、PF、SAF和表层流场等都与海洋学结果相吻合,且局部特征更加清晰。表明卫-卫跟踪重力卫星计划确定的地球重力场模型较之以前存在的重力场模型在中长波部分精度有较大提高,从大地测量(从空间)角度来探测南极绕极流已达到较高的精度。     

全球降雨计划GSMaP与IMERG卫星降雨产品在陕西地区的精度评估

以中国气象局提供的地面降雨观测资料为参考,利用相关系数、均方根误差和相对偏差,降雨误报率、命中率和关键成功指数等6种不同的统计分析指标,从年、季、月、日4种不同的时间尺度对两种高分辨率卫星降雨产品(IMERG和GSMaP)在陕西省的精度进行了对比和评价,并对二者在监测强降雨过程中的表现进行了对比分析.结果表明:(1)在...     

Improved magnetic anomalies of the Antarctic lithosphere from satellite and near-surface data

The Antarctic magnetic anomaly map compiled marine and airborne surveys collected south of 60°S through 1999 and used Magsat data to help fill in the regional gaps between the surveys. Ørsted and CHAMP satellite magnetic observations with greatly improved measurement accuracies and temporal and spatial coverage of the Antarctic, have now supplanted the Magsat data. We combined the new satellite observations with the near-surface survey data for an improved magnetic anomaly map of the Antarctic lithosphere. Specifically, we separated the crustal from the core and external field components in the satellite data using crustal thickness variations estimated from the terrain and the satellite-derived free-air gravity observations. Regional gaps in the near-surface surveys were then filled with predictions from crustal magnetization models that jointly satisfied the near-surface and satellite crustal anomalies. Comparisons in some of the regional gaps that also considered newly acquired aeromagnetic data demonstrated the enhanced anomaly estimation capabilities of the predictions over those from conventional minimum curvature and spherical harmonic geomagnetic field models. We also noted that the growing number of regional and world magnetic survey compilations involve coverage gaps where these procedures can contribute effective near-surface crustal anomaly estimates.     

Vegetation-related dry deposition of global PM2.5 from satellite observations

Journal of Geographical Sciences - Vegetation plays an important role in the dry deposition of particles with significant spatial variability, but the magnitude remains unclear at the global scale....     

On the construction of geomagnetic timescales from non-prejudicial treatment of magnetic anomaly data from multiple ridges

When marine magnetic-anomaly data are used to construct geomagnetic polarity timescales, the usual assumption of a smooth spreading-rate function at one seafloor spreading ridge forces much more erratic rate functions at other ridges. To eliminate this problem, we propose a formalism for the timescale problem that penalizes non-smooth spreading behaviour equally for all ridges. Specifically, we establish a non-linear Lagrange multiplier optimization problem for finding the timescale that (1) agrees with known chron ages and with anomaly-interval distance data from multiple ridges and (2) allows the rate functions for each ridge to be as nearly constant as possible, according to a cumulative penalty function. The method is applied to a synthetic data set reconstructed from the timescale and rate functions for seven ridges, derived by Cande & Kent (1992) under the assumption of smooth spreading in the South Atlantic. We find that only modest changes in the timescale (less than 5 per cent for each reversal) are needed if no one ridge is singled out for the preferential assumption of smoothness. Future implementation of this non-prejudicial treatment of spreading-rate data from multiple ridges to large anomaly-distance data sets should lead to the next incremental improvement to the pre-Quaternary geomagnetic polarity timescale, as well as allow a more accurate assessment of global and local changes in seafloor spreading rates over time.     

1992–2015全球耕地时空变化（英文）

人口和人均食物需求的增加对全球耕地产生了显著的影响。利用欧空局提供的精度为300m的最新土地覆被产品,文章分析了1992-2015年全球耕地的时空变化趋势和耕地转化特征。结果显示:1)在1992-2004年间全球耕地面积增长迅速,而在2004-2012年间耕地增长缓慢,2012年后耕地有缓慢减少的趋势。2)在洲尺度上,非洲耕地有一直增长的态势,而其他洲耕地都经历了耕地转型,有先增长后下降的趋势;在收入较高的国家,耕地多有下降的趋势。3)全球耕地增长的热点区域主要分布在亚马逊林地、欧亚大草原和撒哈拉沙漠边缘。全球耕地减少的中心从欧洲转移到亚洲。由于迅速的城市化,亚洲耕地扩张侵占了大量农田。     

Spatial and Temporal Variation of Cropland at the Global Level from 1992 to 2015

Correlated increases in population and demand for food over recent decades have caused remarkable changes in cropland area globally. Utilizing the latest data product provided by the European Space Agency (ESA), this paper analyzes annual trends and spatiotemporal variations in cropland area and discusses cropland conversion, losses, and increases globally between 1992 and 2015 at a 300 m resolution. The results show that the global area of cropland increased rapidly between 1992 and 2004, more slowly between 2004 and 2012, and began to decrease gradually since 2012. First, an increasing trend in cropland area has been maintained solely in Africa; all other regions are characterized by decreasing trends in later periods despite different transition points and change rates. A reduction in cropland area frequently emerged earlier in high-income countries. Second, increase rates in cropland area have largely decreased in recent years while the overall rate of loss has remained almost the same. Hotspot areas of global cropland increases are mainly concentrated around the edge of the Amazon forest, Eurasian Steppe, and Sahara Desert. These hotspot areas of global cropland loss shifted from Europe to Asia while built-up areas have expanded at the expense of cropland.     

上海城市功能的转型:从全球生产系统角度的透视

本文首先回顾了价值链与全球生产系统的研究进展，结合跨国公司最新发展趋势对全球生产系统特征进行了总结；其次以改革开放以来上海对外贸易结构的改变，描述了90年代以来上海在全球生产系统分工中的变化，认为出口产品结构的转型并没有真正改变上海在国际劳动分工中处于低端环节的地位；最后就上海建设世界城市的差距与优势进行了总结，提出全球化发展的外部力量和地方化响应的内部力量将是上海发展世界城市的两个重要驱动力。