Efficient generation of conditional simulations by chebyshev matrix polynomial approximations to the symmetric square root of the covariance matrix

Consider the problem of generating a realization y₁ of a Gaussian random field on a dense grid of points ₁ conditioned on field observations y₂ collected on a sparse grid of points ₂. An approach to this is to generate first an unconditional realization y over the grid =_{1 2}, and then to produce y₁ by conditioning y on the data y₂. As standard methods for generating y, such as the turning bands, spectral or Cholesky approaches can have various limitations, it has been proposed by M. W. Davis to generate realizations from a matrix polynomial approximations to the square root of the covariance matrix. In this paper we describe how to generate a direct approximation to the conditional realization y₁, on ₁ using a variant of Davis' approach based on approximation by Chebyshev polynomials. The resulting algorithm is simple to implement, numerically stable, and bounds on the approximation error are readily available. Furthermore we show that the conditional realization y₁ can be generated directly with a lower order polynomial than the unconditional realization y, and that further reductions can be achieved by exploiting a nugget effect if one is present. A pseudocode version of the algorithm is provided that can be implemented using the fast Fourier transform if the field is stationary and the grid ₁ is rectangular. Finally, numerical illustrations are given of the algorithm's performance in generating various 2-D realizations of conditional processes on large sampling grids.     

On the physical model of earthquake precursor fields and the mechanism of precursors＇timespace distribution──origin and evidences of the strong body earthquake──generating model

Ｏｎｔｈｅｐｈｙｓｉｃａｌｍｏｄｅｌｏｆｅａｒｔｈｑｕａｋｅｐｒｅｃｕｒｓｏｒｆｉｅｌｄｓａｎｄｔｈｅｍｅｃｈａｎｉｓｍｏｆｐｒｅｃｕｒｓｏｒｓ＇ｔｉｍｅｓｐａｃｅｄｉｓｔｒｉｂｕｔｉｏｎ──ｏｒｉｇｉｎａｎｄｅｖｉｄｅｎｃｅｓｏｆ　ｔｈｅｓｔｒｏｎ...     

Strain analysis in Jurassic argillites of the Monte Sirino area (Lagonegro Zone,southern Apennines,Italy) and implications for deformation paths in pelitic rocks

Analysis of strain in Jurassic argillites forming part of the folded and thrusted sedimentary succession of the Lagonegro basin (southern Italian Apennines) has been carried out using ellipsoid-shaped reduction spots as strain markers. Most of the determined finite strain ellipsoids are of oblate type and show a peculiar distribution of the maximum extension direction (X), with maxima either subparallel or subperpendicular to the local fold axes. Using the strain matrix method, two different deformation histories have been considered to assist the interpretation of the observed finite strain pattern. A first deformation history involved vertical compaction followed by horizontal shortening (occurring by a combination of true tectonic strain and volume loss), whereby all strain is coaxial and there is no change in the intermediate axis of the strain ellipsoid. By this type of deformation sequence, which produces a deformation path where total strain moves from the oblate to the prolate strain field and back to the oblate field, prolate strain ellipsoids can be generated and may be recorded where tectonic deformation has not been large enough to reverse pretectonic compaction. This type of deformation history may be of local importance within the study area (i.e. it may characterize some fold hinge regions) and, more generally, is probably of limited occurrence in deformed pelitic rocks. A second deformation sequence considered the superposition of pre-tectonic compaction and tectonic strain consisting of initial layer-parallel shortening followed by layer-parallel shear (related to flexural folding). Also in this instance, volume change during tectonic deformation and tectonic plane strain have been assumed. For geologically reasonable amounts of volume loss due to compaction and of initial layer-parallel shortening, this type of deformation history is capable of producing a deformation path entirely lying within the oblate strain field, but still characterized by a changeover, during deformation, of the maximum extension axis (X) from a position parallel to the fold axis to one perpendicular to it. This type of deformation sequence may explain the main strain features observed in the study area, where most of the measured finite strain ellipsoids, determined from the limb regions of flexural folds, display an oblate shape, irrespective of the orientation of their maximum extension direction (X) with respect to the local structural trends. More generally, this type of deformation history provides a mechanism to account for the predominance of oblate strains in deformed pelitic rocks.     

中国三个大地震的震源参数及讨论

