The ionospheric eclipse factor method (IEFM) and its application to determining the ionospheric delay for GPS

A new method for modeling the ionospheric delay using global positioning system (GPS) data is proposed, called the ionospheric eclipse factor method (IEFM). It is based on establishing a concept referred to as the ionospheric eclipse factor (IEF) λ of the ionospheric pierce point (IPP) and the IEF’s influence factor (IFF) . The IEF can be used to make a relatively precise distinction between ionospheric daytime and nighttime, whereas the IFF is advantageous for describing the IEF’s variations with day, month, season and year, associated with seasonal variations of total electron content (TEC) of the ionosphere. By combining λ and with the local time t of IPP, the IEFM has the ability to precisely distinguish between ionospheric daytime and nighttime, as well as efficiently combine them during different seasons or months over a year at the IPP. The IEFM-based ionospheric delay estimates are validated by combining an absolute positioning mode with several ionospheric delay correction models or algorithms, using GPS data at an international Global Navigation Satellite System (GNSS) service (IGS) station (WTZR). Our results indicate that the IEFM may further improve ionospheric delay modeling using GPS data.     

Temperature fluctuations in interstellar dust grains

Temperature fluctuations in interstellar dust grains caused by absorption of ultraviolet photons are discussed. Temperature probability distributions are presented for various assumptions about the grain specific heat and far infrared emissivity. Grains smaller than or about 0.01 may suffer large temperature fluctuations and would spend most of their time at very low temperatures. The equilibrium temperature is not a good measure for most average temperatures of such grains. It is shown that, although small grains spend only a small fraction of their time at the higher temperatures, the emitted far infrared spectrum is significantly shifted towards shorter wavelengths than predicted for grains assumed to be at the equilibrium temperature.Invited contribution to the Proceedings of a Workshop onThermodynamics and Kinetics of Dust Formation in the Space Medium held at the Lunar and Planetary Institute, Houston, 6–8 September, 1978.     

Interstellar clouds and the formation of stars

Part I gives a survey of the drastic revision of cosmic plasma physics which is precipitated by the exploration of the magnetosphere throughin situ measurements. The pseudo-plasma formalism, which until now has almost completely dominated theoretical astrophysics, must be replaced by an experimentally based approach involving the introduction of a number of neglected plasma phenomena, such as electric double layers, critical velocity, and pinch effect. The general belief that star light is the main ionizer is shown to be doubtful; hydromagnetic conversion of gravitational and kinetic energy may often be much more important.In Part II the revised plasma physics is applied to dark clouds and star formation. Magnetic fields do not necessarily counteract the contraction of a cloud; they may just as well pinch the cloud. Magnetic compression may be the main mechanism for forming interstellar clouds and keeping them together.Part III treats the formation of stars in a dusty cosmic plasma cloud. Star formation is due to an instability, but it is very unlikely that it has anything to do with the Jeans instability. A reasonable mechanism is that the sedimentation of dust (including solid bodies of different size) is triggering off a gravitationally assisted accretion. A stellesimal accretion analogous to the planetesimal accretion leads to the formation of a star surrounded by a very low density hollow in the cloud. Matter falling in from the cloud towards the star is the raw material for the formation of planets and satellites.The study of the evolution of a dark cloud leads to a scenario of planet formation which is reconcilable with the results obtained from studies based on solar system data. This means that the new approach to cosmical plasma physics discussed in Part I logically leads to a consistent picture of the evolution of dark clouds and the formation of solar systems.     

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.     

Methodology for construction, calibration and validation of a national hydrological model for Denmark

An integrated groundwater/surface water hydrological model with a 1 km² grid has been constructed for Denmark covering 43,000 km². The model is composed of a relatively simple root zone component for estimating the net precipitation, a comprehensive three-dimensional groundwater component for estimating recharge to and hydraulic heads in different geological layers, and a river component for streamflow routing and calculating stream–aquifer interaction. The model was constructed on the basis of the MIKE SHE code and by utilising comprehensive national databases on geology, soil, topography, river systems, climate and hydrology. The present paper describes the modelling process for the 7330 km² island of Sjælland with emphasis on the problems experienced in combining the classical paradigms of groundwater modelling, such as inverse modelling of steady-state conditions, and catchment modelling, focussing on dynamic conditions and discharge simulation. Three model versions with different assumptions on input data and parameter values were required until the performance of the final, according to pre-defined accuracy criteria, model was evaluated as being satisfactory. The paper highlights the methodological issues related to establishment of performance criteria, parameterisation and assessment of parameter values from field data, calibration and validation test schemes. Most of the parameter values were assessed directly from field data, while about 10 ‘free’ parameters were subject to calibration using a combination of inverse steady-state groundwater modelling and manual trial-and-error dynamic groundwater/surface water modelling. Emphasising the importance of tests against independent data, the validation schemes included combinations of split-sample tests (another period) and proxy-basin tests (another area).     

