Annama H chondrite—Mineralogy,physical properties,cosmic ray exposure,and parent body history

The fall of the Annama meteorite occurred early morning (local time) on April 19, 2014 on the Kola Peninsula (Russia). Based on mineralogy and physical properties, Annama is a typical H chondrite. It has a high Ar‐Ar age of 4.4 Ga. Its cosmic ray exposure history is atypical as it is not part of the large group of H chondrites with a prominent 7–8 Ma peak in the exposure age histograms. Instead, its exposure age is within uncertainty of a smaller peak at 30 ± 4 Ma. The results from short‐lived radionuclides are compatible with an atmospheric pre‐entry radius of 30–40 cm. However, based on noble gas and cosmogenic radionuclide data, Annama must have been part of a larger body (radius >65 cm) for a large part of its cosmic ray exposure history. The ¹⁰Be concentration indicates a recent (3–5 Ma) breakup which may be responsible for the Annama parent body size reduction to 30–35 cm pre‐entry radius.     

Tutorial: unified 1D inversion of the acoustic reflection response

Acoustic inversion in one-dimension gives impedance as a function of travel time. Inverting the reflection response is a linear problem. Recursive methods, from top to bottom or vice versa, are known and use a fundamental wave field that is computed from the reflection response. An integral over the solution to the Marchenko equation, on the other hand, retrieves the impedance at any vertical travel time instant. It is a non-recursive method, but requires the zero-frequency value of the reflection response. These methods use the same fundamental wave field in different ways. Combining the two methods leads to a non-recursive scheme that works with finite-frequency bandwidth. This can be used for target-oriented inversion. When a reflection response is available along a line over a horizontally layered medium, the thickness and wave velocity of any layer can be obtained together with the velocity of an adjacent layer and the density ratio of the two layers. Statistical analysis over 1000 noise realizations shows that the forward recursive method and the Marchenko-type method perform well on computed noisy data.     

Orbit and origin of the LL7 chondrite Dishchii'bikoh (Arizona)

The trajectory and orbit of the LL7 ordinary chondrite Dishchii'bikoh are derived from low‐light video observations of a fireball first detected at 10:56:26 UTC on June 2, 2016. Results show a relatively steep ~21° inclined orbit and a short 1.13 AU semimajor axis. Following entry in Earth's atmosphere, the meteor luminosity oscillated corresponding to a meteoroid spin rate of 2.28 ± 0.02 rotations per second. A large fragment broke off at 44 km altitude. Further down, mass was lost to dust during flares at altitudes of 34, 29, and 25 km. Surviving meteorites were detected by Doppler weather radar and several small 0.9–29 g meteorites were recovered under the radar reflection footprint. Based on cosmogenic radionuclides and ground‐based radiometric observations, the Dishchii'bikoh meteoroid was 80 ± 20 cm in diameter assuming the density was 3.5 g/cm³. The meteoroid's collisional history confirms that the unusual petrologic class of LL7 does not require a different parent body than three previously observed LL chondrite falls. Dishchii'bikoh was ejected 11 Ma ago from parent body material that has a 4471 ± 6 Ma U‐Pb age, the same as that of Chelyabinsk (4452 ± 21 Ma). The distribution of the four known pre‐impact LL chondrite orbits is best matched by dynamical modeling if the source of LL chondrites is in the inner asteroid belt in a low inclined orbit, with the highly inclined Dishchii'bikoh being the result of interactions with Earth before impacting.     

3D Marchenko applications: implementation and examples

We implement the 3D Marchenko equations to retrieve responses to virtual sources inside the subsurface. For this, we require reflection data at the surface of the Earth that contain no free-surface multiples and are densely sampled in space. The required 3D reflection data volume is very large and solving the Marchenko equations requires a significant amount of computational cost. To limit the cost, we apply floating point compression to the reflection data to reduce their volume and the loading time from disk. We apply the Marchenko implementation to numerical reflection data to retrieve accurate Green's functions inside the medium and use these reflection data to apply imaging. This requires the simulation of many virtual source points, which we circumvent using virtual plane-wave sources instead of virtual point sources. Through this method, we retrieve the angle-dependent response of a source from a depth level rather than of a point. We use these responses to obtain angle-dependent structural images of the subsurface, free of contamination from wrongly imaged internal multiples. These images have less lateral resolution than those obtained using virtual point sources, but are more efficiently retrieved.     

