Debris flow magnitudes: a global catalog

The episodic occurrence of debris flow events in response to stochastic precipitation and wildfire events makes hazard prediction challenging. Previous work has shown that frequency-magnitude distributions of non-fire-related debris flows follow a power law, but less is known about the distribution of post-fire debris flows. As a first step in parameterizing hazard models, we use frequency-magnitude distributions and cumulative distribution functions to compare volumes of post-fire debris flows to non-fire-related debris flows. Due to the large number of events required to parameterize frequency-magnitude distributions, and the relatively small number of post-fire event magnitudes recorded in the literature, we collected data on 73 recent post-fire events in the field. The resulting catalog of 988 debris flow events is presented as an appendix to this article. We found that the empirical cumulative distribution function of post-fire debris flow volumes is composed of smaller events than that of non-fire-related debris flows. In addition, the slope of the frequency-magnitude distribution of post-fire debris flows is steeper than that of non-fire-related debris flows, evidence that differences in the post-fire environment tend to produce a higher proportion of small events. We propose two possible explanations: 1) post-fire events occur on shorter return intervals than debris flows in similar basins that do not experience fire, causing their distribution to shift toward smaller events due to limitations in sediment supply, or 2) fire causes changes in resisting and driving forces on a package of sediment, such that a smaller perturbation of the system is required in order for a debris flow to occur, resulting in smaller event volumes.

Stable isotope composition and ice rafted debris in early Late Pliocene sediments from Maud Rise, South Atlantic

We have conducted an integrated study of ice-rafted debris (IRD) and oxygen isotopes (measured on Cibicides, Globigerina bulloides, and Neogloboquadrina pachyderma, using identical samples). We used samples from the early Late Pliocene Gauss Chron from ODP Site 114-704 on the Meteor Rise in the subantarctic South Atlantic. During the early Gauss Chron, the oxygen isotopic ratios are generally up to 0.5?-0.6? less than their respective Holocene values. The lowest values in this record can accommodate a warming of about 2.5°C or a sea-level rise of about 50 m, but not both, and probably result from some warming and a small reduction in global ice volume. Starting with isotope stage MG2 [ 3.23 Ma on the Berggren et al. ( 1985) time scale; 3.38 on the Shackleton et al. ( 1995b) time scale] oxygen-isotopic values generally increase (and oscillate about a Holocene mean). The first significant IRD appears at the same time. There is a subsequent increase in IRD amounts upsection. In order to reach the site, this material must have been transported by large, tabular icebergs derived from Antarctic ice shelves or ice tongues, similar to occasional, large modern icebergs. This combined record suggests strongly that the Antarctic ice sheet was essentially intact; some warming at the drill site is indicated, but not a major reduction in ice-volume on Antarctica.

Ice-rafted debris in ODP Holes 177-1090D and 188-1165B

