The crystal structure of (Fe<Subscript>4</Subscript>Cr<Subscript>4</Subscript>Ni)<Subscript>9</Subscript>C<Subscript>4</Subscript>

(Fe4Cr4Ni)9C4 is a metal carbide mineral formed by combination of Fe, Cr and Ni with C. It occurs in a chromite deposit in the Luobusha ophiolite, Tibet. Based on the determination of its crystal structure, the empirical formula is (Fe4.12Cr3.84Ni0.96)8.92C3.70 and the simplified formula is (Fe4Cr4Ni)4C9. The mineral is hexagonal with a = 1.38392(2) nm, c = 0.44690(9) nm, pace group P63 m c, Z=6 and the calculated specific gravity Dx = 7.089 g/cm3. Fe, Cr and Ni occupy different crystallographic sites and their coordination numbers are approximately 12, forming an alternate stacking sequence of flat and puckered layers along the c axis. Some metallic atoms have a defect structure. The interatomic distances of Fe, Cr and Ni are 0.2525-0.2666 nm, and the distances between Fe, Cr, Ni and C are 0.1893-0.2169 nm. The coordination number of carbon is 6. It occurs in interstices of the metallic atoms Fe, Cr and Ni to form trigonalprismatically coordinated polyhedra. These coordination polyhedra are linked with each other via shared corners or shared edges into a new type of metal carbide structure.     

The crystal structure of (Fe_4Cr_4Ni)_9C_4

The Luobusha podiform chromite deposit occursin the Luobusha ophiolite exposed at the In-dus-Yarlung Zangbo suture[1]. The wall rocks of chro-mitites are harzburgites. Many mantle minerals werefound in chromitites and harzburgites. Among of them,there are 120 grains of diamond[2], many elementalmetals (mainly Au, Cu, Fe, Ni, Cr, Al, W, Zn, Pb, Sn,Os, Ir, Ru, C, etc.) and their intermetallic compounds.Besides diamond and graphite, the element C isalso combined with Fe, Cr, Ni, Ti and …     

Structure refinement of astrophyllite

The crystal structure of astrophyllite K₂Na(Fe, Mn, Mg,□)₇[Ti₂(Si₄O₁₂)₂|O₃](OH, F)₄ has been refined. The dimensions of the triclinic unit cell are: a = 0.5359(2) nm,b = 1.1614(4) nm, c = 1.1861(4) nm, α= 113.16(2)°, β= 103.04(2)°,γ= 94.56(2)°,V = 0.6495(5) nm³, Z= 1, space group P1, R=0.057 for 5308 reflections |F_o|>3σ|F_o|. According to structural and compositional differences the monoclinic astrophyllite K₂NaNa(Fe, Mn)₄Mg₂Ti2[Si₄O₁₂]₂(OH)₄(OH, F)₂ and astrophyllite should be considered as two different mineral species. Astrophyllite, monoclinic astrophyllite, bafertisite and lamprophyllite contain heteropolyhedral sheets which topologically are related with Si, O sheets of mica where one or several SiO₄ tetrahedra are replaced by TiO _n polyhedra. Therefore this heterophyllotitanosilicate series represents a kind of functional substitution in inorganic crystals.     

小球藻和铜绿微囊藻的高浓度Chl-a高光谱定量模型

