Evaluation of a combined wavelet and a combined principal component analysis classification system for BCG diagnostic problem

Heart disease is one of the main factor causing death in the developed countries. Over several decades, variety of electronic and computer technology have been developed to assist clinical practices for cardiac performance monitoring and heart disease diagnosis. Among these methods, Ballistocardiography (BCG) has an interesting feature that no electrodes are needed to be attached to the body during the measurement. Thus, it is provides a potential application to asses the patients heart condition in the home. In this paper, a comparison is made for two neural networks based BCG signal classification models. One system uses a principal component analysis (PCA) method, and the other a discrete wavelet transform, to reduce the input dimensionality. It is indicated that the combined wavelet transform and neural network has a more reliable performance than the combined PCA and neural network system. Moreover, the wavelet transform requires no prior knowledge of the statistical distribution of data samples and the computation complexity and training time are reduced.

Reconstruction of vertical thermal structure from several subsurface temperatures in the China Seas and adjacent waters

Empirical Orthogonal Function (EOF) analysis is used in this study to generate main eigenvector fields of historical temperature for the China Seas (here referring to Chinese marine territories) and adjacent waters from 1930 to 2002 (510 143 profiles). A good temperature profile is reconstructed based on several subsurface in situ temperature observations and the thermocline was estimated using the model. The results show that: 1) For the study area, the former four principal components can explain 95% of the overall variance, and the vertical distribution of temperature is most stable using the in situ temperature observations near the surface. 2) The model verifications based on the observed CTD data from the East China Sea (ECS), South China Sea (SCS) and the areas around Taiwan Island show that the reconstructed profiles have high correlation with the observed ones with the confidence level > 95%, especially to describe the characteristics of the thermocline well. The average errors between the reconstructed and observed profiles in these three areas are 0.69A degrees C, 0.52A degrees C and 1.18A degrees C respectively. It also shows the model RMS error is less than or close to the climatological error. The statistical model can be used to well estimate the temperature profile vertical structure. 3) Comparing the thermocline characteristics between the reconstructed and observed profiles, the results in the ECS show that the average absolute errors are 1.5m, 1.4 m and 0.17A degrees C/m, and the average relative errors are 24.7%, 8.9% and 22.6% for the upper, lower thermocline boundaries and the gradient, respectively. Although the relative errors are obvious, the absolute error is small. In the SCS, the average absolute errors are 4.1 m, 27.7 m and 0.007A degrees C/m, and the average relative errors are 16.1%, 16.8% and 9.5% for the upper, lower thermocline boundaries and the gradient, respectively. The average relative errors are all < 20%. Although the average absolute error of the lower thermocline boundary is considerable, but contrast to the spatial scale of average depth of the lower thermocline boundary (165 m), the average relative error is small (16.8%). Therefore the model can be used to well estimate the thermocline.

Two modes of dipole events in tropical Indian Ocean

By analyzing the distributions of subsurface temperature and the surface wind stress anomalies in the tropical Pacific and Indian Oceans during the Indian Ocean Dipole (IOD) events, two major modes of the IOD and their formation mechanisms are revealed. (1) The subsurface temperature anomaly (STA) in the tropical Indian Ocean during the IOD events can be described as a "<" -shaped and west-east-oriented dipole pattern; in the east side of the "<" pattern, a notable tongue-like STA extends westward along the equator in the tropical eastern Indian Ocean; while in the west side of the "<" pattern, the STA has opposite sign with two centers (the southern one is stronger than the northern one in intensity) being of rough symmetry about the equator in the tropical mid-western Indian Ocean. (2) The IOD events are composed of two modes, which have similar spatial pattern but different temporal variabilities due to the large scale air-sea interactions within two independent systems. The first mode of the IOD event originates from the air-sea interaction on a scale of the tropical Pacific-Indian Ocean and coexists with ENSO. The second mode originates from the air-sea interaction on a scale of the tropical Indian Ocean and is closely associated with changes in the position and intensity of the Mascarene high pressure. The strong IOD event occurs when the two modes are in phase, and the IOD event weakens or disappears when the two modes are out of phase. Besides, the IOD events are normally strong when either of the two modes is strong. (3) The IOD event is caused by the abnormal wind stress forcing over the tropical Indian Ocean, which results in vertical transports, leading to the upwelling and pileup of seawater. This is the main dynamic processes resulting in the STA. When the anomalous easterly exists over the equatorial Indian Ocean, the cold waters upwell in the tropical eastern Indian Ocean while the warm waters pileup in the tropical western Indian Ocean, hence the thermocline in the tropical Indian Ocean is shallowed in the east and deepened in the west. The off-equator component due to the Coriolis force in the equatorial area causes the upwelling of cold waters and the shallowing of the equatorial India Ocean thermocline. On the other hand, the anomalous anticyclonic circulations and their curl fields located on both sides of the equator, cause the pileup of warm waters in the central area of their curl fields and the deepening of the equatorial Indian Ocean thermocline off the equator. The above three factors lead to the occurrence of positive phase IOD events. When anomalous westerly dominates over the tropical Indian Ocean, the dynamic processes are reversed, and the negative-phase IOD event occurs.