Constraining the nature of Antarctic Ice Sheet (AIS) response to major past climate changes may provide a window onto future ice response and rates of sea level rise. One approach to tracking AIS dynamics, and differentiating whole system versus potentially heterogeneous ice sheet sector changes, is to integrate multiple climate proxies for a specific time slice across widely distributed locations. This study presents new iceberg-rafted debris (IRD) data across the interval that includes Marine Isotope Stage 31 (MIS 31: 1.081-1.062 Ma, a span of ~19 kyr; Lisiecki and Raymo, 2005), which lies on the cusp of the mid-Brunhes climate transition (as glacial cycles shifted from ~41,000 yr to ~100,000 yr duration). Two sites are studied - distal Ocean Drilling Program (ODP) Leg 177 Site 1090 (Site 1090) in the eastern subantarctic sector of the South Atlantic Ocean, and proximal ODP Leg 188 Site 1165 (Site 1165), near Prydz Bay, in the Indian Ocean sector of the Antarctic margin. At each of these sites, MIS 31 is marked by the presence of the Jaramillo Subchron (0.988-1.072 Ma; Lourens et al., 2004) which provides a time-marker to correlate these two sites with relative precision. At both sites, records of multiple climate proxies are available to aid in interpretation. The presence of IRD in sediments from our study areas, which include garnets indicating a likely East Antarctic Ice Sheet (EAIS) origin, supports the conclusion that although the EAIS apparently withdrew significantly over MIS 31 in the Prydz Bay region and other sectors, some sectors of the EAIS must still have maintained marine margins capable of launching icebergs even through the warmest intervals. Thus, the EAIS did not respond in complete synchrony even to major climate changes such as MIS 31. Further, the record at Site 1090 (supported by records from other subantarctic locations) indicates that the glacial MIS 32 should be reduced to no more than a stadial, and the warm interval of Antarctic ice retreat that includes MIS 31 should be expanded to MIS 33-31. This revised warm interval lasted about 52 kyr, in line with several other interglacials in the benthic d18O records stack of Lisiecki and Raymo (2005), including the super-interglacials MIS 11 (duration of 50 kyr) and MIS 5 (duration of 59 kyr). The record from Antarctica-proximal Site 1165, when interpreted in accord with the record from ANDRILL-1B, indicates that in these southern high latitude sectors, ice sheet retreat and the effects of warming lasted longer than at Site 1090, perhaps until MIS 27. In the current interpretations of the age models of the proximal sites, ice sheet retreat began relatively slowly, and was not really evident until the start of MIS 31. In another somewhat more speculative interpretation, ice sheet retreat began noticeably with MIS 33, and accelerated during MIS 31. Ice sheet inertia (the lag-times in the large-scale responses of major ice sheets to a forcing) likely plays an important part in the timing and scale of these events in vulnerable sectors of the AIS.

The VLT/NaCo large program to probe the occurrence of exoplanets and brown dwarfs at wide orbits II. Survey description, results, and performances

Context. Young, nearby stars are ideal targets for direct imaging searches for giant planets and brown dwarf companions. After the first-imaged planet discoveries, vast efforts have been devoted to the statistical analysis of the occurence and orbital distributions of giant planets and brown dwarf companions at wide (>= 5-6 AU) orbits. Aims. In anticipation of the VLT/SPHERE planet-imager, guaranteed-time programs, we have conducted a preparatory survey of 86 stars between 2009 and 2013 to identify new faint comoving companions to ultimately analyze the occurence of giant planets and brown dwarf companions at wide (10-2000 AU) orbits around young, solar-type stars. Methods. We used NaCo at VLT to explore the occurrence rate of giant planets and brown dwarfs between typically 0.1 and 8 ''. Diffraction-limited observations in H-band combined with angular differential imaging enabled us to reach primary star-companion brightness ratios as small as 10(-6) at 1.5 ''. Repeated observations at several epochs enabled us to discriminate comoving companions from background objects. Results. During our survey, twelve systems were resolved as new binaries, including the discovery of a new white dwarf companion to the star HD8049. Around 34 stars, at least one companion candidate was detected in the observed field of view. More than 400 faint sources were detected; 90% of them were in four crowded fields. With the exception of HD8049 B, we did not identify any new comoving companions. The survey also led to spatially resolved images of the thin debris disk around HD61005 that have been published earlier. Finally, considering the survey detection limits, we derive a preliminary upper limit on the frequency of giant planets for the semi-major axes of [10, 2000] AU: typically less than 15% between 100 and 500 AU and less than 10% between 50 and 500 AU for exoplanets that are more massive than 5 M-Jup and 10 M-Jup respectively, if we consider a uniform input distribution and a confidence level of 95%. Conclusions. The results from this survey agree with earlier programs emphasizing that massive, gas giant companions on wide orbits around solar-type stars are rare. These results will be part of a broader analysis of a total of similar to 210 young, solar-type stars to bring further statistical constraints for theoretical models of planetary formation and evolution.

Orbital debris mitigation through deorbiting with passive electrodynamic drag

