Regional variations in upper mantle structure beneath North America

Several types of seismological data, including surface wave group and phase velocities, travel times from large explosions, and teleseismic travel time anomalies, have indicated that there are significant regional variations in the upper few hundred kilometers of the mantle beneath continental areas. Body wave travel times and amplitudes from large chemical and nuclear explosions are used in this study to delineate the details of these variations beneath North America.

As a preliminary step in this study, theoretical P wave travel times, apparent velocities, and amplitudes have been calculated for a number of proposed upper mantle models, those of Gutenberg, Jeffreys, Lehman, and Lukk and Nersesov. These quantities have been calculated for both P and S waves for model CIT11GB, which is derived from surface wave dispersion data. First arrival times for all the models except that of Lukk and Nersesov are in close agreement, but the travel time curves for later arrivals are both qualitatively and quantitatively very different. For model CIT11GB, there are two large, overlapping regions of triplication of the travel time curve, produced by regions of rapid velocity increase near depths of 400 and 600 km. Throughout the distance range from 10 to 40 degrees, the later arrivals produced by these discontinuities have larger amplitudes than the first arrivals. The amplitudes of body waves, in fact, are extremely sensitive to small variations in the velocity structure, and provide a powerful tool for studying structural details.

Most of eastern North America, including the Canadian Shield has a Pn velocity of about 8.1 km/sec, with a nearly abrupt increase in compressional velocity by ~ 0.3 km/sec near at a depth varying regionally between 60 and 90 km. Variations in the structure of this part of the mantle are significant even within the Canadian Shield. The low-velocity zone is a minor feature in eastern North America and is subject to pronounced regional variations. It is 30 to 50 km thick, and occurs somewhere in the depth range from 80 to 160 km. The velocity decrease is less than 0.2 km/sec.

Consideration of the absolute amplitudes indicates that the attenuation due to anelasticity is negligible for 2 hz waves in the upper 200 km along the southeastern and southwestern margins of the Canadian Shield. For compressional waves the average Q for this region is > 3000. The amplitudes also indicate that the velocity gradient is at least 2 x 10^-3 both above and below the low-velocity zone, implying that the temperature gradient is < 4.8°C/km if the regions are chemically homogeneous.

In western North America, the low-velocity zone is a pronounced feature, extending to the base of the crust and having minimum velocities of 7.7 to 7.8 km/sec. Beneath the Colorado Plateau and Southern Rocky Mountains provinces, there is a rapid velocity increase of about 0.3 km/sec, similar to that observed in eastern North America, but near a depth of 100 km.

Complicated travel time curves observed on profiles with stations in both eastern and western North America can be explained in detail by a model taking into account the lateral variations in the structure of the low-velocity zone. These variations involve primarily the velocity within the zone and the depth to the top of the zone; the depth to the bottom is, for both regions, between 140 and 160 km.

The depth to the transition zone near 400 km also varies regionally, by about 30-40 km. These differences imply variations of 250 °C in the temperature or 6 % in the iron content of the mantle, if the phase transformation of olivine to the spinel structure is assumed responsible. The structural variations at this depth are not correlated with those at shallower depths, and follow no obvious simple pattern.

The computer programs used in this study are described in the Appendices. The program TTINV (Appendix IV) fits spherically symmetric earth models to observed travel time data. The method, described in Appendix III, resembles conventional least-square fitting, using partial derivatives of the travel time with respect to the model parameters to perturb an initial model. The usual ill-conditioned nature of least-squares techniques is avoided by a technique which minimizes both the travel time residuals and the model perturbations.

Spherically symmetric earth models, however, have been found inadequate to explain most of the observed travel times in this study. TVT4, a computer program that performs ray theory calculations for a laterally inhomogeneous earth model, is described in Appendix II. Appendix I gives a derivation of seismic ray theory for an arbitrarily inhomogeneous earth model.

Using Earth Deformation Caused by Surface Mass Loading to Constrain the Elastic Structure of the Crust and Mantle