Modeling vertical heat transport in the wave affected surface layer of the ocean

In considering the vertical heat transport problems in the upper ocean, the flat upper boundary approximation for the free surface and the horizontal homogenous hypothesis are usually applied. However, due to the existence of the wave motion, the application of this approximation may result in some errors to the solar irradiation since it decays quickly in respect to the actual thickness of the water layer below the surface; on the other hand, due to the fluctuation of the water layer depth, it is improper to neglect the effects of the horizontal advection and turbulent diffusion since they also contribute to the vertical heat transport. A new model is constructed in this study to reflect these effects. The corresponding numerical simulations show that the wave motion may remarkably accelerate the vertical heat transferring process and the variation of the temperature in the wave affected layer appears in an oscillating manner.

Ship-board measurements and estimations of air-sea fluxes in the western tropical Pacific during TOGA COARE

Direct air-sea flux measurements were made on RN Kexue #1 at 40 degrees S, 156 degrees E during the Tropical Ocean Global Atmosphere (TOGA) Coupled Ocean-Atmospheric Response Experiment (COARE) Intensive Observation Period (IOP). An array of six accelerometers was used to measure the motion of the anchored ship, and a sonic anemometer and Lyman-alpha hygrometer were used to measure the turbulent wind vector and specific humidity. The contamination of the turbulent wind components by ship motion was largely removed by an improvement of a procedure due to Shao based on the acceleration signals. The scheme of the wind correction for ship motion is briefly outlined. Results are presented from data for the best wind direction relative to the ship to minimize flow distortion effects. Both the time series and the power spectra of the sonic-measured wind components show swell-induced ship motion contamination, which is largely removed by the accelerometer correction scheme, There was less contamination in the longitudinal wind component than in the vertical and transverse components. The spectral characteristics of the surface-layer turbulence properties are compared with those from previous land and ocean results, Momentum and latent heat fluxes were calculated by eddy correlation and compared to those estimated by the inertial dissipation method and the TOGA COARE bulk formula. The estimations of wind stress determined by eddy correlation are smaller than those from the TOGA COARE bulk formula, especially for higher wind speeds, while those from the bulk formula and inertial dissipation technique are generally in agreement. The estimations of latent heal flux from the three different methods are in reasonable agreement. The effect of the correction for ship motion on latent heat fluxes is not as large as on momentum fluxes.

Vertical structure and evolution of the Luzon Warm Eddy

