Middle Eocene Fe and Ca counts, Ca/Fe ratios, revised splices and astronomical age models form different Sites of IODP Exp342

Accurate age models are a tool of utmost important in paleoclimatology. Constraining the rate and pace of past climate change are at the core of paleoclimate research, as such knowledge is crucial to our understanding of the climate system. Indeed, it allows for the disentanglement of the various drivers of climate change. The scarcity of highly resolved sedimentary records from the middle Eocene (Bartonian - Lutetian Stages; 47.8 - 37.8 Ma) has led to the existence of the "Eocene astronomical time scale gap" and hindered the establishment of a comprehensive astronomical time scale (ATS) for the entire Cenozoic. Sediments from the Newfoundland Ridge drilled during Integrated Ocean Drilling Program (IODP) Expedition 342 span the Eocene gap at an unprecedented stratigraphic resolution with carbonate bearing sediments. Moreover, these sediments exhibit cyclic lithological changes that allow for an astronomical calibration of geologic time. In this study, we use the dominant obliquity imprint in XRF-derived calcium-iron ratio series (Ca/Fe) from three sites drilled during IODP Expedition 342 (U1408, U1409, U1410) to construct a floating astrochronology. We then anchor this chronology to numerical geological time by tuning 173-kyr cycles in the amplitude modulation pattern of obliquity to an astronomical solution. This study is one of the first to use the 173-kyr obliquity amplitude cycle for astrochronologic purposes, as previous studies primarily use the 405-kyr long eccentricity cycle as a tuning target to calibrate the Paleogene geologic time scale. We demonstrate that the 173-kyr cycles in obliquity's amplitude are stable between 40 and 50 Ma, which means that one can use the 173-kyr cycle for astrochronologic calibration in the Eocene. Our tuning provides new age estimates for magnetochron reversals C18n.1n - C21r and a stratigraphic framework for key sites from Expedition 342 for the Eocene. Some disagreements emerge when we compare our tuning for the interval between C19r and C20r with previous tuning attempts from the South Atlantic. We therefore present a revision of the original astronomical interpretations for the latter records, so that the various astrochronologic age models for the middle Eocene in the North- and South-Atlantic are consistent.

Compilation of organic carbon distribution and sedimentology in the surface sediments on the continental margin offshore southwestern Africa

In this study we demonstrate the relevance of lateral particle transport in nepheloid layers for organic carbon (OC) accumulation and burial across high-productive continental margins. We present geochemical data from surface sediments and suspended particles in the bottom nepheloid layer (BNL) from the most productive coastal upwelling area of the modern ocean, the Benguela upwelling system offshore southwest Africa. Interpretation of depositional patterns and comparison of downslope trends in OC content, organic matter composition, and 14C age between suspended particles and surface sediments indicate that lateral particle transport is the primary mechanism controlling supply and burial of OC. We propose that effective seaward particle transport primarily along the BNL is a key process that promotes and maintains local high sedimentation rates, ultimately causing high preservation of OC in a depocenter on the upper slope offshore Namibia. As lateral transport efficiently displaces areas of enhanced OC burial from maximum production at highly productive continental margins, vertical particle flux models do not sufficiently explain the relationship between primary production and shallow-marine OC burial. On geologic time scales, the widest distribution and strongest intensity of lateral particle transport is expected during periods of rapid sea-level change. At times in the geologic past, widespread downslope lateral transport of OC thus may have been a primary driver of enhanced OC burial at deeper continental slopes and abyssal basins.

Amino acid racemization for several DSDP samples (Table 1)

Kinetic parameters for the epimerization of isoleucine in multispecific foraminiferal asemblages were used to establish the effects of burial depth and the geothermal gradient on the extent of reaction. It was observed that with a little as thirty meters of burial in a normal thermal regime there were differences between the extent of epimerization measured and that which would have been predicted for thermal equilibrium with bottom water temperatures. As would be expected, these differences are greatest when the heat flow (the geothermal gradient) and/or the sedimentation rates are highest. These effects were observed in most of the DSDP samples studied, and have been used to estimate the average heat flux since the time of sample deposition. Occasional anomalous effects were observed which could not be related to past or present heat flux. These were determined to be due to such geologic occurrences as slumping and reworking or to recent sample contamination. Other problems emerged related to bottom water temperatures including changes over geologic time which are unknown and could not be deduced. Thus, the presence of epimerization anomalies in DSDP cores as noted above limits the effectiveness of amino acid geochronology in such cores, unless these anomalies can be recognized as ab initio.

