Microtextures and crystal chemistry of pigeonite in the ureilites ALHA77257, RKPA80239, Y‐791538, and ALHA81101

Abstract— The microtextures of pigeonite in four ureilites, Allan Hills (ALH) 77257, Reckling Peak (RKP) A80239, Yamato (Y‐) 791538, and Allan Hills A81101, chosen to span a range of composition and shock level, were investigated by transmission electron microscopy (TEM); two of the samples were also investigated by single crystal X‐ray diffraction to determine Fe²⁺‐Mg cation site partitioning. The low‐shock and compositionally homogeneous pigeonites in ALHA77257 and RKPA80329 (Wo 6.4 for both, mg 86.3 and 84.3 respectively) display irregularly spaced, shock‐induced stacking faults oriented parallel to (100), and large antiphase domains (50–100 nm). Antiphase domains have no preferential orientation. No evidence of exsolution was observed. The low‐shock Y‐791538 pigeonite is homogeneous and has higher Ca and mg (Wo 9.4, mg 91.2). TEM investigation showed spinodal decomposition, indicative of incipient exsolution; small antiphase domains were observed (～5 nm). Single crystal refinement yielded R_4s? = 5.71%, with Fe²⁺‐Mg partitioning coefficient k_d = 0.077(8) and T_c = 658(35) °C. ALHA81101 has compositionally heterogeneous pyroxenes, with large local variations in Wo and mg (Wo = 4–13, mg = 86–68). No compositional gradients from core to rim of grains were observed, and the heterogeneity is interpreted as related to cation migration during shock. In one relatively Ca‐rich region (Wo～12), TEM analysis showed augite‐pigeonite exsolution lamellae, with spacing 145(20) nm. Results for ALHA77257, RKPA80239, and Y‐791538 support a model of rapid cooling following breakup of the ureilite parent body. The presence of exsolution lamellae in ALHA81101 can be related to a local shock‐induced Ca enrichment and provides no constraint on the late cooling history.     

The main belt as a source of near-Earth asteroids

We investigate the flux of main-belt asteroid fragments into resonant orbits converting them into near-Earth asteroids (NEAs), and the variability of this flux due to chance interasteroidal collisions. A numerical model is used, based on collisional physics consistent with the results of laboratory impact experiments. The assumed main-belt asteroid size distribution is derived from that of known asteroids extrapolated down to sizes of ≈ 40 cm, modified in such a way to yield a quasi-stationary fragment production rate over times ≈ 100 Myr. The results show that the asteroid belt can supply a few hundred km-sized NEAs per year, well enough to sustain the current population of such bodies. On the other hand, if our collisional physics is correct, the number of existing 10-km objects implies that these objects either have very long-lived orbits, or must come from a different source (i.e., comets). Our model predicts that the fragments supplied from the asteroid belt have initially a power-law size distribution somewhat steeper than the observed one, suggesting preferential removal of small objects. The component of the NEA population with dynamical lifetimes shorter than or of the order of 1 Myr can vary by a factor reaching up to a few tens, due to single large-scale collisions in the main belt; these fluctuations are enhanced for smaller bodies and faster evolutionary time scales. As a consequence, the Earth's cratering rate can also change by about an order of magnitude over the 0.1 to 1 Myr time scales. Despite these sporadic spikes, when averaged over times of 10 Myr or longer the fluctuations are unlikely to exceed a factor two.     

The luminosity oscillation of X-ray bursters and recurrent novae

Thermo-mechanical oscillations of a radiating spherically-symmetric shell containing a fermion gas (i.e., oscillation where only mechanical, thermodynamical, and eventually radiation phenomena are taken into account) are studied. With a reasonable choice of the three relevant parameters all other observational data are similar to the ones of X-ray bursters and recurrent novae.Therefore, thermo-mechanical oscillations could play an important role in the oscillation spectrum of these astronomical objects.     

Evaluation of approaches to calculate debris-flow parameters for hazard assessment