Eddies are frequently observed in the northeastern South China Sea (SCS). However, there have been few studies on vertical structure and temporal-spatial evolution of these eddies. We analyzed the seasonal Luzon Warm Eddy (LWE) based on Argo float data and the merged data products of satellite altimeters of Topex/Poseidon, Jason-1 and European Research Satellites. The analysis shows that the LWE extends vertically to more than 500 m water depth, with a higher temperature anomaly of 5A degrees C and lower salinity anomaly of 0.5 near the thermocline. The current speeds of the LWE are stronger in its uppermost 200 m, with a maximum speed of 0.6 m/s. Sometimes the LWE incorporates mixed waters from the Kuroshio Current and the SCS, and thus has higher thermohaline characteristics than local marine waters. Time series of eddy kinematic parameters show that the radii and shape of the LWE vary during propagation, and its eddy kinetic energy follows a normal distribution. In addition, we used the empirical orthogonal function (EOF) here to analyze seasonal characteristics of the LWE. The results suggest that the LWE generally forms in July, intensifies in August and September, separates from the coast of Luzon in October and propagates westward, and weakens in December and disappears in February. The LWE's westward migration is approximately along 19A degrees N latitude from northwest of Luzon to southeast of Hainan, with a mean speed of 6.6 cm/s.

Diel vertical migration of zooplankton following optimal food intake under predation

A model is developed to investigate the trade-offs between benefits and costs involved in zooplanktonic diel vertical migration (DVM) strategies. The 'venturous revenue' (VR) is used as the criterion for optimal trade-offs. It is a function of environmental factors and the age of zooplankter. During vertical migration, animals are assumed to check instantaneously the variations of environmental parameters and thereby select the optimal behavioral strategy to maximize the value of VR, i.e. taking up as much food as possible with a certain risk of mortality. The model is run on a diel time scale (24 h) in four possible scenarios during the animal's life history. The results show that zooplankton can perform normal DVM balancing optimal food intake against predation risk, with the profile of DVM largely modified by the age of zooplankter.

Control mechanisms of diel vertical migration: Theoretical assumptions

We explore control mechanisms underlying the vertical migration of zooplankton in the water column under the predator-avoidance hypothesis. Two groups of assumptions in which the organisms are assumed to migrate vertically in order to minimize realized or effective predation pressure (type-I) and to minimize changes in realized or effective predation pressure (type-II), respectively, are investigated. Realized predation pressure is defined as the product of light intensity and relative predation abundance and the part of realized predation pressure that really affects organisms is termed as effective predation pressure. Although both types of assumptions can lead to the migration of zooplankton to avoid the mortality from predators, only the mechanisms based on type-II assumptions permit zooplankton to undergo a normal diel vertical migration (morning descent and evening ascent). The assumption of minimizing changes in realized predation pressure is based on consideration of DVM induction only by light intensity and predators. The assumption of minimizing changes in effective predation pressure takes into account, apart from light and predators also the effects of food and temperature. The latter assumption results in the same expression of migration velocity as the former one when both food and temperature are constant over water depth. A significant characteristic of the two type-II assumptions is that the relative change in light intensity plays a primary role in determining the migration velocity. The photoresponse is modified by other environmental variables: predation pressure, food and temperature. Both light and predation pressure are necessary for organisms to undertake DVM. We analyse the effect of each single variable. The modification of the phototaxis of migratory organisms depends on the vertical distribution of these variables. (C) 2001 Academic Press.

Vertical transferring process of rare elements in coral reef lagoons of Nansha Islands, South China Sea

The vertical fluxes and vertical transferring forms of 18 rare elements were studied for the first time in the coral reef ecosystem of Nansha Islands, South China Sea, by deploying sediment traps, The results showed that the vertical transferring flux of most of the measured rare elements in Yongshu lagoon were higher than that in Zhubi lagoon. The vertical transferring forms of rare elements were mainly in the carbonate form, but Ta, As, Th mainly in the ion-exchange form, Ag in iron-manganese oxide form and Sb in the organic matter + sulphide form. None of the 18 rare elements was transferred mainly in the form of detritus silicate to sea floor. This proved that rare elements originating from the earth's crust were redistributed in sinking particulates after they were brought into ocean. The relation between the fluxes and surface seawater temperature (STT) was also studied. The sensitivity of rare elements to SST was in order: Rb>V>As>Ti>U>Zn>Sb>Hf>Ag>Cs.

Spatial variations of size-fractionated chlorophyll, cyanobacteria and heterotrophic bacteria in the Central and Western Pacific

