Sediment particle selectivity and its response to overland flow hydraulics within grass strips

Particle selectivity plays an important role in clarifying sediment transport processes in vegetative filter strips (VFS). 10-m long grass strips at slopes of 5° and 15° were subjected to a series of sediment-laden inflows experiments with different particle sizes to investigate the sediment transport and its response to overland flow hydraulics. The inflow sediments came from local soil, river-bed sand, and mixed, with median particle size d₅₀ of 39.9, 207.9 and 77.4 μm, respectively. Three independent repeated experiments were carried for each treatment. The results show that when the sediment trapping lasted for a certain length of time, the re-entrainment of some small-sized particles was greater than the deposition; that is, net loss occurred, which was not erosion of the original soil. Net loss of particles is mainly determined by the particle diameter. The coarser the inflow sediment particles and/or the steeper the slope, the coarser the particles can be net lost. Deposited sediment causes the VFS bed surface to become smooth and hydraulic resistance decrease exponentially. Unit stream power P is more suitable than shear stress τ of overland flow to be used to describe the process of sediment particle transport in VFS. The relationship between P and d₅₀ of outflow sediment is very consistent with the form of power function with a constant term. These results are helpful to understand the physical process of sediment transport on vegetation hillslopes.     

Organic litter: dominance over stones as a source of interrill flow roughness on low‐gradient desert slopes at Fowlers Gap,arid western NSW,Australia

Thirty‐six runoff plot experiments provide data on flow depths, speeds, and Darcy–Weisbach friction coefficients (f) on bare soil surfaces, and surfaces to which were added sufficient extra plant litter or surface stones to provide projected cover of 5, 10 and 20 per cent. Precision flow depth data were derived with a computer‐controlled gantry and needle gauge for two different discharges for each plot treatment. Taking a fixed flow intensity (Reynolds number, Re = 150) for purposes of comparison shows means of f = 17·7 for bare soil surfaces, f = 11·4 for added stone treatments, and f = 23·8 for added litter treatments. Many individual values of f for stone treatments are lower than for the bare soil surface, but all litter treatments show increases in fcompared to bare soil. The lowering of f in stone treatments relates to the submerged volume that the stones occupied, and the associated concentration of flow onto a smaller part of the plot surface. This leads to locally higher flow intensities and lower frictional drag along threads of flow that the obstacles create. Litter causes higher frictional drag because the particles are smaller, and, for the same cover fraction, are 100 times more numerous and provide 20 times the edge or perimeter length. Along these edges, which in total exceed 2·5 m g^?1 (equivalent to 500 m m^?2 for a loading of 2 t ha^?1), surface tension draws up water from between the litter particles. This reduces flow depth there, and as a consequence of the lower flow intensity, frictional drag rises. Furthermore, no clear passage remains for the establishment of flow threads. These findings apply to shallow interrill flows in which litter is largely immobile. The key new result from these experiments is that under these conditions, a 20 per cent cover of organic litter can generate interrill frictional retardation that exceeds by nearly 41 per cent that of a bare soil surface, and twice that contributed by the same cover fraction of surface stones. Even greater dominance by litter can be anticipated at the many dryland sites where litter covers exceed those tested here. Copyright © 2002 John Wiley & Sons, Ltd.     

Transport and deposition of fine sediment in open channels with different aspect ratios

This paper investigates by means of several large eddy simulations how the channel aspect ratio affects the transport and settling of suspended sediments. The numerical method is successfully validated using data of a physical experiment of fine sediment net deposition in an open channel flow. The channel aspect ratio, A, is known to be the determining factor for the development, strength and distribution of the turbulence‐driven secondary flow, and it is demonstrated that A influences the primary flow, turbulence quantities and the transport and fate of fine sediments. The secondary flow locally supports or hinders the falling of fine sediment particles in a turbulent flow, which results in a non‐uniform deposition of fine sediments over the cross‐section. While the channel aspect ratio has a large influence on the distribution of suspended sediments within the cross‐section, its effect on the cross‐sectional averaged deposition is negligibly small. Copyright © 2012 John Wiley & Sons, Ltd.     

Turbulent “Dynamo” in the Passive Scalar Problem

