Single frequency processing of Ørsted GPS radio occultation measurements

The Global Positioning System (GPS) radio occultation measurements obtained using the TurboRogue GPS receiver on the Danish satellite Ørsted have been processed using the single frequency method. Atmospheric profiles of refractivity and temperature are derived and validated against numerical weather prediction data from the European Centre for Medium-Range Weather Forecast (ECMWF). Results from the Ørsted GPS measurement campaign in February 2000 indicate that the single frequency method can provide retrievals with accuracy comparable to that of using two frequencies. From comparisons between measured dry temperature profiles and corresponding dry temperature profiles derived from ECMWF analysis fields, we find a mean difference of less than 0.5 K and a standard deviation of 2–4 K between 500 and 30 hPa in height. Above 30 hPa the impact of the ionosphere becomes more dominant and more difficult to eliminate using the single frequency method, and the results show degraded accuracy when compared to previous analysis results of occultation data from other missions using the dual frequency method. At latitudes less than 40° (denoted low latitudes), the standard deviation is generally smaller than at latitudes higher than 40° (denoted high latitudes). A small temperature bias is observed centered at 200 hPa for low latitudes and at 300 hPa for high latitudes. This indicates that the ECMWF analyses do not adequately resolve the tropopause temperature minimum. In the lowest part of the troposphere an observed warm bias is thought to be due to erroneous tracking of the GPS signal in cases of atmospheric multipath propagation.     

Modelling Korean extreme rainfall using a Kappa distribution and maximum likelihood estimate

Summary Attempts to use the 4-parameter Kappa distribution (K4D) with the maximum likelihood estimates (MLE) on the summer extreme daily rainfall data at 61 gauging stations over South Korea have been made to obtain reliable quantile estimates for several return periods. A numerical algorithm for searching MLE of K4D by minimizing the negative log-likelihood function with penalty method has been described. The isopluvial maps of estimated design values corresponding to selected return periods have been presented. The highest return values are centered at sites in the south-western part of the Korean peninsula. The distribution of return values for annual maxima of 2-day precipitation (AMP2) is more similar to the climatological features of annual total precipitation of Korea than that of annual maxima of daily precipitation (AMP1). Our results of return values delineate well the horizontal patterns of the heavy precipitation over the Korean peninsula. Received January 15, 2001 Revised October 8, 2001     

Advances in lithium analysis in solids by means of laser-induced breakdown spectroscopy: an exploratory study

Lithium is an important geochemical tracer for fluids or solids. However, because the electron microprobe cannot detect Li, variations of Li abundance at the micrometric scale are most often estimated from bulk analyses. In this study, the Li intense emission line at 670.706 nm in optical emission spectroscopy was used to perfect the analysis of Li at the micrometric scale by means of laser-induced breakdown spectroscopy (LIBS). To estimate lithium content for different geological materials, LIBS calibration of the emission line at 670.706 nm was achieved by use of synthetic glasses and natural minerals. The detection limit for this method is ∼5 ppm Li. Three applications to geological materials show the potential of LIBS for lithium determination, namely for Li-bearing minerals, melt inclusions, quartz, and associated fluid inclusions.For spodumene and petalite from granite pegmatite dikes (Portugal), the Li₂O concentrations are 7.6 ± 1.6 wt% and 6.3 ± 1.3 wt%, respectively, by use of LIBS. These values agree with ion microprobe analyses, bulk analyses, or both. For eucryptite crystals, the Li concentrations are scattered because grain size is smaller than the LIBS spatial resolution (6 to 8 μm). Lithium concentrations of melt inclusions from the Streltsovka U deposit (Siberia) are in the range of 2 to 6.2 wt% (Li₂O) for Li-rich daughter minerals. Lithium estimations on silicate glasses display values between 90 and 400 ppm.Lithium was also analyzed as a trace element in quartz. Transverse profiles were performed in hydrothermal barren quartz veins from the Spanish Central System (Sierra de Guadarrama). The highest Li concentrations (250 to 370 ppm) were found in specific growth bands in conjunction with the observed variation in optical cathodoluminescence intensity. Considering the fluid inclusion analysis, the source of fluid responsible to the Li enrichment in quartz is probably high-salinity fluids derived from sedimentary basins.     

