Advances in Numerical Modeling of Astrophysical and Space Plasmas

Advances in the simulation of astrophysical and cosmic plasmas are the direct result of advances in computational capabilities, today consisting of new techniques such as multilevel concurrent simulation, multi-teraflop computational platforms and experimental facilities for producing and diagnosing plasmas under extreme conditions for the benchmarking of simulations. Examples of these are the treatment of mesoscalic plasma and the scaling to astrophysical and cosmic dimensions and the Accelerated Strategic Computing Initiative whose goal is to construct petaflop (10¹⁵ floating operations per second) computers, and pulsed power and laser inertial confinement plasmas where megajoules of energy are delivered to highly-diagnosed plasmas. This paper concentrates on the achievements to date in simulating and experimentally producing plasmas scaled to both astrophysical and cosmic plasma dimensions. A previous paper (Part I, Peratt, 1997) outlines the algorithms and computational growth.     

Measurement of seabed sound speeds from head waves in shallow water

Acoustic signals from small explosive charges have been measured with sonobuoys on twelve tracks in Australian northern shallow waters with the aim of assessing whether useful geoacoustic information could be obtained. Using the frequency band from 14 to 70 Hz, travel times of head waves were monitored, and the sound speeds and depths of corresponding interfaces in the seabed were derived. The water sound speed varied a little with range, and its depth dependence was allowed for by using its average value. Head waves from interfaces indistinguishable from the seafloor (the water/seabed interface) were detected on only three of the tracks, with derived sound speeds of 2100 to 2300 m/s. The first sub-bottom interfaces were from 50 to 600 m beneath the seafloor, and their sound speeds ranged from around 2000 m/s to 6400 m/s. Thus the head waves were from chalk or limestone, cemented sediments in which sound-speed gradients would be small. The amount of data obtained for the seafloor was limited by incoherence of the signals and, for some tracks, by excessive spacing between shots. The incoherence is generally attributed to multiple head waves that are individually unresolvable, while on two tracks there were indications of medal ground waves. Occasional anomalous data were obtained, but generally the assumptions of the simple interpretation method were found to be valid. Since no curvature in the range-time lines was observed, there was no evidence of sub-bottom sound-speed gradients being significant     

Glacial Lake Victoria,a high‐level Antarctic Lake inferred from lacustrine deposits in Victoria Valley

We present evidence of a large lake (Glacial Lake Victoria) that existed in Victoria Valley in the dry valleys region of Antarctica between at least 20 000 and 8600 ¹⁴C yr BP. At its highstands, Glacial Lake Victoria covered 100 km² and was ca. 200 m deep. The chronology for lake‐level changes comes from 87 AMS radiocarbon dates of lacustrine algae preserved in deltas and glaciolacustrine deposits that extend up to 185 m above present‐day lakes on the valley floor. The existence of Glacial Lake Victoria, as well as other large lakes in the dry valleys, indicates a climate regime significantly different from that of today at the last glacial maximum and in the early Holocene. Copyright © 2002 John Wiley & Sons, Ltd.     

Exchange of ejecta between Telesto and Calypso: Tadpoles, horseshoes, and passing orbits

We have numerically integrated the orbits of ejecta from Telesto and Calypso, the two small Trojan companions of Saturn’s major satellite Tethys. Ejecta were launched with speeds comparable to or exceeding their parent’s escape velocity, consistent with impacts into regolith surfaces. We find that the fates of ejecta fall into several distinct categories, depending on both the speed and direction of launch.The slowest ejecta follow suborbital trajectories and re-impact their source moon in less than one day. Slightly faster debris barely escape their parent’s Hill sphere and are confined to tadpole orbits, librating about Tethys’ triangular Lagrange points L₄ (leading, near Telesto) or L₅ (trailing, near Calypso) with nearly the same orbital semi-major axis as Tethys, Telesto, and Calypso. These ejecta too eventually re-impact their source moon, but with a median lifetime of a few dozen years. Those which re-impact within the first 10 years or so have lifetimes near integer multiples of 348.6 days (half the tadpole period).Still faster debris with azimuthal velocity components ?10 m/s enter horseshoe orbits which enclose both L₄ and L₅ as well as L₃, but which avoid Tethys and its Hill sphere. These ejecta impact either Telesto or Calypso at comparable rates, with median lifetimes of several thousand years. However, they cannot reach Tethys itself; only the fastest ejecta, with azimuthal velocities ?40 m/s, achieve “passing orbits” which are able to encounter Tethys. Tethys accretes most of these ejecta within several years, but some 1% of them are scattered either inward to hit Enceladus or outward to strike Dione, over timescales on the order of a few hundred years.     

Saturn eddy momentum fluxes and convection: First estimates from Cassini images

