晚中新世以来中国北方风成沉积的磁性地层学和沉积学研究

The past two decades have witnessed an unprecedented growth of interest in the palaeoenvironmental significance of the Pleistocene loess deposits in northern China. However, it is only several years ago that the Tertiary red clay sequence underlying Pleistocene loess attracted much attention. One of the major advances in recent studies of eolian deposits on the Loess Plateau is the verification of the eolian origin for the Tertiary red clay sediments. The evidence of the eolian origin for the red clay is mainly from geochemical and sedimentological studies. However, sedimentological studies of the red clay deposits are still few compared with those of the overlying loess sediments. To date, the red clay sections located near Xifeng, Baoji, Lantian, Jiaxian, and Lingtai have been studied, with an emphasis on magnetostratigraphy. These sections have a basal age ranging from ～4.3 Ma to ～7.0 Ma. The thickness of the sections varies significantly, depending perhaps on the development of local geomorphological conditions and the drainage system. Although the stratigraphy of the red clay sections has been recorded in some detail, correlation of the red clay sequences has not yet been undertaken. Geological records (Sun J. et al., 1998) have shown that during glacial periods of the Quaternary the deserts in northem China were greatly expanded compared with modern desert distribution. During interglacial periods, desert areas contracted and retreated mostly to northwestern China because of the increase in inland penetration of monsoonal precipitation. According to pedogenic characteristics of the red clay deposits, the climatic conditions of the Loess Plateau is warmer and wetter generally in the Neogene than in the late Pleistocene. Panicle analyses show that grain size distribution of the red clay sequence is similar to that of the paleosols in the Pleistocene loess record, thus implying a relatively remote provenance of the red clay materials. However, the quantitative or semiquantitative estimates of the distance from the source region to the Loess Plateau during the red clay development remains to be investigated. In this study, magnetostratigraphic and sedimentological studies are conducted at two thick red clay sequences-Jingchuan and Lingtai section. The objectives of these studies are focused on further sedimentological evidence for the eolian origin of the red clay, correlation of red clay sequences, provenance of the red clay, and the palaeoclimate reconstruction in the Neogene. Paleomagnetic studies show that the Jingchuan red clay has a basal age of 8.0 Ma, which is 1 million years older than the previously studied Lingtai section. The Lingtai red clay sequence was divided into five units on the basis of pedogenica characteristics (Ding et al., 1999a). The Jingchuan red clay sequence, however, can be lithologically divided into six units according to field observations. The upper five units of the Jingchuan red clay can generally correlate well with the five units of the Lingtai red clay. Comparison of magnetic susceptibility and color reflectance records of four red clay sections suggests that the Lingtai red clay sequence can be the type-section of the Neogene red clay deposits in northern China. Pleistocene loess and modem dust deposits have a unimodal grain-size distribution. The red clay sediments at Jingchuan and Lingtai also have a unimodal grain-size distribution especially similar to the paleosols in the Pleistocene loess record. Sedimentological studies of a north-south transect of loess deposits above S2 on the Loess Plateau show that loess deposits had distinct temporal and spatial sedimentary differentiation. The characteristics of such sedimentary differentiation can be well presented in a triangular diagram of normalized median grain size, normalized skewness, and normalized kurtosis. The triangular diagrams of the red clay-loess sequence at Lingtai and Jingchuan indicate that loess-paleosol-red clay may be transported and sorted by the same agent wind, thus extending the eolian record in the Loess Plateau from 2.6 Ma back to about 8.0 Ma. It has been recognized that during the last glacial maximum (LGM) the deserts in northern China had a distribution similar to the present, whereas during the Holocene Optimum the deserts retreated to the area west of the Helan Mountains. Advance-retreat cycles of the deserts will lead to changes in the distance of the Loess Plateau to the dust source regions, thereby controlling changes in grain size of the loess deposited in a specific site. To observe spatial changes in sedimentological characteristics of loess during the last glacial-interglacial cycle, the texture of loess was measured along the north-south transect of the Loess Plateau. Since the southern margin of the Mu Us desert during the LGM is already known, several models of grain size parameters versus the minimum distance from the source region to depositional areas were developed. According to these semiquantitative models, the minimum distance from the source region to Lingtai and Jingchuan areas is about 600 km during the Neogene. Therefore the estimated provenance of the Tertiary red clay deposits is the areas now occupied by the Badain Jaran desert and arid regions west of it. The ratio of the free iron to total iron concentration attests to being a good proxy indicator for the summer monsoon evolution. The Lingtai Fe_20_3 ratio record shows high values over three time intervals: 4.8-4.1 Ma, 3.4-2.6 Ma, and during the interglacial periods of the past 0.5 Ma. The increase in summer monsoon intensity over the three intervals also coincides with the well-developed soil characteristics. It is therefore concluded that the East-Asia summer monsoon has experienced a non-linear evolution since the late Miocene. In general, the East Asia summer monsoon was stronger in Neogene than in Quaternary and the strongest East Asia summer monsoon may occur between 4.1 and 4.8 Ma. The relatively small ice volume and high global temperature may be responsible for the strong summer monsoon during the early Pliocene.

