Using a shifted Grünwald-Letnikov scheme for the Caputo derivative to study anomalous solute transport in porous medium     

Pedro Victor Serra Mascarenhas  Ricardo Mendonça de Moraes  André Luís Brasil Cavalcante 《国际地质力学数值与分析法杂志》2019,43(11):1956-1977

We propose an extension of the shifted Grünwald-Letnikov method to solve fractional partial differential equations in the Caputo sense with arbitrary fractional order derivative α and with an advective term. The method uses the relation between Caputo and Riemann-Liouville definitions, the shifted Grünwald-Letnikov, and the traditional backward and forward finite difference method. The stability of the method is investigated for the implicit and explicit scheme with homogeneous boundary conditions, and a stability criterion is found for the advective-dispersive equation. An application of the method is used to solve contaminant diffusion and advective-dispersive problems. The numerical solution for the fractional diffusion and fractional advection-dispersion is compared with their respective analytical solutions for different time and space grid refinements. The diffusion simulation exhibited a good fit between the analytical and numerical solutions, with the explicit scheme going from stable to unstable as the time and space refinement changes. The fractional advection-dispersion application produced small deviations from the analytical solution. These deviations, however, are analogous to the numerical dispersions encountered in conventional finite difference solutions of the advection-dispersion equation. The new method is also compared with the traditional L2 method. Notably, an example that involves asymmetrical fractional conditions, a fractional diffusivity that depends on time, and a source term show how the methods compare. Overall, this study assesses the quality and easiness of use of the numerical method.  相似文献   

Slip rate of trench‐parallel normal faulting along the Mejillones Fault (Atacama Fault System): Relationships with the northern Chile subduction and implications for seismic hazards     

Jean‐Franois Ritz  Rgis Braucher  Didier Bourls  Bertrand Delouis  Carlos Marquardt  Alain Lavenu  Herv Philip  Luc Ortlieb 《地学学报》2019,31(4):390-404

The recent tectonics of the arid northern Chile Andean western forearc is characterized by trench‐parallel normal faults within the Atacama Fault System (AFS). Since the 1995‐Mw 8.1 Antofagasta earthquake, the mechanism driving this recent and localized extension is considered to be associated with the seismic cycle within the subduction zone. Analyzing morphotectonic patterns along these faults allows examining the seismic potential associated with the subduction zone. Using field Digital Elevation Models and in situ‐produced cosmogenic ¹⁰Be, we determined a 0.2 mm/a long‐term slip rate along the Mejillones Fault, one of the most prominent structures within the AFS. This result suggests that the AFS corresponds to slow slip rate faults despite the rapid subduction context. However, the size of coseismic slips observed along the AFS faults suggests that larger subduction earthquakes (Mw > 8.1) may occur episodically in the area.  相似文献   

Indirect dating of an olive tree planting event using luminescence of the sediments lying beneath the roots of the tree: a pilot study in the south-western part of Anatolia,Turkey     

Eren Şahiner  Georgios S. Polymeris  M. Altay Atlihan  Selcuk Aktürk  Niyazi Meriç 《第四纪科学杂志》2020,35(5):706-715

The aim of the present study is to attempt assessing the age of a monumental olive tree located between the Antique Cities of Militus (Didim-Aydın region) and Iasos (Milas-Muğla region). Wood from the trunk of an olive tree is not appropriate for conventional dating approaches such as dendrochronology or ¹⁴C. The sediments closely located surrounding and beneath the roots of the olive tree are considered indicative of the age of the planting event; therefore these sediments were dated using both quartz and feldspar luminescence signal protocols. Methodological aspects including the preheating plateaus, equivalent dose statistical approaches and dose rate using gamma spectrometry are also discussed, as dating of the associated palaeochannel sediments of the area are presented for the first time in the dating literature. The optically stimulated luminescence and/or infrared stimulated luminescence ages are extrapolated to date the event of the tree planting; it is the first time in the literature that an age is reported for an olive tree in the eastern Mediterranean region. The present study stands as the first experimental evidence that olive trees have been cultivated in the area since the Iron Age.  相似文献   

