太阳活动区磁场与CMEs和太阳质子事件

文中选了５ 个典型活动区, 分析了这些活动区的磁场, 与活动区相应的ＣＭＥｓ, 太阳爆发事件和太阳质子事件我们发现, 对于Ｅ ≥１０ｍｅＶ 的太阳质子事件有相应的源活动区, 源耀斑和ＣＭＥ; 活动区矢量磁场有剪切, 磁场剪切越强质子事件越强; 多数在质子耀斑发生前出现磁流浮现; 太阳１０ｃｍ 射电爆发持续时间长文中结果还佐证了Ｓｈｅａｌｙ 等的结果： Ｘ 射线耀斑的长持续时间与ＣＭＥ 的发生正相关另外,在５ 个活动区中, 有三个大耀斑发生前没有明显的磁剪切作为它们的先兆, 它们是非质子源耀斑这是Ｍｏｏｒｅ, Ｈａｇｙａｒｄ 和Ｄａｖｉｓ 的磁场强剪切是耀斑产生的必要条件的反例     

太阳质子事件、双带耀斑及日冕物质抛射

本文统计了第２２ 太阳活动周期间（１９９１ ～１９９５ 年） 发生的２５ 个太阳质子事件与太阳耀斑及日冕物质抛射（ＣＭＥ） 事件的关系　　统计结果表明, 所有的太阳质子事件都与耀斑发生相关, 除２ 个质子事件（１９９４１０２０ 和１９９５１０２０ 日发生的太阳质子事件） 与ＣＭＥ发生无关, 其余质子事件也都与ＣＭＥ 相关　　值得注意的是, 与质子事件相关的耀斑有１６ 个是双带耀斑, 其中包括与ＣＭＥ无关的２ 个事件的耀斑, 占总数的６４ ％ 　　上述统计结果证实了无论是太阳耀斑, 还是物质抛射, 它们对太阳质子事件的发生同样起着非常重要的作用     

太阳耀斑电子束在日冕区产生的静电不稳定性

依据太阳耀斑爆发特征,建立了双耀斑电子束与日冕背景相作用的模型,数值结果表明,该等离子体系统将激发静电不稳定性,其时间增长率ω_i受耀斑热束密度与日冕背景密度比值（n_h/n_o）以及耀斑冷束相对论电子密度与日冕密度值（n_c/n_o）影响较大,并随它们增大而增大,其实频大小在耀斑热束等离子体频率附近。因此,此系统可激发大于1GHz的高频静电辐射,这些结果对揭示耀斑粒子在日冕空间传播行为有一定作用,并可用于探讨高频Ⅲ型射电机制。     

日冕物质抛射的能量来源

观测表明, 日冕物质抛射（ＣＭＥ） 伴随着闭合磁场形成开放的过程理论研究指出, 磁场在ＣＭＥ 过程中起主导作用, 而且对应于相同的边界条件, 开放场中储存的磁能最多一个系统如何能在向外输送能量的同时增加自身的能量？ 本文就该问题进行了讨论     

2000年7月14日太阳高能粒子事件引起地球同步轨道区相对论电子通量巨幅增加

2000年7月14日10:24UT一个X5.6级的耀斑暴发生在太阳中心子午线附近(AR 9077), 同时伴随着一个朝向地球的CME事件及太阳高能粒子(Solar Energetic Particle, SEP)事件. 这次耀斑暴发及CME事件引起了地球磁层、电离层及高层大气的强烈扰动. 中国“风云二号(B)”卫星上的高能粒子探测器(EPD)观测到SEP事件期间, 同步轨道区高能质子、相对论电子有非常剧烈的增加. SEP期间, 高能质子对相对论电子通量的探测造成严重的污染. 结合“风云二号(B)”卫星上的高能粒子探测器(EPD)的特性, 建立了一种从相对论电子通量探测中“清除”高能质子“污染”的方法, 并对相对论电子通量的探测数据实施“清洁”处理. “纯净的”相对论电子通量探测结果显示, 当行星际磁场南向时, 上游太阳风中的高能电子使同步轨道区相对论电子通量有大幅度的增加.     