Diatom‐based reconstruction of summer sea‐surface salinity in the South China Sea over the last 15 000 years

We present a new reconstruction of summer sea‐surface salinity (SSS) over the past 15 000 years based on a diatom record from piston core 17940, located on the northern slope of the South China Sea (SCS). The reconstructed diatom‐based summer SSS values for the modern period are in accord with instrumental observations of summer SSS in the area. Here, the modern summer SSS is primarily controlled by river runoff, in particular from the Pearl River. The reconstruction presented in this study shows that the summer SSS varied between 33.3 and 34.2 psu over the past 15 000 years. The long‐term summer SSS trend closely followed the trend of the orbitally controlled solar insolation at 20°N, suggesting that orbital forcing was the dominant driver of changes in summer SSS in this area. Comparisons to speleothem δ¹⁸O data and studies of surface hydrography in the region suggest that changes in solar insolation affected the summer SSS through changes in the East Asian Monsoon and sea‐level changes associated with the last deglaciation. Univariate spectral analyses indicate that centennial‐scale oscillatory variations in summer SSS were superimposed on the long‐term trend. During the deglacial period (c. 12 000–9000 cal. a BP), the dominant periodicity was centred around 230–250 years, whereas a ～350‐year oscillation dominated in the period 2200–4500 cal. a BP. The balance of evidence suggests that these centennial‐scale changes in summer SSS may have been driven by solar‐induced changes in the East Asian Monsoon, but further evidence is needed to firmly establish this relationship.     

Simulation of nonlinear wave run-up with a high-order Boussinesq model

This paper considers the numerical simulation of nonlinear wave run-up within a highly accurate Boussinesq-type model. Moving wet–dry boundary algorithms based on so-called extrapolating boundary techniques are utilized, and a new variant of this approach is proposed in two horizontal dimensions. As validation, computed results involving the nonlinear run-up of periodic as well as transient waves on a sloping beach are considered in a single horizontal dimension, demonstrating excellent agreement with analytical solutions for both the free surface and horizontal velocity. In two horizontal dimensions cases involving long wave resonance in a parabolic basin, solitary wave evolution in a triangular channel, and solitary wave run-up on a circular conical island are considered. In each case the computed results compare well against available analytical solutions or experimental measurements. The ability to accurately simulate a moving wet–dry boundary is of considerable practical importance within coastal engineering, and the extension described in this work significantly improves the nearshore versatility of the present high-order Boussinesq approach.     

Velocity potential formulations of highly accurate Boussinesq-type models

The highly accurate Boussinesq-type equations of Madsen et al. (Madsen, P.A., Bingham, H.B., Schäffer, H.A., 2003. Boussinesq-type formulations for fully nonlinear and extremely dispersive water waves: Derivation and analysis. Proc. R. Soc. Lond. A 459, 1075–1104; Madsen, P.A., Fuhrman, D.R., Wang, B., 2006. A Boussinesq-type method for fully nonlinear waves interacting with a rapidly varying bathymetry. Coast. Eng. 53, 487–504); Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945) are re-derived in a more general framework which establishes the correct relationship between the model in a velocity formulation and a velocity potential formulation. Although most work with this model has used the velocity formulation, the potential formulation is of interest because it reduces the computational effort by approximately a factor of two and facilitates a coupling to other potential flow solvers. A new shoaling enhancement operator is introduced to derive new models (in both formulations) with a velocity profile which is always consistent with the kinematic bottom boundary condition. The true behaviour of the velocity potential formulation with respect to linear shoaling is given for the first time, correcting errors made by Jamois et al. (Jamois, E., Fuhrman, D.R., Bingham, H.B., Molin, B., 2006. Wave-structure interactions and nonlinear wave processes on the weather side of reflective structures. Coast. Eng. 53, 929–945). An exact infinite series solution for the potential is obtained via a Taylor expansion about an arbitrary vertical position z = zˆ. For practical implementation however, the solution is expanded based on a slow variation of zˆ and terms are retained to first-order. With shoaling enhancement, the new models obtain a comparable accuracy in linear shoaling to the original velocity formulation. General consistency relations are also derived which are convenient for verifying that the differential operators satisfy a potential flow and/or conserve mass up to the order of truncation of the model. The performance of the new formulation is validated using computations of linear and nonlinear shoaling problems. The behaviour on a rapidly varying bathymetry is also checked using linear wave reflection from a shelf and Bragg scattering from an undulating bottom. Although the new models perform equally well for Bragg scattering they fail earlier than the existing model for reflection/transmission problems in very deep water.