We apply an automated cloud feature tracking algorithm to estimate eddy momentum fluxes in Saturn's southern hemisphere from Cassini Imaging Science Subsystem near-infrared continuum image sequences. Voyager Saturn manually tracked images had suggested no conversion of eddy to mean flow kinetic energy, but this was based on a small sample of <1000 wind vectors. The automated procedure we use for the Cassini data produces an order of magnitude more usable wind vectors with relatively unbiased sampling. Automated tracking is successful in and around the westward jet latitudes on Saturn but not in the vicinity of most eastward jets, where the linearity and non-discrete nature of cloud features produces ambiguous results. For the regions we are able to track, we find peak eddy fluxes and a clear positive correlation between eddy momentum fluxes and meridional shear of the mean zonal wind, implying that eddies supply momentum to eastward jets and remove momentum from westward jets at a rate . The behavior we observe is similar to that seen on Jupiter, though with smaller eddy-mean kinetic energy conversion rates per unit mass of atmosphere (). We also use the appearance and rapid evolution of small bright features at continuum wavelengths, in combination with evidence from weak methane band images where possible, to diagnose the occurrence of moist convective storms on Saturn. Areal expansion rates imply updraft speeds of over the convective anvil cloud area. As on Jupiter, convection preferentially occurs in cyclonic shear regions on Saturn, but unlike Jupiter, convection is also observed in eastward jet regions. With one possible exception, the large eddy fluxes seen in the cyclonic shear latitudes do not seem to be associated with convective events.     

Relative sea‐level observations in western Scotland since the Last Glacial Maximum for testing models of glacial isostatic land movements and ice‐sheet reconstructions

Observations of relative sea‐level change and local deglaciation in western Scotland provide critical constraints for modelling glacio‐isostatic rebound in northern Britain over the last 18 000 years. The longest records come from Skye, Arisaig and Knapdale with a shorter, Holocene, record from Kintail. Biostratigraphic (diatom, pollen, dinoflagellate, foraminifera and thecamoebian), lithological and radiocarbon analyses provide age and elevation parameters for each sea‐level index point. All four sites reveal relative sea‐level change that is highly non‐monotonic in time as the local vertical component of glacio‐isostatic rebound and eustasy (or global meltwater influx) dominate at different periods. Copyright © 2006 John Wiley & Sons, Ltd.     

The Nippon/Norway Svalbard Meteor Radar: First results of small-scale structure observations

1　ＩｎｔｒｏｄｕｃｔｉｏｎＡｔｌａｔｉｔｕｄｅｓａｒｏｕｎｄ 80°Ｎａｎｄｉｎｔｈｅｍｅｓｏｓｐｈｅｒｅ /ｌｏｗｅｒｔｈｅｒｍｏｓｐｈｅｒｅ (ＭＬＴ) ,ｔｈｅｒｅｈａｖｅｂｅｅｎｆｅｗｍｅａｓｕｒｅｍｅｎｔｓｏｆｎｅｕｔｒａｌｄｙｎａｍｉｃｓ.ＡｒｅｑｕｉｒｅｍｅｎｔｗａｓｓｅｅｎｆｏｒｌｏｎｇｔｅｒｍｃｏｎｔｉｎｕａｌｍｏｎｉｔｏｒｉｎｇａｎｄａＶＨＦｍｅｔｅｏｒｒａｄａｒｗａｓｉｄｅｎｔｉｆｉｅｄａｓｂｅｉｎｇａｍｏｓｔｓｕｉｔａｂｌｅｉｎｓｔｒｕｍｅｎｔ.ＲａｄａｒｓｌｉｋｅｔｈｅＥＩＳＣＡ…     

Rainfall interception by trees of Pinus radiata,and Eucalyptus viminalis,in a 1300 mm rainfall area of southeastern New South Wales: I. Gross losses and their variability

Interception loss, I, was determined by continuous concurrent measurements of the canopy precipitation balances of a mature seed orchard tree of Pinus radiata, and a dominant tree of Eucalyptus viminalis at a mountainous high rainfall site (900 m a.s.l.) in Tallaganda State Forest of the Upper Shoalhaven Catchment. Approximate canopy storage capacity (Sc) of the pine was 54 l, and that of the eucalypt was 11·3 l. Gross pine I was 26·5 per cent and eucalypt I was 8·3 per cent of total incident rainfall over a period of 18 months, from June 1975 to December 1976. The exponential model that provided the best fit to overall data relating I to gross rainfall (P_g) was of good precision for the pine (r² = 0·73) but rather poor precision for the eucalypt (r² = 0·27). A consistent pattern in interception data of the two canopy types suggested that variation in I was related to change in pervasive conditions influencing rates of evaporation from wet canopies during rainfall. Multiple regression analyses confirmed that factors such as rainfall intensity and windspeed explained some of the variation in eucalypt I but little in pine I. Negative eucalypt I and corresponding low values of pine I over a wide range of P_g (up to 20 mm) suggest that capture of wind-borne precipitation (cloud, mist, or fog) had also complicated the canopy precipitation balances.