Analytical and numerical analysis of tides and salinities in estuaries; part II: salinity distributions in prismatic and convergent tidal channels

Estuaries, commonly, are densely populated areas serving the needs of the inhabitants in multiple ways. Often the interests are conflicting and decisions need to be made by the local managers. Intake of fresh water for consumption, agricultural purposes or use by industries may take place within a region not far landward of the limit of salt intrusion. Human interventions (e.g. deepening of the navigation channels) or climate changes (sea level rise, reduction of the river discharge) can bring these intake locations within the reach of saline or brackish water and consequently endanger their function. To support policy and managerial decisions, a profound knowledge of processes associated with the salinity structure in estuaries is required. Although nowadays advanced numerical three-dimensional models are available that are able to cope with the complexity of the physics there is still a need for relatively simple tools for quick-scan actions in a pre-phase of a project or for instructive purposes. The analytical model described in this paper may serve these needs. It computes the maximum salinity distribution using the dispersion coefficient in the mouth as the only model parameter. The model has been calibrated using observational data in a large number of estuaries and experimental data in a tidal flume. The dispersion coefficient was successfully related to geometric and hydrodynamic parameters resulting in an expression that can be used for convergent estuaries as well as prismatic channels, see Eqs. 25a and 25b. Application of the model in a predictive mode showed its promising capabilities. Comparison with three-dimensional numerical models indicates that the channel geometry in the estuary mouth largely influences dispersive processes. The analytical model for salt intrusion may be used in combination with the analytical model for tidal propagation in convergent estuaries and tidal channels by Van Rijn (part I). In this way, input is obtained on the tidal velocity amplitude and the Chézy roughness following calibration of this model on tidal amplitudes along the estuary.     

A noble gas and cosmogenic radionuclide analysis of two ordinary chondrites from Almahata Sitta

Abstract– We present the results of a noble gas (He, Ne, Ar) and cosmogenic radionuclide (¹⁰Be, ²⁶Al, ³⁶Cl) analysis of two chondritic fragments (#A100, L4 and #25, H5) found in the Almahata Sitta strewn field in Sudan. We confirm their earlier attribution to the same fall as the ureilites dominating the strewn field, based on the following findings: (1) both chondrite samples indicate a preatmospheric radius of approximately 300 g cm^?2, consistent with the preatmospheric size of asteroid 2008 TC₃ that produced the Almahata Sitta strewn field; (2) both have, within error, a ²¹Ne/²⁶Al‐based cosmic ray exposure age of approximately 20 Ma, identical to the reported ages of Almahata Sitta ureilites; (3) both exhibit hints of ureilitic Ar in the trapped component. We discuss a possible earlier irradiation phase for the two fragments of approximately 10–20 Ma, visible only in cosmogenic ³⁸Ar. We also discuss the approximately 3.8 Ga (⁴He) and approximately 4.6 Ga (⁴⁰Ar) gas retention ages, measured in both chondritic fragments. These imply that the two chondrite fragments were incorporated into the ureilite host early in solar system evolution, and that the parent asteroid from which 2008 TC₃ is derived has not experienced a large break‐up event in the last 3.8 Ga.     

Finite element methods for variable density flow and solute transport

Saltwater intrusion into coastal freshwater aquifers is an ongoing problem that will continue to impact coastal freshwater resources as coastal populations increase. To effectively model saltwater intrusion, the impacts of increased salt content on fluid density must be accounted for to properly model saltwater/freshwater transition zones and sharp interfaces. We present a model for variable density fluid flow and solute transport where a conforming finite element method discretization with a locally conservative velocity post-processing method is used for the flow model and the transport equation is discretized using a variational multiscale stabilized conforming finite element method. This formulation provides a consistent velocity and performs well even in advection-dominated problems that can occur in saltwater intrusion modeling. The physical model is presented as well as the formulation of the numerical model and solution methods. The model is tested against several 2-D and 3-D numerical and experimental benchmark problems, and the results are presented to verify the code.