Wiener optimal filtering of GRACE data

We present a spatial averaging method for Gravity Recovery and Climate Experiment (GRACE) gravity-field solutions based on the Wiener optimal filtering. The optimal filter is designed from the least-square minimization of the difference between the desired and filtered signals. It requires information about the power spectra of the desired gravitational signal and the contaminating noise, which is inferred from the average GRACE degree-power spectrum. We show that the signal decreases with increasing spherical harmonic degree j with approximately j^−b, where b = 1.5 for GRACE data investigations. This is termed the Second Kaula rule of thumb for temporal variations of the Earth’s gravity field. The degree power of the noise increases, in the logarithmic scale, linearly with increasing j. The Wiener optimal filter obtained for the signal model with b = 1.5 closely corresponds to a Gaussian filter with a spatial half width of 4° (∼440 km). We find that the filtered GRACE gravity signal is relatively insensitive to the exponent b of the signal model, which indicates the robustness of Wiener optimal filtering. This is demonstrated using the GFZ-GRACE gravity-field solution for April 2004.     

Thermal and exhumation history of the central Rwenzori Mountains, Western Rift of the East African Rift System, Uganda

The Rwenzori Mountains (Mtns) in west Uganda are the highest rift mountains on Earth and rise to more than 5,000 m. We apply low-temperature thermochronology (apatite fission-track (AFT) and apatite (U–Th–Sm)/He (AHe) analysis) for tracking the cooling history of the Rwenzori Mtns. Samples from the central and northern Rwenzoris reveal AFT ages between 195.0 (±8.4) Ma and 85.3 (±5.3) Ma, and AHe ages between 210.0 (±6.0) Ma to 24.9 (±0.5) Ma. Modelled time–temperature paths reflect a protracted cooling history with accelerated cooling in Permo-Triassic and Jurassic times, followed by a long period of constant and slow cooling, than succeeded by a renewed accelerated cooling in the Neogene. During the last 10 Ma, differentiated erosion and surface uplift affected the Rwenzori Mtns, with more pronounced uplift along the western flank. The final rock uplift of the Rwenzori Mtns that partly led to the formation of the recent topography must have been fast and in the near past (Pliocene to Pleistocene). Erosion could not compensate for the latest rock uplift, resulting in Oligocene to Miocene AHe ages.     

Geological, mineralogical and geochemical characteristics of the Radzimowice Au–As–Cu deposit from the Kaczawa Mountains (Western Sudetes, Poland): an example of the transition of porphyry and epithermal style

The sheeted quartz–sulfide veins of the Radzimowice Au–As–Cu deposit in the Kaczawa Mountains are related to Upper Carboniferous post-collisional potassic magmatism of the composite Zelezniak porphyry intrusion. Multiple intrusive activity ranges from early calc-alkaline to sub-alkaline and alkaline rocks and is followed by multiple hydrothermal events. Early crustally derived dacitic magma has low mg# (<63) and very low concentrations of mantle-compatible trace elements, high large-ion lithophile elements (LILE), moderate light rare-earth elements (LREE), and low high-field-strength elements (HFSE). Later phases of more alkaline rocks have higher mg# (60–70), and LILE, LREE, and HFSE characteristics that indicate mafic magma contributions in a felsic magma chamber. The last episode of the magmatic evolution is represented by lamprophyre dikes which pre-date ore mineralization and are spatially related to quartz–sulfide–carbonate veins. The dikes consist of kersantite and spessartite of calc-alkaline affinity with K₂O/Na₂O ratios of 1.1–1.9, mg# of 77–79, and high abundances of mantle-compatible trace elements such as Cr, Ni, and V. They have high LILE, low LREE, and low HFSE contents suggesting a subduction-related post-collisional arc-setting. The mineralization started with arsenopyrite that was strongly brecciated and overprinted by multiple quartz–carbonate phases associated with base-metal sulfides and Au–Ag–Bi–Te–Pb±S minerals. The sulfur isotope composition of sulfides ranges from –1.1 to 2.8 ³⁴S and suggests a magmatic source. At least two generations of gold deposition are recognized: (1) early refractory, and (2) subsequent non-refractory gold mineralization of epithermal style. Co-rich arsenopyrite with refractory gold and pyrite are the most abundant minerals of the early stage of sulfide precipitation. Early arsenopyrite formed at 535–345°C along the arsenopyrite–pyrrhotite–loellingite buffer and late arsenopyrite crystallized below 370°C along the arsenopyrite–pyrite buffer. Non-refractory gold associated with base-metal sulfides and with Bi–Te–Ag–Pb–S mineral assemblages has an average fineness of about 685, and is represented by electrum of two generations, and minor maldonite (Au₂Bi). Fluid inclusions from various quartz generations co-genetic with base-metal sulfides and associated with carbonates, tellurides and non-refractory gold indicate fluids with moderate salinity (9–15 wt% NaCl equiv.) and a temperature and pressure drop from 350 to 190°C and 1.2 to 0.8 kbar, respectively. According to the result of the sulfur isotope fractionation geothermometer the temperature of base-metal crystallization was in the range from 322 to 289°C. Preliminary results of oxygen isotope studies of quartz from veins indicate a gradual increase in the proportion of meteoric water in the epithermal stage. The gold to silver ratio in ore samples with >3 ppm Au is about 1:5 (geometric mean). Hydrothermal alteration started with sericitization, pyritization, and kaolinitization in vein selvages followed by alkaline hydrothermal alteration of propylitic character (illitization and chloritization), albitization and carbonatization. The mineralization of the Radzimowice deposit is considered as related to alkaline magmatism and is characterized by the superposition of low-sulfidation epithermal mineralization on higher-temperature and deeper-seated mesothermal/porphyry style.Editorial handling: B. Lehmann     