CME和冕洞边缘结构的可能关系

本文分析了１９８０ 、１９８４ 和１９８９ 年ＳＭＭ 卫星观测到的１４０ 次日冕物质抛射（ＣＭＥ） 事件在时空分布上与“冕洞边缘结构”、耀斑爆发和爆发日珥等事件的相关关系结果表明, ＣＭＥ 事件与日冕边缘结构的关系最密切此外, ＣＭＥ与赤道冕洞具有同步的长期演化关系由此认为, 冕洞边缘结构对ＣＭＥ的可能贡献是不可忽视的     

2001年4月2日太阳耀斑及其太阳质子事件的观测结果研究

2001年4月2日, 太阳爆发了一个近年来X射线通量最大的一次耀斑并伴有质子事件, 利用“资源一号”卫星星内粒子探测器和神舟二号飞船X射线探测器的观测资料, 对这一事件的高能粒子响应进行了特例研究. “资源一号”卫星运行于太阳同步轨道, 高度约800km, 和宁静时期的统计结果对比, 这次耀斑后, 星内粒子探测器在地球极盖区（地球开磁场区）观测到耀斑粒子的出现, 这是宁静时期没有的; 神舟二号飞船轨道高度400km, 倾角为42°, X射线探测器在42°中高纬地区也观测到高能电子通量比宁静时明显的增加, 这表明, 太阳耀斑引起的近地空间辐射环境的变化遍及纬度约40°以上的区域, 甚至在40°N附近400 km左右的高度上仍然有响应. 但是, 中高纬度、极光带和极盖区的粒子来源, 加速机制和响应方式却不一定相同, 需要分别讨论. 资料分析和对比还表明, 质子事件的强度并不一定和耀斑的X射线通量成正比, 因此, 近地空间高能粒子对耀斑的响应也不是完全决定于X射线强度.     

地球同步轨道附近哨声湍流对“种子电子”的加速

本文在等离子体准线性理论下研究了地球同步轨道附近哨声湍流对亚暴“种子电子”的波-电子共振相互作用. 当发生这种共振时，“种子电子”的动量分布函数经动量扩散而随时间演化，部分低能电子数减少了，而高能尾部分的相对论电子（能量大于1MeV）数增加了，说明“种子电子”得到了哨声湍流的有效加速，且哨声湍流的能量越高，其加速效率越高. 另外，哨声湍流的频率越低（或波数越小），共振电子的能量越高（或单位质量的动量越大）；频率范围越宽，共振电子的能量范围越宽，被加速的电子数也越多. 磁层哨声湍流加速“种子电子” 大约在30h内就可以造成相对论电子数显著增加，这正好和大多数磁暴期间观测到的相对论电子通量的增长时间相当.     

高能电子爆发与绕月卫星表面电位大幅下降的联动效应

"嫦娥"一号、二号绕月飞行经历地球磁尾边界层区域时,分别在2007年11月26日—2008年2月5日和2010年10月3日—2011年2月28日,发现了15次月球轨道0.1~2 MeV电子急剧增加(Bursts of 0.1~2 MeV Energetic Electrons,BEE),卫星周围等离子体离子加速的现象.统计研究表明,这类现象发生在稳定太阳风和弱行星际磁场条件下,且无显著空间环境扰动事件发生时,离子的加速滞后于高能电子爆发,离子能量的变化与高能电子通量的时间演化正相关,地球磁鞘内侧或边界层过渡区域是该类现象的高发区,离子能量增加时卫星表面电位大幅下降可达负几千伏.为了研究高能电子爆发与绕月卫星表面电位变化的关系及其对月球表面电位的影响,本文用电流平衡法建立绕月卫星和月球表面充电模型,并假设能量电子(2eV~2 MeV)满足幂律谱的分布,模拟急剧增加的能量电子对卫星和月球表面电位的影响.模拟结果表明,能量电子急剧增加使得绕月卫星和月球表面电位大幅下降;能量电子总流量1011 cm-2时,绕月卫星和月球表面充电电位可达负上千伏;月球充电到大的负电位的时间仅为卫星充电时间的1/10.鉴于高能电子急剧增加事件的高发生率(~125次/年),能量电子急剧增加使得绕月卫星表面电位大幅下降的发生率应大于实测等离子体离子加速现象的发生率(~25次/年).     