本研究将地物高光谱遥感技术应用于室内实验,从而得到小球藻和铜绿微囊藻的高光谱特征.通过多种半经验方法,如单波段、波段比值和微分法,建立了两藻种最优的Chl-a高光谱定量模型,并与室外情况进行了对比.结果表明:小球藻的最优定量模型为Chl-a=174.6 1138292(R703) 2.3(109[(R703)]2(p＜0.01),相应的方法适宜性为:一阶微分法＞单波段法＞波段比值法;铜绿微囊藻的最优定量模型为Chl-a=5299164(R757)1.9773(p＜0.01),相应的方法适宜性为:单波段法＞波段比值法＞一阶微分法;从高光谱特征来看,小球藻在540 nm和700 nm附近存在明显的特征波峰,其位置随Chl-a浓度增大而向长波方向偏离,铜绿微囊藻在530 nm、660 nm和700 nm附近存在3个较强的特征波峰,在610 nm和680 nm附近存在明显的波谷;与以往室外研究不同的是铜绿微囊藻的反射率在400-500 nm之间的R值并不低,是因为没有非藻类颗粒物的影响,总吸收明显降低.     

适用于多种卫星数据的太湖水体漫衰减系数估算算法

下行漫衰减系数(K_d)是描述水下光场的重要参数,决定水体真光层深度,影响着浮游藻类初级生产力及其分布特征.基于2008—2013年太湖4次大规模野外试验数据,分析太湖水体漫衰减系数特征及其影响因素,建立适用于多种卫星数据且较高精度的太湖水体490 nm处下行漫衰减系数估算模型.结果表明:无机悬浮物是太湖水体漫衰减系数的主要影响因素;红绿波段比值(674 nm/555 nm)最适合于太湖K_d(490)估算,模型反演精度较高(N=72,R~2=0.72,RMSE=0.89 m~(-1),MAPE=21.58%);利用实测光谱数据,模拟得到MODIS/EOS、OLCI/Sentinel-3、GOCI/COMS和MSI/Sentinel-2等主要传感器波段的信号,构建适用于多种卫星传感器K_d(490)估算的红绿波段模型,建模精度较高(N=72,R~20.7,RMSE0.9 m~(-1),MAPE22.0%),且进行了验证(N=37,R~20.7,RMSE0.9 m~(-1),MAPE22.0%).     

湖北道观河水库细鳞鲴的个体生殖力

研究了道观河水库细鳞鲴的个体生殖力。细鳞鲴样本由2~ -5~ 龄鱼组成,其个体绝对生殖力与体长、体重和年龄分别呈直线、对数函数和幂函数相关,体长相对生殖力与体长、体重和年龄分别呈双曲线、对数函数和幂函数相关,体重相对生殖力与它们相关不显著。多元逐步回归分析其绝对生殖力和体长相对生殖力均与体重和成熟系数密切相关,而体重相对生殖力仅与成熟系数相关。     

Remote sensing algorithms for monitoring eutrophication in Ishizuchi storm water reservoir in Kochi Prefecture,Japan / Algorithmes de télédétection pour la surveillance de l'eutrophisation du réservoir d'orage d'Ishizuchi dans la Préfecture de Kochi,Japon

Abstract

Remote sensing technology is being conceived as a cost-effective tool to monitor the eutrophication problem in shallow inland lakes and reservoirs. The objective of this study is to develop hyperspectral remote sensing algorithms for chlorophyll and accessory pigments retrieval in the Ishizuchi Reservoir in Japan. The laboratory reflectance data of 207 field samples with manipulated pigment and suspended sediment concentrations were measured between 400 and 700 nm with a hyperspectral spectroradiometer, and integrated into 60 contiguous narrow bands. Band ratio and spectral absorption algorithms were developed by regression analysis with the observed chlorophyll measurements, and their estimation accuracies were evaluated using the RMSE of the estimates. Algorithms developed from spectral absorption ratios a(675 nm)/a*(675 nm) for five different classes of total pigment concentrations were more accurate than the algorithms based on band ratios R_rs (538 nm)/R_rs (488 nm), R_rs (581 nm)/R_rs (463 nm) and R_rs (674 nm)/R_rs (622 nm) as a result of removing the influence of overlapping absorption and scattering by water molecules and suspended solids from total absorption spectra.     

A three-dimensional,iterative mapping procedure for the implementation of an ionosphere-magnetosphere anisotropic Ohm’s law boundary condition in global magnetohydrodynamic simulations

The mathematical formulation of an iterative procedure for the numerical implementation of an ionosphere-magnetosphere (IM) anisotropic Ohm’s law boundary condition is presented. The procedure may be used in global magnetohydrodynamic (MHD) simulations of the magnetosphere. The basic form of the boundary condition is well known, but a well-defined, simple, explicit method for implementing it in an MHD code has not been presented previously. The boundary condition relates the ionospheric electric field to the magnetic field-aligned current density driven through the ionosphere by the magnetospheric convection electric field, which is orthogonal to the magnetic field B, and maps down into the ionosphere along equipotential magnetic field lines. The source of this electric field is the flow of the solar wind orthogonal to B. The electric field and current density in the ionosphere are connected through an anisotropic conductivity tensor which involves the Hall, Pedersen, and parallel conductivities. Only the height-integrated Hall and Pedersen conductivities (conductances) appear in the final form of the boundary condition, and are assumed to be known functions of position on the spherical surface R=R₁ representing the boundary between the ionosphere and magnetosphere. The implementation presented consists of an iterative mapping of the electrostatic potential , the gradient of which gives the electric field, and the field-aligned current density between the IM boundary at R=R₁ and the inner boundary of an MHD code which is taken to be at R₂>R₁. Given the field-aligned current density on R=R₂, as computed by the MHD simulation, it is mapped down to R=R₁ where it is used to compute by solving the equation that is the IM Ohm’s law boundary condition. Then is mapped out to R=R₂, where it is used to update the electric field and the component of velocity perpendicular to B. The updated electric field and perpendicular velocity serve as new boundary conditions for the MHD simulation which is then used to compute a new field-aligned current density. This process is iterated at each time step. The required Hall and Pedersen conductances may be determined by any method of choice, and may be specified anew at each time step. In this sense the coupling between the ionosphere and magnetosphere may be taken into account in a self-consistent manner.     

Kinematic stationary geodynamo models with separated toroidal and poloidal motions

Abstract

This paper is concerned with a three-dimensional spherical model of a stationary dynamo that consists of a convective layer with a simple poloidal flow of the S^2c ₂ kind between a rotating inner body core and solid outer shell. The rotation of the inner core and the outer shell means that there are regions of concentrated shear or differential rotation at the convective layer boundaries. The induction equation for the inside of the convective layer was solved numerically by the Bullard-Gellman method, the eigenvalue of the problem being the magnetic Reynolds number of the poloidal flow (R _M2) and it was assumed that the magnetic Reynolds number of the core (R _M1) and of the shell (R _M3) were prescribed parameters. Hence R _M2 was studied as a function of R _M1 and R _M3, along with the orientation of the rotation axis, the radial dependence of the poloidal velocity and the relative thickness of the layers for the three different situations, (i) the core alone rotating, (ii) the shell alone rotating and (iii) the core and the shell rotating together. In all three cases it was found that, at definite orientations of the rotation axis, there is a good convergence of both the eigenvalues and the eigenfunctions of the problem as the number of spherical harmonics used to represent the problem increases. For R _M1 =R _M3= 10³, corresponding to the westward drift velocity and the parameters of the Earth's core, the critical values of R _M2 are found to be three orders of magnitude lower than R _M1, R _M3 so that the poloidal flow velocity sufficient for maintaining the dynamo process is 10-20 m/yr. With only the core or the shell rotating, the velocity field generally differs little from the axially symmetric case. However, for R _M2 (or R _M3) lying in the range 10² to 10⁵, the self-excitation condition is found to be of the form R _M2˙R ^½ _M1=constant (or R _M2˙R^½ _M3=constant) and the solution does not possess the properties of the Braginsky near-axisymmetric dynamo. We should expect this, in particular, in the Braginsky limit R _M2˙R^?½; _M1=constant.

An analysis of known three-dimensional dynamo models indicates the importance of the absence of mirror symmetry planes for the efficient generation of magnetic fields.     

南方水库热分层消亡期混合层深度及缺氧区时空变化特征——以南宁市天雹水库为例

湖库热分层消亡引起的突发性水质恶化现象引起了广泛的关注,我国南方水库大多是暖单次混合型湖泊,每年混合一次,导致水库水质周期性下降,但目前对南方水库热分层消亡过程的高频监测研究较少。为探究我国南方水库热分层消亡期水体混合过程的时空变化规律及驱动因素,以广西南宁天雹水库为例,于冬季热分层消亡期(2019年11月—2020年2月)对水库多点位水体理化指标开展原位监测,并利用自建气象站获取气象水文数据。结果表明:(1)水库热分层消亡期间,过渡区水深较浅可在短期内达到完全混合状态且缺氧区同步消失;而湖泊区混合过程整体滞后于过渡区,混合层深度由6.85 m增加至13.65 m,缺氧区逐渐减小,缺氧指数(AI)由0.40减小至0.07,直至水体完全混合后缺氧区消失;水库过渡区较湖泊区提前约40 d达到完全混合状态。(2)气象因子是引起热分层结构变化的主要因素,气温(T)、辐射(R)与混合层深度(MLD)呈现显著负相关(R_T=0.927、R_R=0.925,P<0.01),风速(WS)与MLD呈现显著正相关(R_WS=0.728,P&...     

High-pressure phases of CaCO3: Crystal structure prediction and experiment

Post-aragonite phase of CaCO₃, experimentally known to be stable above 40 GPa [S. Ono, T. Kikegawa, Y. Ohishi, J. Tsuchiya, Post-aragonite phase transformation in CaCO₃ at 40 GPa, Am. Mineral. 90 (2005) 667–671], is believed to be a major carbon-containing mineral in the Earth's mantle. Crystal structure of this mineral phase could not be solved using experimental data or traditional theoretical simulation methods and remained a controversial issue. Using a combination of advanced ab initio simulation techniques and high-pressure experiment, we have been able to determine the crystal structure of CaCO₃ post-aragonite. Here, we performed simulations with the USPEX code [C.W. Glass, A.R. Oganov, and N. Hansen, (in preparation). USPEX: a universal structure prediction program], which is based on an evolutionary algorithm using ab initio free energy as the fitness function. This novel methodology for crystal structure prediction, which uses only the chemical composition as input, is described in detail. For CaCO₃, we identify a number of energetically competitive structures, the most stable of which closely matches the experimental powder diffraction pattern and, in agreement with experiment, becomes more stable than aragonite above 42 GPa. This structure belongs to a new structure type, which is also adopted by the high-pressure post-aragonite phases of SrCO₃ and BaCO₃. It has 2 formula units in the orthorhombic unit cell (space group Pmmn) and contains triangular CO₃²⁻ ions and Ca²⁺ ions in the 12-fold coordination. Above 137 GPa, a pyroxene-type structure (space group C222₁) with chains of CO₄⁴⁻ tetrahedra becomes more stable than post-aragonite. For MgCO₃, this structure becomes more stable than magnesite above 106 GPa and is a good candidate structure for MgCO₃ post-magnesite.     

Plasma pressure distribution in the equatorial plane of the Earth’s magnetosphere at geocentric distances of 6–10<Emphasis Type="Italic">R</Emphasis><Subscript><Emphasis Type="Italic">E</Emphasis></Subscript> according to the international THEMIS mission data

The structure of the averaged plasma pressure distribution in the plasma ring around the Earth at geocentric distances of ∼6–10R _E has been determined. The distribution function moments measured on the international THEMIS mission satellites have been used. The plasma pressure distribution in the equatorial plane at 15R _E > X_SM > −15R _E and 15R _E > Y_SM > −15R _E has been statistically studied. The radial dependence of the plasma pressure at the day-night and morning-evening meridians has been analyzed. It has been indicated that the plasma ring around the Earth has a structure, which is close to being azimuthally symmetric. The achieved results have been compared with the pressure distributions obtained previously. It has been indicated that in the overlapping regions, the achieved results agree with the previously obtained data within the pressure determination errors.     

Experimental evidence for flexibility of a building foundation supported by concrete friction piles

This article explores the possibility to measure deformations of building foundations from measurements of ambient noise and strong motion recordings. The case under study is a seven-storey hotel building in Van Nuys, California. It has been instrumented by strong motion accelerographs, and has recorded several earthquakes, including the 1971 San Fernando (M_L=6.6, R=22 km), 1987 Whittier–Narrows (M_L=5.9, R=41 km), 1992 Landers (M_L=7.5, R=186 km), 1992 Big Bear (M_L=6.5, R=149 km), and 1994 Northridge (M_L=6.4, R=1.5 km) earthquake and its aftershocks (20 March: M_L=5.2, R=1.2 km; 6 December, 1994: M_L=4.3, R=11 km). It suffered minor structural damage in 1971 earthquake and extensive damage in 1994. Two detailed ambient vibration tests were performed following the Northridge earthquake, one before and the other one after the 20 March aftershock. These included measurements at a grid of points on the ground floor and in the parking lot surrounding the building, presented and analyzed in this article. The analysis shows that the foundation system, consisting of grade beams on friction piles, does not act as a “rigid body” but deforms during the passage of microtremor and therefore earthquake waves. For this geometrically and by design essentially symmetric building, the center of stiffness of the foundation system appears to have large eccentricity (this is seen both from the microtremor measurements and from the earthquake recordings). This eccentricity may have contributed to strong coupling of transverse and torsional responses, and to larger than expected torsional response, contributing to damage during the 1994 Northridge, earthquake.     

The crystal structure and thermal history of orthopyroxene from lunar anorthosite 15415

A single crystal of untwinned orthopyroxene from lunar anorthosite sample 15415, with composition (Mg_1.14Fe_0.80Mn_0.02Ca_0.04)(Si_1.97Al_0.03)O₆, has a unit cell in space groupPbca witha = 18.310(15)Å,b = 8.904(10)Å,c = 5.214(7)Å, containing 2 formula units. A set of 742 counter-measured intensity data made with MoK_α radiation has been used to refine the crystal structure in isotropic thermal mode toR = 0.116. Anisotropic refinement led toR = 0.092, but thermal parameters are distorted by non-random errors resulting from poor crystal texture. The resulting structure is in close agreement with that obtained by Ghose [9] for a hypersthene from Greenland. A parameterq, which gives (Mg_qFe_1?q) for cation siteM(1) and (Mg_1.14?qFe_q?0.18Ca_0.04) for siteM(2), was included in the least-squares analysis, yieldingq = 0.90(1).This orthopyroxene has the high degree of cation order expected of pyroxenes subjected to Apollonian metamorphism at lower than 500–600°C. No evidence exists for a subsequent thermal event of sufficient intensity to disorder the pyroxene. On the basis of previous laboratory studies of argon-release patterns of lunar plagioclase and order-disorder kinetics of terrestrial pyroxenes, we attribute the reported isotopic age (3.9–4.1 AE) to cessation of metamorphism, perhaps caused by impact excavation.     

Shaking table tests of a reinforced concrete waffle–flat plate structure designed following modern codes: seismic performance and damage evaluation

Shaking table tests were conducted on a scaled reinforced concrete waffle–flat plate structure. It represented a conventional construction design under current building codes in the Mediterranean area. The test structure was subjected to a sequence of four seismic simulations of increasing magnitude. Each simulation was associated with a seismic hazard level characterized by the mean return period P_R. The test structure performed well for the simulations associated with P_R = 95, 475 and 975 years but collapsed under the maximum considered earthquake of P_R = 2475 years. Damage concentrated at column bases, where the maximum chord rotation reached 93% of the ultimate capacity, and at the transverse beams of the exterior plate‐to‐column connection that failed in torsion. It is shown that most (from 85% to 90%) of the energy input by the earthquake that contributes to damage is dissipated by the plate. The capacity curve of the tested structure estimated from the experimental base shear vs. top displacement relationship allowed us to compute the overstrength (1.4). It is close to the maximum established by European code EN 1998‐1 (1.5). Based on a detailed study of the test results, potential updates to current codes and design recommendations are suggested. Copyright © 2015 John Wiley & Sons, Ltd.     

Source parameters of bohai earthquake, July 18, 1969 and yongshan earthquake, May 11, 1974 determined by synthetic seismograms of teleseismic P waves

The source parameters of the Bohai Sea earthquake, July 18, 1969 and Yongshan, Yunnan earthquake, May 11, 1974 were determined by full — wave theory synthetic seismograms of teleseismic P waves. P+pP+sP wereform were calculated with WKBJ approximation and real integral paths. One — dimensional unilateral, finite propagation source was also considered. By trail — and — error in comparing the theoretical seismograms with the observational ones of WWSSN stations, the source parameters were obtained as follow: for Bohai earthquake, φ=195°, δ=85°, λ=65°,M _o=0.9×10¹⁹Nm,L=59.9km.V _R=3.5km/s, ∧ _R =160°; for Yongshan earthquake, φ=240°, δ=80°, ∧=150°,M _o=1.3×10¹⁸Nm,L=48.8km,V _R=3km/s, ∧ _R =−10°, where φ is strike, δ dip angle, λ slip angle,M _o seismic moment,L rupture length,V _R rupture propagation speed. As III type fractures the faulting propagated along the fault planes, and ∧ _R is the angle from the strike to the propagation direction. Yongshan earthquake showed complexity in its focal process, having four sub—ruptures during the first 60 seconds. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 1–8, 1991.     

Transitions to chaotic thermal convection in a rapidly rotating spherical fluid shell

Abstract

Numerical simulations of thermal convection in a rapidly rotating spherical fluid shell heated from below and within have been carried out with a nonlinear, three-dimensional, time-dependent pseudospectral code. The investigated phenomena include the sequence of transitions to chaos and the differential mean zonal rotation. At the fixed Taylor number T _a =10⁶ and Prandtl number Pr=1 and with increasing Rayleigh number R, convection undergoes a series of bifurcations from onset of steadily propagating motions SP at R=R _c = 13050, to a periodic state P, and thence to a quasi-periodic state QP and a non-periodic or chaotic state NP. Examples of SP, P, QP, and NP solutions are obtained at R = 1.3R _c , R = 1.7 R _c , R = 2R _c , and R = 5 R _c , respectively. In the SP state, convection rolls propagate at a constant longitudinal phase velocity that is slower than that obtained from the linear calculation at the onset of instability. The P state, characterized by a single frequency and its harmonics, has a two-layer cellular structure in radius. Convection rolls near the upper and lower surfaces of the spherical shell both propagate in a prograde sense with respect to the rotation of the reference frame. The outer convection rolls propagate faster than those near the inner shell. The physical mechanism responsible for the time-periodic oscillations is the differential shear of the convection cells due to the mean zonal flow. Meridional transport of zonal momentum by the convection cells in turn supports the mean zonal differential rotation. In the QP state, the longitudinal wave number m of the convection pattern oscillates among m = 3,4,5, and 6; the convection pattern near the outer shell has larger m than that near the inner shell. Radial motions are very weak in the polar regions. The convection pattern also shifts in m for the NP state at R = 5R _c , whose power spectrum is characterized by broadened peaks and broadband background noise. The convection pattern near the outer shell propagates prograde, while the pattern near the inner shell propagates retrograde with respect to the basic rotation. Convection cells exist in polar regions. There is a large variation in the vigor of individual convection cells. An example of a more vigorously convecting chaotic state is obtained at R = 50R _c . At this Rayleigh number some of the convection rolls have axes perpendicular to the axis of the basic rotation, indicating a partial relaxation of the rotational constraint. There are strong convective motions in the polar regions. The longitudinally averaged mean zonal flow has an equatorial superrotation and a high latitude subrotation for all cases except R = 50R _c , at this highest Rayleigh number, the mean zonal flow pattern is completely reversed, opposite to the solar differential rotation pattern.     

一种基于MODIS影像可业务化运行的巢湖水体叶绿素a估算算法

现有水色卫星主要是针对大洋清洁水体设计,内陆浑浊水体多数波段经常饱和;而发展可以业务化运行的内陆水体叶绿素a算法,为生产实践服务,一直是水色遥感的重点和难点之一.利用2013年巢湖星地同步数据(N=55),通过经验正交函数(empirical orthogonal function,EOF)分析方法,选用MODIS唯一不饱和的4个波段(469、555、645、859 nm)数据进行分解,然后回归建模;并使用第三方独立的巢湖实测数据(N=40)进行验证(R2=0.63,URMSE=85.46%).结果表明:该算法用于MODIS影像上,空间分布合理,季节差异明显,且在高悬浮物水体、不同气溶胶条件下均有很好的抗扰动性.实践证明EOF算法可以应用于业务化运行的内陆水体叶绿素a浓度估算,并对其他水色参数反演具有一定的借鉴意义.     

Establishment of a sediment transport capacity equation on loessal slope via experimental investigation

Areliable sediment transport capacity function provides response against challenges of soil erosion prediction on the Loess Plateau of China. The popular sediment transport capacity functions are questionable on loess slopes because the experimental conditions from which they were derived, like bed materials, gentle slopes,and surface roughness, are different from soil erosion processes on the loess slopes. Due to the foregoing uncertainty, a suitable sediment transport capacity function was dev...