To what extent can atmospheric deposition influence the natural background of metals in ground waters? A case study in the Czech Republic

In harmony with requirements of Water Framework Directive the natural backgrounds of selected metals were established on the territory of Czech Republic. This paper is intended to provide characteristics and results of testing methods focused on quantification of the impacts of anthropogenic activities based on data of regional monitoring and their separation from acid atmospheric deposition including an appropriate discussion regarding the natural backgrounds of relevant metals. The results indicate that even in regions extremely affected by acid atmospheric deposition as is the territory of the Czech Republic this acid deposition will not have direct impact on the quality of ground waters on a regional scale. Consequently, the proved changes in quality of ground waters seem to have been caused by secondary processes, specifically by mobilization of metals in a rock environment due to increased acidity.     

Linear mechanism of generation and intencification of internal gravity waves in the ionosphere at their interaction with a nonuniform zonal wind: 1. Model of the medium and initial dynamic equations

The specific features of the generation and intensification of internal gravity wave structures in different atmospheric-ionospheric regions, caused by zonal local nonuniform winds (shear flows), are studied. The model of the medium has been explained and an initial closed system of equations has been obtained in order to study the linear and nonlinear dynamics of internal gravity waves (IGWs) when they interact with the geomagnetic field in a dissipative ionosphere (for the D, E, and F regions).     

Estimation of parameters in a split functional model of geodetic observations (<Emphasis Type="Italic">M</Emphasis><Subscript>split</Subscript> estimation)

The paper presents a method of estimating parameters in two competitive functional models. The models considered here are concerned with the same observation set and are based on the assumption that an observation may result from a realization of either of two different random variables. These variables differ from one another at least in the main characteristic (for example, outliers can be realizations of one variable). A quantity that describes the opportunity of identifying a single observation with one random variable is assumed to be known. That quantity, called the elementary split potential, is strictly referred to the amount of information that an observation can provide about two competitive assumptions concerning the observation distribution. Parameter assessments that maximize the global elementary split potential (concerning all observations), are called M _split estimators. A generalization of M _split estimation presented in the paper refers to the theoretical foundation of M-estimation. An erratum to this article can be found at     

基于SAR干涉点目标分析技术的城市地表形变监测

通过深入研究干涉点目标的相位模型,提出基于空间搜索的邻近点目标干涉相位差解缠方法,用以计算点目标的地形残差和线性形变,以及分离点目标大气延迟相位和非线性形变相位的时空域滤波方法,解决干涉点目标分析中的关键问题.最后,以苏州地区地表沉降监测为应用试验,利用形成的SAR干涉点目标形变信息提取技术,获取苏州市区1992-2002年间的地表沉降信息.研究结果与已有文献记录保持比较好的一致性,证明SAR干涉点目标技术完全可以发展成为应用于城市地表形变监测的实用化技术.

Abstract：