Surface mass loads come in many different varieties, including the oceans, atmosphere, rivers, lakes, glaciers, ice caps, and snow fields. The loads migrate over Earth's surface on time scales that range from less than a day to many thousand years. The weights of the shifting loads exert normal forces on Earth's surface. Since the Earth is not perfectly rigid, the applied pressure deforms the shape of the solid Earth in a manner controlled by the material properties of Earth's interior. One of the most prominent types of surface mass loading, ocean tidal loading (OTL), comes from the periodic rise and fall in sea-surface height due to the gravitational influence of celestial objects, such as the moon and sun. Depending on geographic location, the surface displacements induced by OTL typically range from millimeters to several centimeters in amplitude, which may be inferred from Global Navigation and Satellite System (GNSS) measurements with sub-millimeter precision. Spatiotemporal characteristics of observed OTL-induced surface displacements may therefore be exploited to probe Earth structure. In this thesis, I present descriptions of contemporary observational and modeling techniques used to explore Earth's deformation response to OTL and other varieties of surface mass loading. With the aim to extract information about Earth's density and elastic structure from observations of the response to OTL, I investigate the sensitivity of OTL-induced surface displacements to perturbations in the material structure. As a case study, I compute and compare the observed and predicted OTL-induced surface displacements for a network of GNSS receivers across South America. The residuals in three distinct and dominant tidal bands are sub-millimeter in amplitude, indicating that modern ocean-tide and elastic-Earth models well predict the observed displacement response in that region. Nevertheless, the sub-millimeter residuals exhibit regional spatial coherency that cannot be explained entirely by random observational uncertainties and that suggests deficiencies in the forward-model assumptions. In particular, the discrepancies may reveal sensitivities to deviations from spherically symmetric, non-rotating, elastic, and isotropic (SNREI) Earth structure due to the presence of the South American craton.

Otolith morphometrics and population structure of Pacific sardine (Sardinops sagax) along the west coast of North America

The broad distribution of Pacific sardine (Sardinops sagax) along the Pacif ic coast of North America makes it difficult for fisheries managers to identify regional stocks of this dominant small pelagic species. An investigation of morphometric characteristics of otoliths of Pacific sardine across most of their range revealed regional differences in populations. In a survey of over 2000 otoliths, all ages (with an emphasis on age-1 recruits) were compared. Principal components analysis, multivariate analysis of variance, and a novel method derived from regression and residuals calculations, termed perimeter-weight profiles (PWPs), revealed otolith similarities and differences. The results of the different approaches to statistical comparisons did not always agree. Sardine otoliths from Mexican waters were generally lighter and more lobate than those from U.S. and Canadian populations. Age-1 otoliths from northern California in 2006–07 tended to be heavier and smoother than those from other areas, including year-class cohorts from southern California. Comparisons of age-groups and year-classes of northern California otoliths with the use of the PWP models indicated signif icant trends in year-to-year patterns. In conjunction with other established indices of population structure, otolith PWPs are a useful tool for identifying local and regional stocks of Pacific sardine and may help distinguish populations of other fish species as well.

Effects of starvation on energy density of juvenile chum salmon (Oncorhynchus keta) captured in marine waters of Southeastern Alaska

We conducted laboratory starvation experiments on juvenile chum salmon (Oncorhynchus keta) captured in the neritic marine waters of northern Southeast Alaska in June and July 2003. Temporal changes in fish energy density (whole body energy content [WBEC], cal/g dry weight), percent moisture content, wet weight (g), length (mm), and size-related condition residuals were measured in the laboratory and were then compared to long-term field data. Laboratory water temperatures and salinities averaged 9°C and 32 psu in both months. Trends in response variables were similar for both experimental groups, although sampling intervals were limited in July because fewer fish were available (n= 54) than in June (n=101). Overall, for June (45-d experimental period, 9 intervals), WBEC, wet weight, and condition residuals decreased and percent moisture content increased, whereas fork length did not change. For July (20-d experimental period, 5 intervals), WBEC and condition residuals decreased, percent moisture content and fork length increased, and wet weight did not change. WBEC, percent moisture content, and condition residuals fell outside the norm of longterm data ranges within 10–15 days of starvation, and may be more useful than fork length and wet weight for detecting fish condition responses to suboptimal environments.

