基于EM算法和单幅雷达图像阴影的控制点坡度校正

对在我国现有的条件下进行控制点坡度校正的必要性进行了阐述，并分析了EM算法。根据EM算法以及基于区域增长的余弦散射模型建立了控制点坡度校正模型，并用河北省张北地区的雷达影像进行了实验，取得了较高的精度初值。     

Interevent Time Distribution in Seismicity: A Theoretical Approach

This paper presents an analysis of the distribution of the time τ between two consecutive events in a stationary point process. The study is motivated by the discovery of unified scaling laws for τ for the case of seismic events. We demonstrate that these laws cannot exist simultaneously in a seismogenic area. Under very natural assumptions we show that if, after rescaling to ensure Eτ =1, the interevent time has a universal distribution F, then F must be exponential. In other words, Corral’s unified scaling law cannot exist in the whole range of time. In the framework of a general cluster model we discuss the parameterization of an empirical unified law and the physical meaning of the parameters involved. An erratum to this article is available at .     

A Possible Detection of the 26 December 2004 Great Sumatra-Andaman Islands Earthquake with Solution Products of the International GNSS Service

The main goal of this work is to critically review the IGS solution products and Precise Point Positioning (PPP) in order to demonstrate their potential to contribute to studies of large earthquakes such as the one that devastated Southeast Asia on December 26th, 2004. In view of a possible detection of the Mw 9.0 Sumatra-Andaman Islands Earthquake of December 26, 2004, position solutions, ranging from intervals of years to one second, of four International GNSS Service (IGS) stations within 3000 km of the epicenter were examined. The IGS combined, cumulative solution product (IGS04P51), consisting of epoch and station velocity solutions and based on data spans of several years prior to the earthquake, was used as a reference. Four IGS combined weekly position solutions (igs04P1301-4), two weeks before and after the earthquake, were utilized for the weekly solution resolution. PPP static and kinematic solutions with IGS Final combined orbits and clocks were used for the mean daily and instantaneous 5-min and 1-sec epoch solutions, respectively. The most significant changes, detected by both weekly and daily solutions occurred in longitude. The nearest IGS station ntus, about 1000 km east of the epicenter, moved westward about 15 mm, while the more distant Indian station iisc (∼ 2300 km NW from the epicenter), shifted about 15 mm eastward. In spite of position errors caused by interpolation of the 5-min IGS clocks, the 1-sec solutions, based on separate data sets, available only for two stations (iisc, dgar), still showed seismic surface waves, in particular at the Indian station iisc. Precise daily IGS combined polar motion and length-of-day products, after correcting for the atmospheric effects, also likely detected, statistically significant, anomalistic excitations on December 26, 2004 that could be caused by this great earthquake.     

鲁东南岚山头超高压变质岩流体包裹体特征与板片折返史

鲁东南岚山头含柯石英榴辉岩主要产于花岗质片麻岩内,是苏鲁超高压变质带主要榴辉岩体密集分布区之一。流体包裹体研究表明,榴辉岩矿物及细脉石英中捕获有四种类型包裹体:在超高压-高压榴辉岩相条件下捕获的 N_2±CO_2包裹体;在高压榴辉宕重结晶阶段被捕获的 CO_2-H_2O 包裹体和含子矿物高盐度 H_2O 溶液包裹体:在超高压岩石折返过程中的最晚(角闪岩相退变质甚至更晚)阶段捕获的低盐度 H_2O 溶液包裹体。利用榴辉岩矿物及脉体石英中捕获的流体包裹体相互期次关系,可以对本区超高压变质作用板片折返过程中的流体演化史进行重建。     

Microthermometric measurement of fluid inclusions and its constraints on genesis of PGE-polymetallic deposits in Lower Cambrian black rock series, southern China

