Effect of snow and soil frost melting on the concentrations of suspended solids and phosphorus in two rural watersheds in Western Finland     

Seppo Rekolainen 《Aquatic Sciences - Research Across Boundaries》1989,51(3):211-223

Hydrological processes are known to have a considerable effect on nutrient transport from agricultural land to watercourses. In cold temperate regions peak discharges are caused not only by storm conditions but also by melting of snow and frost. The objective of this work was to investigate the effects of snow and frost melt on concentrations of phosphorus and suspended solids. The samples were taken using flowweighted automatic sampling techniques from two agricultural drainage basins. During the beginning of the snowmelt period the concentration of suspended solids was rather low by comparison with the total phosphorus concentration and the discharge. The different behaviour compared with the relationships found during storm conditions was probably caused by continuous extraction of the soil surface by low ionic concentration melt water, and to some extent by leaching from plant residues. The gradual decrease of total phosphorus concentration during the frost thawing period could be attributed to the gradually increasing proportion of the subsurface and ground water discharge in the runoff.  相似文献   

Long-term evolution of the Neptune Trojan 2001 QR322     

R. Brasser  S. Mikkola  T.-Y. Huang  P. Wiegert  K. Innanen 《Monthly notices of the Royal Astronomical Society》2004,347(3):833-836

  相似文献   

An implementation ofN-body chain regularization     

Seppo Mikkola  Sverre J. Aarseth 《Celestial Mechanics and Dynamical Astronomy》1993,57(3):439-459

The chain regularization method (Mikkola and Aarseth 1990) for high accuracy computation of particle motions in smallN-body systems has been reformulated. We discuss the transformation formulae, equations of motion and selection of a chain of interparticle vectors such that the critical interactions requiring regularization are included in the chain. The Kustaanheimo-Stiefel (KS) coordinate transformation and a time transformation is used to regularize the dominant terms of the equations of motion. The method has been implemented for an arbitrary number of bodies, with the option of external perturbations. This formulation has been succesfully tested in a generalN-body program for strongly interacting subsystems. An easy to use computer program, written inFortran, is available on request.  相似文献   

A Time-Transformed Leapfrog Scheme     

Seppo Mikkola  Sverre Aarseth 《Celestial Mechanics and Dynamical Astronomy》2002,84(4):343-354

We present a time-transformed leapfrog scheme combined with the extrapolation method to construct an integrator for orbits in N-body systems with large mass ratios. The basic idea can be used to transform any second-order differential equation into a form which may allow more efficient numerical integration. When applied to gravitating few-body systems this formulation permits extremely close two-body encounters to be considered without significant loss of accuracy. The new scheme has been implemented in a direct N-body code for simulations of super-massive binaries in galactic nuclei. In this context relativistic effects may also be included.  相似文献   

Orbits of short period comets captured from the Oort cloud     

J. Q. Zheng  M. J. Valtonen  S. Mikkola  M. Korpi  H. Rickman 《Earth, Moon, and Planets》1995,71(3):45-50

Oort cloud comets occasionally obtain orbits which take them through the planetary region. The perturbations by the planets are likely to change the orbit of the comet. We model this process by using a Monte Carlo method and cross sections for orbital changes, i.e. changes in energy, inclination and perihelion distance, in a single planet-comet encounter. The influence of all major planets is considered. We study the distributions of orbital parameters of observable comets, i.e. those which have perihelion distance smaller than a given value. We find that enough comets are captured from the Oort cloud in order to explain the present populations of short period comets. The median value of cos i for the Jupiter family is 0.985 while it is 0.27 for the Halley types. The results may explain the orbital features of short period comets, assuming that the active lifetime of a comet is not much greater than 400 orbital revolutions.  相似文献   

Radio-continuum study of the Nearby sculptor group galaxies. Part 1: NGC 300 at λ=20 cm     

Timothy J. Galvin  Miroslav D. Filipović  Evan J. Crawford  Graeme Wong  Jeff L. Payne  Ain De Horta  Graeme L. White  Nick Tothill  Danica Drašković  Thomas G. Pannuti  Caleb K. Grimes  Benjamin J. Cahall  William C. Millar  Seppo Laine 《Astrophysics and Space Science》2012,340(1):133-142

A series of new radio-continuum (λ=20 cm) mosaic images focused on the NGC?300 galactic system were produced using archived observational data from the VLA and/or ATCA. These new images are both very sensitive (rms?=60 μJy) and feature high angular resolution (<10?″). The most prominent new feature is the galaxy’s extended radio-continuum emission, which does not match its optical appearance. Using these newly created images a number of previously unidentified discrete sources have been discovered. Furthermore, we demonstrate that a joint deconvolution approach to imaging this complete data-set is inferior when compared to an immerge approach.  相似文献   