Anatomy and variability of "Cuvierichelys parisiensis", a geoemydid turtle that crosses the Eocene-Oligocene boundary in Belgium

Abundant material of turtles from the early Oligocene site of Boutersem-TGV (Boutersem, Belgium), is presented here. No information on the turtles found there was so far available. All the turtle specimens presented here are attributable to a single freshwater taxon that is identified as a member of Geoemydidae, Cuvierichelys. It is the first representative of the ‘Palaeochelys s. l.–Mauremys’ group recognized in the Belgian Paleogene record. This material, which allows to know all the elements of both the carapace and the plastron of the taxon, cannot be attributed to the only species of the genus Cuvierichelys so far identified in the Oligocene, the Spanish form Cuvierichelys iberica. The taxon from Boutersem is recognized as Cuvierichelys parisiensis. Thus, both the paleobiogeographic and the biostratigraphic distributions of Cuvierichelys parisiensis are extended, its presence being confirmed for the first time outside the French Eocene record. The validity of some European forms is refuted, and several characters previously proposed as different between Cuvierichelys iberica and Cuvierichelys parisiensis are recognized as subjected to intraspecific variability.

Revisiting the Late Jurassic-Early Cretaceous of the NW South Iberian Basin: new ages and sedimentary environments

The study of the Upper Jurassic-Lower Cretaceous deposits (Higueruelas, Villar del Arzobispo and Aldea de Cortés Formations) of the South Iberian Basin (NW Valencia, Spain) reveals new stratigraphic and sedimentological data, which have significant implications on the stratigraphic framework, depositional environments and age of these units. The Higueruelas Fm was deposited in a mid-inner carbonate platform where oncolitic bars migrated by the action of storms and where oncoid production progressively decreased towards the uppermost part of the unit. The overlying Villar del Arzobispo Fm has been traditionally interpreted as an inner platform-lagoon evolving into a tidal-flat. Here it is interpreted as an inner-carbonate platform affected by storms, where oolitic shoals protected a lagoon, which had siliciclastic inputs from the continent. The Aldea de Cortés Fm has been previously interpreted as a lagoon surrounded by tidal-flats and fluvial-deltaic plains. Here it is reinterpreted as a coastal wetland where siliciclastic muddy deposits interacted with shallow fresh to marine water bodies, aeolian dunes and continental siliciclastic inputs. The contact between the Higueruelas and Villar del Arzobispo Fms, classically defined as gradual, is also interpreted here as rapid. More importantly, the contact between the Villar del Arzobispo and Aldea de Cortés Fms, previously considered as unconformable, is here interpreted as gradual. The presence of Alveosepta in the Villar del Arzobispo Fm suggests that at least part of this unit is Kimmeridgian, unlike the previously assigned Late Tithonian-Middle Berriasian age. Consequently, the underlying Higueruelas Fm, previously considered Tithonian, should not be younger than Kimmeridgian. Accordingly, sedimentation of the Aldea de Cortés Fm, previously considered Valangian-Hauterivian, probably started during the Tithonian and it may be considered part of the regressive trend of the Late Jurassic-Early Cretaceous cycle. This is consistent with the dinosaur faunas, typically Jurassic, described in the Villar del Arzobispo and Aldea de Cortés Fms.

Using a Combination of Phylogenetic Comparative and Paleobotanical Methods to Elucidate Patterns of Lineage Selection in Rosales (Plantae: Anthophyta)