Análise das dosagens do fator de crescimento endotelial vascular no plasma e fluidos peritoneais de pacientes com câncer epitelial de ovário

A carcinogênese epitelial ovariana tem sido foco de estudos científicos em todo o mundo desenvolvido. A angiogênese tumoral ovariana é um processo multifatorial que resulta em vários produtos pró-angiogênicos. Entre eles, o fator de crescimento vascular endotelial (VEGF) é predominante. Os objetivos deste estudo foram relacionar as dosagens do VEGF dos fluidos peritoneais, do plasma periférico e do infundíbulo pélvico aos níveis de citorredução em pacientes operadas de adenocarcinoma epitelial de ovário (CEO); formular um modelo probabilístico de citorredução e utilizar estas dosagens para estimar a probabilidade do desfecho de citorredução. Além disto, foi criada uma nova variável chamada carga de VEGF. Pelo procedimento step-wise a citorredução foi melhor descrita pela carga de VEGF, mas faltou Normalidade aos resíduos, não sendo possível a adequação de um modelo matemático. Porém, a curva ROC, forneceu uma área sob a curva de 0,84, com sensibilidade de 71,4 % e especificidade variando de 69,5 a 73,9%. O ponto de corte ótimo foi 15,52 log de picograma de carga de VEGF. A odds-ratio (OR) calculada para citorredução ótima descrita pela carga de VEGF foi de 11 (IC= 2,59 ; 46,78). No grupo com estágio avançado (III e IV), a OR foi de 6 (IC= 1,15; 31,22). Apesar do pequeno número de casos, esta nova variável pode vir a ser um auxiliar na determinação de situações onde cirurgia citorredutora deixa de ser a pedra fundamental do tratamento primário do CEO e a indução quimioterápica passe a ter o principal papel na citorredução química antes da cirugia nestes casos.

Growth and survival of sockeye salmon (Oncorhynchus nerka) from Karluk Lake and River, Alaska, in relation to climatic and oceanic regimes and indices, 1922–2000

We examined whether the relationship between climate and salmon production was linked through the effect of climate on the growth of sockeye salmon (Oncorhynchus nerka) at sea. Smolt length and juvenile, immature, and maturing growth rates were estimated from increments on scales of adult sockeye salmon that returned to the Karluk River and Lake system on Kodiak Island, Alaska, over 77 years, 1924–2000. Survival was higher during the warm climate regimes and lower during the cool regime. Growth was not correlated with survival, as estimated from the residuals of the Ricker stock-recruitment model. Juvenile growth was correlated with an atmospheric forcing index and immature growth was correlated with the amount of coastal precipitation, but the magnitude of winter and spring coastal downwelling in the Gulf of Alaska and the Pacific Northwest atmospheric patterns that influence the directional bifurcation of the Pacific Current were not related to the growth of Karluk sockeye salmon. However, indices of sea surface temperature, coastal precipitation, and atmospheric circulation in the eastern North Pacific were correlated with the survival of Karluk sockeye salmon. Winter and spring precipitation and atmospheric circulation are possible processes linking survival to climate variation in Karluk sockeye salmon.

Age validation of Dover sole (Microstomus pacificus) by means of bomb radiocarbon