The increase of orbital debris and the consequent proliferation of smaller objects through fragmentation are driving the need for mitigation strategies. The issue is how to deorbit the satellite with an efficient system that does not impair drastically the propellant budget of the satellite and, consequently, reduces its operating life. We have been investigating, in the framework of a European-Community-funded project, a passive system that makes use of an electrodynamics tether to deorbit a satellite through Lorentz forces. The deorbiting system will be carried by the satellite itself at launch and deployed from the satellite at the end of its life. From that moment onward the system operates passively without requiring any intervention from the satellite itself. The paper summarizes the results of the analysis carried out to show the deorbiting performance of the system starting from different orbital altitudes and inclinations for a reference satellite mass. Results can be easily scaled to other satellite masses. The results have been obtained by using a high-fidelity computer model that uses the latest environmental routines for magnetic field, ionospheric density, atmospheric density and a gravity field model. The tether dynamics is modelled by considering all the main aspects of a real system as the tether flexibility and its temperature-dependent electrical conductivity. Temperature variations are computed by including all the major external and internal input fluxes and the thermal flux emitted from the tether. The results shows that a relatively compact and light system can carry out the complete deorbit of a relatively large satellite in a time ranging from a month to less than a year starting from high LEO with the best performance occurring at low orbital inclinations.

Analysis of tape tether survival in LEO against orbital debris

The low earth orbit (LEO) environment contains a large number of artificial debris, of which a significant portion is due to dead satellites and fragments of satellites resulted from explosions and in-orbit collisions. Deorbiting defunct satellites at the end of their life can be achieved by a successful operation of an Electrodynamic Tether (EDT) system. The effectiveness of an EDT greatly depends on the survivability of the tether, which can become debris itself if cut by debris particles; a tether can be completely cut by debris having some minimal diameter. The objective of this paper is to develop an accurate model using power laws for debris-size ranges, in both ORDEM2000 and MASTER2009 debris flux models, to calculate tape tether survivability. The analytical model, which depends on tape dimensions (width, thickness) and orbital parameters (inclinations, altitudes) is then verified with fully numerical results to compare for different orbit inclinations, altitudes and tape width for both ORDEM2000 and MASTER2009 flux data.

Inception of debris avalanches: remarks on geomechanical modelling

Debris avalanches are complex phenomena due to the variety of mechanisms that control the failure stage and the avalanche formation. Regarding these issues, in the literature, either field evidence or qualitative interpretations can be found while few experimental laboratory tests and rare examples of geomechanical modelling are available for technical and/or scientific purposes. As a contribution to the topic, the paper firstly highlights as the problem can be analysed referring to a unique mathematical framework from which different modelling approaches can be derived based on limit equilibrium method (LEM), finite element method (FEM), or smooth particle hydrodynamics (SPH). Potentialities and limitations of these approaches are then tested for a large study area where huge debris avalanches affected shallow deposits of pyroclastic soils (Sarno-Quindici, Southern Italy). The numerical results show that LEM as well as uncoupled and coupled stress–strain FEM analyses are able to individuate the major triggering mechanisms. On the other hand, coupled SPH analyses outline the relevance of erosion phenomena, which can modify the kinematic features of debris avalanches in their source areas, i.e. velocity, propagation patterns and later spreading of the unstable mass. As a whole, the obtained results encourage the application of the introduced approaches to further analyse real cases in order to enhance the current capability to forecast the inception of these dangerous phenomena.

Bare-tape scaling laws for de-orbit missions in a space debris environment