An increasing focus in evolutionary biology is on the interplay between mesoscale ecological and evolutionary processes such as population demographics, habitat tolerance, and especially geographic distribution, as potential drivers responsible for patterns of diversification and extinction over geologic time. However, few studies to date connect organismal processes such as survival and reproduction through mesoscale patterns to long-term macroevolutionary trends. In my dissertation, I investigate how mechanism of seed dispersal, mediated through geographic range size, influences diversification rates in the Rosales (Plantae: Anthophyta). In my first chapter, I validate the phylogenetic comparative methods that I use in my second and third chapters. Available state speciation and extinction (SSE) models assumptions about evolution known to be false through fossil data. I show, however, that as long as net diversification rates remain positive – a condition likely true for the Rosales – these violations of SSE’s assumptions do not cause significantly biased results. With SSE methods validated, my second chapter reconstructs three associations that appear to increase diversification rate for Rosalean genera: (1) herbaceous habit; (2) a three-way interaction combining animal dispersal, high within-genus species richness, and geographic range on multiple continents; (3) a four-way interaction combining woody habit with the other three characteristics of (2). I suggest that the three- and four-way interactions represent colonization ability and resulting extinction resistance in the face of late Cenozoic climate change; however, there are other possibilities as well that I hope to investigate in future research. My third chapter reconstructs the phylogeographic history of the Rosales using both non-fossil-assisted SSE methods as well as fossil-informed traditional phylogeographic analysis. Ancestral state reconstructions indicate that the Rosaceae diversified in North America while the other Rosalean families diversified elsewhere, possibly in Eurasia. SSE is able to successfully identify groups of genera that were likely to have been ancestrally widespread, but has poorer taxonomic resolution than methods that use fossil data. In conclusion, these chapters together suggest several potential causal links between organismal, mesoscale, and geologic scale processes, but further work will be needed to test the hypotheses that I raise here.

Biogeography in deep time - What do phylogenetics, geology, and paleoclimate tell us about early platyrrhine evolution?

Molecular data have converged on a consensus about the genus-level phylogeny of extant platyrrhine monkeys, but for most extinct taxa and certainly for those older than the Pleistocene we must rely upon morphological evidence from fossils. This raises the question as to how well anatomical data mirror molecular phylogenies and how best to deal with discrepancies between the molecular and morphological data as we seek to extend our phylogenies to the placement of fossil taxa. Here I present parsimony-based phylogenetic analyses of extant and fossil platyrrhines based on an anatomical dataset of 399 dental characters and osteological features of the cranium and postcranium. I sample 16 extant taxa (one from each platyrrhine genus) and 20 extinct taxa of platyrrhines. The tree structure is constrained with a "molecular scaffold" of extant species as implemented in maximum parsimony using PAUP with the molecular-based 'backbone' approach. The data set encompasses most of the known extinct species of platyrrhines, ranging in age from latest Oligocene (∼26 Ma) to the Recent. The tree is rooted with extant catarrhines, and Late Eocene and Early Oligocene African anthropoids. Among the more interesting patterns to emerge are: (1) known early platyrrhines from the Late Oligocene through Early Miocene (26-16.5Ma) represent only stem platyrrhine taxa; (2) representatives of the three living platyrrhine families first occur between 15.7 Ma and 13.5 Ma; and (3) recently extinct primates from the Greater Antilles (Cuba, Jamaica, Hispaniola) are sister to the clade of extant platyrrhines and may have diverged in the Early Miocene. It is probable that the crown platyrrhine clade did not originate before about 20-24 Ma, a conclusion consistent with the phylogenetic analysis of fossil taxa presented here and with recent molecular clock estimates. The following biogeographic scenario is consistent with the phylogenetic findings and climatic and geologic evidence: Tropical South America has been a center for platyrrhine diversification since platyrrhines arrived on the continent in the middle Cenozoic. Platyrrhines dispersed from tropical South America to Patagonia at ∼25-24 Ma via a "Paraná Portal" through eastern South America across a retreating Paranense Sea. Phylogenetic bracketing suggests Antillean primates arrived via a sweepstakes route or island chain from northern South America in the Early Miocene, not via a proposed land bridge or island chain (GAARlandia) in the Early Oligocene (∼34 Ma). Patagonian and Antillean platyrrhines went extinct without leaving living descendants, the former at the end of the Early Miocene and the latter within the past six thousand years. Molecular evidence suggests crown platyrrhines arrived in Central America by crossing an intermittent connection through the Isthmus of Panama at or after 3.5Ma. Any more ancient Central American primates, should they be discovered, are unlikely to have given rise to the extant Central American taxa in situ.

