Electromagnetic 2.5-dimensional forward modelling with a boundary integral formulation

An effective and accurate technique for the numerical solution of 2-D electromagnetic scattering problems with 3-D sources is presented. This solution introduces a set of the usual boundary integral equations and uses a scalar Green's function. In this scalar version, the unknowns of the problem are the boundary values of the longitudinal fields and their normal derivatives in the Fourier domain. A generalization of the usual boundary integral formulation enables us to handle a large class of models composed of piecewise homogeneous domains, including contiguous domains, multiply-connected domains and unbounded domains. This formulation involves the solution of a system of linear equations, and results in a significant saving in computation time in comparison with other rigorous methods.
  The requirements for the numerical implementation of this solution are described in detail. Numerical tests were carried out using the important example of electromagnetic tomography. The specific symmetry properties of the response function in this case are illustrated. Numerical accuracy is verified over a large frequency range, up to 1  MHz.     

西南太平洋地区地震序列的特征

共搜集到１９８４￣１９９０年西南太平洋地区１２个板缘地震序列。多数地震序列的特征是：震中分布区域的长轴较长并且随主震震级和序列中强震次数而增加;震中分布区域的长、短轴长度的比值较高;地震序列的余震震源机制和主震的差异不大;震源深度下限超过地壳,可达７０ｋｍ以上。走滑型主震占的比例低,高倾角滑动面的走向既有与俯冲带走向平行的也有横切的,个别逆冲型地震的断层面走向横切俯冲带。它们显示出与板块俯冲带主体     

CUTLASS Finland radar observations of the ionospheric signatures of flux transfer events and the resulting plasma flows

The CUTLASS Finland radar has been run in a two-beam special scan mode, which offered excellent temporal and spatial information on the flows in the high-latitude ionosphere. A detailed study of one day of this data revealed a convection reversal boundary (CRB) in the CUTLASS field of view (f.o.v) on the dayside, the direction of plasma flow either side of the boundary being typical of a dawn-cell convection pattern. Poleward of the CRB a number of pulsed transients are observed, seemingly moving away from the radar. These transients are identified here as the ionospheric signature of flux transfer events (FTEs). Equatorward of the CRB continuous backscatter was observed, believed to be due to the return flow on closed field lines. The two-beam scan offered a new and innovative opportunity to determine the size and velocity of the ionospheric signatures associated with flux transfer events and the related plasma flow pattern. The transient signature was found to have an azimuthal extent of 1900 ± 900 km and an poleward extent of 250 km. The motion of the transient features was in a predominantly westward azimuthal direction, at a velocity of 7.5 ± 3 km.     

Estimating open boundary conditions from coastal tidal observations by adjoint approach

A numerical study for estimating the tidal open boundary conditions of a shelf current modrl from tb coastal tidal observations is presented. The method is based on the optimal control/adjoint method. A lrast square fitting of the model state to simulated data is used. Two ideal domains and coastlines are considered. Using the IAP shallow. water model and its adjoint model, some identical twin experiments are carried out to test efficiency and lirnilsd of the method. The results show that the adjoint method can efficiently estimate the open boundary conditions well for gulf/bay like domains. The adjoint method seems to have great potential to improve the accuracy of tide and shelf current modeling in coastal regions. Project supported hy the National Natural Science Fuundation of China (Grant No. 49376256)     

The introduction of neural network system and its applications in rock engineering

The neural network system has been developing very fast recently. It has been widely used in many industries such as automation, nuclear power plant, chemical industry, etc. Neural network systems have a great advantage in dealing with problems in which many factors influence the process and result, and the understanding of this process is poor, and there are experimental data or field data. In rock engineering, many problems are of this nature. In this paper, a brief introduction to neural network systems is given. Problems such as what is a neural network, how it works and what kind of advantages it has are discussed. After this, several applications in rock engineering, made by us, are presented. Case 1 is ore boundary delineation. In this case, the rock are divided into three classes, i.e.: (1) waste rock; (2) semi-ore; and (3) ore for mining purposes. The neural network system built can judge whether it is ore, semi-ore or waste rock along the borehole according its corresponding geophysical logging data, such as gamma-ray, gamma-gamma, neutron and resistivity. Case 2 is aggregate quality prediction. In this case, the quality parameters: (1) impact value; (2) abrasion value I; and (3) abrasion value II are predicted by using a neural network system based on density, point load, content of quarts and content of brittle minerals. Case 3 is rock indentation depth prediction. In this case, the rock indentation depth under indentation load is predicted by the established neural network system based on the indentation load on rock, indenter type and rock mechanical properties, such as uniaxial compressive strength, Young's modulus. Poisson's ratio, critical energy release rate and density. In all these cases, the neural network systems have been applied successfully. The testing results are satisfactory and better than the existing techniques.     