El entorno espacial actual hay un gran numero de micro-meteoritos y basura espacial generada por el hombre, lo cual plantea un riesgo para la seguridad de las operaciones en el espacio. La situación se agrava continuamente a causa de las colisiones de basura espacial en órbita, y los nuevos lanzamientos de satélites. Una parte significativa de esta basura son satélites muertos, y fragmentos de satélites resultantes de explosiones y colisiones de objetos en órbita. La mitigación de este problema se ha convertido en un tema de preocupación prioritario para todas las instituciones que participan en operaciones espaciales. Entre las soluciones existentes, las amarras electrodinámicas (EDT) proporcionan un eficiente dispositivo para el rápido de-orbitado de los satélites en órbita terrestre baja (LEO), al final de su vida útil. El campo de investigación de las amarras electrodinámicas (EDT) ha sido muy fructífero desde los años 70. Gracias a estudios teóricos, y a misiones para la demostración del funcionamiento de las amarras en órbita, esta tecnología se ha desarrollado muy rápidamente en las últimas décadas. Durante este período de investigación, se han identificado y superado múltiples problemas técnicos de diversa índole. Gran parte del funcionamiento básico del sistema EDT depende de su capacidad de supervivencia ante los micro-meteoritos y la basura espacial. Una amarra puede ser cortada completamente por una partícula cuando ésta tiene un diámetro mínimo. En caso de corte debido al impacto de partículas, una amarra en sí misma, podría ser un riesgo para otros satélites en funcionamiento. Por desgracia, tras varias demostraciones en órbita, no se ha podido concluir que este problema sea importante para el funcionamiento del sistema. En esta tesis, se presenta un análisis teórico de la capacidad de supervivencia de las amarras en el espacio. Este estudio demuestra las ventajas de las amarras de sección rectangular (cinta), en cuanto a la probabilidad de supervivencia durante la misión, frente a las amarras convencionales (cables de sección circular). Debido a su particular geometría (longitud mucho mayor que la sección transversal), una amarra puede tener un riesgo relativamente alto de ser cortado por un único impacto con una partícula de pequeñas dimensiones. Un cálculo analítico de la tasa de impactos fatales para una amarra cilindrica y de tipo cinta de igual longitud y masa, considerando el flujo de partículas de basura espacial del modelo ORDEM2000 de la NASA, muestra mayor probabilidad de supervivencia para las cintas. Dicho análisis ha sido comparado con un cálculo numérico empleando los modelos de flujo el ORDEM2000 y el MASTER2005 de ESA. Además se muestra que, para igual tiempo en órbita, una cinta tiene una probabilidad de supervivencia un orden y medio de magnitud mayor que una amarra cilindrica con igual masa y longitud. Por otra parte, de-orbitar una cinta desde una cierta altitud, es mucho más rápido, debido a su mayor perímetro que le permite capturar más corriente. Este es un factor adicional que incrementa la probabilidad de supervivencia de la cinta, al estar menos tiempo expuesta a los posibles impactos de basura espacial. Por este motivo, se puede afirmar finalmente y en sentido práctico, que la capacidad de supervivencia de la cinta es bastante alta, en comparación con la de la amarra cilindrica. El segundo objetivo de este trabajo, consiste en la elaboración de un modelo analítico, mejorando la aproximación del flujo de ORDEM2000 y MASTER2009, que permite calcular con precisión, la tasa de impacto fatal al año para una cinta en un rango de altitudes e inclinaciones, en lugar de unas condiciones particulares. Se obtiene el numero de corte por un cierto tiempo en función de la geometría de la cinta y propiedades de la órbita. Para las mismas condiciones, el modelo analítico, se compara con los resultados obtenidos del análisis numérico. Este modelo escalable ha sido esencial para la optimización del diseño de la amarra para las misiones de de-orbitado de los satélites, variando la masa del satélite y la altitud inicial de la órbita. El modelo de supervivencia se ha utilizado para construir una función objetivo con el fin de optimizar el diseño de amarras. La función objectivo es el producto del cociente entre la masa de la amarra y la del satélite y el numero de corte por un cierto tiempo. Combinando el modelo de supervivencia con una ecuación dinámica de la amarra donde aparece la fuerza de Lorentz, se elimina el tiempo y se escribe la función objetivo como función de la geometría de la cinta y las propietades de la órbita. Este modelo de optimización, condujo al desarrollo de un software, que esta en proceso de registro por parte de la UPM. La etapa final de este estudio, consiste en la estimación del número de impactos fatales, en una cinta, utilizando por primera vez una ecuación de límite balístico experimental. Esta ecuación ha sido desarollada para cintas, y permite representar los efectos tanto de la velocidad de impacto como el ángulo de impacto. Los resultados obtenidos demuestran que la cinta es altamente resistente a los impactos de basura espacial, y para una cinta con una sección transversal definida, el número de impactos críticos debidos a partículas no rastreables es significativamente menor. ABSTRACT The current space environment, consisting of man-made debris and tiny meteoroids, poses a risk to safe operations in space, and the situation is continuously deteriorating due to in-orbit debris collisions and to new satellite launches. Among these debris a significant portion is due to dead satellites and fragments of satellites resulted from explosions and in-orbit collisions. Mitigation of space debris has become an issue of first concern for all the institutions involved in space operations. Bare electrodynamic tethers (EDT) can provide an efficient mechanism for rapid de-orbiting of defunct satellites from low Earth orbit (LEO) at end of life. The research on EDT has been a fruitful field since the 70’s. Thanks to both theoretical studies and in orbit demonstration missions, this technology has been developed very fast in the following decades. During this period, several technical issues were identified and overcome. The core functionality of EDT system greatly depends on their survivability to the micrometeoroids and orbital debris, and a tether can become itself a kind of debris for other operating satellites in case of cutoff due to particle impact; however, this very issue is still inconclusive and conflicting after having a number of space demonstrations. A tether can be completely cut by debris having some minimal diameter. This thesis presents a theoretical analysis of the survivability of tethers in space. The study demonstrates the advantages of tape tethers over conventional round wires particularly on the survivability during the mission. Because of its particular geometry (length very much larger than cross-sectional dimensions), a tether may have a relatively high risk of being severed by the single impact of small debris. As a first approach to the problem, survival probability has been compared for a round and a tape tether of equal mass and length. The rates of fatal impact of orbital debris on round and tape tether, evaluated with an analytical approximation to debris flux modeled by NASA’s ORDEM2000, shows much higher survival probability for tapes. A comparative numerical analysis using debris flux model ORDEM2000 and ESA’s MASTER2005 shows good agreement with the analytical result. It also shows that, for a given time in orbit, a tape has a probability of survival of about one and a half orders of magnitude higher than a round tether of equal mass and length. Because de-orbiting from a given altitude is much faster for the tape due to its larger perimeter, its probability of survival in a practical sense is quite high. As the next step, an analytical model derived in this work allows to calculate accurately the fatal impact rate per year for a tape tether. The model uses power laws for debris-size ranges, in both ORDEM2000 and MASTER2009 debris flux models, to calculate tape tether survivability at different LEO altitudes. The analytical model, which depends on tape dimensions (width, thickness) and orbital parameters (inclinations, altitudes) is then compared with fully numerical results for different orbit inclinations, altitudes and tape width for both ORDEM2000 and MASTER2009 flux data. This scalable model not only estimates the fatal impact count but has proved essential in optimizing tether design for satellite de-orbit missions varying satellite mass and initial orbital altitude and inclination. Within the frame of this dissertation, a simple analysis has been finally presented, showing the scalable property of tape tether, thanks to the survivability model developed, that allows analyze and compare de-orbit performance for a large range of satellite mass and orbit properties. The work explicitly shows the product of tether-to-satellite mass-ratio and fatal impact count as a function of tether geometry and orbital parameters. Combining the tether dynamic equation involving Lorentz drag with space debris impact survivability model, eliminates time from the expression. Hence the product, is independent of tether de-orbit history and just depends on mission constraints and tether length, width and thickness. This optimization model finally led to the development of a friendly software tool named BETsMA, currently in process of registration by UPM. For the final step, an estimation of fatal impact rate on a tape tether has been done, using for the first time an experimental ballistic limit equation that was derived for tapes and accounts for the effects of both the impact velocity and impact angle. It is shown that tape tethers are highly resistant to space debris impacts and considering a tape tether with a defined cross section, the number of critical events due to impact with non-trackable debris is always significantly low.

