Modification of premare impact craters by volcanism and tectonism

Many lunar craters greater than 10 km in diam exhibit a variety of morphological characteristics which are not produced by meteorite impact or meteorite erosion. Most such craters are located in or near the margins of the maria. Although some could have resulted from processes such as cauldron resurgence, caldera formation, or ring dike emplacement, most have formed by modification of impact craters by endogenic processes including erosion by flowing lava, fissure volcanism, plutonism and uplift of crater floors along ring fractures of impact origin.     

Transport and emplacement of crater and basin deposits

Material is ejected from impact craters in ballastic trajectories; it impacts first near the crater rim and then at progressively greater ranges. Ejecta from craters smaller than approximately 1 km is laid predominantly on top of the surrounding surface. With increasing crater size, however, more and more surrounding surface will be penetrated by secondary cratering action and these preexisting materials will be mixed with primary crater ejecta. Ejecta from large craters and especially basin forming events not only excavate preexisting, local materials, but also are capable of moving large amounts of material away from the crater. Thus mixing and lateral transport give rise to continuous deposits that contain materials from within and outside the primary crater. As a consequence ejecta of basins and large highland craters have eroded and mixed highland materials throughout geologic time and deposited them in depressions inside and between older crater structures.Because lunar mare surfaces contain few large craters, the mare regolith is built up by successive layers of predominantly primary ejecta. In contrast, the lunar highlands are dominated by the effects of large scale craters formed early in lunar history. These effects lead to thick fragmental deposits which are a mixture of primary crater material and local components. These deposits may also properly be named regolith though the term has been traditionally applied only to the relatively thin fine grained surficial deposit on mare and highland terranes generated during the past few billion year. We believe that the surficial highland regolith - generated over long periods of time - rests on massive fragmental units that have been produced during the early lunar history.     

Low crater frequencies and low model ages in lunar maria: Recent endogenic activity or degradation effects?

Recently a number of studies have identified small lunar geologic structures to be <100 Ma in age using standard remote sensing techniques. Here we present new crater size frequency distributions (CSFDs) and model ages using craters D > 10 m for five small target units: one irregular mare patch (IMP) in Mare Nubium and four regions located on lunar wrinkle ridges in Mare Humorum. For comparison we also date another IMP found in a recent study in Mare Tranquillitatis (Braden et al. 2014 ). Absolute model age (AMA) derivation corresponds to 46 ± 5 Ma and 22 ± 1 Ma for Nubium and Sosigenes IMP, respectively. We show that for IMPs and in nearby control mare regions, similar production-like cumulative log–log SFD slopes of −3 are observed. In contrast, control mare regions in Mare Humorum exhibit shallower equilibrium slopes from −1.83 to −2. Three out of four wrinkle ridges appear to be in equilibrium but with crater lifetimes lower than on the corresponding maria. Low crater frequencies on one wrinkle ridge result in an age of 8.6 ± 1 Ma. This study region contains 80% fresh craters, which suggests that the crater population is still in production indicative of a recent resurfacing event.     

Structural uplift and ejecta thickness of lunar mare craters: New insights into the formation of complex crater rims

We investigate the elevated crater rims of lunar craters. The two main contributors to this elevation are a structural uplift of the preimpact bedrock and the emplacement of ejecta on top of the crater rim. Here, we focus on five lunar complex mare craters with diameters ranging between 16 and 45 km: Bessel, Euler, Kepler, Harpalus, and Bürg. We performed 5281 measurements to calculate precise values for the structural rim uplift and the ejecta thickness at the elevated crater rim. The average structural rim uplift for these five craters amounts to S_RU = 70.6 ± 1.8%, whereas the ejecta thickness amounts to E_T = 29.4 ± 1.8% of the total crater rim elevation. Erosion is capable of modifying the ratio of ejecta thickness to structural rim uplift. However, to minimize the impact of erosion, the five investigated craters are young, pristine craters with mostly preserved ejecta blankets. To quantify how strongly craters were enlarged by crater modification processes, we reconstructed the dimensions of the transient crater. The difference between the transient crater diameter and the final crater diameter can extend up to 11 km. We propose reverse faulting and thrusting at the final crater rim to be one of the main contributing factors of forming the elevated crater rim.     

Individual lava flow thicknesses in Oceanus Procellarum and Mare Serenitatis determined from Clementine multispectral data

