IMPACT OF DIABATIC PROCESSES IN AGCM ON 4-DIMENSIONAL VARIATIONAL DATA ASSIMILATION*

The impact of diabatic processes on 4-dimensional variational data assimilation (4D-Var) was studied using the 1995 version of NCEP's global spectral model with and without full physics.The adjoint was coded manually.A cost function measuring spectral errors of 6-hour forecasts to "observation" (the NCEP reanalysis data) was minimized using the L-BFGS (the limited memory quasi-Newton algorithm developed by Broyden,Fletcher,Goldfard and Shanno) for optimizing parameters and initial conditions.Minimization of the cost function constrained by an adiabatic version of the NCEP global model converged to a minimum with a significant amount of decrease in the value of the cost function.Minimization of the cost function using the diabatic model, however,failed after a few iterations due to discontinuities introduced by physical parameterizations.Examination of the convergence of the cost function in different spectral domains reveals that the large-scale flow is adjusted during the first 10 iterations,in which discontinuous diabatic parameterizations play very little role.The adjustment produced by the minimization gradually moves to relatively smaller scales between 10-20th iterations.During this transition period,discontinuities in the cost function produced by "on-off" switches in the physical parameterizations caused the cost function to stay in a shallow local minimum instead of continuously decreasing toward a deeper minimum. Next,a mixed 4D-Var scheme is tested in which large-scale flows are first adiabatically adjusted to a sufficient level,followed by a diabatic adjustment introduced after 10 to 20 iterations. The mixed 4D-Var produced a closer fit of analysis to observations,with 38% and 41% more decrease in the values of the cost function and the norm of gradient,respectively,than the standard diabatic 4D-Var,while the CPU time is reduced by 21%.The resulting optimal initial conditions improve the short-range forecast skills of 48-hour statistics.The detrimental effect of parameterization discontinuities on minimization was also reduced.     

Noril’sk–Talnakh Cu–Ni–PGE deposits: a revised tectonic model

The Norilsk mining district is located at the northwest margin of the Tunguska basin, in the centre of the 3,000×4,000 km Siberian continental flood basalt (CFB) province. This CFB province was formed at the Permo-Triassic boundary from a superplume that ascended into the geometric centre of the Laurasian continent, which was surrounded by subducting slabs of oceanic crust. We suggest that these slabs could have reached the core–mantle boundary, and they may have controlled the geometric focus of the superplume. The resulting voluminous magma intruded and erupted in continental rifts and related extensive flood basalt events over a 2–4 Ma period. Cu–Ni–PGE sulfide mineralization is found in olivine-bearing differentiated mafic intrusions beneath the flood basalts at the northwestern margin of the Siberian craton and also in the Taimyr Peninsula, some 300 km east of a triple junction of continental rifts, now buried beneath the Mesozoic–Cenozoic sedimentary basin of western Siberia. The Norilsk-I and Talnakh-Oktyabrsky deposits occur in the Norilsk–Kharaelakh trough of the Tunguska CFB basin. The Cu–Ni–PGE-bearing mineralized intrusions are 2–3 km-wide and 20 km-long differentiated chonoliths. Previous studies suggested that parts of the magma remained in intermediate-level crustal chambers where sulfide saturation and accumulation took place before emplacement. The 5–7-km-thick Neoproterozoic to Palaeozoic country rocks, containing sedimentary Cu mineralization and evaporites, may have contributed additional metal and sulfur to this magma. Classic tectonomagmatic models for these deposits proposed that subvertical crustal faults, such as the northeast-trending Norilsk–Kharaelakh fault, were major trough-parallel conduits providing access for magmas to the final chambers. However, geological maps of the Norilsk region show that the Norilsk–Kharaelakh fault offsets the mineralization, which was deformed into folds and offset by related reverse faults, indicating compressional deformation after mineralization in the Late Triassic to Early Jurassic. In addition, most of the intrusions are sills, not dykes as should be expected if the vertical faults were major conduits. A revised tectonic model for the Norilsk region takes into account the fold structure and sill morphology of the dominant intrusions, indicating a lateral rather than vertical emplacement direction for the magma into final chambers. Taking into account the fold structure of the country rocks, the present distribution of the differentiated intrusions hosting the Norilsk-I and Talnakh–Oktyabrsky deposits may represent the remnants of a single, >60 km long, deformed and eroded palm-shaped cluster of mineralized intrusions, which are perceived as separate intrusions at the present erosional level. The original direction of sill emplacement may have been controlled by a northeast-trending paleo-rise, which we suggest is present at the southeastern border of the Norilsk–Kharaelakh trough based on analysis of the unconformity at the base of the CFB. The mineralized intrusions extend along this rise, which we interpret as a structure that formed above the extensionally tilted block in the metamorphic basement. Geophysical data indicate the presence of an intermediate magma chamber that could be linked with the Talnakh intrusion. In turn, this T-shaped flat chamber may link with the Yenisei–Khatanga rift along the northwest-trending Pyasina transform fault, which may have served as the principal magma conduit to the intermediate chamber. It then produced the differentiated mineralized intrusions that melted through the evaporites with in situ precipitation of massive, disseminated, and copper sulfide ore. The Norilsk–Kharaelakh crustal fault may not relate to mineralization and possibly formed in response to late Mesozoic spreading in the Arctic Ocean.Editorial handling: P. Lightfoot     