太阳风湍流和磁层亚暴的一种机制

太阳风的动量涨落将通过磁层边界在磁尾激发磁流体力学波。快磁声波携带扰动能量传到等离子体片中,发展为激波,或者通过激波的相互作用而耗散能量,使等离子体加热。等离子体片中的随机费米加速机制,使麦克斯韦分布尾巴部分的高能量粒子被加速到更高能。在宁静态时,加热、加速与耗散过程平衡。当太阳风的动量或者其涨落较大时,整个加热和加速过程加剧,更多的高能粒子产生,并从等离子体片中逃逸,形成高速的等离子体流注入近地轨道和极区,表现为磁层亚暴过程。利用这种机制,可以解释地球磁层亚暴的定性特征。     

A Monte Carlo study of rainfall forecasting with a stochastic model

A procedure for short-term rainfall forecasting in real-time is developed and a study of the role of sampling on forecast ability is conducted. Ground level rainfall fields are forecasted using a stochastic space-time rainfall model in state-space form. Updating of the rainfall field in real-time is accomplished using a distributed parameter Kalman filter to optimally combine measurement information and forecast model estimates. The influence of sampling density on forecast accuracy is evaluated using a series of a simulated rainfall events generated with the same stochastic rainfall model. Sampling was conducted at five different network spatial densities. The results quantify the influence of sampling network density on real-time rainfall field forecasting. Statistical analyses of the rainfall field residuals illustrate improvement in one hour lead time forecasts at higher measurement densities.     

Earthquake Engineering and Engineering Vibration Aims And Scope

正This journal is established by the Institute of Engineering Mechanics(IEM),China Earthquake Administration,to promote scientific exchange between Chinese and foreign scientists and engineers so as to improve the theory and practice of earthquake hazards mitigation,preparedness,and recovery.To accomplish this purpose,the journal aims to attract a balanced number of papers between Chinese and     

Solute concentrations in the lower niger river and the source rock contribution

Water quality analyses for the Niger River for the 1980/81 hydrological year are presented. The samples were collected from the main river at Lokoja, and from two main tributaries, the Kaduna and the Benue Rivers. Different water types were distinguished by the concentrations of major ions. The type Ca > Na > Mg > K - HCO₃ > SO₄ > Cl was represented at all stations during at least part of the year. Chloride was found to dominate the sulphate ion in the Kaduna and Niger, while the Benue maintained a higher concentration of sulphate relative to chloride all year round. Distinct patterns of seasonal variation in the ion concentrations were observed, particularly for the samples collected at Lokoja. Low ion concentrations were prominent during periods of high discharge, while low flow periods coincided with high dissolved ion concentrations. The contribution of rainwater to the total dissolved solids in the river waters was assessed indirectly using rainwater chemistry data from the Gulf of Guinea. The estimated rainwater contribution to the Lower Niger amounts to 5.15 mg 1^?1. Geochemical weathering calculations involving reactions of the four major minerals of granitic rocks - anorthite, biotite, albite, and K-feldspar - with carbon dioxide and water, can account for the average water composition of the Lower Niger. The proportion of the ionic components was also related to the occurrence of the respective element in the minerals.     

Scenario predictions for potential near-field and far-field earthquakes affecting Hong Kong