We use multispectral reflectance data from the lunar Clementine mission to investigate the impact ejecta deposits of simple craters in two separate lunar mare basalt regions, one in Oceanus Procellarum and one in Mare Serenitatis. Over 100 impact craters are studied, and for a number of these we observe differences between the TiO₂ (and FeO) contents of their ejecta deposits and the lava flow units in which they are located. We demonstrate that, in the majority of cases, these differences cannot plausibly be attributed to uncorrected maturity effects. These observations, coupled with morphometric crater relationships that provide maximum crater excavation depths, allow the investigation of sub-surface lava flow stratigraphy. We provide estimated average thicknesses for a number of lava flow units in the two study regions, ranging from ∼80 m to ∼600 m. In the case of the Serenitatis study area, our results are consistent with the presence of sub-surface horizons inferred from recent radar sounding measurements from the JAXA Kaguya spacecraft. The average lava flow thicknesses we obtain are used to make estimates of the average flux of volcanic material in these regions. These are in broad agreement with previous studies, suggesting that the variation in mare basalt types we observe with depth is similar to the lateral variations identified at the surface.     

Transitional impact craters on the Moon: Insight into the effect of target lithology on the impact cratering process

We studied a data set of 28 well‐preserved lunar craters in the transitional (simple‐to‐complex) regime with the aim of investigating the underlying cause(s) for morphological differences of these craters in mare versus highland terrains. These transitional craters range from 15 to 42 km in diameter, demonstrating that the transition from simple to complex craters is not abrupt and occurs over a broad diameter range. We examined and measured the following crater attributes: depth (d), diameter (D), floor diameter (D_f), rim height (h), and wall width (w), as well as the number and onset of terraces and rock slides. The number of terraces increases with increasing crater size and, in general, mare craters possess more terraces than highland craters of the same diameter. There are also clear differences in the d/D ratio of mare versus highland craters, with transitional craters in mare targets being noticeably shallower than similarly sized highland craters. We propose that layering in mare targets is a major driver for these differences. Layering provides pre‐existing planes of weakness that facilitate crater collapse, thus explaining the overall shallower depths of mare craters and the onset of crater collapse (i.e., the transition from simple to complex crater morphology) at smaller diameters as compared to highland craters. This suggests that layering and its interplay with target strength and porosity may play a more significant role than previously considered.     

Chandrayaan-1 observation of distant secondary craters of Copernicus exhibiting central mound morphology: Evidence for low velocity clustered impacts on the Moon

Analysis of the Chandrayaan-1 Terrain Mapping Camera image of a 20 km×27 km area in the Mare Imbrium region revealed a cluster of thousands of fresh and buried impact craters in the size range of 20-1300 m. A majority of the large fresh craters with diameter ranging from 160 to 1270 m exhibit near-circular mounds (30-335 m diameter and 10-40 m height) in the crater floor, and their size depends on the host crater size. The origin of this cluster of secondary craters may be traced to Copernicus crater, based on global lunar image and the analysis of Chandrayaan-1 Hyper Spectral Imager data. Our findings provide further evidence for secondary crater formation by low-velocity impact of a cloud of clustered fragments. The presence of central mounds can also distinguish the secondary craters from the primary craters and refine the chronology of lunar surface based on counting of small craters.     

Floor-fractured lunar craters

Numerous lunar craters (206 examples, mean diameter = 40km) contain pronounced floor rilles (fractures) and evidence for volcanic processes. Seven morphologic classes have been defined according to floor depth and the appearance of the floor, wall, and rim zones. Such craters containing central peaks exhibit peak heights (approximately 1km) comparable to those within well-preserved impact craters but exhibit smaller rim-peak elevation differences (generally 0–1.5km) than those (2.4km) within impact craters. In addition, the morphology, spatial distribution, and floor elevation data reveal a probable genetic association with the maria and suggest that a large number of floor-fractured craters represent pre-mare impact craters whose floors have been lifted tectonically and modified volcanically during the epochs of mare flooding. Floor uplift is envisioned as floating on an intruded sill, and estimates of the buoyed floor thickness are consistent with the inferred depth of brecciation beneath impact craters, a zone interpreted as a trap for the intruding magma. The derived model of crater modification accounts for (1) the large differences in affected crater size and age; (2) the small peak-rim elevation differences; (3) remnant central peaks within mare-flooded craters and ringed plains; (4) ridged and flat-topped rim profiles of heavily modified craters and ringed plains; and (5) the absence of positive gravity anomalies in most floor-fractured craters and some large mare-filled craters. One of the seven morphologic classes, however, displays a significantly smaller mean size, larger distances from the maria, and distinctive morphology relative to the other six classes. The distinctive morphology is attributed, in part, to the relatively small size of the affected crater, but certain members of this class represent a style of volcanism unrelated to the maria - perhaps triggered by the last major basin-forming impacts.     

