Analytical investigation of the exact groundwater divide between rivers beyond the Dupuit–Forchheimer approximation

The groundwater divide is a key feature of river basins and significantly influenced by subsurface hydrological processes. For an unconfined aquifer between two parallel rivers or ditches, it has long been defined as the top of the water table based on the Dupuit–Forchheimer approximation. However, the exact groundwater divide is subject to the interface between two local flow systems transporting groundwater to rivers from the infiltration recharge. This study contributes a new analytical model for two-dimensional groundwater flow between rivers of different water levels. The flownet is delineated in the model to identify groundwater flow systems and the exact groundwater divide. Formulas with two dimensionless parameters are derived to determine the distributed hydraulic head, the top of the water table and the groundwater divide. The locations of the groundwater divide and the top of the water table are not the same. The distance between them in horizontal can reach up to 8.9% of the distance between rivers. Numerical verifications indicate that simplifications in the analytical model do not significantly cause misestimates in the location of the groundwater divide. In contrast, the Dupuit–Forchheimer approximation yields an incorrect water table shape. The new analytical model is applied to investigate groundwater divides in the Loess Plateau, China, with a Monte Carlo simulation process taking into account the uncertainties in the parameters.     

Petrochronology of oxidized granulites from southern Peru

Sapphirine–quartz granulites from the Cocachacra region of the Arequipa Massif in southern Peru record early Neoproterozoic ultrahigh‐temperature metamorphism. Phase equilibrium modelling and zircon petrochronology are used to quantify timing and pressure–temperature (P–T) conditions of metamorphism. Modelling of three magnetite‐bearing sapphirine–quartz samples indicates peak temperatures of >950°C at ~0.7 GPa and a clockwise P–T evolution. Elevated concentrations of Al in orthopyroxene are also consistent with ultrahigh‐temperature conditions. Neoblastic zircon records ages of c. 1.0–0.9 Ga that are interpreted to record protracted ultrahigh‐temperature metamorphism. Th/U ratios of zircon of up to 100 reflect U‐depleted whole‐rock compositions. Concentrations of heavy rare earth elements in zircon do not show systematic trends with U–Pb age but do correlate with variable whole‐rock compositions. Very large positive Ce anomalies in zircon from two samples probably relate to strongly oxidizing conditions during neoblastic zircon crystallization. Low concentrations of Ti‐in‐zircon (<10 ppm) are interpreted to result from reduced titania activities due to the strongly oxidized nature of the granulites and the sequestration of titanium‐rich minerals away from the reaction volume. Whole‐rock compositions and oxidation state have a strong influence on the trace element composition of metamorphic zircon, which has implications for interpreting the geological significance of ages retrieved from zircon in oxidized metamorphic rocks.     

Determination of the regional deformation rates of Shanghai and Kashima VLBI stations based on ITRF97

The vertical deformation rates (VDRs) and horizontal deformation rates (HDRs) of Shanghai VLBI station in China and Kashima and Kashima34 VLBI stations in Japan were re-analysed using the baseline length change rates from Shanghai to 13 global VLBI stations, and from Kashima to 27 stations and from Kashima34 to 12 stations, based on the NASA VLBI global solution glb1123 (Ma, 1999). The velocity vectors of the global VLBI stations were referred to the ITRF97 reference frame, and the Eulerian vectors of different models of plate motion were used for comparative solutions. The VDR of Shanghai station is estimated to be −1.91±0.56 mm/yr, and those of Kashima and Kashima34 stations, −3.72±0.74mm/yr and −8.81±0.84mm/yr, respectively. The difference between the last two was verified by further analysis. Similar estimates were also made for the Kokee, Kauai and MK_VLBA VLBI stations in mid-Pacific.     

合并会计报表的电算编制方法探讨

利用用友软件对如何合并会计报表进行了初步研究，并以案例的形式阐明了采用用友软件实现合并会计报表的步骤，提出了合并会计报表的编制方法。     

青海可可西里地区综合科学考察初报

本区至今还是综合科学考察的空白区。为查明基本情况,1989年夏以乘汽车路线考察方式穿行该区而作了预查,搜集了当地地质、地理、生物等基本资料。经考察发现,区内野生动物资源丰富,上新世以来火山活动普遍,湖泊众多,矿产资源可能有良好前景;所在区又为自然地域和地质构造的重要界线,对青藏高原研究具有重大科学价值。     

The formation and stability of the configuration of the planetary system HD 126611

Recent Doppler velocity measurements have revealed the existence of two planets orbiting the star HD 12661 on medium-eccentricity orbits. The inner planet has a period of 263.6 d and a mass of 2.3 Jupiter masses, and the outer planet has a period of 1444.5d and a mass of 1.57 Jupiter masses. The stability of this system requires the two planets to be in a state of mean-motion orbit resonances. By numerical method we have studied the orbit migration and stability of the system in its early ages under the action of the proto-stellar disk, and calculated the probabilities of the planets being captured into the mean -motion resonances during their migrations. It is found that at present the two planets are probably situated at the edge of the 11:2 mean-motion resonance and are in chaotic motions. This result may be helpful to clarify the arguments on the present configuration. Besides, it is indicated that very probably, after the formation of the system, the gaseous disk has almost disappeared before the planets migrated to the present configuration.     

对"数字地球"的几点认识

自“数字地球”的概念被提出以来，已引起全球的广泛关注和积极响应。然而，数字地球在目前还仅仅是一个概念。为了逐步建立“数字地球”，需要就“数字地球”的理论、技术以及实现进行深入地探讨。从“数字地球”的产生背景、基本概念、研究内容、构成框架、关键技术、实施难点和发展战略等方面谈了对“数字地球”的一些认识。     

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.