Transition from symmetry to asymmetry during continental rifting: an example from the Bight Basin–Terre Adélie (Australian and Antarctic conjugate margins)     

Nicolas Espurt  Jean‐Paul Callot  François Roure  Jennifer M. Totterdell  Heike I. M. Struckmeyer  Roland Vially 《地学学报》2012,24(3):167-180

  相似文献   

Quantifying peakflow attenuation/amplification in a karst river using the diffusive wave model with lateral flow     

Jean‐Baptiste Charlier  Roger Moussa  Pierre‐Yann David  Jean‐Franois Desprats 《水文研究》2019,33(17):2337-2354

The aim of this paper is to quantify peakflow attenuation and/or amplification in a river, investigating lateral flow from the intermediate catchment during floods. This is a challenge for the study of the hydrological response of permeable/intermittent streams, and our contribution refers to a modelling framework based on the inverse problem for the diffusive wave model applied in a karst catchment. Knowing the upstream and downstream hydrographs on a reach between two stations, we can model the lateral one, given information on the hydrological processes involved in the intermediate catchment. The model is applied to 33 flood events in the karst reach of the Iton River in French Normandy where peakflow attenuation is observed. The monitored zone consists of a succession of losing and gaining reaches controlled by strong surface‐water/groundwater (SW/GW) interactions. Our results show that despite a high baseflow increase in the reach, peakflow is attenuated. Model application shows that the intensity of lateral outflow for the flood component is linked to upstream discharge. A combination of river loss and overbank flow for highest floods is proposed for explaining the relationships. Our approach differentiates the role of outflow (river loss and overbank flow) and that of wave diffusion on peakflow attenuation. Based on several sets of model parameterization, diffusion is the main attenuation process for most cases, despite high river losses of up to several m³/s (half of peakflow for some parameterization strategies). Finally, this framework gives new insight into the SW/GW interactions during floods in karst basins, and more globally in basins characterized by disconnected river‐aquifer systems.  相似文献   

Elevated dust deposition in Tierra del Fuego (Chile) resulting from Neoglacial Darwin Cordillera glacier fluctuations          下载免费PDF全文

Heleen Vanneste  François De Vleeschouwer  Sébastien Bertrand  Antonio Martínez‐Cortizas  Aubry Vanderstraeten  Nadine Mattielli  Andrea Coronato  Natalia Piotrowska  Catherine Jeandel  Gaël Le Roux 《第四纪科学杂志》2016,31(7):713-722

  相似文献   

de Sitter expansion with anisotropic fluid in Bianchi type-I space-time     

Özgür Akarsu  Can Battal Kılınç 《Astrophysics and Space Science》2010,326(2):315-322

Some features of the Bianchi type-I universes in the presence of a fluid that wields an anisotropic equation of state (EoS) parameter are discussed in the context of general relativity. The models that exhibit de Sitter volumetric expansion due to the constant effective energy density (the sum of the energy density of the fluid and the anisotropy energy density) are of particular interest. We also introduce two locally rotationally symmetric models, which exhibit de Sitter volumetric expansion in the presence of a hypothetical fluid that has been obtained by minimally altering the conventional vacuum energy. In the first model, the directional EoS parameter on the x axis is assumed to be −1, while the ones on the other axes and the energy density of the fluid are allowed to be functions of time. In the second model, the energy density of the fluid is assumed to be constant, while the directional EoS parameters are allowed to be functions of time.  相似文献   

High-order <Emphasis Type="Italic">h</Emphasis>-adaptive discontinuous Galerkin methods for ocean modelling     

Paul-Emile Bernard  Nicolas Chevaugeon  Vincent Legat  Eric Deleersnijder  Jean-François Remacle 《Ocean Dynamics》2007,57(2):109-121