Interferometric point target analysis (IPTA) is one of the latest developments in radar interferometric processing. It is achieved by analysis of the interferometric phases of some individual point targets, which are discrete and present temporarily stable backscattering characteristics, in long temporal series of interferometric SAR images. This paper analyzes the interferometric phase model of point targets, and then addresses two key issues within IPTA process. Firstly, a spatial searching method is proposed to unwrap the interferometric phase difference between two neighboring point targets. The height residual error and linear deformation rate of each point target can then be calculated, when a global reference point with known height correction and deformation history is chosen. Secondly, a spatial-temporal filtering scheme is proposed to further separate the atmosphere phase and nonlinear deformation phase from the residual interferometric phase. Finally, an experiment of the developed IPTA methodology is conducted over Suzhou urban area. Totally 38 ERS-1/2 SAR scenes are analyzed, and the deformation information over 3 546 point targets in the time span of 1992-2002 are generated. The IPTA-derivecl deformation shows very good agreement with the published result, which demonstrates that the IPTA technique can be developed into an operational tool to map the ground subsidence over urban area.     

Variation in tree mortality and regeneration affect forest carbon recovery following fuel treatments and wildfire in the Lake Tahoe Basin, California, USA

Background

Forest fuel treatments have been proposed as tools to stabilize carbon stocks in fire-prone forests in the Western U.S.A. Although fuel treatments such as thinning and burning are known to immediately reduce forest carbon stocks, there are suggestions that these losses may be paid back over the long-term if treatments sufficiently reduce future wildfire severity, or prevent deforestation. Although fire severity and post-fire tree regeneration have been indicated as important influences on long-term carbon dynamics, it remains unclear how natural variability in these processes might affect the ability of fuel treatments to protect forest carbon resources. We surveyed a wildfire where fuel treatments were put in place before fire and estimated the short-term impact of treatment and wildfire on aboveground carbon stocks at our study site. We then used a common vegetation growth simulator in conjunction with sensitivity analysis techniques to assess how predicted timescales of carbon recovery after fire are sensitive to variation in rates of fire-related tree mortality, and post-fire tree regeneration.

Results

We found that fuel reduction treatments were successful at ameliorating fire severity at our study site by removing an estimated 36% of aboveground biomass. Treated and untreated stands stored similar amounts of carbon three years after wildfire, but differences in fire severity were such that untreated stands maintained only 7% of aboveground carbon as live trees, versus 51% in treated stands. Over the long-term, our simulations suggest that treated stands in our study area will recover baseline carbon storage 10?C35?years more quickly than untreated stands. Our sensitivity analysis found that rates of fire-related tree mortality strongly influence estimates of post-fire carbon recovery. Rates of regeneration were less influential on recovery timing, except when fire severity was high.

Conclusions

Our ability to predict the response of forest carbon resources to anthropogenic and natural disturbances requires models that incorporate uncertainty in processes important to long-term forest carbon dynamics. To the extent that fuel treatments are able to ameliorate tree mortality rates or prevent deforestation resulting from wildfire, our results suggest that treatments may be a viable strategy to stabilize existing forest carbon stocks.     

Use of the L2C signal for inversions of GPS radio occultation data in the neutral atmosphere

Results from processing FORMOSAT-3/COSMIC radio occultations (RO) with the new GPS L2C signal acquired both in phase locked loop (PLL) and open loop (OL) modes are presented. Analysis of L2P, L2C, and L1CA signals acquired in PLL mode shows that in the presence of strong ionospheric scintillation not only L2P tracking, but also L1CA tracking often fails, while L2C tracking is most stable. The use of L2C improves current RO processing in the neutral atmosphere mainly by increasing the number of processed occultations (due to significant reduction in the number of L2 tracking failures) and marginally by a reduction in noise in statistics. The latter is due to the combination of reduced L2C noise (compared to L2P) and increased L1CA noise in those occultations where L2P would have failed. This result suggests application of OL tracking for L1CA and L2C signals throughout an entire occultation to optimally acquire RO data. Two methods of concurrent processing of L1CA and L2C RO signals are considered. Based on testing of individual occultations, these methods allow: (1) reduction in uncertainty of bending angles retrieved by wave optics in the lower troposphere and (2) reduction in small-scale residual errors of the ionospheric correction in the stratosphere.     

EXPERIENCE IN COMPILING LARGE-SCALE MAPS OF EROSION-PRONE AREAS

Specific criteria used in the compilation of large-scale maps of erosion-prone areas are defined and evaluated for tracts of land representing a number of different agricultural cover types. Formulas are derived which facilitate the mapping of such areas based on such readily observable morphometric characteristics as slope length and steepness. Translated by Edward Torrey, Alexandria, VA 22308 from: Vestnik Moskovskogo Universiteta, geografiya, 1988, No. 5, pp. 25-31.