Space-time relation of post-collisional granitic magmatism in Santa Catarina, southern Brazil: U-Pb LA-MC-ICP-MS zircon geochronology of coeval mafic-felsic magmatism related to the Major Gercino Shear Zone

LA-MC-ICP-MS U-Pb zircon dating was performed on syntectonic, early post-collisional granitic and associated mafic rocks that are intrusive in the Brusque Metamorphic Complex and in the Florianopolis Batholith, major tectonic domains separated by the Neoproterozoic Major Gercino Shear Zone (MGSZ) in south Brazil. The inferred ages of magmatic crystallization are consistent with field relationships, and show that the syntectonic granites from both domains are similar, with ages around 630-620 Ma for high-K calc-alkaline metaluminous granites and ca. 610 Ma for slightly peraluminous granites. Although ca. 650 Ma inherited zircon components are identified in granites from both domains, important contrasts on the crustal architecture in each domain are revealed by the patterns of zircon inheritance, indicating different crustal sources for the granites in each domain. The granites from the southern domain (Floriandpolis Batholith) have essentially Neoproterozoic (650-700 Ma and 900-950 Ma) inheritance; with a single 2.0-2.2 Ga inherited age obtained in the peraluminous Mariscal Granite. In the northern Brusque Metamorphic Complex, the metaluminous Rio Pequeno Granite and associated mafic rocks have scarce inherited cores with ages around 1.65 Ga, whereas the slightly peraluminous Serra dos Macacos Granite has abundant Paleoproterozoic (1.8-2.2 Ga) and Archean (2.9-3.4 Ga) inherited zircons. Our results are consistent with the hypothesis that the MGSZ separates domains with distinct geologic evolution; however, the contemporaneity of 630-610 Ma granitic magmatism with similar structural and geochemical patterns on both sides of this major shear zone indicates that these domains were already part of a single continental mass at 630 Ma, reinforcing the post-collisional character of these granites. (C) 2012 Elsevier B.V. All rights reserved.

Analysis and developments of uncertainty processors for real time flood forecasting

The hydrologic risk (and the hydro-geologic one, closely related to it) is, and has always been, a very relevant issue, due to the severe consequences that may be provoked by a flooding or by waters in general in terms of human and economic losses. Floods are natural phenomena, often catastrophic, and cannot be avoided, but their damages can be reduced if they are predicted sufficiently in advance. For this reason, the flood forecasting plays an essential role in the hydro-geological and hydrological risk prevention. Thanks to the development of sophisticated meteorological, hydrologic and hydraulic models, in recent decades the flood forecasting has made a significant progress, nonetheless, models are imperfect, which means that we are still left with a residual uncertainty on what will actually happen. In this thesis, this type of uncertainty is what will be discussed and analyzed. In operational problems, it is possible to affirm that the ultimate aim of forecasting systems is not to reproduce the river behavior, but this is only a means through which reducing the uncertainty associated to what will happen as a consequence of a precipitation event. In other words, the main objective is to assess whether or not preventive interventions should be adopted and which operational strategy may represent the best option. The main problem for a decision maker is to interpret model results and translate them into an effective intervention strategy. To make this possible, it is necessary to clearly define what is meant by uncertainty, since in the literature confusion is often made on this issue. Therefore, the first objective of this thesis is to clarify this concept, starting with a key question: should be the choice of the intervention strategy to adopt based on the evaluation of the model prediction based on its ability to represent the reality or on the evaluation of what actually will happen on the basis of the information given by the model forecast? Once the previous idea is made unambiguous, the other main concern of this work is to develope a tool that can provide an effective decision support, making possible doing objective and realistic risk evaluations. In particular, such tool should be able to provide an uncertainty assessment as accurate as possible. This means primarily three things: it must be able to correctly combine all the available deterministic forecasts, it must assess the probability distribution of the predicted quantity and it must quantify the flooding probability. Furthermore, given that the time to implement prevention strategies is often limited, the flooding probability will have to be linked to the time of occurrence. For this reason, it is necessary to quantify the flooding probability within a horizon time related to that required to implement the intervention strategy and it is also necessary to assess the probability of the flooding time.

