Observations and Simulations of the Polar Field Reversals in Cycle 23

We have used observations obtained by the National Solar Observatory at Kitt Peak to study the reversals of the polar magnetic fields in Cycle 23. We have compared them with corresponding data obtained by the Mt. Wilson Observatory, when these are available, testing both data sets against the locations of H filaments. Because of the unreliability of the data at extreme latitudes and because the apparent time of reversal varies with the degree of smoothing applied to the data, it is difficult to determine precise reversal time in each hemisphere from direct observations. However, we show that it is possible to obtain a better-defined and more precise reversal time using polar maps derived from simulations of the synoptic fields. These indirect values, however, depend critically on the diffusivity used in the simulations. We applied various tests to confirm an empirical value for the diffusivity parameter of about 600 km ² s ^–1 and hence determined empirical reversal times of CR 1976 in the northern hemisphere and CR 1981 in the south.     

Dust Trails of 8P/Tuttle and the Unusual Outbursts of the Ursid Shower

We calculate the position of dust trails from comet 8P/Tuttle, in an effort to explain unusual Ursid meteor shower outbursts that were seen when the comet was near aphelion. Comet 8P/Tuttle is a Halley-type comet in a 13.6-year orbit, passing just outside of Earth's orbit. We find that the meteoroids tend to be trapped in the 12:14 mean motion resonance with Jupiter, while the comet librates in a slightly shorter period orbit around the 13:15 resonance. It takes 6 centuries to decrease the perihelion of the meteoroid orbits enough to intersect Earth's orbit, during which time the meteoroids and comet separate in mean anomaly by 6 years, thus explaining the 6-year lag between the comet's return and Ursid outbursts. The resonances also prevent dispersion along the comet orbit and limit viewing to only one year in each return. We identified past dust trail encounters with dust trails from 1392 (Dec. 1945) and 1378 (Dec. 1986) and predicted another outburst on 2000 December 22 at around 7:29 and 8:35 UT, respectively, from dust trails dating to the 1405 and 1392 returns. This event was observed from California using video and photographic techniques. At the same time, five Global-MS-Net stations in Finland, Japan, and Belgium counted meteors using forward meteor scatter. The outburst peaked at 8:06±07 UT, December 22, at zenith hourly rate ∼90 per hour, and the Ursid rates were above half peak intensity during 4.2 h. We find that most Ursid orbits do scatter around the anticipated positions, confirming the link with comet 8P/Tuttle and the epoch of ejection. The 1405 and 1392 dust trails appear to have contributed similar amounts to the activity profile. Some orbits provide a hint of much older debris being present as well. This work is the strongest evidence yet for the relevance of mean motion resonances in Halley-type comet dust trail evolution.     

High-resolution Studies of the Polar Magnetic Fields during Cycle 23

High-resolution mosaics of the solar polar magnetic fields have been constructed using individual magnetograms obtained with the video magnetograph of the Big Bear Solar Observatory, and the properties of these mosaics are demonstrated in this paper. The mosaics show selected regions of the polar fields on several days during the rising phase of Cycle 23, and are related to the global polar fields (i) by superposing the mosaic for a given day on to a full-disk SOHO-MDI magnetogram obtained on the same day, (ii) by plotting the mosaics in polar projection and using these to identify the approximate regions reported by the mosaics on the NSOKP polar synoptic plots, and (iii) by imposing the locations of the H filaments on to the mosaics in order to infer the neutral lines of the large-scale fields. We have studied the fine structure of the large-scale unipolar fields near the poles and, in particular, have constructed histograms of the magnetic field intensities within particular regions of the mosaics and, in this way, have estimated the ratios of the number of magnetic knots of opposite polarities within the unipolar plumes. We have also generated enlargements of the polar regions of the NSOKP daily magnetograms. These and statistical studies have shown that on days for which the BBSO mosaics are not available, the NSOKP enlargements may be used to study the high-resolution polar fields. Time-series of mosaics obtained over four-hour periods on September 6 and November 18 show that considerable evolution in the structure of existing flux knots and the formation of several new knots has taken place during these periods.     

Observations of the Polar Magnetic Fields During the Polarity Reversals of Cycle 22