A turbulent magnetic dynamo can be considered as the evolution of a vector field in a turbulent fluid flow. The problem of evolution of scalar fields (e.g., number density of small particles) in a turbulent fluid flow is similar to the turbulent magnetic dynamo. The dynamo instability results in generation of magnetic field. The most important effect which can cause a generation of mean magnetic field in a turbulent fluid flow is the -effect: = – (1/3) u · ( × u), where u is the turbulent velocity field with the correlation time . A similar instability in the passive scalar problem results in formation of large-scale inhomogeneous structures in a spatial distribution of particles due to the -effect: = u_p ( · u_p), where u _p is the random velocity field of the particles which they acquire in a turbulent fluid velocity field. The effect is caused by inertia of particles which results in divergent velocity field of the particles. This results in additional turbulent nondiffusive flux of particles. The mean-field dynamics of inertial particles are studied by considering the stability of the equilibrium solution of the derived evolution equation for the mean number density of the particles in the limit of large Péclet numbers. The resulting equation is reduced to an eigenvalue problem for a Schrödinger equation with a variable mass, and a modified Rayleigh-Ritz variational method is used to estimate the lowest eigenvalue (corresponding to the growth rate of the instability). This estimate is in good agreement with obtained numerical solution of the Schrödinger equation. Similar effects arise during turbulent transport of gaseous admixtures (or light noninertial particles) in a low-Mach-number compressible fluid flow. The discussed effects are important in planetary and atmospheric physics (cloud formation, pollutant dynamics, preferential concentration of particles in protoplanetary disks and also planetesimals in them).     

Minimal advective travel time along arbitrary streamlines of porous media flows: The Fermat–Leibnitz–Bernoulli problem revisited

Travel time of marked fluid particles along arbitrary streamlines in arbitrary porous streamtubes is estimated from below based on the Cauchy–Bunyakovskii (Schwartz) and Jensen inequalities. In homogeneous media the estimate is strict and expressed through the length of the streamline, hydraulic conductivity, porosity and the head fall. The minimum is attained at streamlines of unidirectional flow. The bounds for heterogeneous soils, non-Darcian flows and unsaturated media are also written. If such bounds are attained the corresponding trajectories become brachistochrones. For example, in a two-layered aquifer and seepage perpendicular to the layers there is a unique conductivity–porosity ratio which makes a broken streamline brachistocronic. Similarly, if conductivities of two layers are fixed there is a unique incident angle between flow in one medium and the interface which makes a refracted streamline brachistocronic.     

The impact of slope length on the discharge of sediment by rain impact induced saltation and suspension

Simulations using a mechanistic model of raindrop driven erosion in rain‐impacted flow were performed with particles travelling by suspension, raindrop induced saltation and flow driven saltation. Results generated by both a high intensity storm, and a less intense one, indicate that, because of the effect of flow depth on the delivery of raindrop energy to the bed, there is a decline in sediment concentration, and hence soil loss per unit area, with slope length when particles are transported by raindrop induced saltation. However, that decline is reversed when the critical velocities that lead to flow driven saltation are episodically exceeded during an event. The simulations were performed on smooth surfaces and a single drop size but the general relationships are likely to apply for rain made up of a wide range of drop size. Although runoff is not always produced uniformly, as a general rule, flow velocities increase with slope length so that, typically, the distance particles travel before being discharged during an event increase with slope length. The effect of slope length on soil loss per unit area is often considered to vary with slope length to a power greater than zero and less that 1·0. The simulations show that effect of slope length on sediment discharge is highly dependent on the variations in runoff response resulting from variations in rainfall duration‐intensity‐infiltration conditions rather than plot length per se. Consequently, predicting soil loss per unit area using slope length with positive powers close to zero when sheet erosion occurs may not be as effective as commonly expected. Erosion by rain‐impacted flow is a complex process and that complexity needs to be considered when analysing the results of experiments associated with rain‐impacted flow under both natural and artificial conditions. Copyright © 2009 John Wiley & Sons, Ltd.     

Coarse bed load transport in an alluvial gravel bed stream,Dupuyer Creek,Montana

