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.     

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     

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).     

Site effect on vulnerability of high-rise shear wall buildings under near and far field earthquakes

Worldwide experience repeatedly shows that damages in structures caused by earthquakes are highly dependent on site condition and epicentral distance. In this paper, a 21-storey shear wall-structure built in the 1960s in Hong Kong is selected as an example to investigate these two effects. Under various design earthquake intensities and for various site conditions, the fragility curves or damage probability matrix of such building is quantified in terms of the ductility factor, which is estimated from the ratio of storey yield shear to the inter-storey seismic shear. For high-rise buildings, a higher probability of damage is obtained for a softer site condition, and damage is more severe for far field earthquakes than for near field earthquakes. For earthquake intensity of VIII, the probability of complete collapse (P) increases from 1 to 24% for near field earthquakes and from 1 to 41% for far field earthquakes if the building is moved form a rock site to a site consisting a 80 m thick soft clay. For intensity IX, P increases from 6 to 69% for near field earthquake and from 14 to 79% for far field earthquake if the building is again moved form rock site to soft soil site. Therefore, site effect is very important and not to be neglected. Similar site and epicentral effects should also be expected for other types of high-rise structures.     

A comparative study of geothermal and meteorological records of climate change in Kamchatka

To reconstruct the recent climate history in Kamchatka, a series of repeated precise temperature logs were performed in a number of boreholes located in a broad east-west strip (between 52 and 54°N) in the central part of Kamchatka west of Petropavlovsk-Kamchatski. Within three years more than 30 temperature logs were performed in 10 holes (one up to six logs per hole) to the depth of up to 400 metres. Measured temperature gradients varied in a broad interval 0 to 60 mK/m and in some holes a sizeable variation in the subsurface temperatures due to advective heat transport by underground water was observed. Measured data were compared with older temperature profiles obtained in the early eighties by Sugrobov and Yanovsky (1993). Even when older data are of poorer precision (accuracy of about 0.1 K), they presented valuable information of the subsurface temperature conditions existing 20–25 years ago. Borehole observations and the inverted ground surface temperature histories (GSTHs) used for the paleoclimate reconstruction were complemented with a detailed survey of meteorological data. Namely, the long-term surface air temperature (SAT) and precipitation records from Petropavlovsk station (in operation since 1890) were used together with similar data from a number of local subsidiary meteo-stations operating in Central Kamchatka since 1950. Regardless of extreme complexity of the local meteorological/climate conditions, diversity of borehole sites and calibration of measuring devices used during the whole campaign, the results of the climate reconstruction supported a general warming of about 1 K characteristic for the 20th century, which followed an inexpressive cooler period typical for the most of the 19th century. In the last three to four decades the warming rate has been locally increasing up to 0.02 K/year. It was also shown that the snow cover played a dominant role in the penetration of the climate “signal” to depth and could considerably smooth down the subsurface response to the changes occurred on the surface.     

Modelling atmospheric and hydrologic processes for assessment of meadow restoration impact on flow and sediment in a sparsely gauged California watershed

The restoration of meadowland using the pond and plug technique of gully elimination was performed in a 9‐mile segment along Last Chance Creek, Feather River Basin, California, in order to rehabilitate floodplain functions such as mitigating floods, retaining groundwater, and reducing sediment yield associated with bank erosion and to significantly alter the hydrologic regime. However, because the atmospheric and hydrological conditions have evolved over the restoration period, it was difficult to obtain a comprehensible evaluation of the impact of restoration activities by means of field measurements. In this paper, a new use of physically based models for environmental assessment is described. The atmospheric conditions over the sparsely gauged Last Chance Creek watershed (which does not have any precipitation or weather stations) during the combined historical critical dry and wet period (1982–1993) were reconstructed over the whole watershed using the atmospheric fifth‐generation mesoscale model driven with the US National Center for Atmospheric Research and US National Center for Environmental Prediction reanalysis data. Using the downscaled atmospheric data as its input, the watershed environmental hydrology (WEHY) model was applied to this watershed. All physical parameters of the WEHY model were derived from the existing geographic information system and satellite‐driven data sets. By comparing the prerestoration and postrestoration simulation results under the identical atmospheric conditions, a more complete environmental assessment of the restoration project was made. Model results indicate that the flood peak may be reduced by 10–20% during the wet year and the baseflow may be enhanced by 10–20% during the following dry seasons (summer to fall) in the postrestoration condition. The model results also showed that the hydrologic impact of the land management associated with the restoration mitigates bank erosion and sediment discharge during winter storm events. Copyright © 2013 John Wiley & Sons, Ltd.     

Hα emission objects in Circinus