The Mount Wilson synoptic magnetic data for the period September 1987 through March 1996 are completely revised and used to provide polar plots of the solar magnetic fields for both hemispheres. This period, from Carrington rotations 1793 to 1906, covers the reversals of the polar magnetic fields in cycle 22. Comparison of our plots with the presently available H filtergrams for this period shows that the polarity boundaries are consistent in these two data sets where they overlap. The Mount Wilson plots show that the polar field reversals involve a complex sequence of events. Although the details differ slightly, the basic patterns are similar in each hemisphere. First the old polarity becomes isolated at the pole, then shortly thereafter, the isolation is broken, and the polar field includes unipolar regions of both polarities. The old polarity then reclaims the polar region, but when the isolation of this field is established for a second time, it declines in both area and strength. We take the reversal to be complete when the old polarity field is no longer observed in the Mount Wilson plots. With this criterion we find that the polar field reversal is completed in the north by CR 1836, i.e., by December 1990, and in the south by CR 1853, i.e., March 1992.     

Effects of variable attribute weights on landform classification

This paper focuses on the attribute weight issue and advocates use of modi?able attribute weights in terrain‐based environmental analysis and classi?cation. A question was asked: ‘How much will the result of a terrain‐based environmental analysis be affected if the weights of used terrain attributes are changed?’ The literature on landform classi?cation and the fuzzy k‐means method was reviewed in particular to help clarify the background and importance of this weight assignment issue. As an example, the effects of modifying attribute weights were evaluated for fuzzy k‐means landform classi?cation in a case study area. A total of 102 classi?cations were compared with each other and with a soil map, and comparison methods were speci?cally designed to evaluate the differences between these classi?cations. The results show that fuzzy k‐means landform classi?cation is sensitive to weight adjustments of adopted terrain attributes. The sensitivity is particularly high when the attribute weights started to be tuned away from the standard (i.e. uniform) weight of one. Better matching between landform classi?cation and a soil map may be produced when attribute weights are tuned. In all, we advocate the widespread adoption of an exploratory attitude in assigning attribute weights for environmental analysis and classi?cation. Copyright © 2006 John Wiley & Sons, Ltd.     

An approximate analytical solution for the deposition of heavy particles released from an elevated line source

This work re-examines and further develops an analytical solution for the deposition swath of heavy particles released in the atmosphere from an elevated source over uniform terrain, correcting the particle diffusivity for the crossing trajectory effect. The revised (approximate) analytical solution proves to be accurate within 20% over a wide range of micrometeorological conditions and particle size, despite its neglect of the turbulence component of the deposition flux. It compares very satisfactorily with experimental data and with the simulations of a Lagrangian stochastic model, provided the variable U(H)/w _g ≤7 (ratio of the mean horizontal wind speed at source height to the particle settling velocity). In this domain of validity, simple formulae relating the statistics of the deposition swath to U(H)/w _g are derived.     

A numerical simulation of residual circulation in Tampa Bay. Part I: Low-frequency temporal variations

The residual (time-average) salinity and circulation in a numerical ocean model of the Tampa Bay estuary are shown to experience significant temporal variation under realistic forcing conditions. A version of the Estuarine Coastal Ocean Model developed for Tampa Bay with 70 by 100 horizontal grid points and 11 sigma levels is examined for the years 2001–2003. Model output variables are averaged over the entire time of the simulation to generate long-term residual fields. The residual axial current is found to be dominated by the buoyancy-driven baroclinic circulation with an outflow (southwestward) at the surface and to the sides of the shipping channel, and an inflow (northeastward) usually occurring subsurface within or above the shipping channel. Averages over 30 d are used to examine variations in the residual fields. During the simulation the average surface salinity near the head of Tampa Bay varies with the freshwater inflow, from 12‰ to 33%. At the bay mouth salinity varies from 30%. to 36%.. A localized measure of the baroclinic circulation in the shipping channel indicates the residual circulation can vary strongly, attaining a magnitude triple the long-term mean value. The baroclinic circulation can be disrupted, going to near zero or even reversing, when the buoyancy-driven flow is weak and the surface winds are to the northeast. Three time periods, representing different environmental conditions, are chosen to examine these results in detail. A scaling argument indicates the relative strength of buoyancy versus wind as ΔρgH²(LC _Dω²)⁻¹, where δρ is head-to-mouth density difference across the bay,g is gravitational acceleration,H is depth,L is bay length,C _D is the surface wind drag coefficient, andw is wind speed. Tampa Bay is usually in the buoyancy dominated regime. The importance of winds in the weak-buoyancy case is demonstrated in an additional simulation without wind stress.