Historical earthquakes noted in the written records of the South China region, including Hong Kong, are not well delineated along identified prominent fault sources. Despite the lack of any definitive, localised trend in the spatial distribution of seismic activity in the region, there does appear to be some major disparity in the seismic activity rates (especially for large magnitude earthquakes) between the near-field and the far-field regions of Hong Kong. Despite this observation, previous studies of the regional seismic activity and seismic ground motion hazard (the latter using a probabilistic seismic hazard assessment, PSHA) have considered very broad source zone regions, in which uniform levels of seismic activity have been assumed. The present paper further scrutinises this broad source zone (BSZ) approach by adopting a novel expanding circular disc (ECD) method to determine the rates of earthquake recurrence. Such a method is intended to counter-check previously developed models by determining earthquake scenario events in terms of magnitude–distance (M–R) pairs or combinations, having defined values of average return period. Unlike the BSZ approach, the ECD method specifically accounts for the supposed variations in the seismic activity rates between events in the near-field and the far-field of Hong Kong.The form of the developed method is particularly suited to the determination of design-level earthquake ground motions for bedrock sites, since it assumes a directionally-independent attenuation model as described in the companion paper. It is found that, whilst the BSZ approach may indicate the overall average levels of hazard that are representative of the South China region as a whole, it does not capture the large disparity in seismic activity rates between near-field and far-field events. This important feature is expected to have a significant impact upon engineering assessments of the seismic safety of structures in Hong Kong and elsewhere in the South China region. For example, it is found that for events with M≥6, the seismic activity rate (normalised by time and area) in the very far-field is around 3.5 times larger than in the near-field and medium-field of Hong Kong. The resulting design M-R combinations, covering a range of return periods from 70 to 2500 years, are limited, for very long return periods and for distant events, by the maximum credible earthquake (MCE) magnitude. Intensive research to determine this seismic hazard parameter is recommended, in order to refine further the results of the ECD analysis, which presently conservatively assumes the MCE to range between M=6 in the near-field of Hong Kong to M=8 in the very far-field, at distances greater than 280 km from Hong Kong.     

Depositional ages and characteristics of Middle–Upper Jurassic and Lower Cretaceous lacustrine deposits in southeastern Mongolia

Lower Cretaceous lacustrine oil shales are widely distributed in southeastern Mongolia. Due to the high organic carbon content of oil shale, many geochemical studies and petroleum exploration have been conducted. Although most of the oil shales are considered to be Early Cretaceous in age, a recent study reveals that some were deposited in the Middle Jurassic. The present study aims at establishing depositional ages and characteristics of the Jurassic and Cretaceous lacustrine deposits in Mongolia. The Lower Cretaceous Shinekhudag Formation is about 250 m thick and composed of alternating beds of shale and dolomite. The Middle Jurassic Eedemt Formation is about 150 m thick and composed of alternating beds of shale, dolomitic marl, and siltstone. The alternations of shale and dolomite in both formations were formed by lake level changes, reflecting precipitation changes. Shales were deposited in the center of a deep lake during highstand, while dolomites were formed by primary precipitation during lowstand. Based on the radiometric age dating, the Shinekhudag Formation was deposited between 123.8 ±2.0 Ma and 118.5 ±0.9 Ma of the early Aptian. The Eedemt Formation was deposited at around 165–158 Ma of Callovian–Oxfordian. The calculated sedimentation rate of the Shinekhudag Formation is between 4.7 ±2.6 cm/ky and 10.0 ±7.6 cm/ky. Shales in the Shinekhudag Formation show micrometer‐scale lamination, consisting of algal organic matter and detrital clay mineral couplets. Given the average thickness of micro‐laminae and calculated sedimentation rate, the micro‐lamination is most likely of varve origin. Both Middle–Upper Jurassic and Lower Cretaceous lacustrine oil shales were deposited in intracontinental basins in the paleo‐Asian continent. Tectonic processes and basin evolution basically controlled the deposition of these oil shales. In addition, enhanced precipitation under humid climate during the early Aptian and the Callovian–Oxfordian was another key factor inducing the widespread oil shale deposition in Mongolia.     