Structural changes and valence states in the MgCr<Subscript>2</Subscript>O<Subscript>4</Subscript>–FeCr<Subscript>2</Subscript>O<Subscript>4</Subscript> solid solution series

The influence on the structure of Fe²⁺ Mg substitution was studied in synthetic single crystals belonging to the MgCr₂O₄–FeCr₂O₄ series produced by flux growth at 900–1200 °C in controlled atmosphere. Samples were analyzed by single-crystal X-ray diffraction, electron microprobe analyses, optical absorption-, infrared- and Mössbauer spectroscopy. The Mössbauer data show that iron occurs almost exclusively as ^IVFe²⁺. Only minor Fe³⁺ (<0.005 apfu) was observed in samples with very low total Fe. Optical absorption spectra show that chromium with few exceptions is present as a trivalent cation at the octahedral site. Additional absorption bands attributable to Cr²⁺ and Cr³⁺ at the tetrahedral site are evident in spectra of end-member magnesiochromite and solid-solution crystals with low ferrous contents. Structural parameters a₀, u and T–O increase with chromite content, while the M–O bond distance remains nearly constant, with an average value equal to 1.995(1) Å corresponding to the Cr³⁺ octahedral bond distance. The ideal trend between cell parameter, T–O bond length and Fe²⁺ content (apfu) is described by the following linear relations: a₀=8.3325(5) + 0.0443(8)Fe²⁺ (Å) and T–O=1.9645(6) + 0.033(1)Fe²⁺ (Å) Consequently, Fe²⁺ and Mg tetrahedral bond lengths are equal to 1.998(1) Å and 1.965(1) Å, respectively.     

Stable carbon isotope ratios in lake and swamp sediments as a proxy for prehistoric forest clearance and crop cultivation in the Neotropics

Close correspondence between stable carbon isotope ratios ( ¹³ C), pollen, and charcoal profiles in sediment cores from Laguna Zoncho and Machita swamp, Costa Rica, shows that prehistoric forest clearance and crop cultivation can be detected in the stable carbon isotope ratios of total organic carbon ( ¹3C _TOC ). Analyses of δ 13C _TOC complement evidence from pollen, charcoal, and phytoliths and provide a proxy that is sensitive to the intensity and/or proximity to core sites of prehistoric forest clearance and agriculture in watersheds. Stable carbon isotope analyses are particularly useful in situations in which other evidence of forest clearance and agriculture is limited.     

Effects of the number of stations and time resolution on derivation

