Correcting systematic effects in a large set of photometric light curves

We suggest a new algorithm to remove systematic effects in a large set of light curves obtained by a photometric survey. The algorithm can remove systematic effects, such as those associated with atmospheric extinction, detector efficiency, or point spread function changes over the detector. The algorithm works without any prior knowledge of the effects, as long as they linearly appear in many stars of the sample. The approach, which was originally developed to remove atmospheric extinction effects, is based on a lower rank approximation of matrices, an approach which has already been suggested and used in chemometrics, for example. The proposed algorithm is especially useful in cases where the uncertainties of the measurements are unequal. For equal uncertainties, the algorithm reduces to the Principal Component Analysis (PCA) algorithm. We present a simulation to demonstrate the effectiveness of the proposed algorithm and we point out its potential, in the search for transit candidates in particular.     

An HLLC Riemann solver for relativistic flows – II. Magnetohydrodynamics

An approximate Riemann solver for the equations of relativistic magnetohydrodynamics (RMHD) is derived. The Harten–Lax–van Leer contact wave (HLLC) solver, originally developed by Toro, Spruce and Spears, generalizes the algorithm described in a previous paper to the case where magnetic fields are present. The solution to the Riemann problem is approximated by two constant states bounded by two fast shocks and separated by a tangential wave. The scheme is Jacobian-free, in the sense that it avoids the expensive characteristic decomposition of the RMHD equations and it improves over the HLL scheme by restoring the missing contact wave.
Multidimensional integration proceeds via the single step, corner transport upwind (CTU) method of Colella, combined with the constrained transport (CT) algorithm to preserve divergence-free magnetic fields. The resulting numerical scheme is simple to implement, efficient and suitable for a general equation of state. The robustness of the new algorithm is validated against one- and two-dimensional numerical test problems.     

关于水资源承载力理论与方法的研究

本文首先明确了水资源承载力的基本概念 ,指出可持续发展理论、水 -生态 -社会经济复合系统理论、二元模式下的水文循环机制和过程是水资源承载力研究的理论基础 ,对水资源承载力研究的基本内容、评价指标体系的建立和评价方法进行了较全面的概括和总结。     

闽西南古田-小陶花岗质杂岩体的形成时代和成因

古田-小陶花岗质杂岩体是由小陶、古田、华家、蛟洋和姑田等单元组成.小陶、古田和华家单元为含黑云母二长花岗岩,蛟洋和姑田单元是含角闪石黑云母石英二长岩.小陶和古田花岗岩以高SiO2（＞71.77%）、K2O（＞4.67%）和低CaO（＜1.71%）为特征,属于弱过铝-强过铝花岗岩;它们含高的Rb、Nb、Ta、Th、U和低的Sr、Ba含量;稀土（REE）含量和轻重稀土分异（（La/Yb）N）变化大,分别变化于143×10-6～548×10-6,和3.3～26.5,Eu亏损中等到强.姑田和蛟洋单元都以低SiO2,（62.88%～66.38%）和高Al2O3（15.44%～16.51%）为特征.姑田单元的REE总量低而蛟洋单元的REE总量高,但两者的Eu负异常都不明显.华家单元花岗岩的地球化学特征介于上述两组花岗岩之间.LA-ICPMS锆石U-Pb定年显示小陶花岗岩形成于222±3Ma,而古田花岗岩形成于161±1Ma.小陶、古田和姑田花岗岩具有相似的87Sr/86Sr初始比值（0.7089～0.7099）和εNd（t）值（-11.1～-11.5）,而且古田和小陶花岗岩的锆石Hf同位素组成也相似,表明它们起源于相同的源区.地球化学和Sr-Nd-Hf同位素指示它们的源区物质是具有1.73～1.5Ga平均地壳存留年龄的变质沉积岩.而位于杂岩体南部的华家花岗岩显示明显高的87Sr/86Sr初始比值（0.7310）和更低的εNd（t）值（-16.4）,指示它的母岩浆起源于古元古代壳源物质,暗示杂岩体南北基底组成可能不同.古田花岗岩形成于燕山早期的板内伸展构造背景,而小陶花岗岩形成于印支晚期,地球化学特征以及同时代岩石组合以及构造应力场分析表明其最可能形成于后造山构造背景.     

渭河中游全新世黄土剖面光释光测年及记录的古洪水事件

通过对渭河流域进行广泛的考察,在中游咸阳附近一处阶地发现保存完好的全新世黄土-古土壤剖面里夹有古洪水滞流沉积层,对其进行了年代学和沉积学研究.对于采集的全新世地层样品,进行粒度、磁化率测量分析,证明所夹沉积物是典型的古洪水滞流沉积物,该层记录了古洪水事件发生的气候水文信息.应用红外后蓝光的SAR光释光测年技术,获得该剖面9个OSL年龄值,确定渭河在距今3.2～2.8 ka之间为一个洪水多发时期,并建立了渭河古洪水事件的年代序列.这一结果揭示了古洪水的频发与全新世中期向晚期过渡的转折时期,气候由温湿向干旱化发展,大气系统失稳,气候变化剧烈,降水量异常变化,变率增大,是导致渭河流域特大古洪水多发的主要原因.     

青藏高原北羌塘他念他翁山吉塘岩群酉西岩组时代的确定