In this paper, we present an h-adaptive discontinuous Galerkin formulation of the shallow water equations. For a discontinuous Galerkin scheme using polynomials up to order , the spatial error of discretization of the method can be shown to be of the order of , where is the mesh spacing. It can be shown by rigorous error analysis that the discontinuous Galerkin method discretization error can be related to the amplitude of the inter-element jumps. Therefore, we use the information contained in jumps to build error metrics and size field. Results are presented for ocean modelling problems. A first experiment shows that the theoretical convergence rate is reached with the discontinuous Galerkin high-order h-adaptive method applied to the Stommel wind-driven gyre. A second experiment shows the propagation of an anticyclonic eddy in the Gulf of Mexico. An erratum to this article can be found at  相似文献   

Phytoplankton community structure defined by key environmental variables in Tagus estuary, Portugal   总被引：4，自引：0，他引：4  

Brogueira MJ  Oliveira Mdo R  Cabeçadas G 《Marine environmental research》2007,64(5):616-628

In this work, we analyze environmental (physical and chemical) and biological (phytoplankton) data obtained along Tagus estuary during three surveys, carried out in productive period (May/June/July) at ebb tide. The main objective of this study was to identify the key environmental factors affecting phytoplankton structure in the estuary. BIOENV analysis revealed that, in study period, temperature, salinity, silicate and total phosphorus were the variables that best explained the phytoplankton spatial pattern in the estuary (Spearman correlation, rho=0.803). A generalized linear model (GLM) also identified salinity, silicate and phosphate as having a high explanatory power (63%) of phytoplankton abundance. These selected nutrients appear to be consistent with the requirements of the dominant phytoplankton group, Baccilariophyceae. Apparently, phytoplankton community is adapted to fluctuations in light intensity, as suspended particulate matter did not come out as a key factor in shaping phytoplankton structure along Tagus estuary.  相似文献   

Analysis of phyllosilicate deposits in the Nili Fossae region of Mars: Comparison of TES and OMEGA data     

Joseph Michalski  François Poulet  Nicolas Mangold 《Icarus》2010,206(1):269-289

Diverse phyllosilicate deposits discovered previously in the Nili Fossae region with near infrared reflectance data are a window into the complex history of aqueous alteration on Mars. In this work, we used thermal infrared data from the Thermal Emission Spectrometer (TES) in combination with near infrared data from the Observatoire pour la Minéralogie, l’Eau, les Glaces, et l’Activité (OMEGA) to better constrain the mineralogy and geologic origin of these deposits. We developed a TES spectral index for identification of clay minerals, which correctly identifies the phyllosilicates in the Nili Fossae area and points to several other interesting deposits in the Syrtis Major region. However, detailed inspection of the TES spectral features of Nili Fossae phyllosilicates shows a feature at low wavenumbers (350-550 cm⁻¹) that is not an exact match to any specific Fe³⁺-, Al-, or Mg-rich phyllosilicate phase. Instead, the feature is more similar to basaltic glass and may indicate that the phyllosilicates in this region are: (1) rich in Fe²⁺ (based on similarity to trends seen in laboratory data of clay minerals), (2) poorly crystalline/extremely disordered, and/or (3) present within a matrix of actual basalt glass. This feature is similar to spectral features seen in altered rocks in the Columbia Hills region of Gusev Crater by previous authors. By calibrating measured spectral index values against mathematical spectral mixtures of typical martian dark surfaces and known abundances of alteration minerals, we are able to estimate an enrichment in abundance of alteration minerals in the altered surfaces. Many dark, Noachian deposits in the Nili Fossae area are enriched phyllosilicates by 20-30% (±10-15%) relative to dark, volcanic surfaces in the same region. The distribution and abundance of these phases indicates that alteration in the region was pervasive, but did not completely erase the original mineralogy of what was likely an Fe-rich basalt protolith. As a group, the Nili Fossae phyllosilicate deposits are fundamentally different from those found in the Mawrth Vallis region. Nili Fossae deposits have strong thermal infrared features related to admixed pyroxene, plagioclase, and occasionally olivine, whereas the Mawrth Vallis deposits contain no mafic minerals. Comparison of TES and OMEGA data also illustrates some more general differences between the datasets, including the impact of physical character of the martian surface on detectability of minerals in each spectral range.  相似文献