Regenerative morphological traits in a woody species community in Tumbesian tropical dry forest.

The study of functional morphological traits enables us to know fundamental aspects of the dynamics of plant communities in local and global habitats. Regenerative morphological traits play an important role in defining plant history and ecological behavior. Seed and fruit characteristics determine to a large extent the patterns for dispersal, germination, establishment and seedling recruitment a given species exhibits on its natural habitat. Despite their prominent role, seed and fruit traits have been poorly studied at the community level of woody plant species in neo-tropical dry forests. In the present study we aimed at i) evaluate the functional role of morphological traits of seeds, fruits and embryo in woody plant species; ii) determine which are the morphological patterns present in seeds collected from the community of woody species that occur in neo-tropical dry forests; and iii) compare woody plant species seed mass values comparatively between neo-tropical dry and tropical forests. To do so, mature seeds were collected from 79 plant species that occur in the Tumbesian forest of Southwest Ecuador. The studied species included the 42 and 37 most representative tree and shrubbery species of the Tumbesian forest respectively. A total of 18 morphological traits (seven quantitative and 11 qualitative) were measured and evaluated in the seeds, fruits and embryos of the selected species, and we compared the seeds mass with other forest types. Our results showed a huge heterogeneity among traits values in the studied species. Seed mass, volume and number were the traits that vary the most at the community level, i.e. seed length ranged from 1.3 to 39 mm, and seed width from 0.6 to 25 mm. Only six embryo types were found among the 79 plant species. In 40 % of the cases, fully developed inverted embryos with large and thick cotyledons to store considerable amount of nutrients were recorded. We concluded that highly variable and functionally complementary morphological traits occur among the studied woody plants of the dry Tumbesian forest. The latter favors a plethora of behavioral mechanisms to coexist among woody species of the dry forest in response to the environmental stress that is typical of arid areas.