关于位于青藏高原北羌塘地区的吉塘岩群是否属于前寒武纪变质岩系历来存在较大争议,对其形成年龄一直缺少有效同位素年代学限定.通过高精度的LA-ICP-MS(激光剥蚀等离子体质潜仪)锆石微区原位U-Pb同位素测年,获得他念他翁山一带吉塘岩群酉西岩组石英绿泥片岩(原岩为中性火山岩)的原岩形成年龄为965±55Ma,绿片岩(原岩为中基性火山岩)的原岩形成年龄为1048.2±3.3Ma,表明吉塘岩群酉西岩组原岩形成时代为蓟县纪末一南华纪初.吉塘岩群酉西岩组可能代表Rodinia超大陆裂解早期伸展期间火山沉积作用的地质记录.     

闽东大湾钼铍矿的辉钼矿Re-Os同位素年龄——兼论福建省钼矿时空分布及构造背景

闽东大湾钼铍矿是海西(海峡西岸)经济区新近勘查取得重要进展的燕山期钼矿之一.铍钼矿体主要赋存于上侏罗统南园组火山岩中,隐伏的黑云母花岗岩体为主要的成矿母岩,属中-高温热液充填-交代型钼铍矿床.辉钼矿的Re-Os同位素测年结果显示,等时线年龄为92.2±1.3 Ma,表明钼矿成矿时代为晚白垩世.结合近年来福建省钼矿成矿年代学资料,发现自闽西(武夷成矿带)到闽东(沿海火山岩成矿带),钼矿成矿时代从165～102Ma过渡到112 ～90 Ma,显示逐渐变新的趋势.综合分析认为,165～ 131Ma的钼矿可能形成于岩石圈伸展背景,110 ～ 90 Ma的钼矿可能形成于古太平洋板块俯冲背景之下.     

A survey of numerical methods for stochastic differential equations

The development of numerical methods for stochastic differential equations has intensified over the past decade. The earliest methods were usually heuristic adaptations of deterministic methods, but were found to have limited accuracy regardless of the order of the original scheme. A stochastic counterpart of the Taylor formula now provides a framework for the systematic investigation of numerical methods for stochastic differential equations. It suggests numerical schemes, which involve multiple stochastic integrals, of higher order of convergence. We shall survey the literature on these and on the earlier schemes in this paper. Our discussion will focus on diffusion processes, but we shall also indicate the extensions needed to handle processes with jump components. In particular, we shall classify the schemes according to strong or weak convergence criteria, depending on whether the approximation of the sample paths or of the probability distribution is of main interest.     

Late quaternary glacial stades in the Cordillera Central,Colombia, based on glacial geomorphology,tephra–soil stratigraphy,palynology, and radiocarbon dating

Using data from glacial geomorphology, tephra–soil stratigraphy and mineralogy, palynology, and radiocarbon dating, a sequence of glacial and bioclimatic stades and interstades has been identified for the last ca. 50000 yr in the Ruiz-Tolima massif, Cordillera Central, Colombia. Six Pleistocene cold stades separated by warmer interstades occurred: before 48000, between 48000 and 33000, between 28000 and 21000, from ≥16000 to ca. 14000, ca. 13000–12400, and ca. 11000–10000 yr BP. Although these radiocarbon ages are minimum-limiting ages obtained from tephra layers on top of tills, the tills are not significantly older because most are bracketed by dated tephra sets in measured stratigraphic sections. Two minor moraine stages likely reflect glacier standstill during cold intervals ca. 7400 yr BP and slightly earlier. Finally, glaciers readvanced between the seventeenth and nineteenth centuries. In contrast to the ice-clad volcanoes of the massif, ca. 34 km² in area above an altitude of ca. 4800 m, the ice cover expanded to 1200 km² during the Last Glacial Maximum (LGM) and was still 800 km² during Late-glacial time (LGT). Glacier reconstructions based on the moraines suggest depression of the equilibrium line altitude (ELA) by ca. 1100 m during the LGM and 500–600 m during LGT relative to the modern ELA, which lies at ca. 5100 m in the Cordillera Central. Glaciers in this region apparently reached their greatest extent when the climate was cold and wet, e.g. during stades corresponding to Oxygen Isotope Stage 3; glaciers were still expanding during the LGM ca. 28000–21000 yr BP, but they shrank considerably after 21000 yr BP because of greatly reduced precipitation. © 1997 John Wiley & Sons, Ltd.     

基于自定义RDD的海量遥感图像并行镶嵌方法

图像镶嵌是遥感图像处理中的重要内容,在跨区域遥感图像分析中发挥重要作用。为了解决传统遥感图像并行算法中存在的计算节点利用率低、频繁数据I/O等问题,本文根据Spark分布式内存计算框架,充分利用Spark利于迭代数据处理的优势,提出了一种基于Spark自定义RDD（弹性分布式数据集）的并行镶嵌方法。该方法首先在集群的多个节点上通过相位相关法执行图像重叠区域估计操作,从而提高了图像重叠区域估计的多节点并行计算;然后,通过重写Spark中RDD的compute和getPartitions方法,自定义针对遥感图像处理的RDD,并将图像镶嵌中的重叠区域估计、图像配准和图像融合3个关键步骤作为自定义RDD的Transformation类型的操作算子;最后,通过隐式转换创建自定义RDD,并调用自定义RDD的操作算子实现图像镶嵌的并行处理。实验结果表明,与传统基于MPI的并行镶嵌算法相比,该方法在保证图像镶嵌效果的基础上,能够有效提高大数据量的图像镶嵌效率。