Preservation potential of implanted solar wind volatiles in lunar paleoregolith deposits buried by lava flows

The lunar surface is bathed in a variety of impacting particles originating from the solar wind, solar flares, and galactic cosmic rays. These particles can become embedded in the regolith and/or produce a range of other molecules as they pass through the target material. The Moon therefore contains a record of the variability of the solar and galactic particle fluxes through time. To obtain useful temporal snapshots of these processes, discrete regolith units must be shielded from continued bombardment that would rewrite the record over time. One mechanism for achieving this preservation is the burial of a regolith deposit by a later lava flow. The archival value of such deposits sandwiched between lava layers is enhanced by the fact that both the under- and over-lying lava can be dated by radiometric techniques, thereby precisely defining the age of the regolith layer and the geologic record contained therein. The implanted volatile species would be vulnerable to outgassing by the heat of the over-lying flow, at temperatures exceeding 300-700 °C. However, the insulating properties of the finely particulate regolith would restrict significant heating to shallow depths. We have therefore modeled the heat transfer between lunar mare basalt lavas and the regolith in order to establish the range of depths below which implanted volatiles would be preserved. We find that the full suite of solar wind volatiles, consisting predominantly of H and He, would survive at depths of ∼13-290 cm (for 1-10 m thick lava flows, respectively). A substantial amount of CO, CO₂, N₂ and Xe would be preserved at depths as shallow as 3.7 cm beneath meter-thick flows. Given typical regolith accumulation rates during mare volcanism, the optimal localities for collecting viable solar wind samples would involve stacks of thin mare lava flows emplaced a few tens to a few hundred Ma apart, in order for sufficient regolith to develop between burial events. Obtaining useful archives of Solar System processes would therefore require extraction of regolith deposits buried at quite shallow depths beneath radiometrically-dated mare lava flows. These results provide a basis for possible lunar exploration activities.     

A preliminary analysis of lunar extra-mare basalts: Distribution,compositions, ages,volumes, and eruption styles

Extra-mare basalts occupy 8.5% of the lunar basalt area and comprise 1% of the total mare basalt volume. They are preferentially located where the crust is thin and topographically low. In terms of age, eruption style, and composition they are as variable as the mare basalts. In some instances extrusion in extra-mare craters was preceded by floor-fracturing whereas in other cases it apparently was not. The volume of lava erupted may have been controlled more by the volume of magma produced than by hydrostatic effects. A minimum of nearly 1300 separate basalt eruptions, is indicated; the true value could be nearer 30 000 separate eruptions.     

A mechanism for the production of lunar crater rays

Observations of high resolution photographs of part of one of the prominent rays of the lunar crater Copernicus show that there is a concentration of small bright rayed and haloed craters within the ray. These craters contribute to the overall ray brightness; they have been measured and their surface distribution has been mapped. Sixty-two percent of the bright craters can be identified from study of high resolution photographs as concentric impact craters. These craters contain in their ejecta blankets, rocks from the lunar substrate that are brighter than the adjacent mare surface. It is concluded that the brightness of the large ray from the crater Copernicus is due to the composite effect of many small concentric impact craters with rocky ejecta blankets. If this is the dominant mechanism for the production of other rays from Copernicus and other large lunar craters, then rays may not contain significant amounts of ejecta from the central crater or from large secondary craters. They may in fact only reflect local excavation of mare substrate material by myriads of small secondary or tertiary impact craters.     

Volatiles,vesicles, and vugs: Unraveling the magmatic and eruptive histories of Steno crater basalts

