Algorithm Studies on How to Obtain a Conditional Nonlinear Optimal Perturbation (CNOP)

The conditional nonlinear optimal perturbation (CNOP), which is a nonlinear generalization of the linear singular vector (LSV), is applied in important problems of atmospheric and oceanic sciences, including ENSO predictability, targeted observations, and ensemble forecast. In this study, we investigate the computational cost of obtaining the CNOP by several methods. Differences and similarities, in terms of the computational error and cost in obtaining the CNOP, are compared among the sequential quadratic programming (SQP) algorithm, the limited memory Broyden-Fletcher-Goldfarb-Shanno (L-BFGS) algorithm, and the spectral projected gradients (SPG2) algorithm. A theoretical grassland ecosystem model and the classical Lorenz model are used as examples. Numerical results demonstrate that the computational error is acceptable with all three algorithms. The computational cost to obtain the CNOP is reduced by using the SQP algorithm. The experimental results also reveal that the L-BFGS algorithm is the most effective algorithm among the three optimization algorithms for obtaining the CNOP. The numerical results suggest a new approach and algorithm for obtaining the CNOP for a large-scale optimization problem.     

L波段探测系统使用心得

呼和浩特观象台自2008年9月1日使用L波段雷达─GTS1型数字式探空仪气象探测系统,综合探测已2个多月了,在工作期间积累了一些经验,与各位同行交流。1     

准线形对流系统中强度—速度负相关性的特征研究

准线形对流系统(如飑线、弓形回波)是常见的强对流性系统,伴随它出现的一般是大风、旋风甚至冰雹等严重的灾害性天气现象.弓形回波有的是由一个单体产生的,也有嵌套在飑线或其他线形系统中,作为整个线形回波的一个部分.弓形回波能在很长的带上产生灾害性的直线风.在我们的研究中发现,准线形对流系统中,反射率因子(R)和径向速度(V)之间存在一定的负相关性,对发生在中国东部地区(24.5°-34.29°N,116.01°-119.52°E)的16个个例进行统计分析,发现在这些系统中都存在这种R-V负相关的特征,负相关的程度随高度有规律性的变化.通过对2006年6月29日个例的风场分析证实这种特征确实存在.对流系统与雷达射线方向接近垂直时,R-V负相关特征明显,这是因为该方向上,观测到的径向速度近乎等于实际风速,数值最大;而因垂直运动造成的水平风速减小并不因探测位置的不同而有很大变化.文中通过分析,推断R-V的这种负相关性与准线形系统中对流的强弱、中层径向辐合的强弱以及后向入流都有很大的关系.     

重庆2008年7月21日强对流天气成因及其特征

对2008年7月21日重庆市一次强对流天气过程进行诊断分析,通过分析天气背景,雷达回波,并应用四维变分同化方法反演风场,研究了强对流天气的发展机制.结果表明:高层冷空气入侵和中低层低涡系统强盛的暖湿气流形成的辐合切变是发生强对流的主要原因.雷达观测具有明显的弱回波区和"逆风区";反演水平风场表现为低层强回波前部有一条辐合上升带,到中层辐合上升带与强回波区重合,高层则对应辐散区;从风场的垂直剖面来看,强对流后部为入流区,到达中部后,形成强烈的上升气流,低层辐合带和强烈的上升气流形成强回波.在前部弱回波区处存在强降水下落形成的出流,该气流与暖湿气流辐合造成了浅薄的出流边界.     

On the origin of the compact H ii regions around massive stellar clusters

The existence of compact H?ii regions around massive stellar clusters with ages exceeding several Myr challenges our understanding of the physical processes occurring inside such clusters, and their impact on the interstellar medium of the host galaxy. Here, we summarize our recent results dealing with the hydrodynamics of matter ejected by massive stars inside the cluster and show that compact H?ii regions found around some massive clusters may indicate that these are evolving in a bimodal hydrodynamic regime. The latter is characterized by the accumulation of the injected matter in the central, thermally unstable zone, and by the ejection of mass supplied by massive stars in the outer regions of the cluster.     