THEORETICAL APPROACH TO THE SOLUTION OF THE INFILTRATION PROBLEM

ABSTRACT

The one-dimensional transient downward entry of water in unsaturated soils is investigated theoretically. The mathematical equation describing the infiltration process is derived by combining Darcy's dynamic equation of motion with the continuity and thermodynamic state equations adjusted for the unsaturated flow conditions. The resulting equation together with the corresponding initial and boundary conditions constitues a mathematical initial boundary value problem requiring the solution of a nonlinear partial differential equation of the parabolic type. The volumetric water content is taken as the dependent variable and the time and the position along the vertical direction are taken as the independent variables. The governing equation is of such nature that a solution exists for t > 0 and is uniquely determined if two relationships are defined, together with the specified state of the system, at the initial time t = 0 and at the two boundaries. The two required relations are those of pressure versus permeability and pressure versus volumetric water content.

Since the partial differential equation has strong non-linear terms, a discrete solution is obtained by approximating the derivatives with finite-differences at discrete mesh points in the solution domain and integrated for the corresponding initial and boundary conditions. The use of an implicit difference scheme is employed in order to generate a system of simultaneous non-linear equations that has to be solved for each time increment. For n mesh points the two boundary conditions provide two equations and the repetition of the recurrence formula provides n—2 equations, the total being n equations for each time increment. The solution of the system is obtained by matrix inversion and particularly with a back-substitution technique. The FORTRAN statements used for obtaining the solution with an electronic digital computer (IBM 704) are presented together with the input data.

Analysis of the errors involved in the numerical solution is made and the stability and convergence of the solution of the approximate difference equation to that of the differential equation is investigated. The method applied is that of making a Fourier series expansion of a whole line of errors and then following the progress of the general term of the series expansion and also the behavior of each constituent harmonic. The errors (forming a continuous function of points in an abstract Banach space) are represented by vectors with the Fourier coefficients constituting a second Banach space. The amplification factor of the difference equation is shown to be always less than unity which guarantees the stability of the employed implicit recurrence scheme.

Experiments conducted on a vertical column packed uniformly with very fine sand, show a satisfactory agreement between the theoretically and experimentally obtained values. Many experimental results are shown in an attempt to explain the infiltration phenomenon with emphasis on the shape and movement of the wet front, and the effects of the degree of compaction, initial water content and deaired water on the infiltration rate.     

Late Permian post‐ophiolitic trondhjemites from Central Iran: a mark of subduction role in growth of Paleozoic continental crust

Late Permian trondhjemites in the Anarak area occur as stocks and dykes, which cross cut the Anarak ophiolite and its overlying metasedimentary rocks, and are exposed along the northern Anarak east–west main faults. These leucocratic intrusive bodies have enclaves of all ophiolitic units and metamorphic rocks. They are composed of amphibole, plagioclase (oligoclase), quartz, zircon and muscovite. Secondary minerals are chlorite (pycnochlorite), epidote, albite, magnetite and calcite. Whole‐rock major‐ and trace‐element analyses reveal that they are characterized by high SiO ₂ (67.8–71.0 wt%), Al ₂ O ₃ (14.9–17.1 wt%) and Na ₂ O (5.3–8.6 wt%), low K ₂ O (0.1–1.5 wt%; average: 0.8 wt%), low Rb/Sr ratio (0.01–0.40; average: 0.09), low Y (3–6 ppm), negative Ti, Nb and Ta anomalies, slightly negative or positive Eu anomaly, LREE enrichment and fractionated HREE. These rocks present 2 to 40 times enrichment in inclined chondrite‐normalized REE patterns. Geochemical characteristics of the Anarak trondhjemites all reflect melting of a mafic protolith at more than 10 kbar. The field evidence and whole‐rock chemistry reveal that these rocks have been crystallized from magmas derived from melting of subducted Anarak oceanic crust. This study reveals that melting of garnet amphibolite was an important element of continent formation in the study area.