Many different runout prediction methods can be applied to estimate the mobility of future debris flows during hazard assessment. The present article reviews the empirical, analytical, simple flow routing and numerical techniques. All these techniques were applied to back-calculate a debris flow, which occurred in 1982 at La Guingueta catchment, in the Eastern Pyrenees. A sensitivity analysis of input parameters was carried out, while special attention was paid to the influence of rheological parameters. We used the Voellmy fluid rheology for our analytical and numerical modelling, since this flow resistance law coincided best with field observations. The simulation results indicated that the “basal” friction coefficients rather affect the runout distance, while the “turbulence” terms mainly influence flow velocity. A comparison of the velocity computed on the fan showed that the analytical model calculated values similar to the numerical ones. The values of our rheological parameters calibrated at La Guingueta agree with data back-calculated for other debris flows. Empirical relationships represent another method to estimate total runout distance. The results confirmed that they contain an important uncertainty and they are strictly valid only for the conditions, which were the basis for their development. With regards to the simple flow routing algorithm, this methods could satisfactorily simulate the total area affected by the 1982 debris flow, but it was not able to directly calculate total runout distance and velocity. Finally, a suggestion on how different runout prediction methods can be applied to generate debris-flow hazard maps is presented. Taking into account the definition of hazard and intensity, the best choice would be to divide the resulting hazard maps into two types: “final hazard maps” and “preliminary hazard maps”. Only the use of numerical models provided final hazard maps, because they could incorporate different event magnitudes and they supplied output-values for intensity calculation. In contrast, empirical relationships and flow routing algorithms, or a combination of both, could be applied to create preliminary hazard maps. The present study only focussed on runout prediction methods. Other necessary tasks to complete the hazard assessment can be looked up in the “Guidelines for landslide susceptibility, hazard and risk zoning” included in this Special Issue.     

Applied hydrodynamic wave-resistance computation by Fourier transform

An applied Fourier transform computation for the hydrodynamic wave-resistance coefficient is shown, oriented to potential flows with a free surface and infinity depth. The presence of a ship-like body is simulated by its equivalent pressure disturbance imposed on the un-perturbed free surface, where a linearized free surface condition is used. The wave-resistance coefficient is obtained from the wave-height downstream. Two examples with closed solutions are considered: a submerged dipole, as a test-case, and a parabolic pressure distribution of compact support. In the three dimensional case, a dispersion relation is included which is a key resource for an inexpensive computation of the wave pattern far downstream like fifteen ship-lengths.     

False waves in wave records

A methodology to define discrete waves from free sea surface elevation time series is presented. The method allows an objective discrimination of false waves among small waves, avoiding the use of arbitrary criteria associated with zero-up-crossing, zero-down-crossing, or other definitions. The method uses the Hilbert Transform and the representation of the time series in the complex plane. A discrete wave corresponds to a 2π phase-advance in the complex plane. The waves between zero crossings which do not show 2π phase advance are considered false waves. Wave rider buoy records, measured off the west coast of Portugal, were employed in assessing some of the statistical implications of this methodology.     

Application of FLATModel, a 2D finite volume code, to debris flows in the northeastern part of the Iberian Peninsula

FLATModel is a two-dimensional shallow-water approximation code with corrections and modifications that create a simulation tool adapted to debris-flows behaviour. FLATModel uses the finite volume method with the numerical implementation of the Godunov scheme and includes correction terms regarding the effect of flow over high slopes and curvature. Additionally, the stop-and-go phenomenon, the basal entrainment and a correction regarding the front inclination of the final deposit are incorporated into FLATModel. In addition, different flow resistance laws were integrated in the numerical code including Bingham, Herschel–Bulkley and Voellmy fluid model. Firstly, our numerical model was validated using analytical solutions of a dam-break scenario and published data on a laboratory experiment. Secondly, three real events, which occurred in the northeastern part of the Iberian Peninsula, were back-calculated. Although field observations of the three events are not very detailed, the back-analyses revealed interesting patterns on the flow dynamics, and the numerical results generally showed good agreement with field data. Comparing the different flow resistance laws, the Voellmy fluid model presents the best behaviour regarding both the flow behaviour and the deposit characteristics. Preliminary simulation runs incorporating the effect of basal entrainment offered satisfactory results, although the final volume is rather sensitive on the selected friction angle of channel-bed material. The outcomes regarding the correction of the calculated front inclination of the final deposit showed that this implementation strongly improves the simulation results and better represents steep fronts of final deposits.     

A calibration system for superconducting gravimeters

A new method for the calibration of a superconducting gravity meter is described, in which a 273 Kg annular mass is placed around the meter and is moved up and down. The geometry of the apparatus is easy to model and the accuracy in the computation of the gravity variation induced by the mass, 6.7µgal, is limited only by the accuracy in the knowledge of value of the gravitational constant. Measurements done in 91 and 92 for the calibration of the instrument GWR-T015 are described. The calibration factor has been determined with a precision of about 0.3%.     

Rb/Sr dating of the Upper Proterozoic basement of Zambesia,Mozambique

