An iterative scheme for use in MT 2D forward modelling — to speed up production of the Jacobian and to improve accuracy of the inverse matrix

An important stage in two-dimensional magnetotelluric modelling is the calculation of the Earth's response functions for an assumed conductivity model and the calculation of the associated Jacobian relating those response functions to the model parameters. The efficiency of the calculation of the Jacobian will affect the efficiency of the inversion modelling. Rodi (1976) produced all the Jacobian elements by inverting a single matrix and using an approximate first-order algorithm. Since only one inverse matrix required calculation the procedure speeded up the inversion. An iterative scheme to improve the approximation to the Jacobian information is presented in this paper. While this scheme takes a little longer than Rodi's algorithm, it enables a more accurate determination of the Jacobian information. It is found that the Jacobian elements can be produced in 10% of the time required to calculate an inverse matrix or to calculate a 2D starting model. A modification of the algorithm can further be used to improve the accuracy of the original inverse matrix calculated in a 2D finite difference program and hence the solution this program produces. The convergence of the iteration scheme is found to be related both to the originally calculated inverse matrix and to the change in the newly formed matrix arising from perturbation of the model parameter. A ridge regression inverse algorithm is used in conjunction with the iterative scheme for forward modelling described in this paper to produce a 2D conductivity section from field data.     

The CCD photometric study of the newly identified RS CVn binary star DV Piscium

This paper presents new CCD BVRI light curves of the newly discovered RS CVn eclipsing binary star DV Psc. From the asymmetric light curves outside the eclipse, we find there are two depressions in the phase ranges 0.3–0.45 and 0.6–0.9, respectively. By analyzing the light curves using the Wilson–Devinney program, the four photometric solutions of the system are obtained and the starspot parameters are also derived. It turns out that the case of two spots being on the primary is most successful in reproducing the light curve distortion of DV Psc. Moreover, analysis the longitudes of spots suggests that there are two active longitude belts (one is about 90°, the other is about 270°). At the same time, on 22 November 2008, the first flare-like event was detected on DV Psc at phase 0.35 whose duration was about 13.5 min. These findings reveal that DV Psc has a high degree of magnetic activity.     

Using a New Sky Brightness Monitor to Observe the Annular Solar Eclipse on 15 January 2010

For the future development of Chinese Giant Solar Telescope (CGST) in Western China, a new sky brightness monitor (SBM) has been produced for the site survey for CGST. To critically examine the performance and sensitivity of SBM, we used it in the observation of the annular solar eclipse in Dali City, Yunnan, on 15 January 2010. The observation met good weather conditions with an almost clear sky during the eclipse. The SBM measurement translates into the solar illuminance changes at a level of 2.4×10^?4 I?s^?1 during the eclipse. The time of the minimal sky brightness in the field of view (FOV) is found consistent with the time of maximum eclipse. Two local sky regions in the FOV are chosen to make a time series of the calibrated skylight profiles. The evolution of the sky brightness thus calibrated also shows good consistency with the eclipse, particularly between the second and the third contacts. The minimal sky brightness in each local sky region took place within half a minute from the corresponding predicted contact time. Such small time delays were mainly caused by occasional cirri. The minimal sky brightness measured during the eclipse is a few millionths of I _?? with standard deviation of 0.11 millionths of I _??. The observation supports that the single-scattering process (optically thin conditions) is the main contributor to the atmospheric scattering. We have demonstrated that many important aerosol optical parameters can be deduced from our data. We conclude that the new SBM is a sensitive sky photometer that can be used for our CGST and coronagraph site surveys.     

A conscientious investigation into the effect of Re and mass-ratio on two-degree-of-freedom vortex-induced vibrations

Through a serious application of an overlapping mesh, vortex-induced vibration (VIV) of low-mass-rationale cylinders was computed in the range of 1.0×10 3 相似文献   

苡茨河泥石流地质灾害及防治

据统计,1992~2013年苡茨河曾多次发生泥石流地质灾害,并造成村民生命财产巨大损失。2013年7月24日再次爆发泥石流,破坏农田71280m2(108亩)。因此在对苡茨河泥石流形成条件、活动特征分析基础上,分析其泥石流的形成机制及其造成巨大危害的原因,预测其发展趋势,在已有治理措施的基础上,提出两种防治对策进行对比,得出一个较为合理的防治方案。     

江西萍乡杨家湾2号洞哺乳动物化石及其时代探讨