We present a new survey for Hα emission objects in the Circinus cloud complex and introduce an efficient photometric method for detecting Hα emission via observations in a narrow‐band filter. The observed flux is compared to a blackbody fit of the continuum. Our search strategy reveals 20 stars with strong Hα emission (EW > 10 Å), eight of them being new detections. All Hα stars display infrared excess corroborating their youth. On the other hand, the region contains a number of infrared excess objects that do not show Hα emission. Our results support the picture that accretion – as witnessed by Hα emission – is a highly variable phenomenon. Therefore, Hα surveys can only trace the temporarily active objects. In contrast, infrared excess is a more robust tracer that reveals most of the population of young stellar objects in a star forming region. Our analysis of the general stellar content of the region yields a reliable distance of 450 pc for the Circinus cloud. Moreover, we find a ratio of total‐to‐selective extinction of R^V = 2.8 suggesting that smaller‐than‐normal dust grains may be present. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

P–T–t–d evolution of orogenic middle crust of the Roc de Frausa Massif (Eastern Pyrenees): a result of horizontal crustal flow and Carboniferous doming?

Structural, petrological and textural studies are combined with phase equilibria modelling of metapelites from different structural levels of the Roc de Frausa Massif in the Eastern Pyrenees. The pre‐Variscan lithological succession is divided into the Upper, Intermediate and Lower series by two orthogneiss sheets and intruded by Variscan igneous rocks. Structural analysis reveals two phases of Variscan deformation. D1 is marked by tight to isoclinal small‐scale folds and an associated flat‐lying foliation (S1) that affects the whole crustal section. D2 structures are characterized by tight upright folds facing to the NW with steep NE–SW axial planes. D2 heterogeneously reworks the D1 fabrics, leading to an almost complete transposition into a sub‐vertical foliation (S2) in the high‐grade metamorphic domain. All structures are affected by late open to tight, steeply inclined south‐verging NW–SE folds (F3) compatible with steep greenschist facies dextral shear zones of probable Alpine age. In the micaschists of the Upper series, andalusite and sillimanite grew during the formation of the S1 foliation indicating heating from 580 to 640 °C associated with an increase in pressure. Subsequent static growth of cordierite points to post‐D1 decompression. In the Intermediate series, a sillimanite–biotite–muscovite‐bearing assemblage that is parallel to the S1 fabric is statically overgrown by cordierite and K‐feldspar. This sequence points to ~1 kbar of post‐D1 decompression at 630–650 °C. The Intermediate series is intruded by a gabbro–diorite stock that has an aureole marked by widespread migmatization. In the aureole, the migmatitic S1 foliation is defined by the assemblage biotite–sillimanite–K‐feldspar–garnet. The microstructural relationships and garnet zoning are compatible with the D1 pressure peak at ~7.5 kbar and ~750 °C. Late‐ to post‐S2 cordierite growth implies that F2 folds and the associated S2 axial planar leucosomes developed during nearly isothermal decompression to <5 kbar. The Lower series migmatites form a composite S1–S2 fabric; the garnet‐bearing assemblage suggests peak P–T conditions of >5 kbar at suprasolidus conditions. Almost complete consumption of garnet and late cordierite growth points to post‐D2 equilibration at <4 kbar and <750 °C. The early metamorphic history associated with the S1 fabric is interpreted as a result of horizontal middle crustal flow associated with progressive heating and possible burial. The upright F2 folding and S2 foliation are associated with a pressure decrease coeval with the intrusion of mafic magma at mid‐crustal levels. The D2 tectono‐metamorphic evolution may be explained by a crustal‐scale doming associated with emplacement of mafic magmas into the core of the dome.     

青藏铁路唐古拉山-拉萨段全新世控震断裂研究

地表调查表明,沿青藏铁路唐古拉山-拉萨段存在5条重要的全新世控震断裂带,从北到南分别是温泉盆地西缘断裂带、安多盆地北缘断裂带、崩错断裂带、谷露西缘断裂带和当雄-羊八井断裂带.构造-地貌和年代学分析结果表明,北部的温泉盆地西缘断裂和安多盆地北缘断裂带的活动强度相对比较小,平均垂直活动速率约为0.2～0.5mm/a.南侧的谷露西缘断裂带和当雄-羊八井断裂带的全新世垂直活动速率为约(15±0.5)mm/a.而中部的崩错走滑断裂带的活动强度最大,晚第四纪期间的走滑速率可达(11±4.5)mm/a.全新世断裂活动和古地震研究表明,其中温泉盆地西缘断裂带、安多盆地北缘断裂带、崩错断裂带的西北分支、当雄-羊八井断裂带的当雄段等区域未来发生强震的概率相对更大.