南极GPS数据质量分析

以多路径误差和电离层延迟为评价指标,利用TEQC软件和QCVIEW软件对南极埃默里冰架、内陆冰盖、纳拉峡海湾、拉斯曼丘陵的实测GPS数据以及部分IGS站点的数据进行了数据质量分析,得到结论： 1)在同一站点上,L2上的多路径误差总比L1上的多路径误差要大;2)在南极不同环境下的多路径误差不同;3)纬度和太阳活动剧烈程度是影响南极地区GPS数据电离层延迟大小的两个重要因素。     

Persistent organic pollutants and histological lesions in Mayan catfish Ariopsis assimilis from the Bay of Chetumal, Mexico

Livers of catfish (Ariopsis assimilis) from the Bay of Chetumal were analyzed for organochlorine compounds and hydrocarbons as part of a study to diagnose the environmental health of the Bay after a catfish mass mortality that occurred in 1996. The presence of histological lesions in several organs of the fish as result of chemical exposure was also evaluated. The concentrations of organic pollutants found in the Bay may be considered high if compared to the levels reported for sites affected by chemical pollution. High prevalences of cellular alteration histopathologies were found in liver, including hepatic tumors. The presence of some lesions may be related statistically to environmental pollution in the Bay, specially with chlorinated compounds.     

The gas/steam ratio as indicator of heat transfer at the Solfatara fumaroles,Phlegraean Fields (Italy)

A simple geochemical model of Solfatara, Phlegraean Fields (Italy), is proposed on the basis of gas composition and temperature at the surface.Data on the Solfatara fumaroles have been collected since 1979 within the framework of a geochemical monitoring for the surveillance of the Phlegraean volcanic system.Surface manifestations of Solfatara are likely to be fed through isoenthalpic expansion of dry steam, which separates from a geothermal liquid in an intensively fractured zone at about 236°C. This value is consistent both with gas composition and surface temperature.The gas/steam ratio appears to be the most effective parameter to detect changes of heat flow at depth.Actually a remarkable decrease in the gas/steam ratio has been observed since 1981, while the gas composition and the temperature did not change significantly. These facts suggest increased heat flow at depth.     

Saturn’s zonal wind profile in 2004-2009 from Cassini ISS images and its long-term variability

Five years of Cassini Imaging Science Subsystem images, from 2004 to 2009, are analyzed in this work to retrieve global zonal wind profiles of Saturn’s northern and southern hemispheres in the methane absorbing bands at 890 and 727 nm and in their respective adjacent continuum wavelengths of 939 and 752 nm. A complete view of Saturn’s global circulation, including the equator, at two pressure levels, in the tropopause (60 mbar to 250 mbar with the MT filters) and in the upper troposphere (from ∼350 mbar to ∼500 mbar with the CB filter set), is presented. Both zonal wind profiles (available at the Supplementary Material Section), show the same structure but with significant differences in the peak of the eastward jets and the equatorial region, including a region of positive vertical shear symmetrically located around the equator between the 10° < |φ_c| < 25° where zonal velocities close to the tropopause are higher than at 500 mbar. A comparison of previously published zonal wind sets obtained by Voyager 1 and 2 (1980-1981), Hubble Space Telescope, and ground-based telescopes (1990-2004) with the present Cassini profiles (2004-2009) covering a full Saturn year shows that the shape of the zonal wind profile and intensity of the jets has remained almost unchanged except at the equator, despite the seasonal insolation cycle and the variability of Saturn’s emitted power. The major wind changes occurred at equatorial latitudes, perhaps following the Great White Spot eruption in 1990. It is not evident from our study if the seasonal insolation cycle and its associated ring shadowing influence the equatorial circulation at cloud level.     

Applying essentially non‐oscillatory interpolation to controlled‐source electromagnetic modelling*