Systematic microthermometric measurements of fluid inclusions in the PGE-polymetallic deposits hosted in the Lower Cambrian black rock series in southern China were performed, and the results suggest: (1) there exist two types of fluid inclusions. TypeⅠis of NaCl-H2O system with low-medium salinity, and its homogenization temperatures (Th) and salinities are 106.9- 286.4℃ and ( 0.8- 21.8) wt%NaCl eq. respectively; TypeⅡ is of CaCl2-NaCl-H2O system with medium-high salinities, and its homogenization temperatures and salinities range from 120.1℃ to 269.6℃ and ( 11.4- 31.4) wt%NaCl eq., respectively. The typeⅡ fluid inclusions have been discovered for the first time in this kind of deposits; (2) two generations of ore-forming fluids were recognized. Characteristics of fluid inclusions in the PGE-polymetallic ores and carbonate-quartz stockworks in the underlying phosphorites are almost of no difference, they may represent ore-forming fluids at the main metallogenic stage. The peak value of homogenization temperature of those fluid inclusions is about 170℃, while their salinities possess a remarkable bimodal distribution pattern with two peak values of (27-31) wt%NaCl eq. and (4-6) wt%NaCl eq. On the contrary, fluid inclusions in the carbonate-quartz veins in the hanging wall may represent ore-forming fluids at the post-metallogenetic stage. The homogenization temperatures and the peak values of salinities are mostly 130-170℃ and (12-14) wt%NaCl eq., respectively; (3) nobel gas isotopic composition analyses in combination with the microthermometric measurements of fluid inclusions suggest that the ore-forming fluids at the main metallogenetic stage were probably derived from mixing of basinal hot brines with the CaCl2-NaCl-H2O system and seawater with the NaCl-H2O system; (4) in the Early Cambrian, the basinal hot brines were trapped in the Caledonian basins, which were distributed along the southern margin of the Yangtze Craton, and where giant thick sediments were accumulated, and expelled and migrated laterally along the strata because of the pressure caused by overlying sediments. The basinal hot brines absorbed Ni, Mo, V, PGE from the surrounding rocks and were transformed into ore-bearing hydrothermal fluids with the CaCl2-NaCl-H2O system and medium-high salinities, then ascended along faults and mixed with seawater of the NaCl-H2O system, and finally PGE-polymetallic deposits or occurrences were formed in the black rock series.     

Primary Ore Mineral Assemblage and Fluid Inclusion Study of the Batu Hijau Porphyry Cu-Au Deposit, Sumbawa, Indonesia

Abstract. The Batu Hijau porphyry Cu‐Au deposit, Sumbawa Island, Indonesia, is associated with a tonalitic intrusive complex. The temperature‐pressure condition of mineralization at the Batu Hijau deposit is discussed on the basis of fluid inclusion microthermometry. Then, the initial Cu‐Fe sulfide mineral assemblage is discussed. Bornite and chalcopyrite are major copper ore minerals associated with quartz veinlets. The quartz veinlets have been classified into ‘A’ veinlets associated with bornite, digenite, chalcocite and chalcopyrite, ‘B’ veinlets having chalcopyrite bornite along vuggy center‐line, rare ‘C’ chalcopyrite‐quartz veinlets, and late ‘D’ veinlets consisting of massive pyrite and quartz (Clode et al., 1999). Copper and gold mineralization is associated with abundant ‘A’ quartz veinlets. Abundant fluid inclusions are found in veinlet quartz consisting mainly of gas‐rich inclusions and polyphase inclusions throughout the veinlet types. The hydrothermal activity occurred in temperature‐pressure conditions of aqueous fluid immiscibility into hypersaline brine and dilute vapor. The halite dissolution (Tm[halite]) and liquid‐vapor homogenization (Th) temperatures of the polyphase inclusions in veinlet quartz range from 270 to 472d?C and from 280 to 454d?C, respectively. The estimated salinity ranges from 36 to 47 wt% (NaCl equiv.). The apparent pressures lower than 300 bars are estimated to have been along the liquid‐vapor‐halite curve for the fluid inclusions having the Th lower than the Tm that trapped the brine saturated with halite, or at slightly higher pressure relative to liquid‐vapor‐halite curve for the fluid inclusions having the Th higher than the Tm that trapped the brine unsaturated with halite. The actual temperature and pressure during the hydrothermal activity at the Batu Hijau deposit are estimated to have been around 300d?C and 50 bars. At such temperature‐pressure conditions, the principal and initial Cu‐Fe sulfide mineral assemblages are thought to be chalcopyrite + bornite solid solution (bnss) for the chalcopyrite‐bearing assemblage, and chalcocite‐digenite solid solution and bnss for the chalcopyrite‐free assemblage.     