Coarse bed load was sampled in a gravel/cobble bed stream during two major floods in the snowmelt runoff season. The channel is characterized by high rates of bank erosion and, therefore, high rates of sediment supply and bed load flux. Peak discharge reached four times bank‐full, and bed load was sampled at flows 0·7–1·7 times bank‐full. A large aperture bed load sampler (1 m by 0·45 m) captured the largest particles in motion, and specifically targeted the coarse bed load size distribution by using a relatively large mesh (32 mm or D₂₅ of streambed surface size distribution). Bed load flux was highly variable, with a peak value of 0·85 kg/s/m for the coarse fraction above 38 mm. Bed load size distribution and maximum particle size was related to flow strength. Entrainment was size selective for particles D₇₀ and larger (88–155 mm), while particles in the range D₃₀–D₇₀ (35–88 mm) ceased to move at essentially the same flow. Bed load flux was size selective in that coarse fractions of the streambed surface were under‐represented in or absent from the bed load. Painted tracer particles revealed that the streambed surface in the riffles could remain stable even during high rates of bed load transport. These observations suggest that a large proportion of bed load sediments was sourced from outside the riffles. Repeat surveys confirmed major scour and fill in pools (up to 0·75 m), and bank erosion (>2 m), which together contributed large volumes of sediment to the bed load. Copyright © 2007 John Wiley & Sons, Ltd.     

The influence of roughness geometry and Shields parameter on flow resistance in gravel-bed channels

The spatial variability of bed particles of a gravel-bed channel is analysed and treated experimentally in order to simulate the effects of the arrangement of coarse bed elements on the flow resistance law. For the studied bed patterns, characterized by the concentration Γ of coarser elements arranged on the bed layer, a particle arrangement parameter α is proposed. The α parameter is useful for estimating the intercept b₀ of the semi-logarithmic flow resistance law deduced by flume measurements carried out for the hydraulic condition of large-scale and transition roughness. The differences between the experimental friction factor parameter values and the ones calculated by the proposed semi-logarithmic relationship are explained by the ratio between the Shields parameter and its critical value. The analysis shows that estimates of the friction factor parameter are not improved by introducing the Froude number into the flow resistance law. © 1997 by John Wiley & Sons, Ltd.     

Field studies on factors affecting very fine and ultra fine particulate organic matter deposition in low‐gradient sand‐bed streams

The knowledge on particle deposition in streams is mainly based on investigations in mountain streams. No data exist from low‐gradient sand‐bed streams that largely differ in the morphological and hydraulic factors proposed to affect deposition. To identify physical control on particle deposition in low‐gradient streams, we assessed deposition of very fine and ultra fine organic particulate matter in 18 sand‐bed stream reaches. We added particles derived from lake sediment and assessed the mean transport distance S_P and the deposition velocity v_dep. Additionally, reach hydraulics were estimated by injections of a conservative solute tracer (NaCl). Among the low‐gradient streams, particle deposition kinetics were variable but similar to deposition in mountain streams. S_P was solely related to the flow velocity. This relation was confirmed when comprising published data on deposition of fine organic particles. An association between particle deposition and transient storage factors was insignificant. We found significance of the transient storage to S_P only for repeated measures within a single reach, when flow velocity and benthic conditions were nearly constant. Measured v_dep/v_fall ratios were much larger than unity in most reaches. Evidence from this relation suggests that the vertical transport of very fine and ultra fine organic particulate matter through the water column was caused mainly by vertical mixing. Copyright © 2006 John Wiley & Sons, Ltd.     

Stochastic modelling of a diffusive wave for flood propagation using the random walk particle tracking method in a hypothetical city

Deterministic numerical schemes have been widely used for the solution of the diffusive wave (DW) equation, however, these schemes are computationally costly and suffer instability issues. This paper presents a stochastic random walk particle tracking (RWPT) method to solve such an equation for a dam‐break flow problem. Three different wave duration scenarios are presented for simulations of the DW for flood flows in a hypothetical city. The hypothetical city is represented by a domain of size 2,000 m by 500 m in x and y directions, respectively. The domain is divided into 25 m by 25 m cells. A dam is located at the upstream of the hypothetical city. Each scenario has a distinct propagation pattern after the dam is breached. Analysed and presented are 18 different simulations, which are composed of three different building configurations, two different bed slopes, and three different shapes of hydrographs. In this method, the flood volume is divided into a large number of particles where each particle carries a fixed amount of the flood volume. These particles undergo convective and diffusive movements, and their superposition represents propagation of the DW in the flow domain. The solution algorithm of the RWPT‐based equations is used to compute flood inundation depths in the hypothetical city. Comparison is drawn among the simulated results from three different shapes of the inflow hydrographs. The proposed stochastic method has two major advantages over traditional deterministic schemes: (a) greater efficiency, thus lesser computational costs, and (b) no instability issues.     