DIFFERENTIATING PATH AND SOURCE PROCESSES IN COMPLEX GEOLOGIC STRUCTURE THROUGH PASSIVE SEISMIC STUDIES: BERING GLACIER, ALASKA AND VILLARRICA VOLCANO, CHILE

Measuring shallow seismic sources provides a way to reveal processes that cannot be directly observed, but the correct interpretation and value of these signals depend on the ability to distinguish source from propagation effects. Furthermore, seismic signals produced by a resonating source can look almost identical to those produced by impulsive sources, but modified along the path. Distinguishing these two phenomena can be accomplished by examining the wavefield with small aperture arrays or by recording seismicity near to the source when possible. We examine source and path effects in two different environments: Bering Glacier, Alaska and Villarrica Volcano, Chile. Using three 3-element seismic arrays near the terminus of the Bering Glacier, we have identified and located both terminus calving and iceberg breakup events. We show that automated array analysis provided a robust way to locate icequake events using P waves. This analysis also showed that arrivals within the long-period codas were incoherent within the small aperture arrays, demonstrating that these codas previously attributed to crack resonance were in fact a result of a complicated path rather than a source effect. At Villarrica Volcano, seismometers deployed from near the vent to ~10 km revealed that a several cycle long-period source signal recorded at the vent appeared elongated in the far-field. We used data collected from the stations nearest to the vent to invert for the repetitive seismic source, and found it corresponded to a shallow force within the lava lake oriented N75°E and dipping 7° from horizontal. We also used this repetitive signal to search the data for additional seismic and infrasonic properties which included calculating seismic-acoustic delay times, volcano acoustic-seismic ratios and energies, event frequency, and real-time seismic amplitude measurements. These calculations revealed lava lake level and activity fluctuations consistent with lava lake level changes inferred from the persistent infrasonic tremor.

Archiving Montana’s Geologic Data

Since the 19th century, enormous amounts of time and money were invested in the exploration and development of Montana’s natural resources. These investments generated tremendous volumes of geologic, geophysical, and environmental data. Over the years, many of these data have been stored, forgotten, lost, or destroyed. In this lecture, Peggy discusses the steps the Bureau is taking to rescue and make the data available to the public.

Geologic-hydrologic setting and surface sedimentology in the Persian Gulf