2015年11月和2016年10月先后对杨家湾2号洞化石点进行发掘。目前已出土哺乳动物化石1000多件,主要为单个牙齿化石,骨骼较少且较为破碎。初步鉴定出31个哺乳动物属种,其中偶蹄类数量最为丰富,野猪所占比例最大;灵长类和食肉类较为丰富;啮齿类、长鼻类及奇蹄类较少;尚未发现食虫类和翼手类。大部分属种为现生类群,少量灭绝属种如巴氏大熊猫、东方剑齿象及华南巨貘等。从动物群组成来看,该动物群与杨家湾1号洞、湖南道县福岩洞等动物群较为相似,属于我国南方典型的晚更新世大熊猫-剑齿象动物群,但是其所包含的动物属种相对后两个化石点都较少,尤其是小哺乳动物属种单调。杨家湾2号洞的堆积物颗粒比1号洞的较粗,且含较多角砾,缺少红色粘土;据初步判断,该洞的堆积物是由洞口进入。动物群组合指示该地区晚更新世气候较为温暖湿润。

     

Impact of climate change on runoff and uncertainty analysis

Natural Hazards - It is necessary to analyze the future runoff changes using a more realistic climate classification scheme. This paper investigates the climate changes and runoff variation by...     

二维有限元大地电磁正演模拟在Matlab上的实现

Matlab高级语言开发平台对矩阵的计算简便快速,并且绘图输出文件方便。以这些特点为契机,在Matlab平台上,讨论利用有限单元法实现二维大地电磁正演的问题,推导出相应的公式,计算地表或海底视电阻率及相位,并与解析解以及主流方法结果对比,证实了本文方法的正确性。     

顾及有效角动量与IGS超快解数据的极移预报方法

极移高精度预报对卫星实时定轨、深空探测器导航等应用至关重要.本文提出了一种联合有效角动量(Effective Angular Momentum, EAM)与国际全球卫星导航系统服务组织(International GNSS Service, IGS)提供的超快解数据进行极移预报的方法.该方法基于IGS超快解数据得到的极移第1天预报值, 对引入EAM得到的极移预报结果进行校正, 获得联合预报值.首先, 基于LS(Least Squares)+AR(Auto-Regressive)模型实现了引入EAM的极移预报, 相对国际地球自转与参考系统服务组织(International Earth Rotation and Reference Systems Service, IERS)提供的公报A数据, 在超短期(第1~10天)预报跨度可以得到更高精度的极移预报结果, 其中大气和海洋角动量发挥了主要作用.其次, 鉴于IGS超快解数据精度高、更新快的特点, 以IGS超快解为基础数据, 基于LS+AR模型可以得到极移第1天预报值, 其精度显著优于IERS公报A的极移第1天预报值.最后, 利用第1天预报值对顾及有效角动量的预报结果进行校正获得了联合预报值, 进一步提高了超短期极移预报精度(尤其是第1~5天).2020年7月24日—2022年1月30日间的联合预报结果表明: 第1~20天的预报值总体优于IERS公报A的预报值.其中, 第1~10天的预报精度显著提升, 在预报第1天, X、Y方向预报值相对公报A预报值的精度提升分别可达39.5%~62.3%和24.5%~51.9%;在预报第10天, 相对公报A预报值的精度提升分别可达28.0%~28.9%和21.9%~23.4%.

     

Validity of the Bouchet’s complementary relationship at 102 observatories across China

Bouchet in 1963 hypothesized that for large homogeneous land surface with minimum advection of heat and moisture, there exists a 1:1 complementary relationship of potential and actual evaporation coupled through land-atmosphere feedbacks. The complementary relationship has been widely used to estimate regional actual evaporation and explain the pan evaporation paradox. We examine the standardized potential evaporation (potential evaporation divided by wet environment evaporation) at 102 observatories at different elevations across China. Generally, the relationship is appropriate at the low elevations (<1000 m). With the increase of elevation, vapor transfer power becomes much less than radiation energy budget because of lower vapor pressure deficit and stronger global solar radiation. As a result, at the high elevations (over 1000 m), the excess energy resulted by limited moisture availability is not enough to be converted into drying power of the air. This result suggests that the complementary relationship is asymmetric at the high elevations. Supported by the Presidential Special Award Foundation, the Chinese Academy of Sciences (Grant No. O7R70020SD) and the National Key Technology R & D Program (Grant No. 2006BAC08B0408)