The yield strength of subliquidus basalts — experimental results

Yield strength is an important property of particle–fluid suspensions. In basaltic lavas that crystallize during flow emplacement, the onset of yield strength may result in threshold transitions in flow behavior and flow surface morphology. However, yield strength–crystallinity relations are poorly known, particularly in geologic suspensions, where difficulties of experimental and field measurements have limited data acquisition in the subliquidus temperature range. Here we describe two complementary experimental approaches designed to examine the effect of particle shape on the low-shear yield strength of subliquidus basalts. The first involves melting cubes of holocrystalline basalt samples with different initial textures to determine the temperature (crystallinity) at which these samples lose their cubic form. These experiments provide information on the minimum crystal volume fractions (0.20<φ<0.35) required to maintain the structual integrity of the cube. The second set of experiments uses suspensions of corn syrup and neutrally buoyant particles to isolate the effect of particle shape on yield strength development. From these experiments, we conclude that the shape is important in determining the volume fraction range over which suspensions exhibit a finite yield strength. As anisotropic particles may orient during flow, the effect of particle shape will be controlled by the orientation distribution of the constituent particles. We find that the so-called ‘excluded volume’ can be used to relate results of experiments on anisotropic particles to those of suspensions of spherical particles. Recent measurements of yield strength onset in basaltic melts at crystal volume fractions near 0.25 are consistent with our observations that crystal frameworks develop at low to moderate crystal volume fractions when crystals are anisotropic (e.g. plagioclase). We further suggest that conditions leading to yield strength onset at low crystallinities include rapid cooling (increased crystal anisotropy), heterogeneous nucleation (which promotes extensive crystal clustering and large cluster anisotropy) and static conditions (random crystal orientations).     

Displacement and rectification of planetary fluids

Abstract

Laboratory experiments on the decay (spin-up) of fluid motion on the β-plane are compared with theory. Under weakly dissipative conditions, some particles conserve potential vorticity during the entire decay. We also study the rectified mean flow which is produced by the lateral Reynolds Stress when a low frequency force is applied to the planetary fluid. The possible connection of effects with oceanic phenomena is briefly discussed.     

Predicting physical equations of soil detachment by simulated concentrated flow in Ultisols (subtropical China)

Soil detachment in concentrated flow is due to the dislodging of soil particles from the soil matrix by surface runoff. Both aggregate stability and shear strength of the topsoil reflect the erosion resistance of soil to concentrated runoff, and are important input parameters in predicting soil detachment models. This study was conducted to develop a formula to predict soil detachment rate in concentrated flow by using the aggregate stability index (As), root density (R_d) and saturated soil strength (σ_s) in the subtropical Ultisols region of China. The detachment rates of undisturbed topsoil samples collected from eight cultivated soil plots were measured in a 3.8 m long, 0.2 m wide hydraulic flume under five different flow shear stresses (τ = 4.54, 9.38, 15.01, 17.49 and 22.54 Pa). The results indicated that the stability index (As) was well related with soil detachment rate, particularly for results obtained with high flow shear stress (22.54 Pa), and the stability index (As) has a good linear relationship with concentrated flow erodibility factors (K_c). There was a positive linear relationship between saturated soil strength (σ_s) and critical flow shear stress (τ_c) for different soils. A significant negative exponential relationship between erodibility factors (K_c) and root density (R_d) was detected. This study yielded two prediction equations that allowed comparison of their efficiency in assessing soil detachment rate in concentrated flow. The equation including the root density (R_d) may have a better correlation coefficient (R² = 0.95). It was concluded that the formula based on the stability index (As), saturated soil strength (σ_s) and root density (R_d) has the potential to improve methodology for assessing soil detachment rate in concentrated flow for the subtropical Chinese Ultisols. Copyright © 2012 John Wiley & Sons, Ltd.     

Experimental analysis on the impact force of viscous debris flow