We used bomb radiocarbon (14C) in this age validation study of Dover sole (Microstomus pacificus). The otoliths of Dover sole, a commercially important fish in the North Pacific, are difficult to age and ages derived from the current break-andburn method were not previously validated. The otoliths used in this study were chosen on the basis of estimated birth year and for the ease of interpreting growth zone patterns. Otolith cores, material representing years 0 through 3, were isolated and analyzed for 14C. Additionally, a small number of otoliths with difficult-to-interpret growth patterns were analyzed for 14C to help determine age interpretation. The measured Dover sole 14C values in easier-to-interpret otoliths were compared with a 14C reference chronology for Pacific halibut (Hippoglossus stenolepis) in the North Pacific. We used an objective statistical analysis where sums of squared residuals between otolith 14C values of Dover sole and the reference chronology were examined. Our statistical analysis also included a procedure where the Dover sole 14C values were standardized to the reference chronology. These procedures allowed an evaluation of aging error. The 14C results indicated that the Dover sole age estimates from the easier-to-interpret otoliths with the break-and-burn method are accurate. This study validated Dover sole ages from 8 to 47 years.

Using an inverse-logistic model to describe growth increments of blacklip abalone (Haliotis rubra) in Tasmania

A new description of growth in blacklip abalone (Haliotis rubra) with the use of an inverse-logistic model is introduced. The inverse-logistic model avoids the disadvantageous assumptions of either rapid or slow growth for small and juvenile individuals implied by the von Bertalanffy and Gompertz growth models, respectively, and allows for indeterminate growth where necessary. An inverse-logistic model was used to estimate the expected mean growth increment for different black-lip abalone populations around southern Tasmania, Australia. Estimates of the time needed for abalone to grow from settlement until recruitment (at 138 mm shell length) into the fishery varied from eight to nine years. The variability of the residuals about the predicted mean growth increments was described with either a second inverse-logistic relationship (standard deviation vs. initial length) or by a power relationship (standard deviation vs. predicted growth increment). The inverse-logistic model can describe linear growth of small and juvenile abalone (as observed in Tasmania), as well as a spectrum of growth possibilities, from determinate to indeterminate growth (a spectrum that would lead to a spread of maximum lengths).

Sobrevida e fatores prognósticos em uma coorte hospitalar de pacientes com câncer de próstata localizado

O câncer de próstata é a neoplasia mais incidente entre os homens brasileiros. Atualmente, grande parte destes tumores é confinada à próstata no momento do diagnóstico. No entanto, muitos tumores clinicamente classificados como localizados não o são de fato, levando a indicações terapêuticas curativas não efetivas. Por outro lado, muitos pacientes com câncer sem significância clínica são tratados desnecessariamente em função da limitação prognóstica do estadiamento clínicos (pré-tratamento) de pacientes com diagnóstico histológico de adenocarcinoma de próstata localizado (estágios I e II), em coorte hospitalar composta por pacientes tratados no Instituto Nacional de Câncer, Rio de Janeiro, matriculados entre 1990 a 1999. As funções de sobrevida foram calculadas empregando-se o estimados de Kaplan-Meier tomando-se como início a data do diagnóstico histológico e como eventos os óbitos cuja causa básica foi o câncer de próstata. Para avaliação dos fatores prognósticos clínicos foram calculadas as hazard ratios (HR), com intervalos de confiança de 95%, seguindo-se o modelo de riscos proporcionais de Cox. Foram analisadas como fatores prognósticos independentes as variáveis: idade, cor, grau de instrução, data do primeiro tratamento, grau de diferenciação celular d o tumor primário biopsiado (Gleason), estadiamento clínico e PSA total pré-tratamento. O pressuposto dos riscos proporcionais foi avaliado pela análise dos resíduos de Schoenfeld e a influência de valores aberrantes pelos resíduos martingale e escore. Foram selecionados 258 pacientes pelos critérios de elegibilidade do estudo, dos quais 46 foram a óbito durante o período de seguimento. A sobrevida global foi de 88% em 5 anos e de 71% em 10 anos. Idade maior que 80 anos, classificação de Gleason maior que 6, PSA maior que 40ng/ml, estádio B2 e cor branca foram marcadores independentes de pior prognóstico. Fatores prognósticos clássicos na literatura foram úteis na estimativa do prognóstico nesta coorte hospitalar. Os resultados mostram que para pacientes diagnosticados em fases iniciais, os fatores sócio-econômico analisados, não influenciaram o prognóstico. Outros estudos devem ser conduzidos no país para investigar as diferenças no prognóstico em relação à etnia.