The Proterozoic basement of the Province of Zambezia (Mozambique) is located near the southern end of the so called Mozambique Belt. Rb-Sr data derived from the analyses of 66 whole rocks and a few mineral samples bear witness of a magmatic and metamorphic episode referable to the Kibaran (ca. 1000 Ma) and helped to model the late-Proterozoic evolution as follows. At 1100–1050 Ma a pre-orogenic, calc-alkaline magmatism was responsible for the emplacement of now deformed granite and granodiorite. This phase was accompanied by rhyolitic volcanism. A subsequent peak of deformation and metamorphism occurred at 1000-950 Ma and was associated with a second generation of substantially post-tectonic granites. A thermal event resulting in the emplacement of a final generation of granites and pegmatites can be dated to around 500-450 Ma.The metamorphic basement has not preserved its protolith age, probably as a result of extensive isotopic homogeneization of Sr in Kibaran times. Its maximum time of differentiation from the mantle should be around 1600 Ma.The very low values of initial Sr isotopic ratios suggest the absence of Archean crustal precursors and support the rather young, yet pre-Pan-African evolution of the southernmost Mozambique Belt.

---

Zusammenfassung Das proterozoische Grundgebirge der Provinz Zambesi (Mosambique) befindet sich in der Nähe des südlichen Endes des sogenannten Mozambique Belt. Rb/Sr Daten aus 66 Ganzgesteins- und einigen Mineralanalysen zeugen von einer magmatischen und metamorphen Episode die sich dem Kibaran (ca. 1000 Ma) zuordnen läßt, und wie sich im folgenden zeigt, hilfreich bei der Modellisierung der spätproterozoischen Entwicklung ist. Um 1100 bis 1050 Ma erfolgte ein präorogener kalkalkalischer Magmatismus, bestehend aus Graniten und Granodioriten, die heute deformiert vorliegen. Diese Phase wurde von rhyolithischem Vulkanismus begleitet. Von 1000 bis 950 Ma folgte der Höhepunkt der Deformation und Metamorphose, begleitet von einer bedeutenden 2. Generation posttektonischer Granite. Bedingt durch ein weiteres thermales Ereignis intrudierte die letzte Granit- und Pegmatitgeneration vor 500 bis 450 Ma.Im metamorphen Grundgebirge ist das Eduktalter vermutlich aufgrund der starken isotopischen Homogenisierung von Sr während der Kibara-Zeit nicht erhalten. Seine maximale Differentiationszeit vom Mantel muß vor 1600 Ma gewesen sein.Die sehr niedrigen Gehalte der initialen Sr-Isotopenverhältnisse lassen die Abwesenheit von archaischen Krustenvorläufern vermuten, und unterstützen die sehr junge, bis jetzt Prä-Pan-Afrikanische Entwicklung des südlichsten Mozambique Belts.

---

Résumé Le socle protérozoïque de la province du Zambèze, au Mozambique, se trouve à proximité de l'extrémité sud du «Mozambique Belt». Des mesures Rb/Sr effectuées sur 66 roches totales et sur quelques minéraux témoignent d'un épisode magmatique et métamorphique attribuable au Kibarien (environ 1.000 Ma), ce qui permet de se représenter comme suit l'histoire de cette région au Protérozoïque supérieur. Il y a 1.100 à 1.050 Ma, la région a connu une phase de magmatisme calco-alcalin qui s'est traduite par la mise en place de granites et granodiorites, aujord'hui déformés. Cette phase s'est accompagnée d'un volcanisme rhyolitique. Elle a été suivie, entre 1.000 et 950 Ma, d'une activité tectono-métamorphique majeure, à laquelle est associée une deuxième et importante génération de granites, post-tectoniques. Ultérieurement, la région a connu encore un épisode thermique, daté à 500-450 Ma, responsable de la mise en place d'une dernière génération de granite et de pegmatite.Le socle métamorphique n'a pas gardé de témoins de l'âge de son protolithe, en raison probablement de la forte homogénisation isotopique du Sr à l'époque kibarienne. Lâge maximal de sa différenciation à partir du manteau doit se situer autour de 1.600 Ma.La valeur très basse des rapports isotopiques initiaux du Sr suggère l'absence de précurseurs crustaux archéens et incite à admettre une évolution relativement jeune — bien que pré-panafricaine — du «Mozambique Belt» méridional.

---

() .. . , Rb/Sr 66 , , Kibara, 1000 , , , . 1100 1050 - - , , . . 1000 950 , - . , 500 450 . , , , Kibara. , , 1600 my . - .