(Table T1) Ice rafted debris characteristics of ODP Site 177-1092

We have carried out a multiphase analysis of samples from ODP Site 177-1092, Meteor Rise, subantarctic South Atlantic. Samples were analyzed for ice-rafted debris (IRD [see Table T1]) and stable isotopes from benthic foraminifera [see Murphy et al., 2002, doi:10.1016/S0031-0182(01)00495-3]. Both analyses were performed on the same samples. Additional work was performed to identify the paleomagnetic stratigraphy. The analyzed samples range in age from about 2.6(?) Ma to 4.6 Ma, a time span that saw considerable global warmth, but witnessed overall global refrigeration and the transition to truly bipolar glaciations. IRD arrived frequently during the Early and early Late Pliocene, but only as 'background rafting' (occasional grains per sample). The first identifiable IRD above background rafting is associated with marine isotope stage (MIS) KM4 (~3.18 Ma). Successive IRD peaks become larger, the same pattern as noted at nearby Site 114-704. A very large peak near the top of the record, approximately 2.8 Ma, is considered to represent a hiatus. Peaks below 51.3 meters composite depth (mcd) coincide with positive excursions of the oxygen isotopic record, and with negative excursions of the carbon isotopic curve, a pattern also noted at Site 114-704. However, the reasonably large IRD peak at 51 mcd (tentatively identified with MIS G11) coincides with a positive excursion on the carbon isotopic curve and negative excursion on the oxygen isotopic curve. This relationship suggests a northern hemisphere interglacial, rising sea level, destabilization of the Antarctic margin, and delivery of Antarctic icebergs to the Southern Ocean. Such a mechanism has recently been suggested by Kanfoush et al. (2000, doi:10.1126/science.288.5472.1815) for latest Pleistocene stadial/interstadial oscillations. Here we suggest that such a mechanism may have been in place on glacial/interglacial time scales as early as the Late Pliocene.

Ice rafted debris distribution in sediment core CRP-1

Lonestone abundances in CRP-1 were investigated using three methods: core examination at Cape Roberts Camp, analysis of digital core images and follow-up core examination. For all images of split-core, we determined size and depth of every detectable lonestone larger than 3 mm. Lonestone abundance decreases exponentially with clast size. Although no significant depth-dependent variations in lonestone size distribution were detected, a strong 0.5-0.7 m abundance periodicity, of unknown origin, is evident within diamicts. Lonestone volume percentage was estimated from size distribution: most size classes contribute approximately the same volume to the total. Sizes >16 mm have rare enough lonestones that their counts are nonrepresentative when based on short intervals of split core. This problem does not affect total counts significantly, but the volume analysis needs to be confined to <= 6 mm lonestones to avoid instability induced by rare and nonrepresentative larger lonestones. If lonestone abundance can be used as an indicator of glacial proximity, then our CRP-1 lonestone abundance logs confirm the overall character of previously inferred variations in relative distance to the ice margin. Large-scale changes in lonestone abundance also reflect the CRP-1 sequence stratigraphy, with individual sequences generally characterised by basal lonestone-rich diamict overlain by lonestone-poor sands and muds. The relationship between glacial proximity and lonestone abundance within diamicts and within sand-mud intervals is, however, less certain. For example, two or three gradual lonestone increases may indicate regressions during glacial advances, in contrast to the more common CRP-l pattern of dominantly transgressive sequences.