Mudanças de Estado e Multiplicadores Fiscais no Brasil entre 1999-2012

Este trabalho avalia o comportamento dos multiplicadores fiscais no Brasil entre 1999-2012. Para tanto, utiliza a metodologia desenvolvida por Sims, Waggoner e Zha (2008), que é um procedimento Bayesiano de estimação no qual os parâmetros do modelo mudam com alterações no estado da economia e os estados (regimes) seguem um processo de mudança de regime markoviano. Ou seja, foi estimado um modelo VAR Estrutural Bayesiano com mudança de regimes Markoviana (Markov Switching Structural Bayesian Vector Autoregression - MS-SBVAR). A base de dados é composta pelo consumo da administração pública, pela formação bruta de capital fixo da administração pública, pela carga tributária líquida e pelo Produto Interno Bruto (PIB), das três esferas do governo (federal, estadual, incluindo o Distrito Federal, e municipal). O software MATLAB/Dynare foi utilizado na estimação dos modelos e os resultados sugerem a ocorrência de 2 ou 3 regimes nos dois modelos que melhor se ajustaram aos dados. Os multiplicadores estimados apresentaram os sinais esperados e os diferentes tipos de multiplicadores fiscais calculados apresentaram valores maiores para a resposta do PIB a choques na formação bruta de capital fixo da administração pública que são eficazes, uma vez que possuem valores maiores do que um e impacto de longo prazo no PIB - quando comparado aos choques no consumo da administração pública, que possuem pouca persistência e são ineficazes (menores do que um), além de uma resposta negativa e persistente do PIB a choques na carga tributária líquida. Os resultados obtidos não indicam, ainda, multiplicadores fiscais maiores em regimes com maior variância nos resíduos do modelo.

Quantifying annual variation in catchability for commercial and research fishing

Biomass indices, from commercial catch per unit of effort (CPUE) or random trawl surveys, are commonly used in fisheries stock assessments. Uncertainty in such indices, often ex-pressed as a coefficient of variation (CV), has two components: observation error, and annual variation in catchability. Only the former can be estimated directly. As a result, the CVs used for these indices either ignore the annual-variation component or assume a value for it (often implicitly). Two types of data for New Zealand stocks were examined: 48 sets of residuals and catchability estimates from stock assessments using either CPUE or trawl survey indices; and biomass estimates from 17 time series of trawl surveys with between 4 and 25 species per time series. These data show clear evidence of significant annual variation in catchability. With the trawl survey data, catchability was detectably extreme for many species in about one year in six. The assessment data suggest that this annual variability typically has a CV of about 0.2. For commercial CPUE the variability is slightly less, and a typical total CV (including both components) of 0.15 to 0.2. This is much less than the values of 0.3 to 0.35 that have commonly been assumed in New Zealand. Some estimates of catchability are shown to be implausible.

La reproducción de la merluza en el Mar Argentino (Merlucciidae, Meluccius merluccius hubbsi) : 1. Descripción histológica del ciclo del ovario de merluza