Geographic and vertical variations of size-fractionated (0.2-1 mu m, 1-10 mu m, and >10 mu m) Chlorophyll a (Chl.a) concentration, cyanobacteria abundance and heterotrophic bacteria abundance were investigated at 13 stations from 4 degrees S, 160 degrees W to 30 degrees N, 140 degrees E in November 1993. The results indicated a geographic distribution pattern of these parameters with instances of high values occurring in the equatorial region and offshore areas, and with instance of low values occurring in the oligotrophic regions where nutrients were almost undetectable. Cyanobacteria showed the highest geographic variation (ranging from 27x10(3) to 16,582x10(3) cell l(-1)), followed by Chl.a (ranging from 0.048 to 0.178 mu g l(-1)), and heterotrophic bacteria (ranging from 2.84x10(3) to 6.50 x 10(5) cell l(-1)). Positive correlations were observed between nutrients and Chl.a abundance. Correspondences of cyanobacteria and heterotrophic bacteria abundances to nutrients were less significant than that of Chl.a. The total Chl.a was accounted for 1.0-30.9%, 35.9-53.7%, and 28.1-57.3% by the >10 mu m, 1-10 mu m and 0.2-1 mu m fractions respectively. Correlation between size-fractionated Chl.a and nutrients suggest that the larger the cell size, the more nutrient-dependent growth and production of the organism. The ratio of pheophytin to chlorophyll implys that more than half of the > 10 mu m and about one third of the 1-10 mu m pigment-containing particles in the oligotrophic region were non-living fragments, while most of the 1-10 mu m fraction was living cells. In the depth profiles, cyanobacteria were distributed mainly in the surface layer, whereas heterotrophic bacteria were abundant from surface to below the euphotic zone. Chl.a peaked at the surface layer (0-20 m) in the equatorial area and at the nitracline (75-100 m) in the oligotrophic regions. Cyanobacteria were not the principle component of the picoplankton. The carbon biomass ratio of heterotroph to phytoplankton was greater than 1 in the eutrophic area and lower than 1 in oligotrophic waters.

Growing season ecosystem respirations and associated component fluxes in two alpine meadows on the Tibetan Plateau

From 30 June to 24 September in 2003 ecosystem respiration (Re) in two alpine meadows on the Tibetan Plateau were measured using static chamber- and gas chromatography- (GC) based techniques. Simultaneously, plant removal treatments were set to partition Re into plant autotrophic respiration (Ra) and microbial heterotrophic respiration (Rh). Results indicated that Re had clear diurnal and seasonal variation patterns in both of the meadows. The seasonal variability of Re at both meadow sites was caused mainly by changes in Ra, rather than Rh. Moreover, at the Kobresia humilis meadow site (K_site), Ra and Rh accounted for 54% and 46% of Re, respectively. While at the Potentilla fruticosa scrub meadow (P_site), the counterparts accounted for 61% and 39%, respectively. T test showed that there was significant difference in Re rates between the two meadows (t = 2.387, P = 0.022). However, no significant difference was found in Rh rates, whereas a significant difference was observed in Ra rates between the two meadows. Thus, the difference in Re rate between the two meadows was mainly attributed to plant autotrophic respirations. During the growing season, the two meadows showed relatively low Q(10) values, suggesting that Re, especially Rh was not sensitive to temperature variation in the growing season. Additionally, Re and Rh at the K_site, as well as Rh at the P_site was negatively correlated with soil moisture, indicating that soil moisture would also play an important role in respirations.

地球介质非线性与各向异性的地震波响应及其成像效应