Evolution of Hydrothermal System at the Dizon Porphyry Cu-Au Deposit, Zambales, Philippines

Abstract. Evolution of hydrothermal system from initial porphyry Cu mineralization to overlapping epithermal system at the Dizon porphyry Cu‐Au deposit in western central Luzon, Zambales, Philippines, is documented in terms of mineral paragen‐esis, fluid inclusion petrography and microthermometry, and sulfur isotope systematics. The paragenetic stages throughout the deposit are summarized as follows; 1) stockwork amethystic quartz veinlets associated with chalcopyrite, bornite, magnetite and Au enveloped by chlorite alteration overprinting biotite alteration, 2) stockwork quartz veinlets with chalcopyrite and pyrite associated with Au and chalcopyrite and pyrite stringers in sericite alteration, 3) stringer quartz veinlets associated with molybdenite in sericite alteration, and 4) WNW‐trending quartz veins associated with sphalerite and galena at deeper part, while enargite and stibnite at shallower levels associated with advanced argillic alteration. Chalcopyrite and bornite associated with magnetite in quartz veinlet stockwork (stage 1) have precipitated initially as intermediate solid solution (iss) and bornite solid solution (bnss), respectively. Fluid inclusions in the stockwork veinlet quartz consist of gas‐rich inclusions and polyphase inclusions. Halite in polyphase inclusions dissolves at temperatures ranging from 360d?C to >500d?C but liquid (brine) and gas (vapor) do not homogenize at <500d?C. The maximum pressure and minimum temperature during the deposition of iss and bnss with stockwork quartz veinlets are estimated to be 460 bars and 500d?C. Fluid inclusions in veinlet stockwork quartz enveloped in sericite alteration (stage 2) consist mainly of gas‐rich inclusions and polyphase inclusions. In addition to the possible presence of saturated NaCl crystals at the time of entrapment of fluid inclusions that exhibit the liquid‐vapor homogenization temperatures lower than the halite dissolution temperatures in some samples, wide range of temperatures of halite dissolution and liquid‐vapor homogenization of polyphase inclusions from 230d?C to >500d?C and from 270d?C to >500d?C, respectively, suggests heterogeneous entrapment of gaseous vapor and hypersaline brine. The minimum pressure and temperature are estimated to be about 25 bars and 245d?C. Fluid inclusions in veinlet quartz associated with molybdenite (stage 3) are dominated by gas‐rich inclusions accompanied with minor liquid‐rich inclusions that homogenize at temperatures between 350d?C and 490d?C. Fluid inclusions in vuggy veinlet quartz associated with stibnite (stage 4) consist mainly of gas‐rich inclusions with subordinate polyphase inclusions that do not homogenize below 500d?C. Fluid inclusions in veinlet quartz associated with galena and sphalerite (stage 4) are composed of liquid‐rich two‐phase inclusions, and they homogenize into liquid phase at temperatures ranging widely from 190d?C to 300d?C (suggesting boiling) and the salinity ranges from 1.0 wt% to 3.4 wt% NaCl equivalent. A pressure of about 15 bars is estimated for the dilute aqueous solution of 190d?C from which veinlet quartz associated with galena and sphalerite precipitated. In addition to a change in temperature‐pressure regime from lithostatic pressure during the deposition of iss and bnss with stockwork quartz veinlets to hydrostatic pressure during fracture‐controlled quartz veinlet associated with galena and sphalerite, a decrease in pressure is supposed to have occurred due to unroofing or removal of the overlying piles during the temperature decrease in the evolution of hydrothermal system. The majority of the sulfur isotopic composition of sulfides ranges from ±0 % to +5 %. Sulfur originated from an iso‐topically uniform and homogeneous source, and the mineralization occurred in a single hydrothermal system.     