The present paper deals on the histological description of the hake ovary made on the basis of gonad observations of 394 females during the period April 1966 March 1967. The material was obtained from weekly sampling of commercial catches carried out at the Institute of Marine Biology (Prov. Buenos Aires, Argentina). The anatomical and histological description of the standard ovary and the adopted terminology are given. The maturation process is divided into five periods, from ovocyte formation to yolked ovocyte formation, with its histological description. Ovary changes are analyzed on detail. The following conclusions were outlined: 1. Analysis demonstrated that although some specimens were totally spawned others, at the end of the spawning period, retaining a great number of ovocytes in different maturity stages. Therefore, postspawners have been classified as follows: Postspawned II : This stages is characterized by the empty ovarian structure, with ovocytes in stage II, which will remain in the resting phase untill next spawning season. Postspawned III and IV: Their main characteristics are: tissue destruction, bloody residuals and remaining ovocytes in stage III and IV, respectively. 2. Some transformations were found in ovaries of postspawned III and IV. They are classified as follows, according to its origin and structure: Developed from follicular cell membrane – a) Glandular formations, b) Epiteloid formations - Originated from remaining ovocytes, c) Ovocyte disintegration, d) Ovocyte with follicular cell infiltration. 3. All those structures derived from postpawners III and IV have a temporary character and will be reabsorbed. Their presence delay the recuperation of the organ and its reproductive functions. Consequently, the possibility of those structures acting as control mechanisms is suggested. 4. Transformations pointed out in paragraph Nº 2 prevent the possibility of consecutive spawning originated from the remaining ovocytes (II and IV). 5. No structures originated from postspawners III and IV were found during summer season. 6. Reproductive cycle of hake has been described monthly. It was observed that maturing ovaries predominate in summer (November-December).

Metrical dental analysis on golden monkey (Rhinopithecus roxellana)

Dental variation in the Chinese golden monkey (Rhinopithecus roxellana) is here evaluated by univariate, bivariate, and multivariate analyses. Allometric analyses indicate that canines and P3s are positively, but other dimensions negatively scaled to mandible and maxilla, and to body size. With the exception of the mesiodistal dimensions of I-1 and M-3, and the buccolingual dimension of Pq, mandibular dental variables show similar scaling relative to body size. Analysis of residuals shows that males have significantly larger canine, P-3 and buccolingual dimensions of the postcanine teeth (M-2 and M-3) than females. A significant difference in shape between the sexes is found in the buccolingual dimension of the upper teeth, but not in the mandible. Unlike the situation in some other species, Female golden monkeys do nor exhibit relatively larger postcanine teeth than males, in fact, the reverse is true, especially for M(2)s and M(3)s. The fact that most of the dental variables show low negative allometry to body size might be related a cold environment that has led to the development of larger body size with I-educed energy loss. When the raw data are examined by Discriminant Function Analysis the sexes are clearly distinguishable.

基于显式分形插值的高保真地震数据重建方法研究