The practice of geophysical prospecting shows us the complex interior earth. The studies of the complexity play an important role and practical guide for the subsurface structure. At present, the complexity of the earth mainly means lateral and vertical homogeneity, anisotropy and non-linear quality. And the anisotropy and non-linear media studies become the frontier in seismology and exploration seismology. This paper summarizes the development of complexities and presents the forward and inverse in the non-linear and anisotropic media. Firstly, the paper introduces the theory of seismic wave propagation in the non-linear and anisotropic media, the theoretical basis for simulation and inversion research. Secondly, high quality numerical simulation method with little dispersion has been developed to investigate the influence of complexity including anisotropy and non-linear multi-component seismograms. Because most real data in seismology have a single component, we developed two aspects work on anisotropic multi-component imaging. One is prestack reflection migration. The result show that distorted images are obtained if data from anisotropic media are migrated using isotropic extrapolation. Moreover, image quality will be improved greatly after considering anisotropy in subsurface layers. The other one is the we take advantage of multi-component data to inversion of the anisotropic parameters jointly seimic reflection travel time and polarization information. Based on these research works, we get the following results: 1.Combing numerical simulation, systematical studies indicate that anisotropy and non-linear seismograms characters are significant to detect cracked belts in the earth and to understand deformation field and mechanism. 2.Based on anisotropic media models, we developed an efficient prestack migration method for subsurface structure and different observation methods seismic data, which improving the imaging quality with VSP, seismograms and real data. 3.Jointly seismic inversion combining seismic anisotropic reflection traveltimes and polarizations data show that the complete wrong inversion and the following explanation will be resulted by ignoring anisotropy.

多波多分量地震资料走时近似方法的研究与应用

Compared with the conventional P wave, multi-component seismic data can markedly provide more information, thus improve the quality of reservoir evaluation like formation evaluation etc. With PS wave, better imaging result can be obtained especially in areas involved with gas chimney and high velocity formation. However, the signal-to-noise of multi-component seismic data is normally lower than that of the conventional P wave seismic data, while the frequency range of converted wave is always close to that of the surface wave which adds to the difficulty of removing surface wave. To realize common reflection point data stacking from extracted common conversion point data is a hard nut to crack. The s wave static correction of common receiver point PS wave data is not easy neither. In a word, the processing of multi-component seismic data is more complicated than P wave data. This paper shows some work that has been done, addressing those problems mentioned above. (1) Based on the AVO feature of converted wave, this paper has realized the velocity spectrum of converted waves by using Sarkar’s generalized semblance method taking into account of AVO factor in velocity analysis. (2)We achieve a method of smoothly offset division normal method.Firstly we scan the stacking velocities in different offset divisions for a t0, secondly obtain some hyperbolas using these stacking velocities, then get the travel time for every trace using these hyperbolas; in the end we interpolate the normal move out between two t0 for every trace. (3) Here realize a method of stepwise offset division normal moveout.It is similar to the method of smoothly offset division normal moveout.The main difference is using quadratic curve, sixth order curve or fraction curve to fit these hyperbolas. (4)In this paper, 4 types of travel time versus distance functions in inhomogeneous media whose velocity or slowness varies with depth and vertical travel time have been discussed and used to approximate reflection travel time. The errors of ray path and travel time based on those functions in four layered models were analyzed, and it has shown that effective results of NMO in synthetic or real data can be obtained. (5) The velocity model of converted PS-wave can be considered as that of P -wave based on the ghost source theory, thus the converted wave travel time can be approximated by calculation from 4 equivalent velocity functions: velocity or slowness vary linearly with depth or vertical travel time. Then combining with P wave velocity analysis, the converted wave data can be corrected directly to the P-wave vertical travel time. The improvements were shown in Normal Move out of converted waves with numerical examples and real data. (6) This paper introduces the methods to compute conversion point location in vertical inhomogeneous media based on linear functions of velocity or slowness versus depth or vertical travel time, and introduce three ways to choose appropriate equivalent velocity methods, which are velocity fitting, travel time approximation and semblance coefficient methods.

尺度匹配方法在井震联合重建和多次波衰减中的应用

