TRMM/PR资料在中尺度模式MM5中的应用(英)

将中尺度非静力模式MM5中的积云参数化Grell方案作了改进,使它含有降水云的云水,即雨水含量,并用该模式对1998年6月29日08时-30日08时(北京时)和1998年7月1日08时-2日08时(北京时)发生在淮河流域的两次特大暴雨进行数值模拟研究,同时通过采用R-qr关系将TRMM/PR(Tropical Rainfall MeasuringMission/Precipition Radar)得到的降水强度资料R计算出比含水量qr,然后用q'v=qv+qr取代原模式中的比湿qv。结果表明将TRMM/PR资料加入模式后,由于PR雷达具有较高的空间分辨率,能够很好地反映中小尺度系统的空间结构,能够使湿度值接近实际,缩短了中尺度系统的发展时间,使得模拟出来的降雨强度、雨量中心位置以及雨带形状更接近实况。     

Impacts of recent El Niño Modoki on dry/wet conditions in the Pacific rim during boreal summer

Present work uses 1979–2005 monthly observational data to study the impacts of El Niño Modoki on dry/wet conditions in the Pacific rim during boreal summer. The El Niño Modoki phenomenon is characterized by the anomalously warm central equatorial Pacific flanked by anomalously cool regions in both west and east. Such zonal SST gradients result in anomalous two-cell Walker Circulation over the tropical Pacific, with a wet region in the central Pacific. There are two mid-tropospheric wave trains passing over the extratropical and subtropical North Pacific. They contain a positive phase of a Pacific-Japan pattern in the northwestern Pacific, and a positive phase of a summertime Pacific-North American pattern in the northeastern Pacific/North America region. The western North Pacific summer monsoon is enhanced, while the East Asian summer monsoon is weakened. In the South Pacific, there is a basin-wide low in the mid-latitude with enhanced Australian high and the eastern South Pacific subtropical high. Such an atmospheric circulation pattern favors a dry rim surrounding the wet central tropical Pacific. The El Niño Modoki and its climate impacts are very different from those of El Niño. Possible geographical regions for dry/wet conditions influenced by El Niño Modoki and El Niño are compared. The two phenomena also have very different temporal features. El Niño Modoki has a large decadal background while El Niño is predominated by interannual variability. Mixing-up the two different phenomena may increase the difficulty in understanding their mechanisms, climate impacts, and uncertainty in their predictions.     

A Method Solving Kepler’s Equation for Hyperbolic Case

We developed a method to solve Keplers equation for the hyperboliccase. The solution interval is separated into three regions; F 1, F 1, F 1. For the region F is large, we transformed the variablefrom F to exp F and applied the Newton method to the transformed equation.For the region F is small, we expanded the equation with respect to F andapplied the Newton method to the approximated equations. For the middleregion, we adopted a discretization of the Newton method and used the Taylorseries expansion for the evaluation of transcendental functions (Fukushima,1997). Numerical measurements showed that, in the case of Intel Pentiumprocessor, the new method is more than 3.7 times as fast as the combinationof a fourth order correction formula (Fukushima, 1997) and Roysstarting procedure (Roy, 1988) and others.     

Impact of Mascarene High variability on the East African ‘short rains’

The interannual variability of East African ‘short rains’ (EASR) and its link with the Mascarene High (MH) variation are explored, using observations and reanalysis data. Correlation and composite analyses for flood and drought events reveal that the EASR variability is strongly linked to the MH zonal displacement, in particular, the zonal movement of the MH eastern ridge. When the MH eastern ridge is anomalously displaced to the west (east) of its normal position, the south east (SE) trade winds over the South Indian Ocean (SIO) anomalously strengthen (weaken). This enhances (reduces) the relatively cool and dry SE trade winds and induces cold (warm) sea surface temperature anomaly in the SIO. As a result, convection over the western equatorial SIO is suppressed (enhanced) and leads to rainfall deficits (excess) over East Africa. Droughts in East Africa are associated with a westward migration of the MH eastern ridge, while the relationship is less clear for flood events and their link to an eastward migration of the MH. Therefore, the zonal migration of the MH eastern ridge provides a novel indicator for the EASR extremes especially droughts. This revelation has immense social application for rainfall forecast over East Africa where rainfall deficits have become more prevalent against the background of deteriorating conventional forecasts for EASR droughts.     

Tropical Indian Ocean variability revealed by self-organizing maps

The tropical Indian Ocean climate variability is investigated using an artificial neural network analysis called self-organizing map (SOM) for both observational data and coupled model outputs. The SOM successfully captures the dipole sea surface temperature anomaly (SSTA) pattern associated with the Indian Ocean Dipole (IOD) and basin-wide warming/cooling associated with ENSO. The dipole SSTA pattern appears only in boreal summer and fall, whereas the basin-wide warming/cooling appears mostly in boreal winter and spring owing to the phase-locking nature of these phenomena. Their occurrence also undergoes significant decadal variation. Composite diagrams constructed for nodes in the SOM array based on the simulated SSTA reveal interesting features. For the nodes with the basin-wide warming, a strong positive SSTA in the eastern equatorial Pacific, a negative Southern Oscillation, and a negative precipitation anomaly in East Africa are found. The nodes with the positive IOD are associated with a weak positive SSTA in the central equatorial Pacific or positive SSTA in the eastern equatorial Pacific, a positive (negative) sea level pressure anomaly in the eastern (western) tropical Indian Ocean, and a positive precipitation anomaly over East Africa. The warming in the central equatorial Pacific appears to correspond to El Niño Modoki discussed recently. These results suggest usefulness of SOM in studying large-scale ocean–atmosphere coupled phenomena.     