In 1972, Apollo 17 astronauts returned 170.4 kg of lunar material. Within 1 month of their return, a subset of those samples was specially curated with the forethought that future analytical techniques would offer new insight into the formation and evolution of the Moon. Of interest in this work is sample 71036, a basalt collected from the rim of Steno crater in the Taurus–Littrow Valley, which was stored frozen and was processed and released for study 50 years later. We report, for the first time, the detailed mineralogy and petrology of 71036 and its companion samples 71035, 71037, and 71055 using a novel combination of 2-D and 3-D methods. We investigate lunar volatiles through in situ measurements of apatite and 3-D measurements of vesicles to understand the degassing histories of the Steno crater basalts. Our coupled 2-D petrography and 3-D tomography data sets support a model of the Steno crater basalts crystallizing in the upper crust of a mare lava flow. Apatite F and OH chemistry and the late-stage deformation of voids and formation of smaller vesicles provide evidence supporting coeval degassing of volatiles and crystallization of mesostasis apatite in Apollo 17 basalts. This work helps to close knowledge gaps surrounding the origin, magmatic evolution, emplacement, and crystallization history of high-titanium basalts.     

Photometric anomalies of the lunar surface studied with SMART-1 AMIE data

We present new results from the mapping of lunar photometric function parameters using images acquired by the spacecraft SMART-1 (European Space Agency). The source data for selected lunar areas imaged by the AMIE camera of SMART-1 and the data processing are described. We interpret the behavior of photometric function in terms of lunar regolith properties. Our study reveals photometric anomalies on both small (sub-kilometer) and large (tens of kilometers) scales. We found the regolith mesoscale roughness of lunar swirls to be similar in Mare Marginis, Mare Ingenii, and the surrounding terrains. Unique photometric properties related to peculiarities of the millimeter-scale regolith structure for the Reiner Gamma swirl are confirmed. We identified several impact craters of subkilometer sizes as the source of photometric anomalies created by an increase in mesoscale roughness within the proximal crater ejecta zones. The extended ray systems reveal differences in the photometric properties between proximal and distant ejecta blankets. Basaltic lava flows within Mare Imbrium and Oceanus Procellarum indicate higher regolith porosity for the redder soils due to differences in the chemical composition of lavas.     

Volcanic and structural history of the rocks exposed at Pickering Crater (Daedalia Planum, Mars)

In the western hemisphere of Mars Amazonian volcanism from Arsia Mons produced the smooth surfaces of Daedalia Planum and masks older rocks. Close to the southern termination of Daedalia Planum basement rocks are exposed in which are preserved craters that escaped or were only partially filled by this most recent volcanism. Pickering Crater is an approximately 130 km diameter crater. The youngest lavas flowed into this crater from Daedalia Planum by way of a NW rim breach, covering its western part. East of a well-defined flow front an older lava sequence with a distinctive platy surface and derived from a more proximal unestablished source to the northeast is exposed. Several units are identified within this sequence on the basis of surface texture, which is more subdued in progressively older rocks. Only local mapping of the flow front boundaries of these units is possible because of incomplete coverage by high resolution imagery. During emplacement of the older lavas a NE-SW striking en echelon graben system and parallel smaller troughs and dikes formed under inferred regional NW-SE extension. A much earlier strike-slip regime pre-dating the lavas exposed in the crater floor is postulated, based on the highly fretted nature of the rim of Pickering Crater and an elongated smaller crater to its northeast, approximately 40 km long in the NE-SW direction. The rims of these craters contrast with that of a smoother rimmed impact crater in the southeast that was excavated subsequent to strike-slip deformation but prior to the emplacement of platy surfaced lavas.     

The origin of lunar crater rays

Lunar rays are filamentous, high-albedo deposits occurring radial or subradial to impact craters. The nature and origin of lunar rays have long been the subjects of major controversies. We have determined the origin of selected lunar ray segments utilizing Earth-based spectral and radar data as well as FeO, TiO₂, and optical maturity maps produced from Clementine UVVIS images. These include rays associated with Tycho, Olbers A, Lichtenberg, and the Messier crater complex. It was found that lunar rays are bright because of compositional contrast with the surrounding terrain, the presence of immature material, or some combination of the two. Mature “compositional” rays such as those exhibited by Lichtenberg crater, are due entirely to the contrast in albedo between ray material containing highlands-rich primary ejecta and the adjacent dark mare surfaces. “Immaturity” rays are bright due to the presence of fresh, high-albedo material. This fresh debris was produced by one or more of the following: (1) the emplacement of immature primary ejecta, (2) the deposition of immature local material from secondary craters, (3) the action of debris surges downrange of secondary clusters, and (4) the presence of immature interior walls of secondary impact craters. Both composition and state-of-maturity play a role in producing a third (“combination”) class of lunar rays. The working distinction between the Eratosthenian and Copernican Systems is that Copernican craters still have visible rays whereas Eratosthenian-aged craters do not. Compositional rays can persist far longer than 1.1 Ga, the currently accepted age of the Copernican-Eratosthenian boundary. Hence, the mere presence of rays is not a reliable indication of crater age. The optical maturity parameter should be used to define the Copernican-Eratosthenian boundary. The time required for an immature surface to reach the optical maturity index saturation point could be defined as the Copernican Period.     