A miniaturized flume experiment was carried out to measure impact forces of viscous debris flow. The flow depth (7.2–11.2 cm), velocity (2.4–5.2 m/s) and impact force were recorded during the experiment. The impact process of debris flow can be divided into three phases by analyzing the variation of impact signals and flow regime. The three phases are the sudden strong impact of the debris flow head, continuous dynamic pressure of the body and slight static pressure of the tail. The variation of impact process is consistent with the change in the flow regime. The head has strong–rapid impact pressure, which is shown as a turbulent‐type flow; the body approximates to steady laminar flow. Accordingly, the process of debris flows hitting structures was simplified to a triangle shape, ignoring the pressure of the tail. In order to study the distribution of the debris flow impact force at different depths and variation of the impact process over time, the impact signals of slurry and coarse particles were separated from the original signals using wavelet analysis. The slurry's dynamic pressure signal appears to be a smooth curve, and the peak pressure is 12–34 kPa when the debris flow head hits the sensors, which is about 1.54 ± 0.36 times the continuous dynamic pressure of the debris flow body. The limit of application of the empirical parameter α in the hydraulic formula was also noted. We introduced the power function relationship of α and the Froude number of debris flows, and proposed a universal model for calculating dynamic pressure. The impact pressure of large particles has the characteristic of randomness. The mean frequency of large particles impacting the sensor is 210 ± 50–287 ± 29 times per second, and it is 336 ± 114–490 ± 69 times per second for the debris flow head, which is greater than that in the debris flow body. Peak impact pressure of particles at different flow depths is 40–160 kPa, which is 3.2 ± 1.5 times the impact pressure of the slurry at the bottom of the flow, 3.1 ± 0.9 times the flow in the middle, and 3.3 ± 0.9 times the flow at the surface. The differences in impact frequency indicate that most of the large particles concentrate in the debris flow head, and the number of particles in the debris flow head increases with height. This research supports the study of solid–liquid two phase flow mechanisms, and helps engineering design and risk assessment in debris flow prone areas. Copyright © 2015 John Wiley & Sons, Ltd.     

The impact of plant morphology on flow structure: comparative analysis of two types of submerged flexible macrophyte

ABSTRACT

The presence of aquatic vegetation in riverine and lacustrine environments alters the mean and turbulent flow structure and thus impacts the fate and transport of sediment and contaminants. Turbulent flows through Vallisneria natans (V. natans) and Potamogeton malaianus (P. malaianus) were investigated in a laboratory flume. The impact of plant morphology on mean velocity profile and turbulence distribution was analysed and discrepancies in flow alteration caused by different types of macrophyte were highlighted. Results show that a dense canopy of submerged macrophyte leads to a velocity profile featuring a counter velocity gradient in the lower part of the canopy. Negative Reynolds stress and its local maximum were observed there. Discrepancies in flow structure caused by different morphologies of both tested plants were further identified. With smaller frontal area in the lower part of the canopy, P. malaianus causes a much bigger gradient and local maximum in the velocity profile, and thus a larger local stress maximum than V. natans. The mean velocity gradient around the top of canopy, the Reynolds stress and the turbulence kinetic energy at the canopy interface are smaller than for the flow through the V. natans canopy. Larger reduction of the mean velocity within the V. natans canopy makes the suspended sediment of fine particles more easily deposited than in the P. malaianus canopy.     

Contaminant transport in a small deformation aquitard affected by the delayed drainage phenomenon

This paper presents a formulation accounting for the effect of delayed drainage phenomenon (DDP) on the breakthrough of contaminant flux in an aquitard, by considering the movement of soil particles, porosity variation, hydraulic head variation, and transient flow during the consolidation. The water flow equation in an aquitard was based on the Terzaghi's consolidation theory, and the contaminant transport equation was derived on the basis of the mass balance law. Two cases were used to illustrate the effect of DDP on the contaminant transport in an aquitard of small deformation. It is found that the breakthrough time of contaminant in an aquitard is very long, which is mainly ascribed to the low permeability of aquitard and sorption of soil particles. It is also found that the increase of depletion, which is in general induced by the increase of thickness and specific storativity and the decrease of hydraulic conductivity, enhances the impact of DDP on the contaminant transport in an aquitard. A larger delay index (τ₀) of DDP gives a greater delay breakthrough time (DBT) of solute transport in an aquitard, which controls the difference of the breakthrough time of contaminant transport in aquitards with and without the occurrence of DDP. For the cases where advection plays a dominant role during the process of solute transport, τ₀ is almost linearly correlated with DBT, and the ratio of DBT over the breakthrough time without consideration of DDP also approximately shows a linear relationship with the ratio of specific storativity to porosity, given a fixed drawdown in the adjacent aquifer with the sorption being ignored.     

Über die Bindung der Aerosolteilchen auf der Oberfläche der Wolkenelemente