Accurate computation of gravitational field of a tesseroid

We developed an accurate method to compute the gravitational field of a tesseroid. The method numerically integrates a surface integral representation of the gravitational potential of the tesseroid by conditionally splitting its line integration intervals and by using the double exponential quadrature rule. Then, it evaluates the gravitational acceleration vector and the gravity gradient tensor by numerically differentiating the numerically integrated potential. The numerical differentiation is conducted by appropriately switching the central and the single-sided second-order difference formulas with a suitable choice of the test argument displacement. If necessary, the new method is extended to the case of a general tesseroid with the variable density profile, the variable surface height functions, and/or the variable intervals in longitude or in latitude. The new method is capable of computing the gravitational field of the tesseroid independently on the location of the evaluation point, namely whether outside, near the surface of, on the surface of, or inside the tesseroid. The achievable precision is 14–15 digits for the potential, 9–11 digits for the acceleration vector, and 6–8 digits for the gradient tensor in the double precision environment. The correct digits are roughly doubled if employing the quadruple precision computation. The new method provides a reliable procedure to compute the topographic gravitational field, especially that near, on, and below the surface. Also, it could potentially serve as a sure reference to complement and elaborate the existing approaches using the Gauss–Legendre quadrature or other standard methods of numerical integration.     

Coexistence of a Fluid Phase with Granitic Magma: Geochemical Characteristics of REE and Cu in Miyako Granite and in Scheelite-bearing Aplitic Veins at the Yamaguchi Cu-W Skarn Deposit, Iwate, Japan

Abstract: The North granitic body of the Miyako pluton is located in the Northern Kitakami belt, Northeast Japan. The formation of the scheelite–chalcopyrite–magnetite–bearing aplitic veins and scheelite–chalcopyrite–magnetite–bearing Yamaguchi skarn deposit was closely associated with the formation of the Miyako plutons. Petrographic facies of the North granitic body vary from quartz diorite in marginal zone (zone A), to tonalite and granodiorite (zone B), and to granite (zone C) in the central. The large numbers of aplitic veins distributed around the Yamaguchi mining area are divided into two groups: barren and scheelite–mag–netite–chalcopyrite–bearing aplitic veins. The latter cut massive clinopyroxene skarns of the Yamaguchi deposit, and are composed of plagioclase, K‐feldspar and titanite. Some plagioclase crystals have dusty cores with irregularly shaped K‐feldspar flakes, and clear rims of albite. Textures of plagioclase in the mineralized aplitic veins are different from the idiomorphic textures with sharp plagioclase crystal boundaries that occur in the North granitic body and barren aplitic veins. These textural data suggest that the mineralized aplitic veins were formed from hydrothermal fluid. Changes in the contents of major and minor (Rb, Sr, Sc, Co, Th, U) elements in the North Miyako granitic body are similar to those of zoned plutons formed by typical magmatic differentiation processes. On the other hand, concentrations of REE, especially middle to heavy REE, of granitic rocks in zone C and barren aplitic veins are significantly lower than those of granitic rocks in zones A and B. The hypothetical chondrite‐normalized REE patterns, calculated assuming fractional crystallization from zone B granitic melt, suggest that REE concentrations of the residual melt increased with the degree of fractional crystallization, and changed into a pattern with enriched LREE and strongly negative Eu anomaly. However, the REE patterns of granitic rocks in zone C are different from the hypothetical patterns. Moreover, the REE patterns of magnetite–scheelite–chalcopyrite aplitic veins are quite different from those of granitic rocks. The Cu contents of granitic rocks in the North Miyako body increase from zone A (5–26 ppm) to zone B (10–26 ppm), and then clearly decrease to zone C (5–7 ppm) and drastically increase to the barren aplitic veins (39–235 ppm). Concentrations of Cu in the mineralized aplitic veins are also higher than those of the granitic rocks in zone C. The decrease in REE and Cu contents of granitic rocks from zone B to zone C is not a result of simple magmatic fractional differentiation. Fluid inclusions in quartz from mineralized aplitic veins contain 3.3 wt% NaCl equivalent and 5.8 wt% CO₂. It was also demonstrated experimentally that the removal of MREE and HREE by fluid from melt enabled the formation of complexes of REE and ligands of OH^‐ and CO₃^2‐. Based on the possibility that the melt of the granitic rocks of zone C and the mineralized aplitic veins coexisted with CO₂‐bearing fluid, it is thought that REE were extracted from the melt to the CO₂‐bearing fluid, and that the REE in the mineralized aplitic veins were transported by the CO₂‐bearing fluid. It is likely that the low HREE and Cu contents of the granitic rocks in zone C could have been caused by the removal of those elements from the granitic melt by the fluid coexisting with the melt. The expelled materials could have been the sources of scheelite–magnetite–chalcopyrite–bearing aplitic veins and copper mineralization of the Yamaguchi Cu‐W skarn deposit.     