大地震的震源机制及地震矩、应力降等参数,不(又对认识地震的破裂过程,而且对预报强地面运动都是非常重要的。一些文章对板内地震应力降的特点、地震矩与高频源谱的关系及区域特征的讨论,使我们可直接由长周期体波或面波得出的地震矩M₀估计出不同地区中大地震的高频源谱,以供强地面运动预报参考。 由于对发生在中国板内的1976年11月6日盐源地震（M_s=6.5）和1976年11月15日宁河地震（M_s=6.3）的震源参数还没有详尽的报道;另外,对1973年7月14日     

Large eddy simulation of hydrocyclone—prediction of air-core diameter and shape

Hydrocyclones are widely used in the mining and chemical industries. An attempt has been made in this study, to develop a CFD (computational fluid dynamics) model, which is capable of predicting the flow patterns inside the hydrocyclone, including accurate prediction of flow split as well as the size of the air-core. The flow velocities and air-core diameters are predicted by DRSM (differential Reynolds stress model) and LES (large eddy simulations) models were compared to experimental results. The predicted water splits and air-core diameter with LES and RSM turbulence models along with VOF (volume of fluid) model for the air phase, through the outlets for various inlet pressures were also analyzed. The LES turbulence model led to an improved turbulence field prediction and thereby to more accurate prediction of pressure and velocity fields. This improvement was distinctive for the axial profile of pressure, indicating that air-core development is principally a transport effect rather than a pressure effect.     

Increasing critical sensitivity of the Load/Unload Response Ratio before large earthquakes with identified stress accumulation pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.     

Variability in Early Amazonian Tharsis stress state based on wrinkle ridges and strike-slip faulting

We present evidence for a decrease in the magnitude of Tharsis-circumferential compressive stress during the Late Hesperian to the Middle Amazonian based on chronologic changes in the predominant style of faulting in southern Amazonis Planitia. Using high-resolution MOLA topography, we identify a population of strike-slip faults that exhibit Middle Amazonian-aged displacements of regional chrono-stratigraphic units. These strike-slip faults are adjacent to an older population of previously documented Late Hesperian-aged thrust faults (wrinkle ridges). Along-strike orientations of these thrust and strike-slip faults reveal the Tharsis-radial stress to be the area's most compressive remote principal stress and that this stress orientation and magnitude persisted throughout the Late Hesperian to the Middle Amazonian. We show that the change in the predominant style of faulting from thrust faulting to strike-slip faulting during this time requires a decrease of the Tharsis-circumferential compressive stress to a magnitude less than lithostatic load, with negligible change in stress orientation.     

Tectonic relation between northeastern China and the Korean peninsula revealed by interpretation of GRACE satellite gravity data

The major continental blocks in northeastern Asia are the North China block and the South China block, which have collided starting from the Korean peninsula. Geologic and geophysical interpretations reveal a well defined suture zone in northeastern China from Qinling through Dabie to Jiaodong. The discovery of high-pressure metamorphic rocks in the Hongseong area of the Korean peninsula, prominent evidence for the collision zone, indicates extension of the collision zone in northeastern China into the Korean peninsula. Interpretation of the GRACE satellite gravity dataset shows two prominent structural boundaries in the Yellow Sea. One extends from the Jiaodong Belt in eastern China to the Imjingang Belt in the Korean peninsula. The other extends from near Nanjing, eastern China, to Hongseong. Tectonic movement in or near the suture zone may be responsible for seismic activity in the western Korean peninsula and the development of the Yellow Sea sedimentary basin.     

Palaeozoic and Mesozoic tectonic evolution and palaeogeography of East Asian crustal fragments: The Korean Peninsula in context

East and Southeast Asia comprises a complex assembly of allochthonous continental lithospheric crustal fragments (terranes) together with volcanic arcs, and other terranes of oceanic and accretionary complex origins located at the zone of convergence between the Eurasian, Indo-Australian and Pacific Plates. The former wide separation of Asian terranes is indicated by contrasting faunas and floras developed on adjacent terranes due to their prior geographic separation, different palaeoclimates, and biogeographic isolation. The boundaries between Asian terranes are marked by major geological discontinuities (suture zones) that represent former ocean basins that once separated them. In some cases, the ocean basins have been completely destroyed, and terrane boundaries are marked by major fault zones. In other cases, remnants of the ocean basins and of subduction/accretion complexes remain and provide valuable information on the tectonic history of the terranes, the oceans that once separated them, and timings of amalgamation and accretion. The various allochthonous crustal fragments of East Asia have been brought into close juxtaposition by geological convergent plate tectonic processes. The Gondwana-derived East Asia crustal fragments successively rifted and separated from the margin of eastern Gondwana as three elongate continental slivers in the Devonian, Early Permian and Late Triassic–Late Jurassic. As these three continental slivers separated from Gondwana, three successive ocean basins, the Palaeo-Tethys,. Meso-Tethys and Ceno-Tethys, opened between these and Gondwana. Asian terranes progressively sutured to one another during the Palaeozoic to Cenozoic. South China and Indochina probably amalgamated in the Early Carboniferous but alternative scenarios with collision in the Permo–Triassic have been suggested. The Tarim terrane accreted to Eurasia in the Early Permian. The Sibumasu and Qiangtang terranes collided and sutured with Simao/Indochina/East Malaya in the Early–Middle Triassic and the West Sumatra terrane was transported westwards to a position outboard of Sibumasu during this collisional process. The Permo–Triassic also saw the progressive collision between South and North China (with possible extension of this collision being recognised in the Korean Peninsula) culminating in the Late Triassic. North China did not finally weld to Asia until the Late Jurassic. The Lhasa and West Burma terranes accreted to Eurasia in the Late Jurassic–Early Cretaceous and proto East and Southeast Asia had formed. Palaeogeographic reconstructions illustrating the evolution and assembly of Asian crustal fragments during the Phanerozoic are presented.     

Mesozoic transtensional basin history of the Eastern Cordillera, Colombian Andes: Inferences from tectonic models

Backstripping analysis and forward modeling of 162 stratigraphic columns and wells of the Eastern Cordillera (EC), Llanos, and Magdalena Valley shows the Mesozoic Colombian Basin is marked by five lithosphere stretching pulses. Three stretching events are suggested during the Triassic–Jurassic, but additional biostratigraphical data are needed to identify them precisely. The spatial distribution of lithosphere stretching values suggests that small, narrow (<150 km), asymmetric graben basins were located on opposite sides of the paleo-Magdalena–La Salina fault system, which probably was active as a master transtensional or strike-slip fault system. Paleomagnetic data suggesting a significant (at least 10°) northward translation of terranes west of the Bucaramanga fault during the Early Jurassic, and the similarity between the early Mesozoic stratigraphy and tectonic setting of the Payandé terrane with the Late Permian transtensional rift of the Eastern Cordillera of Peru and Bolivia indicate that the areas were adjacent in early Mesozoic times. New geochronological, petrological, stratigraphic, and structural research is necessary to test this hypothesis, including additional paleomagnetic investigations to determine the paleolatitudinal position of the Central Cordillera and adjacent tectonic terranes during the Triassic–Jurassic. Two stretching events are suggested for the Cretaceous: Berriasian–Hauterivian (144–127 Ma) and Aptian–Albian (121–102 Ma). During the Early Cretaceous, marine facies accumulated on an extensional basin system. Shallow-marine sedimentation ended at the end of the Cretaceous due to the accretion of oceanic terranes of the Western Cordillera. In Berriasian–Hauterivian subsidence curves, isopach maps and paleomagnetic data imply a (>180 km) wide, asymmetrical, transtensional half-rift basin existed, divided by the Santander Floresta horst or high. The location of small mafic intrusions coincides with areas of thin crust (crustal stretching factors >1.4) and maximum stretching of the subcrustal lithosphere. During the Aptian–early Albian, the basin extended toward the south in the Upper Magdalena Valley. Differences between crustal and subcrustal stretching values suggest some lowermost crustal decoupling between the crust and subcrustal lithosphere or that increased thermal thinning affected the mantle lithosphere. Late Cretaceous subsidence was mainly driven by lithospheric cooling, water loading, and horizontal compressional stresses generated by collision of oceanic terranes in western Colombia. Triassic transtensional basins were narrow and increased in width during the Triassic and Jurassic. Cretaceous transtensional basins were wider than Triassic–Jurassic basins. During the Mesozoic, the strike-slip component gradually decreased at the expense of the increase of the extensional component, as suggested by paleomagnetic data and lithosphere stretching values. During the Berriasian–Hauterivian, the eastern side of the extensional basin may have developed by reactivation of an older Paleozoic rift system associated with the Guaicáramo fault system. The western side probably developed through reactivation of an earlier normal fault system developed during Triassic–Jurassic transtension. Alternatively, the eastern and western margins of the graben may have developed along older strike-slip faults, which were the boundaries of the accretion of terranes west of the Guaicáramo fault during the Late Triassic and Jurassic. The increasing width of the graben system likely was the result of progressive tensional reactivation of preexisting upper crustal weakness zones. Lateral changes in Mesozoic sediment thickness suggest the reverse or thrust faults that now define the eastern and western borders of the EC were originally normal faults with a strike-slip component that inverted during the Cenozoic Andean orogeny. Thus, the Guaicáramo, La Salina, Bitúima, Magdalena, and Boyacá originally were transtensional faults. Their oblique orientation relative to the Mesozoic magmatic arc of the Central Cordillera may be the result of oblique slip extension during the Cretaceous or inherited from the pre-Mesozoic structural grains. However, not all Mesozoic transtensional faults were inverted.