Uncertainties in quantitative time-lapse seismic analysis

Most seismic time-lapse studies so far have been of a qualitative nature. Identification of areas with minor or no seismic changes has been used to plan new infill drilling targets. Increased accuracy in seismic acquisition methods, in both conventional streamer surveys and newer methods such as multicomponent sea-bed seismic and permanent sensors, opens possibilities for the next step: quantitative time-lapse analysis. Quantitative methods here mean the estimation of, for instance, a change in fluid saturation from 20% water to 90% water or the estimation of a pore pressure change of 5 MPa. Explicit expressions for the uncertainties associated with estimated changes in, for instance, reservoir pressure and fluid saturation are derived. These formulae can be used to compare relative uncertainties between estimated parameters as well as to identify the critical factors in various estimation techniques. The importance of accurate rock physics input, as well as that of highly repeatable time-lapse seismic data, is emphasized. Furthermore, uncertainty analysis can be used to find optimal weight factors when the same parameter (e.g. saturation change) is estimated by two or three different techniques.     

电离层水平不均匀性对无线电掩星资料反演的影响

电离层水平不均匀结构是引起掩星资料反演的主要误差源.本文模拟不均匀结构对无线电掩星资料反演的影响.首先利用3-D全球射线追踪技术进行了模拟.结果指出:电离层电子密度20%的变化可能引起GPS的时延增加约12m.其次,将90%的掩星时延对应的路径长度,定义为电离层掩星测量的水平分辨率,其值为1200-2000 km.最后,计算了不同纬度电离层当地正午、子夜和日出时的电离层水平不对称指数.它表明:低纬夏季日出时,不对称指数达-78%,反演结果可能有较大误差.     

Melting of β-quartz up to 2.0 GPa and thermodynamic optimization of the silica liquidus up to 6.0 GPa

Melting relations of β-quartz were experimentally determined at 1.0 GPa (1900±20 °C), 1.5 GPa (2033±20 °C), and 2.0 GPa (2145±20 °C) using a new high-pressure assembly in a piston–cylinder apparatus and substantial differences were found with data previously reported. The new melting data of β-quartz were combined and optimized with all available thermodynamic, volumetric, and phase equilibria data for β-cristobalite, β-quartz and coesite to produce a P–T liquidus diagram for silica valid up to 6.0 GPa. Using the new optimized thermodynamic parameters, the invariant point β-cristobalite+β-quartz+liquid and β-quartz+coesite+liquid were determined to lie at 1687±17 °C and 0.457 GPa, and 2425±25 °C and 5.00 GPa, respectively.     

The molecular universe

Helen J Fraser , Martin R S McCoustra and David A Williams present a simple guide to astrochemistry.
Molecules play a fundamental role in many regions of our universe. The science where chemistry and astronomy overlap is known as astrochemistry, a branch of astronomy that has risen in importance over recent years. In this article we review the significance of chemistry in several astronomical years. IN this article we review the significance of chemistry in several astronomical environments including the early universe, interstellar clouds, starforming regions and protoplanetary disks. We discuss theoretical models, laboratory experiments and observational data, and present several recent and exciting results that challenge our perception of the "molecular universe".     

Multiparametric Traveltime Inversion

In conventional seismic processing, the classical algorithm of Hubral and Krey is routinely applied to extract an initial macrovelocity model that consists of a stack of homogeneous layers bounded by curved interfaces. Input for the algorithm are identified primary reflections together with normal moveout (NMO) velocities, as derived from a previous velocity analysis conducted on common midpoint (CMP) data. This work presents a modified version of the Hubral and Krey algorithm that is designed to extend the original version in two ways, namely (a) it makes an advantageous use of previously obtained common-reflection-surface (CRS) attributes as its input and (b) it also allows for gradient layer velocities in depth. A new strategy to recover interfaces as optimized cubic splines is also proposed. Some synthetic examples are provided to illustrate and explain the implementation of the method.