Comparison of potentially toxic metals leaching from weathered rocks at a closed mine site between laboratory columns and field observation

Potentially toxic metals, such as Cu, Pb and Zn, are leached from weathered rocks at many closed mine sites due to the acidic environments and mineralogical modifications. The mobilized toxic metals may cause contamination of surrounding water bodies. In this study, both laboratory column experiments and field observations at a former mine located in the north of Japan were carried out to compare the leaching behavior of Cu, Pb and Zn. The thickness of the surface weathered rock was varied (10, 20 and 30 cm) for the columns experiments while porous cups for porewater sampling were set up at different depths (GL-15, -45, -70, and -95 cm) for the field observations. Deionized water was added once a week over 75 weeks to the columns to simulate rainfall while porewater was extracted by a vacuum pump in several sampling campaigns (over 18 months). Similar low pH and leaching behavior of potentially toxic metals were observed for column experiments and field observations. A moderate increase in concentration with depth was observed for Cu and Zn. However, no increase in concentration was observed for Pb. This suggests that the leaching of Cu and Zn is enhanced by the length of the flow pathway through the weathered rock layer while Pb concentration is restricted by the precipitation of insoluble Pb sulfate. Thus, the thickness of the weathered rock layer is an important parameter that should be taken into consideration to estimate the loads of some potentially toxic metals, which is important when designing remediation schemes.     

Hydrothermal Alteration of Oman Ophiolite Extrusives in Ghuzayn Area

Abstract. The hydrothermal alteration in Ghuzayn Volcanics was associated with mineralization and accumulation of three massive sulfide deposits. The Ghuzayn Volcanics were discriminated into basaltic and andesitic lavas. The crossplots of Zr versus Nb, Y, Hf, La and Lu show that they fall in the same linear fractionation trends with more evolved affinities in the andesitic lavas compared to the less evolved affinities in the basaltic lavas. The immobile trace element contents of the Ghuzayn Volcanics show that they are of oceanic to continental tholeiitic affinity which fall in the fields of back‐arc basin basalt, N‐MORB and island arc tholeiite all meet and match with the field of the Lau‐Tonga back‐arc basin lavas. The ore‐body No. 2 is underlain by a zone of intense argillization and silicification which are so‐called alteration pipe as in some of the Cyprus‐type massive sulfide deposits. The alteration in the proximity of the ore‐body No. 2 is characterized by a zone of epidotization and slight silicification and sulfide disseminations, and surrounded by another zone of slight silicification and sulfide dissemination in the outermost rim around the ore‐body No. 2. The alteration zones in the proximity of the ore‐body No. 2 were divided into Zones I and II based on the abundance of the secondary minerals in the altered basaltic lavas. The Zone I is located in the most northern upper part of the ore‐body No. 2 and characterized by corrensite, saponite and prehnite. However, the Zone II is located in the southern part of the ore‐body No. 2 and characterized by chlorite and epidote. Both alteration zones were formed by different alteration stages. Stage 1 was formed by non to partly reacted fluids, to crystallize Mg‐chlorite and albite at temperatures ranging from 150 to 250d?C. Further, these fluids evolved with time and became Mg‐depleted, Si‐ and metal‐enriched to crystallize Fe²⁺‐chlorite, epidote, quartz and sulfides at temperatures ranging from 250 to 350d?C in stage 2. Later, these fluids were evolved again to be Mg‐ and Ca‐enriched in stage 3 to crystallize prehnite, laumontite, corrensite and saponite at temperatures ranging from 150 to 220d?C. Late zeolite and calcite have overgrown the previously crystallized phases and crosscutting veins along the altered basaltic lavas in stage 4. Finally, the lavas were cracked and refractured to facilitate penetration of seawater into deeper parts to heat up again and re‐leach the silica and metals in a new mineralizing event.     

Rectangular rotation of spherical harmonic expansion of arbitrary high degree and order

In order to move the polar singularity of arbitrary spherical harmonic expansion to a point on the equator, we rotate the expansion around the y-axis by \(90^{\circ }\) such that the x-axis becomes a new pole. The expansion coefficients are transformed by multiplying a special value of Wigner D-matrix and a normalization factor. The transformation matrix is unchanged whether the coefficients are \(4 \pi \) fully normalized or Schmidt quasi-normalized. The matrix is recursively computed by the so-called X-number formulation (Fukushima in J Geodesy 86: 271–285, 2012a). As an example, we obtained \(2190\times 2190\) coefficients of the rectangular rotated spherical harmonic expansion of EGM2008. A proper combination of the original and the rotated expansions will be useful in (i) integrating the polar orbits of artificial satellites precisely and (ii) synthesizing/analyzing the gravitational/geomagnetic potentials and their derivatives accurately in the high latitude regions including the arctic and antarctic area.