Explicit solution of a free boundary problem for a nonlinear absorption model of mixed saturated-unsaturated flow

In wet soils, zones of saturation naturally develop in the vicinity of impermeable strata, surface ponds and subterranean cavities. Hydrology must be then concerned with transient flow through coexisting unsaturated and saturated zones. The models of advancing saturated zones necessarily involve a nonlinear free boundary problem.A closed-form analytic solution is presented for a nonlinear diffusion model under conditions of ponding at the surface. The soil water diffusivity is restricted to the special functional form D(θ) = a/(b − θ)², where θ is the water content field to be determined and a, b are positive constants. The explicit solution depends on a parameter C (determined by the data of the problem), according to two cases: 1 < C < C₁ or C ≥ C₁, where C₁ is a constant which is obtained as the unique solution of an equation. This result complements the study given in P. Broadbridge, Water Resources Research, 1990, 26, 2435–2443, in order to established when the explicit solution is available. The behavior of the bifurcation parameter C₁ as a function of the driving potential is studied with the corresponding limits for small and large values. Moreover, the sorptivity is proven to be continuously differentiable function of the variable C.     

Biostratigraphy and depositional facies of the Cretaceous-Tertiary boundary strata in Amakusa-Shimojima,Kyushu, Western Japan

The Upper Himenoura Subgroup exposed in the island of Amakusa-Shimojima, Kyushu, Japan shows an example of the terminal Cretaceous stratigraphic record in the circum Pacific region. This sequence is a part of the Upper Cretaceous intra-arc basins of southwest Japan. Four cycles of upward coarse-graded facies are recognized. Each cycle consists of a basinal mud facies in the lower part and a tide-dominated shallow marine to brackish coarse clastic facies in the upper part. Biostratigraphic correlation chiefly based on ammonites, inocerami and trigoniids indicates that this sequence is Campanian to Maastrichtian in age. The occurrence of the above three fossils decreases upward and is terminated at the top of the sequence, being replaced by a molluscan assemblage similar to the Danian. This suggests that the sedimentation may have continued to the very end of the Cretaceous period and possibly to the beginning of the Tertiary.     

Relative sea-level, tectono-eustasy, geoidal-eustasy and geodynamics during the Cretaceous

A priori, the recorded relative sea-level changes during the Cretaceous must be the combined effect of tectono-eustasy, geoidal-eustasy and various crustal level changes. To this we must add the human factor of differences and errors in interpretations.A posteriori, it is claimed that geoidal-eustasy dominated during the Hauterivian, Barremian, Turonian, Santonian and Maastrichtian, that tectono-eustasy dominated during the Albian, Cenomanian, Campanian and at the Maastrichtian/Danian boundary, and that local influences of sea-floor spreading are identified from the Albian/Cenomanian boundary onwards. To this we must add the local differential crustal movements modulating the global and regional ocean level changes. Geoidal-eustasy is mainly expressed as a latitudinal differentiation of the sea-level with out-of-phase changes between the hemispheres or the both high latitude regions. Furthermore, sedimentological records seem to record short-period geoidal-eustatic cycles.     

A study of a fair weather boundary layer in TOGA-COARE: Parameterization of surface fluxes in large scale and regional models for light wind conditions

The study focuses on a way to parameterize the effect of subgrid scale convective motions on surface fluxes in large scale and regional models for the case of light surface winds. As previously proposed, these subgrid effects are assumed to scale with the convection intensity through the relationship: where is the mean velocity of the wind, U₀ the velocity of the mean wind, w_* the free convection velocity, and an empirical coefficient to be determined. Both observations and numerical simulation are presently used to determine the free convection coefficient .Large eddy simulation of a fair weather convective boundary layer case observed during TOGA-COARE is performed. Comparisons between observations and the simulation of surface properties and vertical profiles in the planetary boundary layer are presented. The simulated vertical turbulent fluxes of heat, moisture and buoyancy range well within estimates from aircraft measurements.The most important result is that the true free convection coefficient , directly estimated from simulation, leads to a value of 0.65, smaller than the ones estimated from temporal and spatial variances. Using observations and simulation, estimates of from temporal and spatial variances are obtained with similar values 0.8. From both theoretical derivations and numerical computations, it is shown that estimates of the true from variances are possible but only after applying a correction factor equal to 0.8. If this correction is not used, is overestimated by about 25%. The time and space sampling problem is also addressed in using numerical simulations.