Model prediction capacity of ephemeral gully evolution in conservation tillage systems

Ephemeral gully (EG) erosion has an important impact on agricultural soil losses and increases field surface hydrology connectivity and transport of pollutants to nearby water bodies. Watershed models including an EG component are scarce and not yet properly evaluated. The objective of this study is to evaluate the capacity of one such tool, AnnAGNPS, to simulate the evolution of two EG formed in a conservation tillage system. The dataset for model testing included runoff measurements and EG morphological characteristics during 3 years. Model evaluation focused on EG evolution of volume, width, and length model outputs, and included calibration and testing phases and a global sensitivity analysis (GSA). While the model did not fully reproduce width and length, the model efficiency to simulate EG volume was satisfactory for both calibration and testing phases, supporting the watershed management objectives of the model. GSA revealed that the most sensitive factors were EG depth, critical shear stress, headcut detachment exponent coefficient b, and headcut detachment leading coefficient a. For EG outputs the model was additive, showing low sensitivity to interactions between the inputs. Prediction of EG spatial evolution on conservation tillage systems requires improved development of gully erosion components, since many of the processes were developed originally for traditional tillage practices or larger channel systems. Our results identify the need for future research when EG form within conservation tillage systems, in particular to study gully headcut, soil erodibility, and width functions specific to these practices.     

甘肃阿尔金南缘白石头沟石墨矿地质特征、成矿条件及找矿标志

甘肃白石头沟石墨矿位于阿尔金南缘敦煌地块,赋存于前寒武纪敦煌岩群大理岩和石英片岩之中,目前圈定矿体33条,主要矿体2条。作为研究白石头沟石墨矿理想地层的敦煌岩群,是一套以变质碎屑沉积岩和变质基性岩为主的变质岩浆杂岩。本文从矿床角度出发,立足于区域地质勘查资料,结合前人研究成果,在对白石头沟石墨矿成矿地质背景、矿床及矿体特征、组合样化学分析、矿石矿物赋存状态及特征研究的基础上,对矿床的成因进行初步探讨并归纳区域找矿标志。研究表明,白石头沟石墨矿可分为大理岩型和石英片岩型两种类型,是古元古代经历区域动力热液变质作用而形成的沉积变质型矿床,成矿物质主要来源于变质碎屑沉积岩类的原岩（海相陆源碎屑岩和碳酸盐岩）本身,成矿环境为古老活动大陆边缘。找矿标志具有特定层位、岩性、矿化及物探特征。     

High resolution facies architecture and digital outcrop modeling of the Sandakan formation sandstone reservoir,Borneo:Implications for reservoir characterization and flow simulation

Advances in photogrammetry have eased the acquisition of high-resolution digital information from outcrops, enabling faster, non-destructive data capturing and improved reservoir modeling. Geocellular models for flow dynamics with in the virtual outcrop in siliciclastic deposits at different sets of sandstone facies architecture remain, however, a challenge. Digital maps of bedding, lithological contrast, spatial-temporal variations of bedding and permeability characteristics make it more easy to understand flow tortuosity in a particular architecture. An ability to precisely model these properties can improve reservoir characterization and flow modeling at different scales. Here we demonstrate the construction of realistic 2 D sandstone facies based models for a pragmatic simulation of flow dynamics using a combination of digital point clouds dataset acquired from LiDAR and field investigation of the Sandakan Formation, Sabah, Borneo.Additionally, we present methods for enhancing the accuracy of outcrop digital datasets for producing high resolution flow simulation. A well-exposed outcrop from the Sandakan Formation, Sabah, northwest Borneo having a lateral extent of 750 m was chosen in order to implement our research approach. Sandstone facies and its connectivity are well constrained by outcrop observations, data from air-permeability measurements, bilinear interpolation of permeability, grid construction and water vector analysis for flow dynamics.These proportions were then enumerated in terms of static digital outcrop model(DOM) and facies model based on sandstone facies bedding characteristics. Flow simulation of water vector analysis through each of the four sandstone facies types show persistent spatial correlation of permeability that align with either cross-bedded orientation or straight with more dispersion high quality sandstone(porosity 21.25%-41.2%and permeability 1265.20-5986.25 mD) and moderate quality sandstone(porosity 10.44%-28.75% and permeability 21.44-1023.33 mD). Whereas, in more heterolithic sandstone(wavy-to flaser-bedded and bioturbated sandstone), lateral variations in permeability show spatially non-correlated patterns over centimeters to tens of meters with mostly of low quality sandstone(porosity 3.4%-12.31% and permeability < 1 mD to 3.21 mD). These variations reflect the lateral juxtaposition in flow dynamics. It has also been resulted that the vertical connectivity and heterogeneities in terms of flow are mostly pragmatic due to the interconnected sandstone rather than the quality of sandstone.     