Results of analyses of unmelted meteoritic debris collected from Eltanin ejecta in three sediment core from POLARSTERN cruise ANT-XII/4

A total of 1.7 g of unmelted meteorite particles have been recovered from FS Polarstern piston cores collected on expedition ANT XII/4 that contain ejecta from the Eltanin impact event. Most of the mass (1.2 g) is a large, single specimen that is a polymict breccia, similar in mineralogy and chemistry to howardites or the silicate fraction of mesosiderites. Most of the remaining mass is in several large individual pieces (20-75mg each) that are polymict breccias, fragments dominated by pyroxene, and an igneous rock fragment. The latter has highly fractionated REE, similar to those reported in mafic clasts from mesosiderites. Other types of specimens identified include fragments dominated by maskelynite or olivine. These pieces of the projectile probably survived impact by being blown off the back surface of the Eltanin asteroid during its impact into the Bellingshausen Sea.

Eulerian-Lagrangian two phase debris flow model

The main objective of this work is to develop a quasi three-dimensional numerical model to simulate stony debris flows, considering a continuum fluid phase, composed by water and fine sediments, and a non-continuum phase including large particles, such as pebbles and boulders. Large particles are treated in a Lagrangian frame of reference using the Discrete Element Method, the fluid phase is based on the Eulerian approach, using the Finite Element Method to solve the depth-averaged Navier-Stokes equations in two horizontal dimensions. The particle’s equations of motion are in three dimensions. The model simulates particle-particle collisions and wall-particle collisions, taking into account that particles are immersed in a fluid. Bingham and Cross rheological models are used for the continuum phase. Both formulations provide very stable results, even in the range of very low shear rates. Bingham formulation is better able to simulate the stopping stage of the fluid when applied shear stresses are low. Results of numerical simulations have been compared with data from laboratory experiments on a flume-fan prototype. Results show that the model is capable of simulating the motion of big particles moving in the fluid flow, handling dense particulate flows and avoiding overlap among particles. An application to simulate debris flow events that occurred in Northern Venezuela in 1999 shows that the model could replicate the main boulder accumulation areas that were surveyed by the USGS. Uniqueness of this research is the integration of mud flow and stony debris movement in a single modeling tool that can be used for planning and management of debris flow prone areas.

Eulerian-Lagrangian Two Phase Debris Flow Model

The main objective of this work is to develop a quasi three-dimensional numerical model to simulate stony debris flows, considering a continuum fluid phase, composed by water and fine sediments, and a non-continuum phase including large particles, such as pebbles and boulders. Large particles are treated in a Lagrangian frame of reference using the Discrete Element Method, the fluid phase is based on the Eulerian approach, using the Finite Element Method to solve the depth-averaged Navier–Stokes equations in two horizontal dimensions. The particle’s equations of motion are in three dimensions. The model simulates particle-particle collisions and wall-particle collisions, taking into account that particles are immersed in a fluid. Bingham and Cross rheological models are used for the continuum phase. Both formulations provide very stable results, even in the range of very low shear rates. Bingham formulation is better able to simulate the stopping stage of the fluid when applied shear stresses are low. Results of numerical simulations have been compared with data from laboratory experiments on a flume-fan prototype. Results show that the model is capable of simulating the motion of big particles moving in the fluid flow, handling dense particulate flows and avoiding overlap among particles. An application to simulate debris flow events that occurred in Northern Venezuela in 1999 shows that the model could replicate the main boulder accumulation areas that were surveyed by the USGS. Uniqueness of this research is the integration of mud flow and stony debris movement in a single modeling tool that can be used for planning and management of debris flow prone areas.