Summary The importance ofStefan's flow by collection of insoluble aerosol particles on the surface of the evaporating droplets by condensation has recently been emphasized mainly in connection with the possible wash-out mechanism in the atmosphere. A simple theory is deduced for the collection efficiency of insoluble particles on droplets under theStefan's force and a comparision is made with the binding of aerosol particles byBrown's motion and microturbulence of air flow. In general we must take in account the influence ofStefan's flow by calculating the wash-out efficiency.     

Deposition and resuspension of fine particles in a riverine ‘dead zone’

Particulate matter suspended in the River Severn (Shropshire, UK) consists chiefly of clay-sized mineral particles, together with living and dead micro-organisms (algae and bacteria). Its concentration depends strongly on discharge, but the particle size distribution shows no systematic variability. For most samples, the particle volume is log-normally distributed with respect to diameter, the mean diameter being ca. 9 μm. The particles are mainly aggregates, including some with linear dimensions of the order of tens or hundreds of micrometres. Particle density depends appreciably on size, decreasing from ca. 2.5 × 10⁶ g m^?3 at a diameter of 2.5 μm to ca. 1.3 × 10⁶ g m ^?3 at 20 μm. The collision efficiency factor for particle aggregation is estimated to be 0.01–0.03. At low discharge, the ‘dead zone’ in the River Severn at Leighton is a well defined region of stagnant water behind a gravel bar. The rate of deposition of fine particles on its bed is of the order of tens of grams per square metre per day. Resuspension requires a critical bed shear velocity of 0.03–0.04 m s^?1, which occurs at main river discharges greater than about 150 m³ s^?1. Under such conditions the gravel bar is underwater and the dead zone is a region of highly turbulent return flow. A simple mechanistic model of particle dynamics in the dead zone accounts reasonably well for particle accumulation rates when run with parameter values based on measured particle and hydraulic properties. Calculations with the model suggest that most of the sedimentation flux to the dead zone bed is due to particles with equivalent sphere diameters in the range 30–240 μm. Simulations indicate that deposition proceeded continuously during spring and summer, whereas repeated deposition and resuspension occurred in autumn and winter.     

Polar patches observed by ESR and their possible origin in the cusp region

Observations by the EISCAT Svalbard radar in summer have revealed electron density enhancements in the magnetic noon sector under conditions of IMF B_z southward. The features were identified as possible candidates for polar-cap patches drifting anti-Sunward with the plasma flow. Supporting measurements by the EISCAT mainland radar, the CUTLASS radar and DMSP satellites, in a multi-instrument study, suggested that the origin of the structures lay upstream at lower latitudes, with the modulation in density being attributed to variability in soft-particle precipitation in the cusp region. It is proposed that the variations in precipitation may be linked to changes in the location of the reconnection site at the magnetopause, which in turn results in changes in the energy distribution of the precipitating particles.     

Energetic Electrons as a Field Line Topology Tracer in the High Latitude Boundary/CUSP Region: Cluster Rapid Observations

Energetic electrons (e.g., 50 keV) travel along field lines with a high speed of around 20 R_Es⁻¹. These swift electrons trace out field lines in the magnetosphere in a rather short time, and therefore can provide nearly instantaneous information about the changes in the field configuration in regions of geospace. The energetic electrons in the high latitude boundary regions (including the cusp) have been examined in detail by using Cluster/RAPID data for four consecutive high latitude/cusp crossings between 16 March and 19 March 2001. Energetic electrons with high and stable fluxes were observed in the time interval when the IMF had a predominately positive B_z component. These electrons appeared to be associated with a lower plasma density exhibiting no obvious tailward plasma flow (<20 keV). On the other hand, no electrons or only spike-like electron events have been observed in the cusp region during southward IMF. At that time, the plasma density was as high as that in the magnetosheath and was associated with a clear tailward flow. The fact that no stable energetic electron fluxes were observed during southward IMF indicates that the cusp has an open field line geometry. The observations indicate that both the South and North high latitude magnetospheric boundary regions (including both North and South cusp) can be energetic particle trapping regions. The energetic electron observations provide new ways to investigate the dynamic cusp processes. Finally, trajectory tracing of test particles has been performed using the Tsyganenko 96 model; this demonstrates that energetic particles (both ions and electrons) may be indeed trapped in the high latitude magnetosphere.