Scale matching method means adjusting information with different scale to the same level. This thesis focuses on scale unification of information with different frequency bandwidth. Well-seismic cooperate inversion is an important component of reservoir geophysics; multiple prediction & subtraction is a development of multiple attenuation in recent years. The common ground of these two methods is that they both related to different frequency bandwidth unification. Well log、cross-hole seismic、VSP、3D seismic and geological information have different spatial resolution, we can decrease multi-solution of reservoir inversion and enhance the vertical and lateral resolution of the geological object by integrate those information together; Compare the predicted multiple generated by SRME with the real multiple, we find the predicted multiple convolutes at least one wavelet more, which brings frequency bandwidth difference between them. So the subtraction method also relates to multi-scale information unification. This thesis gives a method of well constrained seismic high resolution processing basing on auto gain control modulation. It uses base function method which utilizes original well-seismic match result as initial condition and processed seismic trace as initial model to extrapolate the high frequency information of the well logs to the seismic profiles. In this way we can broaden the bandwidth of the seismic and make the high frequency gain geological meaning. In this thesis we introduce the revised base function method to adaptive subtraction and verify the validity of the method using models. Key words: high frequency reconstruction, scale matching, base function, multiple, SRME prediction & subtraction

电离层垂直TEC映射函数的实验观测与统计特性

When used in the determining the total electron content (TEC), which may be the most important ionospheric parameter, the worldwide GPS observation brings a revolutionary change in the ionospheric science. There are three steps in the data processing to retrieve GPS TEC: (1) to estimate slant TEC from the measurements of GPS signals; (2) to map the slant TEC into vertical; and (3) to interpolate the vertical TEC into grid points. In this scientific dissertation we focus our attention on the second step, the mapping theory and method to convert slant TEC into vertical. This is conventionally done by multiplying on the slant TEC a mapping function which is usually determined by certain models of electron density profile. Study of the vertical TEC mapping function is of significance in GPS TEC measurement. This paper first reviews briefly the three steps in GPS TEC mapping process. Then we compare the vertical TEC mapping function which were respectively calculated from the electron density profiles of the ionospheric model and retrieved from the observation of worldwide GPS TEC. We also perform the statistical analysis on the observational mapping functions. The main works and results are as follows: 1. We calculated the vertical TEC mapping functions for both SLM and Chapman models, and discussed the modulation of the ionosphere height to the mapping functions. We use two simple models, single layer model (SLM) and Chapman models, of the ionospheric electron density profiles to calculate the vertical TEC mapping function. In the case of the SLM, we discuss the control of the ionospheric altitude, i.e., the layer height hipp, to the mapping function. We find that the mapping function decreases rapidly as hipp increases. For the Chapman model we study also the control mapping function by both ionospheric altitude indicated by the peak electron density height hmF2, and the scale height, H, which present the thickness of the ionosphere. It is also found that the mapping function decreases rapidly as hmF2 increases. and it also decreases as H increases. 2. Then we estimate the mapping functions from the GPS observations and compare them with those calculated from the electron density models. We first, proposed a new method to estimate the mapping functions from GPS TEC data. This method is then used to retrieve the observational mapping function from both the slant TEC (TECS) provided by International GPS Service (IGS）and vertical TEC provide by JPL Global Ionospheric Maps (GIMs). Then we compare the observational mapping function with those calculated from the electron density models, SLM and Chapman. We find that the values of the observational mapping functions are much smaller than that from the model mapping functions, when the zenith angle is large enough. We attribute this to the effect of the plasmasphere which is above about 1000 km. 3. We statistically analyze the observational mapping functions and reveal their climatological changes. Observational mapping functions during 1999-2007 are used in our statistics. The main results are as follows. (1) The observational mapping functions decrease obviously with the decrement of the solar activity which is represented by the F10.7 index; (2) In annual variations of the observational mapping functions, the semiannual component is found at low-latitudes, and the remarkable seasonal variations at mid- and high-latitudes. (3) The diurnal variation of the observational mapping functions is that they are large in daytime and small at night, they become extremely small in the early morning before sunrise. (4) The observational mapping functions change with latitudes that they are smaller at lower latitudes and larger at higher. All of the above variations of the observational mapping functions are explained by the existence of the plasmasphere, which changes more slowly with time and more rapidly with latitude than the ionosphere does . In summary, our study on the vertical TEC mapping function imply that the ionosphere height has a modulative effect on the mapping function. We first propose the concept of the 'observational mapping functions' , and provide a new method to calculate them. This is important in improving the TEC mapping. It may also possible to retrieving the plasmaspheric information from GPS observations.