Based on the fractal theories, contractive mapping principles as well as the fixed point theory, by means of affine transform, this dissertation develops a novel Explicit Fractal Interpolation Function（EFIF）which can be used to reconstruct the seismic data with high fidelity and precision. Spatial trace interpolation is one of the important issues in seismic data processing. Under the ideal circumstances, seismic data should be sampled with a uniform spatial coverage. However, practical constraints such as the complex surface conditions indicate that the sampling density may be sparse or for other reasons some traces may be lost. The wide spacing between receivers can result in sparse sampling along traverse lines, thus result in a spatial aliasing of short-wavelength features. Hence, the method of interpolation is of very importance. It not only needs to make the amplitude information obvious but the phase information, especially that of the point that the phase changes acutely. Many people put forward several interpolation methods, yet this dissertation focuses attention on a special class of fractal interpolation function, referred to as explicit fractal interpolation function to improve the accuracy of the interpolation reconstruction and to make the local information obvious. The traditional fractal interpolation method mainly based on the randomly Fractional Brown Motion (FBM) model, furthermore, the vertical scaling factor which plays a critical role in the implementation of fractal interpolation is assigned the same value during the whole interpolating process, so it can not make the local information obvious. In addition, the maximal defect of the traditional fractal interpolation method is that it cannot obtain the function values on each interpolating nodes, thereby it cannot analyze the node error quantitatively and cannot evaluate the feasibility of this method. Detailed discussions about the applications of fractal interpolation in seismology have not been given by the pioneers, let alone the interpolating processing of the single trace seismogram. On the basis of the previous work and fractal theory this dissertation discusses the fractal interpolation thoroughly and the stability of this special kind of interpolating function is discussed, at the same time the explicit presentation of the vertical scaling factor which controls the precision of the interpolation has been proposed. This novel method develops the traditional fractal interpolation method and converts the fractal interpolation with random algorithms into the interpolation with determined algorithms. The data structure of binary tree method has been applied during the process of interpolation, and it avoids the process of iteration that is inevitable in traditional fractal interpolation and improves the computation efficiency. To illustrate the validity of the novel method, this dissertation develops several theoretical models and synthesizes the common shot gathers and seismograms and reconstructs the traces that were erased from the initial section using the explicit fractal interpolation method. In order to compare the differences between the theoretical traces that were erased in the initial section and the resulting traces after reconstruction on waveform and amplitudes quantitatively, each missing traces are reconstructed and the residuals are analyzed. The numerical experiments demonstrate that the novel fractal interpolation method is not only applicable to reconstruct the seismograms with small offset but to the seismograms with large offset. The seismograms reconstructed by explicit fractal interpolation method resemble the original ones well. The waveform of the missing traces could be estimated very well and also the amplitudes of the interpolated traces are a good approximation of the original ones. The high precision and computational efficiency of the explicit fractal interpolation make it a useful tool to reconstruct the seismic data; it can not only make the local information obvious but preserve the overall characteristics of the object investigated. To illustrate the influence of the explicit fractal interpolation method to the accuracy of the imaging of the structure in the earth’s interior, this dissertation applies the method mentioned above to the reverse-time migration. The imaging sections obtained by using the fractal interpolated reflected data resemble the original ones very well. The numerical experiments demonstrate that even with the sparse sampling we can still obtain the high accurate imaging of the earth’s interior’s structure by means of the explicit fractal interpolation method. So we can obtain the imaging results of the earth’s interior with fine quality by using relatively small number of seismic stations. With the fractal interpolation method we will improve the efficiency and the accuracy of the reverse-time migration under economic conditions. To verify the application effect to real data of the method presented in this paper, we tested the method by using the real data provided by the Broadband Seismic Array Laboratory, IGGCAS. The results demonstrate that the accuracy of explicit fractal interpolation is still very high even with the real data with large epicenter and large offset. The amplitudes and the phase of the reconstructed station data resemble the original ones that were erased in the initial section very well. Altogether, the novel fractal interpolation function provides a new and useful tool to reconstruct the seismic data with high precision and efficiency, and presents an alternative to image the deep structure of the earth accurately.

太阳辐射对热层和电离层变化性的影响