For at least 30 years now it has been well known that the Dst index can be modelled using the solar wind as input. Since then, many attempts have been made to improve the predictability of Dst using different approaches. These attempts are useful, for instance, to understand which features of the solar wind–magnetosphere interactions are most important in producing magnetospheric activity and how the Dst index would improve the space weather forecast. The Dst index is by far the most reliable and simple indication that a magnetic storm is in progress. In this work, the effect of using more than four magnetic stations and shorter time intervals than the hourly averages used in Sugiura's procedures is evaluated. The discussion is based on the results presented by Burton in 1975 and Feldstein in 1984 considering 4 or 12 magnetic stations and time averages of 2.5 min for a magnetic disturbed period that occurred from February 7–28, 1967, including two geomagnetic storms. The analysis has shown that the global representation of a magnetic storm by the standard Dst (Sugiura) is well preserved either using 4, 6, 12 magnetic stations or using 1 h, 2.5 min 1 min averages. A brief review of the current understanding of Dst has been included to support the discussions. The analysis performed has shown that a more refined Dst index (time and number of stations>4) would be useful to investigate the intrinsic processes and the different current systems involved in the ring current development during magnetic storms; the standard Dst, as it is conceived, is quite adequate to monitor geomagnetic storms and identify their overall features; concerning the magnetic stations normally considered, the inclusion of higher magnetic latitude stations (>35) may underestimate the observed Dst.     

Unstable kinetic Alfvén wave in partially ionized plasma

The stability of kinetic Alfven waves is discussed for a partially ionized plasma with a flux of ionizing electrons which balance the plasma particle losses. Accidental electromagnetic perturbations are shown to be unstable due to the energy change of ionizing electrons.     

4D一体化多媒体电子地图集理论方法和技术研究与开发--以《云南省生态环境多媒体电子地图集》为例

本文首先研究了4D数据一体化多媒体电子地图集的理论和信息机理,包括4D数据在电子地图集中的集成、融合和应用的信息机理,多源多类信息交互认知和传输机理,4D信息的查询、分析、模拟、调控和表达的信息机理,以及4D数据一体化电子地图集的理论框架。其次,文章深入研究和解决了建立4D一体化多媒体电子地图集过程中的诸多技术方法问题,例如多源多类数据的集成和融合、基于主导数据库的多重表达、多媒体地图动画、基于4D数据的信息分析和量算、生态环境信息动态仿真和虚拟表达等。在论文第三部分,文章介绍了作者研发的《云南省生态环境多媒体电子地图集》成果,包括其结构、功能、内容、超媒体链接等方面的设计,查询功能和分析量算功能的开发,以及本电子地图集的分析应用成果。     

POS辅助航空摄影测量精度分析

通过陕西勉县测区的POS辅助航空摄影测量试验,本文对SWDC-4搭载POS/AV 510 V5系统获取的航测数据使用直接定向法和POS辅助空中三角测量2种方法进行1∶10 000成图精度验证。结果证明:国产SWDC-4航摄仪与POS/AV系统集成进行POS辅助空三的加密成果能满足该测区1∶10 000地形图测图的平面和高程精度要求,直接定向也是可行的。     

利用遥感技术建立墨西哥热带地上生物量模型(英文)

Spatially-explicit estimation of aboveground biomass(AGB) plays an important role to generate action policies focused in climate change mitigation,since carbon(C) retained in the biomass is vital for regulating Earth’s temperature.This work estimates AGB using both chlorophyll(red,near infrared) and moisture(middle infrared) based normalized vegetation indices constructed with MCD43A4 MODerate-resolution Imaging Spectroradiometer(MODIS) and MOD44B vegetation continuous fields(VCF) data.The study area is located in San Luis Potosí,Mexico,a region that comprises a part of the upper limit of the intertropical zone.AGB estimations were made using both individual tree data from the National Forest Inventory of Mexico and allometric equations reported in scientific literature.Linear and nonlinear(expo-nential) models were fitted to find their predictive potential when using satellite spectral data as explanatory variables.Highly-significant correlations(p = 0.01) were found between all the explaining variables tested.NDVI62,linked to chlorophyll content and moisture stress,showed the highest correlation.The best model(nonlinear) showed an index of fit(Pseudo-r2) equal to 0.77 and a root mean square error equal to 26.00 Mg/ha using NDVI62 and VCF as explanatory variables.Validation correlation coefficients were similar for both models:lin-ear(r = 0.87**) and nonlinear(r = 0.86**).     

一种平面四参数法坐标转换方法的实现

简述了GPS坐标转换技术,详细介绍了一种平面四参数坐标转换的方法,然后采用VB6.0编写程序实现此坐标转换方法,并对这种坐标转换方法进行了验证。程序界面简捷,直观,实用,程序精度较高,在平时的工作中的到较好的应用。