Satellite attitude dynamics and estimation with the implicit midpoint method     

Christian Hellström  Seppo Mikkola 《New Astronomy》2009,14(5):467-477

We describe the application of the implicit midpoint integrator to the problem of attitude dynamics for low-altitude satellites without the use of quaternions. Initially, we consider the satellite to rotate without external torques applied to it. We compare the numerical solution with the exact solution in terms of Jacobi’s elliptic functions. Then, we include the gravity-gradient torque, where the implicit midpoint integrator proves to be a fast, simple and accurate method. Higher-order versions of the implicit midpoint scheme are compared to Gauss–Legendre Runge–Kutta methods in terms of accuracy and processing time. Finally, we investigate the performance of a parameter-adaptive Kalman filter based on the implicit midpoint integrator for the determination of the principal moments of inertia through observations.  相似文献   

Explicit algorithmic regularization in the few-body problem for velocity-dependent perturbations     

Christian Hellström  Seppo Mikkola 《Celestial Mechanics and Dynamical Astronomy》2010,106(2):143-156

A new algorithm is presented for the numerical integration of second-order ordinary differential equations with perturbations that depend on the first derivative of the dependent variables with respect to the independent variable; it is especially designed for few-body problems with velocity-dependent perturbations. The algorithm can be used within extrapolation methods for enhanced accuracy, and it is fully explicit, which makes it a competitive alternative to standard discretization methods.  相似文献   

One to One Resonance at High Inclination     

R. Brasser  D. C. Heggie  S. Mikkola 《Celestial Mechanics and Dynamical Astronomy》2004,88(2):123-152

We report results from long term numerical integrations and analytical studies of particular orbits in the circular restricted three-body problem. These are mostly high-inclination trajectories in 1 : 1 resonance starting at or near the triangular Lagrangian L₅ point. In some intervals of inclination these orbits have short Lyapunov times, from 100 to a few hundred periods, yet the osculating semi-major axis shows only relatively small fluctuations and there are no escapes from the 1 : 1 resonance. The eccentricity of these chaotic orbits varies in an erratic manner, some of the orbits becoming temporarily almost rectilinear. Similarly the inclination experiences large variations due to the conservation of the Jacobi constant. We studied such orbits for up to 10⁹ periods in two cases and for 10⁶ periods in all others, for inclinations varying from 0° to 180°. Thus our integrations extend from thousands to 10 million Lyapunov times without escapes of the massless particle. Since there are no zero-velocity curves restricting the motion this opens questions concerning the reason for the persistence of the 1 : 1 resonant motion. In the theory sections we consider the mechanism responsible for the observed behavior. We derive the averaged 1 : 1 resonance disturbing function, to second order in eccentricity, to calculate a critical inclination found in the numerical experiment, and examine motion close to this inclination. The cause of the chaos observed in the numerical experiments appears to be the emergence of saddle points in the averaged disturbing potential. We determine the location of several such saddle points in the (, ) plane, with being the mean longitude difference and the argument of pericentre. Some of the saddle points are illustrated with the aid of contour plots of the disturbing function. Motion close to these saddles is sensitive to initial conditions, thus causing chaos.This revised version was published online in October 2005 with corrections to the Cover Date.  相似文献   

Satellite image based vegetation classification of a large area using limited ground reference data: a case study in the Usa Basin, north-east European Russia     

Tarmo Virtanen  Kari Mikkola  Ari Nikula 《Polar research》2004,23(1):51-66

Predicted global changes can be studied effectively by combining spatially explicit data sets on vegetation and other landscape properties with process models. However, detailed knowledge of the vegetation distribution of remote Arctic areas is relatively scarce. This paper shows how a mesoscale vegetation and land cover classification of a large, remote Arctic area can be conducted at a fine spatial resolution (30 m cell size) using a limited ground reference data set. The study area is the catchment of the River Usa (93 500 km²) in north-eastern European Russia. Vegetation zones in the Usa Basin range from taiga in the south to forest-tundra and tundra in the north, and to alpine in the Ural mountains in the east. Classification was done using a mosaic of spectrally adjusted Landsat TM5 images from five different dates and a semi-supervised method. Ground reference data were collected during the summers of 1998, 1999 and 2000. Accuracy of the 21-class vegetation type/land cover classification produced was tested against test points interpreted from oblique aerial photographs taken from a helicopter (logistic limitations prohibited the collection of representative ground reference data). The main vegetation types (forests, willow dominated stands and meadows, peatlands, tundra heaths, mainly unvegetated areas, and water bodies) were distinguished with relatively high accuracy: 84% of the test points were classified correctly. Spatially detailed land cover data sets like the one described here allow detailed landscape-level analysis and process modelling on many different subjects.  相似文献