Implicit integration under mixed controls of a breakage model for unsaturated crushable soils

This paper discusses a series of stress point algorithms for a breakage model for unsaturated granular soils. Such model is characterized by highly nonlinear coupling terms introduced by breakage‐dependent hydro‐mechanical energy potentials. To integrate accurately and efficiently its constitutive equations, specific algorithms have been formulated using a backward Euler scheme. In particular, because implementation and verification of unsaturated soil models often require the use of mixed controls, the incorporation of various hydro‐mechanical conditions has been tackled. First, it is shown that the degree of saturation can be replaced with suction in the constitutive equations through a partial Legendre transformation of the energy potentials, thus changing the thermomechanical state variables and enabling a straightforward implementation of a different control mode. Then, to accommodate more complex control scenarios without redefining the energy potentials, a hybrid strategy has been used, combining the return mapping scheme with linearized constraints. It is shown that this linearization strategy guarantees similar levels of accuracy compared with a conventional strain–suction‐controlled implicit integration. In addition, it is shown that the use of linearized constraints offers the possibility to use the same framework to integrate a variety of control conditions (e.g., net stress and/or water‐content control). The convergence profiles indicate that both schemes preserve the advantages of implicit integration, that is, asymptotic quadratic convergence and unconditional stability. Finally, the performance of the two implicit schemes has been compared with that of an explicit algorithm with automatic sub‐stepping and error control, showing that for the selected breakage model, implicit integration leads to a significant reduction of the computational cost. Such features support the use of the proposed hybrid scheme also in other modeling contexts, especially when strongly nonlinear models have to be implemented and/or validated by using non‐standard hydro‐mechanical control conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

A new magnitude category disaggregation approach for temporal high‐resolution rainfall intensities

Simulation of quick runoff components such as surface runoff and associated soil erosion requires temporal high‐resolution rainfall intensities. However, these data are often not available because such measurements are costly and time consuming. Current rainfall disaggregation methods have shortcomings, especially in generating the distribution of storm events. The objectives of this study were to improve point rainfall disaggregation using a new magnitude category rainfall disaggregation approach. The procedure is introduced using a coupled disaggregation approach (Hyetos and cascade) for multisite rainfall disaggregation. The new procedure was tested with ten long‐term precipitation data sets of central Germany using summer and winter precipitation to determine seasonal variability. Results showed that dividing the rainfall amount into four daily rainfall magnitude categories (1–10, 11–25, 26–50, >50 mm) improves the simulation of high rainfall intensity (convective rainfall). The Hyetos model category approach (Hyetos_Cat) with seasonal variation performs representative to observed hourly rainfall compared with without categories on each month. The mean absolute percentage accuracy of standard deviation for hourly rainfall is 89.7% in winter and 95.6% in summer. The proposed magnitude category method applied with the coupled Hyetos_Cat–cascade approach reproduces successfully the statistical behaviour of local 10‐min rainfall intensities in terms of intermittency as well as variability. The root mean square error performance statistics for disaggregated 10‐min rainfall depth ranges from 0.20 to 2.38 mm for summer and from 0.12 to 2.82 mm for the winter season in all categories. The coupled stochastic approach preserves the statistical self‐similarity and intermittency at each magnitude category with a relatively low computational burden. Copyright © 2014 John Wiley & Sons, Ltd.