How soil texture,channel shape and cross-sectional area affect moisture dynamics and water loss in irrigation channels

To enhance the utilization efficiency of farmland irrigation water and reduce the leakage of water conveyance channels, the leakage process of channels was simulated dynamically. The simulated results were compared with data measured in laboratory experiments, and the performance of the model was evaluated. The results indicated that the simulated values of the model were consistent with the observation values, and the R² values varied between 0.91 and 0.99. In addition, based on the laboratory experiments, a water supply system (Mariotte bottles) and soil box were built using plexiglass. Three influencing factors, namely, the channel form, soil texture and channel cross-sectional area, were varied to observe and calculate the resulting cumulative infiltration amount, infiltration rate and wetting front migration distance. HYDRUS-3D software was used to solve the three-dimensional soil water movement equation under different initial conditions. The results demonstrated that the U-shaped channel was more effective than the trapezoidal channel in increasing the utilization efficiency of the water resources. A U-shaped channel with a small channel cross-sectional area should be adopted and the soil particle size should be prioritized in the construction of water conveyance channels for farmlands. The simulation results were in agreement with the observed results, which indicates that HYDRUS-3D is a reliable tool that can accurately simulate the soil moisture movement in water conveyance channels. The research results can provide a reference for the design and operation of farmland irrigation systems.     

Changes in run‐off and sediment load in the three parts of the Yellow River basin,in response to climate change and human activities

Hydrological regimes in the Yellow River have changed significantly because of climate change and intensive human interventions. These changes present severe challenges to water resource utilization and ecological development. Variation of run‐off, suspended sediment load (SSL), and eight precipitation indices (P1: 0–12 mm·day^?1, P12: 12–25 mm·day^?1, P25: 25–50 mm·day^?1, P50: P ≥ 50 mm·day^?1 and corresponding rainfall day: Pd1, Pd12, Pd25, Pd50 day year^?1) in three critical parts of the Yellow River basin (source region: SRYRB, upper reaches: URYRB, middle reaches: MRYRB) were investigated for the period from 1960 to 2015. The results show that run‐off and SSL significantly decreased (P < 0.01) in the URYRB and the MRYRB, whereas their decline in the SRYRB was insignificant (P > 0.05). Moreover, run‐off in the URYRB had one change point in 1987, and SSL in the URYRB as well as run‐off and SSL in the MRYRB had two change points (in the 1970s and the 1990s). Over the same period, only Pd1 and Pd12 in the SRYRB showed significant increasing trends, and an abrupt change appeared in 1981. The optimal precipitation indices for assessing the effects of precipitation on run‐off and SSL in the URYRB and MRYRB were Pd50 and P12, respectively. A double‐mass curve analysis showed that precipitation and human activities contributed to approximately 20% and 80% of the reduction in run‐off, respectively, for both the SRYRB and the MRYRB. However, the contribution rate of precipitation and human activities on SSL reduction was approximately 40% and 60% in the URYRB and 5% and 95% in the MRYRB, respectively. Human activities, primarily soil and water conservation measures and water extraction (diversion), were the main factors (>50%) that reduced the run‐off. However, the dominant driving factors for SSL reduction were soil and water conservation measures and reservoir interception, for which the contribution rate was higher than 70% in the MRYRB. This work strengthens the understanding of hydrological responses to precipitation change and provides a useful reference for regional water resource utilization.