Modelling and inversion of controlled‐source electromagnetic (CSEM) fields requires accurate interpolation of modelled results near strong resistivity contrasts. There, simple linear interpolation may produce large errors, whereas higher‐order interpolation may lead to oscillatory behaviour in the interpolated result. We propose to use the essentially non‐oscillatory, piecewise polynomial interpolation scheme designed for piecewise smooth functions that contains discontinuities in the function itself or in its first or higher derivatives. The scheme uses a non‐linear adaptive algorithm to select a set of interpolation points that represent the smoothest part of the function among the sets of neighbouring points. We present numerical examples to demonstrate the usefulness of the scheme. The first example shows that the essentially non‐oscillatory interpolation (ENO) scheme better captures an isolated discontinuity. In the second example, we consider the case of sampling the electric field computed by a finite‐volume CSEM code at a receiver location. In this example, the ENO interpolation performs quite well. However, the overall error is dominated by the discretization error. The other examples consider the comparison between sampling with essentially non‐oscillatory interpolation and existing interpolation schemes. In these examples, essentially non‐oscillatory interpolation provides more accurate results than standard interpolation, especially near discontinuities.     

Estimation of the effect of protective piers on wind wave parameters in the Sevastopol Bay

An effectiveness of the storm wave attenuation by protective piers in the Sevastopol Bay of the Black Sea is studied on the basis of numerical simulation using the SWAN spectral model. Analyzed are the parameters of waves generated by winds of four main directions as well as by the southern cyclone during the storm on November 11, 2007. It is obtained that waves from the northwest part of the Black Sea penetrate most intensively into the Sevastopol Bay in case of western wind and, to a lesser degree, in case of northern and southern winds. A protective effect of the piers is observed in the west part of the bay only and the wave attenuation near the southern coast is more significant than near the northern one. The area of the southern coast directly behind the southern pier is completely protected from the storm waves and, as moving away from the pier, the danger of intensive wave effect on the coast is kept.     

The biogeochemical sulfur cycle in the marine boundary layer over the Northeast Pacific Ocean

The major components of the marine boundary layer biogeochemical sulfur cycle were measured simultaneously onshore and off the coast of Washington State, U.S.A. during May 1987. Seawater dimethylsulfide (DMS) concentrations on the continental shelf were strongly influenced by coastal upwelling. Concentration further offshore were typical of summer values (2.2 nmol/L) at this latitude. Although seawater DMS concentrations were high on the biologically productive continental shelf (2–12 nmol/L), this region had no measurable effect on atmospheric DMS concentrations. Atmospheric DMS concentrations (0.1–12 nmol/m³), however, were extremely dependent upon wind speed and boundary layer height. Although there appeared to be an appreciable input of non-sea-salt sulfate to the marine boundary layer from the free troposphere, the local flux of DMS from the ocean to the atmosphere was sufficient to balance the remainder of the sulfur budget.     

The gravitational field of topographic-isostatic masses and the hypothesis of mass condensation II-the topographic-isostatic geoid

In Part I we focussed on a convergent representation of the gravitational potential generated bytopographic masses on top of the equipotential surface atMean Sea Level, thegeoid, and by those masses which compensate topography. Topographic masses have also been condensated, namely represented by a single layer. Part II extends the computation of the gravitational field of topographic-isostatic masses by a detailed analysis of itsforce field in terms ofvector-spherical harmonic functions. In addition, the discontinuous mass-condensated topographic gravitational force vector (head force) is given. Once we identify theMoho discontinuity asone interface of isostatically compensated topographical masses, we have computed the topographic potential and the gravitational potential which is generated by isostatically compensated masses atMean Sea Level, the geoid, and illustrated by various figures of geoidal undulations. In comparison to a data oriented global geoid computation ofJ. Engels (1991) the conclusion can be made that the assumption of aconstant crustal mass density, the basic condition for isostatic modeling, does not apply. Insteaddensity variations in the crust, e.g. betweenoceanic and continental crust densities, have to be introduced in order to match the global real geoid and its topographic-isostatic model. The performed analysis documents that thestandard isostatic models based upon aconstant crustal density areunreal.