Monte Carlo calculations of lunar Regolith thickness distributions

Regolith thickness distributions associated with crater populations observed on selected maria surfaces have been calculated using a Monte Carlo computer technique. The calculations assume that the crater type produced and the volume of debris ejected and added to the growing regolith depends on the ratio of crater diameter and regolith thickness present at the time and place of formation of each crater. Calculated thickness distributions obtained are in agreement with those estimated using a previously described statistical method based on the morphology of small lunar craters. Additionally, the Monte Carlo calculations accurately predict the size frequency distributions of the same types of small, fresh lunar craters used in the statistical method. The model employed is therefore realistic. Furthermore, the model calculations presented are shown to have value (a) in predicting the thickness of the regolith from crater populations at various lunar sites, (b) relative dating applications in which crater populations are compared, and (c) in interpreting the origin and history of regolith deposits at specific locations.     

The geology of 433 Eros

Abstract— The global high‐resolution imaging of asteroid 433 Eros by the Near‐Earth Asteroid Rendezvous (NEAR) Shoemaker spacecraft has made it possible to develop the first comprehensive picture of the geology of a small S‐type asteroid. Eros displays a variety of surface features, and evidence of a substantial regolith. Large scale facets, grooves, and ridges indicate the presence of at least one global planar structure. Directional and superposition relations of smaller structural features suggest that fracturing has occurred throughout the object. As with other small objects, impact craters dominate the overall shape as well as the small‐scale topography of Eros. Depth/diameter ratios of craters on Eros average ～0.13, but the freshest craters approach lunar values of ～0.2. Ejecta block production from craters is highly variable; the majority of large blocks appear to have originated from one 7.6 km crater (Shoemaker). The interior morphology of craters does not reveal the influence of discrete mechanical boundaries at depth in the manner of craters formed on lunar mare regolith and on some parts of Phobos. This lack of mechanical boundaries, and the abundant evidence of regolith in nearly every high‐resolution image, suggests a gradation in the porosity and fracturing with depth. The density of small craters is deficient at sizes below ～200 m relative to predicted slopes of empirical saturation. This characteristic, which is also found on parts of Phobos and lunar highland areas, probably results from the efficient obliteration of small craters on a body with significant topographic slopes and a thick regolith. Eros displays a variety of regolith features, such as debris aprons, fine‐grained “ponded” deposits, talus cones, and bright and dark streamers on steep slopes indicative of efficient downslope movement of regolith. These processes serve to mix materials in the upper loose fragmental portion of the asteroid (regolith). In the instance of “ponded” materials and crater wall deposits, there is evidence of processes that segregate finer materials into discrete deposits. The NEAR observations have shown us that surface processes on small asteroids can be very complex and result in a wide variety of morphologic features and landforms that today seem exotic. Future missions to comets and asteroids will surely reveal still as yet unseen processes as well as give context to those discovered by the NEAR Shoemaker spacecraft.     

Analysis of crater distributions in mare units on the lunar far side

Mare material is asymmetrically distributed on the Moon. The Earth-facing hemisphere, where the crust is believed to be 26 km thinner than on the farside, contains substantially more basaltic mare material. Using Lunar Topographic Orthophoto Maps, we calculated the thickness of the mare material in three farside craters, Aitken (0.59 km), Isaev (1.0 km), and Tsiolkovskiy (1.75 km). We also studied crater frequency distribution in five farside mare units (Aitken, Isaev, Lacus Solitudinis, Langemak, and Tsiolkovskiy) and one light plains unit (in Mendeleev). Nearly 10 000 farside craters were counted. Analysis of the crater frequency on the light plains unit gives an age of 4.3 billion yr. Crater frequency distributions on the mare units indicate ages of 3.7 and 3.8 billion yr, suggesting that the units are distributed over a narrow time period of approximately 100 million yr. Returned lunar samples from nearside maria give dates as young as 3.1 billion yr. The results of this study suggest that mare basalt emplacement on the far side ceased before it did on the near side.     