山东莱芜中生代高镁侵入岩及其超镁铁质岩捕掳体的成因研究

In order to know better about the Phanerozoic lithosphere thinning process of Sino-Korea Plate, four Cretaceous intrusion complexes and their ultramafic xenoliths are investigated by this thesis, which are located in Laiwu, Shandong Province, Eastern China, a region far away from plate margin. The four complexes, Kuanshan, Jiaoyu, Jingniushan and Tietonggou, intruded into Archaeozoic granite gneiss and Paleozoic carbonate rocks with scam iron ore at their contact zone. The four complexes can be divided into two magma series, abyssal rocks for the early and hypabyssal rocks for the later. K-Ar dating show that the abyssal rocks intrusion began with 120 ±2 Ma and the hypabyssal rocks intruded about 113 Ma. Abyssal rocks, mainly made up of augite diorites, amphibole diorites and gabbros for the lesser, are chemically characterized with high-Mg (Mg#>0.5) high-K calcalklic rock, which are depleted with Nb, Ta and Ti related to LILE and extremely enriched with Sr and Pb. Comparatively, augite diorites are the most LREE enriched in abyssal rocks, and they show no Eu abnorrnity or weak positive Eu abnormity. Gabbros show the least LREE enrichment with a strong Eu abnormity relatively. In (~(87)Sr/~(86)Sr)_1 -ε Nd(T) diagram, the abyssal rocks show a mixing trend , a rapid change in ε Nd(T) with a relatively small change in (~(87)Sr/~(86)Sr)_1. Low radiogenic Sr and Pb composition with high radiogeic Nd composition indicate that the mixing processes happened in lower crust Melt-rock interactions in lower crust might be the most possible processes to produce these high-Mg and high-K calcalklic magmas. Hypabyssal rocks, mainly made up of granite porphyry and dioritic porphrite, show much higher ε Nd(T) than abyssal rocks. Granite porphyry are distinct with an adakite geochemical characteristics, high (La/Yb)_N, Sr/Y and low Rb/Sr ratio. The adakitic granite porphyry indicates a new lower crust produced by underplating within plate. Ultramafic xenoliths had been found only in augite diorites and amphibole diorites. Field investigations show that ultramafic xenoliths in augite diorites had been inherited from amphibole diorites, so ultramafic xenoliths had been only entrained by hydrous dioritic magma. Ultramafic xenoliths are mainly made up of dunite and harzburgite, orthopyroxenite and bistangite are the lessor. Coarse olivines in dunite often show many chromite exsolution lamellae. Opx in orthopyroxenite often show dense chromite exsolution lamellae. The presence of exsolution minerals indicates that ultramafic xenoliths had cooled before entraining. Metasomatism phenomenons are popular in dunite and harzburgite xenoliths, including two kinds of assemblage, cpx+phlogapite and opx+pl. The first metasomatism assemblage indicates an ancient enrichment. Rb-Sr dating of xenoliths shows that the ancient enrichment happened in 223 ± 7Ma. The second metasomatism assemblage indicates a recent, relatively not deep melt-rock interaction, which might be related with the genesis of the high-Mg high-K calcalklic rocks. Mineralogy and geochemistry indicate that these ultramafic xenoliths might sample the crust-mante transition zone (or the base of lower crust, moho). Investigation of high-Mg intrusions and their ultramafic xenoliths in Laiwu indicate that the thinning processes of Sino-Korea Plate can be divided into two stages. The first stage is lithosphere mantle thinning with crust thickening by underplating in lower crust. The second stage is that the thickened lower crust delaminated into the underlying mantle.

降雨触发滑坡的非线性模型与斜坡演化的混沌效应

This paper analyzes landsliding process by nonlinear theories, especially the influence mechanism of external factors (such as rainfall and groundwater) on slope evolution. The author investigates landslide as a consequence of the catastrophic slide of initially stationary or creeping slope triggered by a small perturbation. A fully catastrophe analysis is done for all possible scenarios when a continuous change is imposed to the control parameters. As the slip surface continues and erosion due to rainfall occurs, control parameters of the slip surface may evolve such that a previously stable slope may become unstable (e.g. catastrophe occurs), when a small perturbation is imposed. Thus the present analysis offers a plausible explanation to why slope failure occurs at a particular rainfall, which is not the largest in the history of the slope. It is found, by analysis on the nonlinear dynamical model of the evolution process of slope built, that the relationship between the action of external environment factors and the response of the slope system is complicatedly nonlinear. When the nonlinear action of slope itself is equivalent to the acting ability of external environment, the chaotic phenomenon appears in the evolution process of slope, and its route leading to chaos is realized with bifurcation of period-doublings. On the basis of displacement time series of the slope, a nonlinear dynamic model is set up by improved Backus generalized linear inversion theory in this paper. Due to the equivalence between autonomous gradient system and catastrophe model, a standard cusp catastrophe model can be obtained through variable substitution. The method is applied to displacement data of Huangci landslide and Wolongsi landslide, to show how slopes evolve before landsliding. There is convincing statistical evidence to believe that the nonlinear dynamic model can make satisfied prediction results. Most important of all, we find that there is a sudden fall of D, which indicates the occurrence of catastrophe (when D=0).

石油二次运移的物理与数值模拟实验及其应用

Facing the problems met in studies on predominant hydrocarbon migration pathways, experiments and numerical simulating were done in this thesis work to discuss the migration mechanisms. The aim is to analyze quantitatively the pathway pattern in basin scale and to estimate the hydrocarbon loss on the pathway that offer useful information for confirming the potential hydrocarbon accumulation. Based on our understandings on hydrocarbon migration and the fluid dynamic theory, a series of migration experiments were designed to observe the phenomena where kerosene is used as draining phase driven only by buoyancy force that expulses pore water. These experiments allow to study the formation of migration pathways, the distribution of non-wetting oil along these pathways, and the re-utilizing of previously existing pathways marked by residual traces etc. The types of pattern for migration pathways may be characterized by a phase diagram using two dimensionless numbers: the capillary number and the Bond number. The NMR technique is used to measure the average saturation of residual oil within the pathways. Based our experiment works and percolation concept, a numerical simulation model were proposed and realized. This model is therefore called as BP (Buoyancy Percolation) simulator, since buoyancy is taken as the main driving force in hydrocarbon migration. To make sure that BP model is applicable to simulate the process of oil secondary migration, the experimental phenomena are compared with those simulated with BP model by fractal method, and the result is positive. After then, we use BP simulator to simulate the process of migration of oil in the porous media saturated with water at different scale. And the results seem similar to those cited in literatures. In addition, our software is applied in Paris basin to predict the pathway of hydrocarbon migration happened in the Middle Jurassic reservoirs. It is found that the results obtained with our BP model are generally agree with Hindle (1997) and Bekeles'(1999), but our simulated migration pathway pattern and migration direction seem more reasonable than theirs.

油气流线模拟方法研究与系统研制

The modeling of petroleum flow path (petroleum charging) and the detail of corresponding software development are presented in this paper, containing principle of petroleum charging, quantitative method, and practical modeling in two oil fields. The Modeling of Petroleum Flow Path is based on the result of basin modeling, according to the principle of petroleum migrating along the shortest path from the source to trap, Petroleum System Dynamics (Prof. Wu Chonglong, 1998), the concept of Petroleum Migration and Dynamic Accumulation (Zhou Donyan, Li Honhui, 2002), etc. The simulation is done combing with all parameters of basin, and considering the flow potential, non-uniformity of source and porous layer. It's the extending of basin modeling, but not belong to it. It is a powerful simulating tool of petroleum system, and can express quantitatively every kind of geology elements of a petroleum basin, and can recuperate dynamically the geology processes with 3D graphics. At result, we can give a result that the petroleum flow shows itself the phenomena of main path, and without using the special theory such as deflection flow in fractures(Tian Kaiming, 1989, 1994, Zhang Fawang, Hou Xingwei, 1998), and flow potential(England, 1987). The contour map of petroleum flow quantitative show clearly where the coteau - dividing slot is, and which convergence region are the main flow path of petroleum, and where is the favorable play of petroleum. The farsighted trap can be determined if there are enough information about structural diagram and can be evaluated, such as the entrapment extent, spill point, area, oil column thickness, etc. Making full use of the result of basin modeling with this new tool, the critical moment and scheme of the petroleum generation and expulsion can be showed clearly. It's powerful analysis tool for geologist.

矾山杂岩体成因及其与磷（铁）矿的成矿关系

The Fanshan complex consists of layered potassic ultramafic-syenite intrusions. The Fanshan apatite (-magnetite) deposit occurs in the Fanshan complex, and is an important style of phosphorus deposit in China. The Fanshan complex consists of three (First- to Third-) Phases of intrusion, and then the dikes. The First-Phase Intrusive contains ten typical layered rocks: clinopyroxenite, biotite clinopyroxenite, coarse-grained biotite clinopyroxenite, pegmatitic orthoclase-biotite clinopyroxenite, variegated orthoclase clinopyroxenite, interstitial orthoclase clinopyroxenite, biotite rock, biotite-apatite rock, biotite rock and magnetite-apatite rock. This layered intrusive consists of nine rhythmic units. Each rhythmic unit essentially comprises a pair of layers: clinopyroxenite at the bottom and biotite clinopyroxenite at the top. The apatite (-magnetite) deposit is situated near the top of rhythmic Unit no. 6 of the First-Phase Intrusive. The Second-Phase Intrusive contains three typical rocks: coarse-grained orthoclase clinopyroxenite, . coarse-grained salite syenite and schorlomite-salite syenite. The Third-Phase Intrusive includes pseudo-trachytic salite syenite, porphyritic augite syenite, fine-grained orthoclase clinopyroxenite and fine-grained salite syenite. The origin of the Fanshan complex is always paid attention to it in China. Because most layered igneous intrusion in the world not only have important deposit in it, but also carry many useful information for studying the formation of the intrusion and the evolvement of magma. Two sketch maps were drawn through orebodies along no. 25 cross-cut on 425 mL and no. 1 cross-cut on 491 mL in the Fanshan mine. Through this mapping, a small-scaled rhythmic layering (called sub-rhythmic layering in the present study) was newly found at the top of the rhythmic Unit no. 6. The concept of sub-rhythmic layering is defined in this article. The sub-rhythmic layering is recognized throughout this apatite-rich part, except for magnetite-apatite rock. Presence of the layered magnetite-apatite rock is one of the characteristics of the Fanshan apatite (-magnetite) deposit. Thus, from this layer downwards six units of sub-rhythmic layering are recognized in the present study. Each unit consists of biotite clinopyroxenite (or biotite rock and biotite-apatite rock) layer at the bottom and apatite rock layer at the top. To study this feature in detail is an important work for understanding the origin of the Fanshan complex and apatite (-magnetite) deposit. The origin of the Fanshan complex and the relation of the formation of the apatite(-magnetite)deposit will be interpreted by the study of sub-rhythmic layering on the basis of previous research works. The magma formed the Fanshan complex was rich in K2O, early crystallized pyroxene, and after this phase more biotite crystallized, but no amphibole appeared. This indicated that the activity of H2O in the magma was low. Major element compositions of biotite and clinopyroxene (on thin sections) in the sub-rhythmic layering were analyzed using electron microprobe analyzer. The analytical results indicate Mg/(Mg+Fe*+Mn) atomic ratios (Fe*, total iron) of these two minerals rhythmically changed in sub-rhythmic layering. The trends of Mg/(Mg+Fe*+Mn) atomic ratio (Fe*, total iron) of biotite and clinopyroxene indicate that the magma evolved markedly from relatively magnesian bottom layer to less magnesian top layer in each sub-rhythmic unit. A general trend through the sub-rhythmic layering sequence is both minerals becoming relatively magnesian upwards. The formation temperatures for sub-rhythmic layering yield values between 600 and 800 ℃, were calculated using the ratio of Mg/(Mg+Fe+Mn) in the salite and biotite assemblage. The equilibrium pressures in the rhythmic layers calculated using the contents of Al in the salite were plotted in the section map, shown a concave curve. This indicates that the magma formed the First-Phase Intrusive crystallized by two vis-a-vis ways, from its bottom and top to its centre, and the magnetite-apatite rock was crytallized in the latest stage. The values of equilibrium pressures in the sub-rhythmic layering were 3.6-6.8(xlO8) Pa with calculated using the contents of Al in the salite. The characteristics of geochemistry in various intrusive rocks and the rocks or apatite of sub-rhythmic layers indicated that the Fanshan complex formed by the comagmatic crystallization. The contents of immiscible elements and REEs of apatite rock at the top of one sub-rhythmic unit are more than biotite clinopyroxenite at the bottom. The contents of immiscible elements and REEs of apatite of biotite clinopyroxenite at the bottom of one sub-rhythmic unit are higher than apatite rock at the top. The curves of rocks (or apatite) in the upper sub-rhythmic units are between two curves of the below sub-rhythmic unit in the primitive mantle-normalized trace element abundance spider diagram and the primitive mantle-normalized REE pattern. The trend for the contents of immiscible elements and REEs inclines to the same contents from the bottom to the top in sub-rhythmic layering. These characteristics of geochemistry of rocks or apatites from sub-rhythmic layering indicate that the latter sub-rhythmic unit was produced by the residual magma after crystallization of the previous sub-rhythmic unit. The characteristics of petrology, petrochemistry, geochemistry in the Fanshan complex and sub-rhythmic layers and the trends of Mg/(Mg+Fe+Mn) atomic ratio of biotite and clinopyroxene in sub-rhytmic layering rejected the hypotheses, such as magma immiscibility, ravitational settling and multiple and pulse supplement of magma. The hypothesis of differentiation by crystallization lacks of evidences of field and excludes by this study. On the base of the trends of formation temperatures and pressures, the characteristics of petrology, petrochemistry, geochemistry for the Fanshan complex and the characteristics of geochemistry for the rocks (or apatites), the trends of Mg/(Mg+Fe+Mn) atomic ratio of biotite and clinopyroxene in sub-rhytmic layering, and the data of oxygen, hydrogen, strontium and neodymium isotopes, this study suggests that the magma formed the Fanshan complex was formed by low degree partial melting of mantle at a low activity of H2O, and went through the differentiation at the depth of mantle, then multiply intruded and crystallized. The rhythmic layers of the First-Phase Intrusive formed by the magma fractional crystallized in two vis-a-vis ways, from the bottom and top to the centre in-situ fractional crystallization. The apatite (-magnetite) deposit of the Fanshan complex occurs in sub-rhythmic layering sequence. The the origin of the sub-rhythmic layering is substantially the origin of the Fanshan apatite (-magnetite) deposit. The magma formed the rhythmic layers of First-Phase Intrusive was rich in H2O, F and P at the later stage of its in-situ fractional crystallization. The Fanshan apatite (-magnetite) deposit was formed by this residual magma in-situ fractional crystallization. The magnetite-apatite rock was crystallized by two vis-a-vis ways at the latest stage in-situ fractional crystallization in the rhythmic layers. The result was light apatite layer below heavy the magnetite-apatite layer, formed an "inversion" phenomenon.

乡镇级人口统计数据空间化研究

Population data which collected and saved according to administrative region is a kind of statistical data. As a traditional method of spatial data expression, average distribution in every administrative region brings population data on a low spatial and temporal precision. Now, an accurate population data with high spatial resolution is becoming more and more important in regional planning, environment protection, policy making and rural-urban development. Spatial distribution of population data is becoming more important in GIS study area. In this article, the author reviewed the progress of research on spatial distribution of population. Under the support of GIS, correlative geographical theories and Grid data model, Remote Sensing data, terrain data, traffic data, river data, resident data, and social economic statistic were applied to calculate the spatial distribution of population in Fujian province, which includes following parts: (1) Simulating of boundary at township level. Based on access cost index, land use data, traffic data, river data, DEM, and correlative social economic statistic data, the access cost surface in study area was constructed. Supported by the lowest cost path query and weighted Voronoi diagram, DVT model (Demarcation of Villages and Towns) was established to simulate the boundary at township level in Fujian province. (2) Modeling of population spatial distribution. Based on the knowledge in geography, seven impact factors, such as land use, altitude, slope, residential area, railway, road, and river were chosen as the parameters in this study. Under the support of GIS, the relations of population distribution to these impact factors were analyzed quantificationally, and the coefficients of population density on pixel scale were calculated. Last, the model of population spatial distribution at township level was established through multiplicative fusion of population density coefficients and simulated boundary of towns. (3) Error test and analysis of population spatial distribution base on modeling. The author not only analyzed the numerical character of modeling error, but also its spatial distribution. The reasons of error were discussed.

安多地区早中侏罗世岩浆岩岩石地球化学和同位素年代学研究及其对构造演化的制约

Anduo area is located in the Central Tibet, the middle segment of the Bangonghu-Nujiang suture. Anduo Block is the northern part of Lhasa terrane. The relationships among the different geological bodies were determined during the 1: 250000 regional geological surveying. Petrography, petrologic geochemistry, isotopic geochemistry and geochronology of igneous rocks from the suture and granitoids from Anduo Block were analyzed systematically as a whole for the first time. Then, their tectonic setting and history are discussed.Anduo ophiolitic melange consists of metamorphic peridotites, cumulates, plagiogranites, sheeted dykes swarm, pillow lava and radiolarian cherts. The concentration of Cr and Ni in the metamorphic peridotites is very high, with Mg# about 0.94 ~ 0.97, higher 87Sr/86Sr and Pb isotopic ratios, and lower 143Nd/i44Nd ratio. LREE is enriched relative to HREE and positive Eu anomaly is very clear. The REE distribution curve is U shape. Nb and Ta anomalies from cumulate gabbro and sheeted dyke swarm are not clear, while that are slightly negative from pillow lava. Plagiogranite belongs to strong calc-alkaline series with high Si, middle Al, low Fe, Mg and low K contents. Eu anomaly (~ 1.23) from plagiogranites is slightly positive. The character of all components of ophiolite is similar to that of the MORB, while to some extent the ophiolite was influenced by crustal material. Anduo ophiolite formed in a mature back-arc basin. Additionally, intermediate acidity volcanic rocks within Anduo phiolite melange are island arc calc-alkline rocks related to ocean subduction.The early-middle Jurassic plutonic rocks are tonalite, granodiorite bearing-phenocryst, magaporphyritic hornblende monzogranite, magaporphyritic monzogranite, monzogranite bearing-phenocryst and syenogranite in turn. They belong to calc-alkaline series which developed from middle K to high K series temporally. REE distribution curves of all plutonic rocks are similar and parallel to each other. SREE and negative Eu anomaly values decrease. In the multi-element spider diagram, the curves of different plutons are similar to each other, but troughs of Nb, Sr, P and Ti from young plutons become more evident. This suggests that thereare some closely petrogenetic affinities among plutonic rocks which make up amagma plutonism cycle of the early-middle Jurassic. Magma source is mainly crustal,but abundant mafic microgranular enclaves within granitoids indicate that crastalmagma should be mixed with mantle-derived magma and the mantle-derived magmadecreased subsequently. Tonalite has features of I-type granite, magaporphyriticmonzogranite is transition type, and monzogranite bearing-phenocryst is S-typegranite. The characteristic of granitoids from Anduo Block suggest that the formingtectonic setting is active continental margin.Reliable zircon U-Pb SHRIMP ages are obtained in the study area firstly. Plagiogranite from the Anduo ophiolite of the Bangonghu-Nujiang suture is 175.1 Ma, and granitoids from Anduo Block is 172.6-185.4 Ma. Additionally, plagioclase from the plagiogranite dates a 40Ar/39Ar age of 144 Ma, while biotite and hornblend from granitoids of Anduo Block give a 163-165 Ma.Similar cooling ages of plagiogranite from the Anduo ophiolitic melange and granitoids from Anduo Block and the spatial distribution of the ophiolitic rocks between Anduo, Naqu, and Shainzha area suggest that bilateral subduction of the Bangonghu-Nujiang oceanic basin took place in the early-middle Jurassic. During this subduction, Anduo ophiolitic rocks were related to north subduction of the Bangonghu-Nujiang oceanic basin and Anduo back-arc basin spreading, while granitoids from Anduo Block were related to south subduction.

Stable Mixing of Complete and Incomplete Information

An increasing number of parameter estimation tasks involve the use of at least two information sources, one complete but limited, the other abundant but incomplete. Standard algorithms such as EM (or em) used in this context are unfortunately not stable in the sense that they can lead to a dramatic loss of accuracy with the inclusion of incomplete observations. We provide a more controlled solution to this problem through differential equations that govern the evolution of locally optimal solutions (fixed points) as a function of the source weighting. This approach permits us to explicitly identify any critical (bifurcation) points leading to choices unsupported by the available complete data. The approach readily applies to any graphical model in O(n^3) time where n is the number of parameters. We use the naive Bayes model to illustrate these ideas and demonstrate the effectiveness of our approach in the context of text classification problems.

Summarizing Qualitative Behavior from Measurements of NonlinearsCircuits

This report describes a program which automatically characterizes the behavior of any driven, nonlinear, electrical circuit. To do this, the program autonomously selects interesting input parameters, drives the circuit, measures its response, performs a set of numeric computations on the measured data, interprets the results, and decomposes the circuit's parameter space into regions of qualitatively distinct behavior. The output is a two-dimensional portrait summarizing the high-level, qualitative behavior of the circuit for every point in the graph, an accompanying textual explanation describing any interesting patterns observed in the diagram, and a symbolic description of the circuit's behavior which can be passed on to other programs for further analysis.

A Study of Qualitative and Geometric Knowledge in Reasoning about Motion

Reasoning about motion is an important part of our commonsense knowledge, involving fluent spatial reasoning. This work studies the qualitative and geometric knowledge required to reason in a world that consists of balls moving through space constrained by collisions with surfaces, including dissipative forces and multiple moving objects. An analog geometry representation serves the program as a diagram, allowing many spatial questions to be answered by numeric calculation. It also provides the foundation for the construction and use of place vocabulary, the symbolic descriptions of space required to do qualitative reasoning about motion in the domain. The actual motion of a ball is described as a network consisting of descriptions of qualitatively distinct types of motion. Implementing the elements of these networks in a constraint language allows the same elements to be used for both analysis and simulation of motion. A qualitative description of the actual motion is also used to check the consistency of assumptions about motion. A process of qualitative simulation is used to describe the kinds of motion possible from some state. The ambiguity inherent in such a description can be reduced by assumptions about physical properties of the ball or assumptions about its motion. Each assumption directly rules out some kinds of motion, but other knowledge is required to determine the indirect consequences of making these assumptions. Some of this knowledge is domain dependent and relies heavily on spatial descriptions.

Theory of Handwriting

Handwriting production is viewed as a constrained modulation of an underlying oscillatory process. Coupled oscillations in horizontal and vertical directions produce letter forms, and when superimposed on a rightward constant velocity horizontal sweep result in spatially separated letters. Modulation of the vertical oscillation is responsible for control of letter height, either through altering the frequency or altering the acceleration amplitude. Modulation of the horizontal oscillation is responsible for control of corner shape through altering phase or amplitude. The vertical velocity zero crossing in the velocity space diagram is important from the standpoint of control. Changing the horizontal velocity value at this zero crossing controls corner shape, and such changes can be effected through modifying the horizontal oscillation amplitude and phase. Changing the slope at this zero crossing controls writing slant; this slope depends on the horizontal and vertical velocity zero amplitudes and on the relative phase difference. Letter height modulation is also best applied at the vertical velocity zero crossing to preserve an even baseline. The corner shape and slant constraints completely determine the amplitude and phase relations between the two oscillations. Under these constraints interletter separation is not an independent parameter. This theory applies generally to a number of acceleration oscillation patterns such as sinusoidal, rectangular and trapezoidal oscillations. The oscillation theory also provides an explanation for how handwriting might degenerate with speed. An implementation of the theory in the context of the spring muscle model is developed. Here sinusoidal oscillations arise from a purely mechanical sources; orthogonal antagonistic spring pairs generate particular cycloids depending on the initial conditions. Modulating between cycloids can be achieved by changing the spring zero settings at the appropriate times. Frequency can be modulated either by shifting between coactivation and alternating activation of the antagonistic springs or by presuming variable spring constant springs. An acceleration and position measuring apparatus was developed for measurements of human handwriting. Measurements of human writing are consistent with the oscillation theory. It is shown that the minimum energy movement for the spring muscle is bang-coast-bang. For certain parameter values a singular arc solution can be shown to be minimizing. Experimental measurements however indicate that handwriting is not a minimum energy movement.

Motion Planning with Six Degrees of Freedom

The motion planning problem is of central importance to the fields of robotics, spatial planning, and automated design. In robotics we are interested in the automatic synthesis of robot motions, given high-level specifications of tasks and geometric models of the robot and obstacles. The Mover's problem is to find a continuous, collision-free path for a moving object through an environment containing obstacles. We present an implemented algorithm for the classical formulation of the three-dimensional Mover's problem: given an arbitrary rigid polyhedral moving object P with three translational and three rotational degrees of freedom, find a continuous, collision-free path taking P from some initial configuration to a desired goal configuration. This thesis describes the first known implementation of a complete algorithm (at a given resolution) for the full six degree of freedom Movers' problem. The algorithm transforms the six degree of freedom planning problem into a point navigation problem in a six-dimensional configuration space (called C-Space). The C-Space obstacles, which characterize the physically unachievable configurations, are directly represented by six-dimensional manifolds whose boundaries are five dimensional C-surfaces. By characterizing these surfaces and their intersections, collision-free paths may be found by the closure of three operators which (i) slide along 5-dimensional intersections of level C-Space obstacles; (ii) slide along 1- to 4-dimensional intersections of level C-surfaces; and (iii) jump between 6 dimensional obstacles. Implementing the point navigation operators requires solving fundamental representational and algorithmic questions: we will derive new structural properties of the C-Space constraints and shoe how to construct and represent C-Surfaces and their intersection manifolds. A definition and new theoretical results are presented for a six-dimensional C-Space extension of the generalized Voronoi diagram, called the C-Voronoi diagram, whose structure we relate to the C-surface intersection manifolds. The representations and algorithms we develop impact many geometric planning problems, and extend to Cartesian manipulators with six degrees of freedom.

The topological behaviour of 3D null points in the Sun's corona

Brown D. S. and Priest E. R. 2001, The topological behaviour of 3D null points in the Sun's corona, Astronomy and Astrophysics, 367, 339-346

A Topological Analysis of the Magnetic Breakout Model for an Eruptive Solar Flare

Maclean, R., Beveridge, C., Longcope, D.W., Brown, D.S. and Priest, E.R., 2005, A topological analysis of the magnetic breakout model for an eruptive solar flare, Proc. Roy. Soc., 461, 2099-2120. Sponsorship: PPARC/STFC