Solar ultraviolet (UV) radiation at wavelengths less than 400 nm is an important source of energy for aeronomic processes throughout the solar system. Solar UV photons are absorbed in planetary atmospheres, as well as throughout the heliosphere, via photodissociation of molecules, photoionization of molecules and atoms, and photoexcitation toexcitation including resonance scattering. In this paper, the solar irradiances data measured by TIMED SEE, as well as the solar proxies such as F10.7 and Mg II, thermosphere neutral density of CHAMP measurements and topside ionospheric plasmas densities from DMSP, are used to analyze solar irradiance effects on the variabilities of the thermosphere and the ionosphere. First, thermosphere densities near 410 km altitude are analyzed for solar irradiance variability effects during the period 2002-2004. Correlations between the densities and the solar irradiances for different spectral lines and wavelength ranges reveal significantly different characteristics. The density correlates remarkably well with all the selected solar irradiances except the lower chromospheric O I (130.4 nm) emission. Among the chosen solar proxies, the Mg II core-to-wing ratio index, EUV (30-120 nm) and F10.7 show the highest correlations with the density for short-term (< ~27 days) variations. For both long- (> ~27 days) and short-term variations, linear correlation coefficients exhibit a decreasing trend from low latitudes towards high latitudes. The density variability can be effectively modeled (capturing 71% of the variance) using multiple solar irradiance indices, including F10.7, SEUV (the EUV 30-120 nm index), and SFUV (the FUV 120-193 nm index), in which a lag time of 1 day was used for both F10.7 and SEUV, and 5 days for SFUV. In our regression formulation SEUV has the largest contribution to the density variation (40%), with the F10.7 having the next largest contribution (32%) and SFUV accounting for the rest (28%). Furthermore, a pronounced period of about 27.2 days (mean period of the Sun's rotation) is present in both density and solar irradiance data of 2003 and 2004, and a pronounced period of about 54.4 days (doubled period of the solar rotation) is also revealed in 2004. However, soft X-ray and FUV irradiances did not present a pronounced 54.4 day period in 2004, in spite of their high correlation with the densities. The Ap index also shows 54-day periodicities in 2004, and magnetic activity, together with solar irradiance, affects the 54-day variation in density significantly. In addition, NRLMSISE00, DTM-2000 and JB2006 model predictions are compared with density measurements from CHAMP to assess their accuracy, and the results show that these models underestimate the response of the thermosphere to variations induced by solar rotation. Next, the equatorial topside ionospheric plasmas densities Ni are analyzed for solar irradiance variability effects during the period 2002-2005. Linear correlations between Ni and the solar irradiances for different wavelength ranges reveal significantly different characteristics. XUV (0-35 nm) and EUV (115-130 nm) show higher correlation with Ni for the long-term variations, whereas EUV (35-115 nm) show higher correlation for the short-term variations. Moreover, partial correlation analysis shows that the long-term variations of Ni are affected by both XUV (0-35 nm) and EUV (35-115 nm), whereas XUV (0-35 nm) play a more important role; the short-term variations of Ni are mostly affected by EUV (35-115 nm). Furthermore, a pronounced period of about 27 days is present in both Ni and solar irradiance data of 2003 and 2004, and a pronounced period of about 54 days is also revealed in 2004. Finally, prompted by previous studies that have suggested solar EUV radiation as a means of driving the semiannual variation, we investigate the intra-annual variation in thermosphere neutral density near 400 km during 2002-2005. The intra-annual variation, commonly referred to as the ‘semiannual variation’, is characterized by significant latitude structure, hemispheric asymmetries, and inter-annual variability. The magnitude of the maximum yearly difference, from the yearly minimum to the yearly maximum, varies by as much as 60% from year to year, and the phases of the minima and maxima also change by 20-40 days from year to year. Each annual harmonic of the intra-annual variation, namely, annual, semiannual, ter-annual and quatra-annual, exhibits a decreasing trend from 2002 through 2005 that is correlated with the decline in solar activity. In addition, some variations in these harmonics are correlated with geomagnetic activity, as represented by the daily mean value of Kp. Recent empirical models of the thermosphere are found to be deficient in capturing most of the latitude dependencies discovered in our data. In addition, the solar flux and geomagnetic activity proxies that we have employed do not capture some latitude and inter-annual variations detected in our data. It is possible that these variations are partly due to other effects, such as seasonal-latitudinal variations in turbopause altitude (and hence O/N2 composition) and ionosphere coupling processes that remain to be discovered in the context of influencing the intra-annual variations depicted here. Our results provide a new dataset to challenge and validate thermosphere-ionosphere general circulation models that seek to delineate the thermosphere intra-annual variation and to understand the various competing mechanisms that may contribute to its existence and variability. We furthermore suggest that the term “intra-annual” variation be adopted to describe the variability in thermosphere and ionosphere parameters that is well-captured through a superposition of annual, semiannual, ter-annual, and quatra-annual harmonic terms, and that “semiannual’ be used strictly in reference to a pure 6-monthly sinusoidal variation. Moreover, we propose the term “intra-seasonal” to refer to those shorter-term variations that arise as residuals from the above Fourier representation.