Formation of lunar mare domes along crustal fractures: Rheologic conditions, dimensions of feeder dikes, and the role of magma evolution

In this study we examine a set of lunar mare domes located in the Hortensius/Milichius/T. Mayer region and in northern Mare Tranquillitatis with respect to their formation along crustal fractures, their rheologic properties, the dimensions of their feeder dikes, and the importance of magma evolution processes during dome formation. Many of these domes display elongated summit vents oriented radially with respect to major impact basins, and several dome locations are also aligned in these preferential directions. Analysis of Clementine UV/VIS and Lunar Prospector gamma ray spectrometer data reveals that the examined mare domes formed from low-Si basaltic lavas of high FeO and low to moderate TiO₂ content. Based on their morphometric properties (diameter, height, volume) obtained by photoclinometric and shape from shading analysis of telescopic CCD images, we derive rheologic quantities (lava viscosity during eruption, effusion rate, duration of the effusion process, magma rise speed) and the dimensions of the feeder dikes. We establish three rheologic groups characterised by specific combinations of rheologic properties and dike dimensions, where the most relevant discriminative parameter is the lava viscosity η. The first group is characterised by and contains the domes with elongated vents in the Milichius/T. Mayer region and two similar domes in northern Mare Tranquillitatis. The second group with comprises the very low aligned domes in northern Mare Tranquillitatis, and the third group with the relatively steep domes near Hortensius and in the T. Mayer region. The inferred dike dimensions in comparison to lunar crustal thickness data indicate that the source regions of the feeder dikes are situated within the upper crust for six of the domes in northern Mare Tranquillitatis, while they are likely to be located in the lower crust and in the upper mantle for the other examined domes. By comparing the time scale of magma ascent with the time scale on which heat is conducted from the magma into the host rock, we find evidence that the importance of magma evolution processes during ascent such as cooling and crystallisation increases with lava viscosity. We conclude that different degrees of evolution of initially fluid basaltic magma are able to explain the broad range of lava viscosities inferred for the examined mare domes. The spectral data reveal that differences in TiO₂ content may additionally account for the systematic difference in lava viscosity between the two examined lunar regions. We show that the described mechanisms are likely to be valid also for other lunar mare domes situated near Cauchy and Arago, regarded for comparison. On the other hand, we find for the Gruithuisen and Mairan highland domes that despite their inferred high lava viscosities of , no significant magma cooling in the dike occurred during ascent, supporting previous findings that the highland domes were formed during a specific phase of non-mare volcanism by highly silicic viscous lavas.     

Lunar Clinopyroxene and Plagioclase: Surface Distribution and Composition

The Clementine UVVIS images and the spectral and chemical (mineral) characteristics of lunar soil samples previously measured by the Lunar Samples Characterization Consortium were used to map the plagioclase and clinopyroxene abundance in the lunar surface material. An excess of plagioclase was found in young highland craters (e.g., in the crater Tycho) and in their ray systems. For clinopyroxenes, analogous behavior was observed in mare craters (e.g., in the crater Aristarchus). The maps for the FeO and Al₂O₃ bulk contents and the contents of these oxides in plagioclase and clinopyroxene were estimated by the same technique. These maps were compared to each other and to the predicted distribution of the lunar regolith maturity. The regolith of highland ray systems (e.g., the Tycho crater system) is characterized not only by low maturity but also by peculiar iron and aluminum contents: the lower the soil maturity degree, the smaller the iron content and the greater the aluminum content. This is confirmed by the data for the lunar soil samples from the Apollo 16 landing site. A cluster analysis of the “clinopyroxene content-maturity” and “plagioclase content-maturity” correlation diagrams allowed the mineral mapping of the lunar surface to be performed.__________Translated from Astronomicheskii Vestnik, Vol. 39, No. 4, 2005, pp. 291–303.Original Russian Text Copyright © 2005 by Shkuratov, Kaydash, Pieters.