1. Morphology and sedimentation The deepest parts of the Persian Gulf lie off the Iranian coast. Several swells separate the Persian Gulf into the Western Basin, the Central Basin and the Strait of Hormuz, which leads without noticeable morphological interruption onto the Biaban Shelf; the latter gradually drops off towards the continental slope, which itself has a strongly subdivided morphology. The sediment distribution in the Western Basin runs parallel to the basin's axis to a depth of 50 -60 m. This is caused by the shallow and uniform slope of the Iranian coast into the Western Basin, by clear exposure of the area to the Shamal-Winds and by tidal currents parallel to the basin's axis. Most other parameters also show isolines parallel to the coast line. Data from the sediment analyses show a net transport which extends out along the Central Swell: coarse fraction > 63 µ, total carbonate content, carbonate in fine fractions < 2 µ, 2-6 µ and 20-63 µ, calcite-aragonite ratios in the fine fractions 2-6 µ and 20-63 µ and quartz-dolomite ratios in fine fraction 2-6 µ. At least the uppermost 10-40 m of this sediment is late Holocene. This implies sedimentation rates of several meters per 1000 years. The slope from the Iranian coast into the Central Basin (max. depth 100 m) is generally steeper, with interspersed islands and flats. Both facts tend to disturb a sediment dustribition parallel to the basin's axis over extensive areas and may preclude any such trend from being detected by the methods and sample net used. The spatial distribution of the coarse fraction, however, seems to indicate sediment transport at greater water depths perpendicular to the basin's long axis and along the steepest gradients well into the Central Basin. The flats of the Central Basin have a sediment cover distinctly different from those of the deeper basin areas. Characteristic parameters are the extremely high percentages of coarse grained sediments, total content of carbonate CO2 over 40, low total organic carbon content, (however values are high if calculated on the basis of the < 63 µ fraction), low total N-content, and low C/N ratios. These characteristics probably result from the absence of any terrigenous material being brought in as well as from exposure to wave action. Finest terrigenous material is deposited in the innermost protected part of the Hormuz Bay. In the deep channel cut into the Biaban Shelf which carries the Persian Gulf out-flow water to the Indian Ocean, no fine grained sediment is deposited as shown by grain size data. 2. Geographic settings and sedimentation Flat lands border the Arabian coast of the Persian Gulf except for the Oman region. The high and steep Zagros Mountains form the Iranian coastline. Flat topography in combination with generally low precipitation precludes fluviatile sediment being added from the South. Inorganic and biogenic carbonates accumulating under low sedimentation rates are dominant on the shallow Arabic Shelf and the slopes into the Western and Central Basins. The fluviatile sediment brought in from the Iranian side, however decisively determine the composition of the Holocene sediment cover in the Persian Gulf and on the Biaban Shelf. Holocene sediments extend 20-30 km seaward into the Western Basin and about 25 km on to the Biaban Shelf. As mentioned before, sedimentation rates are of several meters/1000 years. The rocks exposed in the hinterland influence the sediments. According to our data the Redbeds of the Zagros Mountains determine the colour of the very fine grained sediments near the Iranian Coast of the Persian Gulf. To the West of Hormuz, addition of carbonate minerals is particularly high. Dolomite and protodolomite, deposited only in this area, as well as palygorskite, have proven to be excellent trace minerals. To the East of Hormuz, the supply of terrigenous carbonates is considerably lower. Clay minerals appear to bring in inorganically bound nitrogen thus lowering the C/N ratio in these sediments especially off river mouths. 3. Climate and sedimentation The Persian Gulf is located in a climatically arid region. This directly affects sedimentation through increased wind action and the infrequent but heavy rainfalls which cause flash floods. Such flash floods could be responsible for transporting sedheats into the Central Basin in a direction perpendicular to the Gulf's axis. Eolian influx is difficult to asses from our data; however, it probably is of minor importance from the Iranian side and may add, at the most, a few centimeters of fine sediment per 1000 years. 4. Hydrology and sedimentation High water temperatures favor inorganic carbonate precipitation in southern margin of the Gulf, and probably on the flats, as well as biogenic carbonate production in general. High evaporation plus low water inflow through rivers and precipitation cause a circulation pattern that is typical for epicontinental seas within the arid climate region. Surface water flows in from the adjoining ocean, in this case the Indian Ocean and sinks to the bottom of the Persian Gulf mainly in the northern part of the Western Basin, on the "Mesopotamischer Flachschelf" ard probably in the area of the "Arabischer Flachschelf". This sinking water continually rejuvenates the bottom out-flow water. The inflowing surface water from the Indian Ocean brings organic matter into the Persian Gulf, additional nutrients are added by the "fresh" upwelling waters of the Gulf of Oman. Both nutrients and organic matter diminish very rapidly as the water moves into the Persian Gulf. This depletion of nutrients and organic matter is the reasonfor generally low organic carbon contents of the Persian Gulf sediments. The Central Swell represents a distinct boundary, to the west of which the organic carbon content are lower than to the east when sediment samples of similar grain size distribution are compared. The outflow carries well oxygenated water over the bottom of the Persian Gulf and the resulting oxidation further decreases the content of organic matter. In the Masandam-Channel and in the Biaban-Shelf channel, the outflowing water prevents deposition of fine material and transports sediment particles well beyond the shelf margin. The outflowing water remains at a depth of 200-300 m depending on its density and releases ist suspending sediment load to the ocean floor, irrespectative of the bottom morphology. This is reflected in several parameters in which the sediments from beneath the outflow differ from nearby sediments not affected by the outflowing water. High carbonate content of total samples and of the individual size fraction as well as high aragonite and dolomite contents of individual size fractions characterize the sediment beneath the outflowing water. The tidal currents, which avt more or less parallel to the Gulf's axis, favor mixing of the water masses, they rework sediments at velocities reported here. This fact enlarges to a certain degree the extent of our interfaces which are based on only a few sample points (Persian Gulf and Biaban Shelf one sample per 620 km**2, continental slope one sample per 1000 km**2). The water on the continental slope shows and oxygen minimum at 200-1200 m which favors preservation of organically-bound carbon in the sediment. The low pH-values may even permit dissolution of carbonate minerals.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.