Search NASA⌕ Search

SEARCH · Search NASA

Results for “Fractional derivative”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

Fractional Modeling in Action: A Survey of Nonlocal Models for Subsurface Transport, Turbulent Flows, and Anomalous Materials

Modeling of phenomena such as anomalous transport via fractional-order differential equations has been established as an effective alternative to partial differential equations, due to the inherent ability to describe large-scale behavior with greater efficiency than fully-resolved classical models. In this review article, we first provide a broad overview of fractional-order derivatives with a clear emphasis on the stochastic processes that underlie their use. We then survey three exemplary application areas – subsurface transport, turbulence, and anomalous materials – in which fractional-order differential equations provide accurate and predictive models. For each area, we report on the evidence of anomalous behavior that justifies the use of fractional-order models, and survey both foundational models as well as more expressive state-of-the-art models. We also propose avenues for future research, including more advanced and physically sound models, as well as tools for calibration and discovery of fractional-order models.

97 MATHEMATICS AND COMPUTING↗

The bar rotation rate as a diagnostic of dark matter content in the centre of disc galaxies

We investigate the link between the bar rotation rate and dark matter content in barred galaxies by concentrating on the cases of the lenticular galaxies NGC 4264 and NGC 4277. These two gas-poor galaxies have similar morphologies, sizes, and luminosities. But, NGC 4264 hosts a fast bar, which extends to nearly the corotation, while the bar embedded in NGC 4277 is slow and falls short of corotation. We derive the fraction of dark matter ƒ DM, bar within the bar region from Jeans axisymmetric dynamical models by matching the stellar kinematics obtained with the MUSE integral-field spectrograph and using SDSS images to recover the stellar mass distribution. We build mass-follows-light models as well as mass models with a spherical halo of dark matter, which is not tied to the stars. We find that the inner regions of NGC 4277 host a larger fraction of dark matter (⁠ƒ DM,bar = 0.53±0.02⁠) with respect to NGC 4264 (⁠(⁠ƒ DM,bar = 0.33±0.04⁠) in agreement with the predictions of theoretical works and the findings of numerical simulations, which have found that fast bars live in baryon-dominated discs, whereas slow bars experienced a strong drag from the dynamical friction due to a dense DM halo. This is the first time that the bar rotation rate is coupled to fDM, bar derived from dynamical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗

The influence of cloud cover on the reliability of satellite-based solar resource data

Satellite-based solar resource data are often developed and validated by using binary cloudiness categories: clear sky or overcast cloudy sky. To investigate the reliability of solar resource data in partially cloudy conditions, we estimate cloud fraction using two distinct algorithms: a physical retrieval model using surface observed global horizontal irradiance (GHI) and direct normal irradiance (DNI) and a temporal average of cloud mask data estimated by the observed DNI. Our analysis reveals a significant presence of scattered clouds, broken clouds, and mismatches between satellite- and surface-based cloud data at 17 surface sites across the contiguous United States, though confidently clear and cloudy conditions collectively account for more than 70 % of the data. Solar radiation is computed using the National Solar Radiation Database (NSRDB) algorithm and validated using surface observations. Here, our findings suggest that, in the presence of scattered clouds, NSRDB data for clear-sky conditions can be subject to significant overestimation. In cloudy-sky conditions classified by satellite data, DNI computed by the Fast All-sky Radiation Model for Solar applications with DNI (FARMS-DNI) can be underestimated when limited clouds are detected by surface observations. The bias observed in several cloudiness categories indicates that the NSRDB is exceptionally accurate in confidently clear conditions. However, clear-sky conditions with scattered clouds and mismatched cloud data contribute significantly to the overall uncertainties in the NSRDB. Therefore, future improvements in solar resource data should involve development and implementation of satellite-derived cloud fraction and should consider a novel radiative transfer model accounting for amplified cloud reflection. The evaluation within cloudiness categories also provides a physical rationale for the superior performance of FARMS-DNI compared to the Direct Insolation Simulation Code (DISC) in both cloudy-sky and all-sky conditions.

14 SOLAR ENERGY↗

Search for heavy neutral Higgs bosons A and H in the $t\bar{t}$Z channel in proton-proton collisions at 13 TeV

A direct search for new heavy neutral Higgs bosons A and H in the $t\bar{t}$Z channel is presented, targeting the process pp → A → ZH with H → $t\bar{t}$. For the first time, the channel with decays of the Z boson to muons or electrons in association with all-hadronic decays of the $t\bar{t}$ system is targeted. The analysis uses proton-proton collision data collected at the CERN LHC with the CMS experiment at $\sqrt{𝑠}$ = 13 TeV, which correspond to an integrated luminosity of 138 fb −1 . No signal is observed. Upper limits on the product of the cross section and branching fractions are derived for narrow resonances A and H with masses up to 2100 and 2000 GeV, respectively, assuming A boson production through gluon fusion. The results are also interpreted within two-Higgs-doublet models, where A and H are CP-odd and CP-even states, respectively, complementing and substantially extending the reach of previous searches.

2HDM↗

Revealing the Influence of Binding Motifs on Electron Transfer and Recombination Kinetics for CdSe Quantum Dots Functionalized with a Modified Viologen

Anchoring of molecules to the surfaces of semiconductor nanocrystals (NCs) presents an opportunity to leverage the precise synthetic tunability of molecular function and the remarkable light harvesting properties of NCs to drive photochemical reactions. However, charge transfer between the two species depends not only on the energy level alignments but also on the details of their binding interactions, which are difficult to probe. Here, in this work, we characterize the binding between CdSe quantum dots (QDs) and a new phosphonated derivative of the electron acceptor methyl viologen, designed to attach to the QD surface via the phosphonate group. We use isothermal titration calorimetry to probe the thermodynamics of the QD–molecule interaction and use the parameters determined therein to analyze transient absorption spectroscopy measurements of forward and back electron transfer from QDs to the viologen. We find that the ligand-like phosphonate binding leads to an electron-transfer rate constant that is 3 orders of magnitude smaller than that for the face-on binding of the bipyridine ring of methyl viologen. Back electron transfer is also significantly slower in the derivative. Interestingly, a minor fraction of the phosphonated derivative also binds in the face-on configuration, with similar forward and back electron transfer kinetics as methyl viologen. Numerical simulations show that the ligand-like binding will lead to significantly improved quantum yields of photocatalysis over a wide range of reaction rates. By independently characterizing binding thermodynamics and charge transfer kinetics, this work reveals how the complexities underlying electron transfer at the NC–molecule interface determine photocatalytic outcomes. This work also represents a step toward controlling forward and back electron transfer kinetics via rational molecular design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Nature of the nonequilibrium phase transition in the non-Markovian driven Dicke model

The Dicke model famously exhibits a phase transition to a superradiant phase with a macroscopic population of photons and is realized in multiple settings in open quantum systems. In this paper, we study a variant of the Dicke model where the cavity mode is lossy due to the coupling to a Markovian environment while the atomic mode is coupled to a colored bath. We analytically investigate this model by inspecting its low-frequency behavior via the Schwinger-Keldysh field theory and carefully examine the nature of the corresponding superradiant phase transition. Integrating out the fast modes, we can identify a simple effective theory allowing us to derive analytical expressions for various critical exponents including the dynamical exponent. We find excellent agreement with previous numerical results when the non-Markovian bath is at zero temperature; however, contrary to these studies, our low-frequency approach reveals that the same exponents govern the critical behavior when the colored bath is at finite temperature unless the chemical potential is zero. Furthermore, we show that the superradiant phase transition is classical in nature, while it is genuinely nonequilibrium. In this work, we derive a fractional Langevin equation and conjecture the associated fractional Fokker-Planck equation that captures the system's long-time memory as well as its nonequilibrium behavior. Finally, we consider finite-size effects at the phase transition and identify the finite-size scaling exponents, unlocking a rich behavior in both statics and dynamics of the photonic and atomic observables.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Velocity Profiling of a Gas–Solid Fluidized Bed Using Electrical Capacitance Volume Tomography

In this study, a method of producing velocity profile maps from electrical capacitance volume tomography (ECVT) measurements by reconstructing displacement from measured changes in capacitance is developed and applied to fluidized bed systems. The mapping of the reconstruction leverages the gradient of the sensitivity distribution of the ECVT sensor to circumvent the need for image cross correlation techniques. Experimental data of both bubbling and slugging fluidized beds are collected in a cold flow model. Adaptation of the technique is discussed in detail, and velocity profiles are obtained for a range of gas flow rates. The produced velocity maps are compared against the established methods of cross correlation and against empirical correlations from the literature and are found to agree well in tracking slug and bubble velocity. The exception is when the tracked object is large relative to the ECVT sensor dimensions, a scenario that can be avoided through proper sensor design. The quantities of average velocity, momentum, and solid and gas volume fraction are derived from the image and velocity profiles. The results demonstrate and extend the power of ECVT as a measurement tool for the study and monitoring of gas–solid fluidized beds by providing a computationally cheaper alternative to 3-D cross correlation for deriving velocity profiles.

47 OTHER INSTRUMENTATION↗

Observations Pertaining to Precipitation within the Northeast Pacific Stratocumulus-to-Cumulus Transition

Three genuine stratocumulus-to-cumulus transitions sampled during the Cloud System Evolution over the Trades (CSET) campaign are documented. The focus is on Lagrangian evolution of in situ precipitation, thought to exceed radar/lidar retrieved values because of Mie scattering. Two of the three initial stratocumulus cases are pristine [cloud droplet number concentrations ($N_d$) of ~22 cm –3 ] but occupied boundary layers of different depths, while the third is polluted ($N_d$ ~ 225 cm –3 ). Hourly satellite-derived cloud fraction along Lagrangian trajectories indicate that more quickly deepening boundary layers tend to transition faster, into more intense but more occasional precipitation. These transitions begin either in the morning or late afternoon, suggesting that preceding night processes can precondition or delay the inevitable transition. The precipitation shifts toward larger drop sizes throughout the transition as the boundary layers deepen, with aerosol concentrations only diminishing in two of the three cases. Ultraclean ($N_d$ < 1 cm –3 ) cumulus clouds evolved from pristine stratocumulus cloud with unusually high precipitation rates occupying a shallow, well-mixed boundary layer. Results from a simple one-dimensional evaporation model and from radar/lidar retrievals suggest subcloud evaporation likely increases throughout the transition. This, coupled with larger drop sizes capable of lowering the latent cooling profile, facilitates the transition to more surface-driven convection. The coassociation between boundary layer depth and precipitation does not provide definitive conclusions on the isolated effect of precipitation on the pace of the transition. Differences between the initial conditions of the three examples provide opportunities for further modeling studies.

54 ENVIRONMENTAL SCIENCES↗

A two-equation soot-in-flamelet modeling approach applied under Spray A conditions

Soot production (including formation and oxidation) is studied in the transient, high-pressure and turbulent n-dodecane Spray A flames from the Engine Combustion Network (ECN) using computational fluid dynamics (CFD) simulations. A two-equation soot-in-flamelet modeling approach is applied within the framework of the Unsteady Flamelet Progress Variable (UFPV) model and results are validated against experimental data. Equations for soot mass fraction and soot number density derived in the mixture fraction space are solved in the context of detailed flamelet calculations. Source terms for the different steps in the soot chemistry are tabulated and incorporated in the flamelet manifold. For the reference condition, the modeling approach based on the tabulated flamelet manifold reduces the computational cost of a CFD calculation by approximately 40 times compared to a non-tabulated well-mixed (WM) modeling approach. The soot-in-flamelet approach is then extended to study the effect of ambient oxygen concentration, ambient mixture composition and ambient temperature on soot production. Results show that the modeling approach is able to capture the experimental trends for the soot volume fraction (SVF) with good quantitative agreement, especially in the soot ramp-up region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Transfer hydrogenation of heavy hydrocarbons with hydrogen-donor solvents

Processes for converting the coal-derived heavy-oil fraction of syncrude to polyols are described. The processes involve mixing a feed stream comprising the coal-derived heavy-oil fraction with an alcohol and aqueous sulfuric acid, heating the mixture, reacting the coal-derived heavy-oil fraction with ozone, and reacting the ozonated heavy-oil fraction with glycerin to form the polyol. In some cases, the ozonated heavy-oil fraction can be neutralized before reacting the ozonated heavy-oil fraction with the glycerin.

Chauhan, Satya P.↗

Generalized fractional approach to solving partial differential equations with arbitrary dispersion relations

The dynamics of physical systems are typically modeled by partial differential equations (PDEs). Because of the complexity in physical systems, PDE models rely on approximations or limiting cases. To model the full complexity of physical systems, it is necessary to use multiscale approaches where appropriate models are used at each level. Alternatively, complexity can be directly tackled by fractional integrodifferential equations that must be derived for each system. Here, we argue that complexity can be restored in PDEs by describing them from a fractional calculus perspective. Rather than deriving a fractional integrodifferential equation, we reinterpret the dispersion relation of the system by use of the Riesz definition, which contains the required information relating the energy and momentum space of the system and thus fully describes their dynamics. The approach is demonstrated by two examples: the Landau–Lifshitz equation in a 1D ferromagnetic chain and a modified KdV equation supporting surface gravity waves or Euler dispersion. The presented approach is applicable to fluids, soft matter, and solid-state matter and can be readily generalized to higher dimensions and more complex systems. While numerical calculations are needed to determine the fractional operator, the approach is analytical and can be utilized to determine analytical solutions and investigate nonlinear problems.

97 MATHEMATICS AND COMPUTING↗

Sequential Precipitation Input Tagging (SPIT) to Estimate Water Transit Times and Hydrologic Tracer Dynamics Within Water‐Tagging Enabled Hydrologic Models

Determining the age distribution of water exiting a catchment is important for understanding groundwater storage and mixing. New water-tagging capabilities within models track precipitation events as they move through simulated storages, yet forward modeling of individual events may not systematically capture the full transit time distribution (TTD). Here, we present a “sequential precipitation input tagging” (SPIT) framework to tag all input precipitation at regular intervals during extended model simulations. Monthly tags over 7 years were applied at six National Ecological Observatory Network sites to calculate TTDs and derive mean virtual tracer age, $\overline{T_{V}}$, fractions of young water, F yw , and hydrologic tracer concentrations (water isotopes δ 18 O and δ 2 H) within a tagging enabled version of the Weather Research and Forecast hydrologic model (WRF-Hydro). Throughout seven simulation years, the fraction of simulated discharge derived from tagged events, F tag , increased each year, with the final year's F tag ranging from 66% to 100% and highlights the need to apply SPIT over many years to understand TTDs. When the F tag was >75%, simulated $\overline{T_{V}}$ ranged 179–923 days and F yw 0.6%–23.9%, with daily values exhibiting a power-law relationship with precipitation, discharge, and groundwater. Through implementation of SPIT, we find this hydrologic model configuration performs poorly in estimation of $\overline{T_{V}}$ and F yw (root mean squared error of 469 days and 14.4% respectively), suggesting it misrepresents subsurface mixing. Thus, the SPIT framework provides a reproducible approach to calculate watershed transit times within tagging enabled models and thereby assess and improve representation of hydrologic processes.

fraction of young water↗

Assessing Arctic low-level clouds and precipitation from above – a radar perspective

Most Arctic clouds occur below 2 km altitude as revealed by CloudSat satellite observations. However, recent studies suggest that the relatively coarse spatial resolution, low sensitivity, and blind zone of the radar installed on CloudSat may not enable it to comprehensively document low-level clouds. We investigate the impact of these limitations on the Arctic low- level cloud fraction, which is the amount of cloudy points with respect to all points as a function of height, derived from CloudSat radar observations. For this purpose, we leverage highly resolved vertical profiles of low-level cloud fraction derived from downlooking Microwave Radar/radiometer for Arctic Clouds (MiRAC) radar reflectivity measurements. MiRAC has been operated during four aircraft campaigns taking place in the vicinity of Svalbard during different times of the year and covering more than 25,000 km. This allows us to study the dependence of CloudSat limitations on different synoptic and surface conditions.

54 ENVIRONMENTAL SCIENCES↗

Composition and Plume Gas Interaction Control Iron Fractional Solubility More Than Particle Size in Volcanic Ash: Implications for Fertilization of the North Atlantic

Deposition of volcanic ash is thought to impact marine biogeochemical cycling by adding soluble iron (Fe) to the surface ocean. The magnitude of this input is a function of the amount of ash deposited, the total Fe content in the ash, and ash‐derived Fe's fractional solubility. However, the relative importance of chemical composition, acidic processing by the volcanic plume, and ash particle size in determining solubility is unclear. We paired an aerosol leach meant to provide an upper limit for fractional Fe solubility with chemical analyses of ash from the Cumbre Vieja (CV) and La Soufrière eruptions, which both impacted the North Atlantic in 2021. Fe in the ash samples is <6% soluble, but Fe fractional solubility in CV ash is approximately triple that of La Soufrière ash. Compared to La Soufrière, a larger proportion of the Fe in CV ash is in silicate rather than oxide minerals, which release more soluble Fe. Elevated levels of surficial fluorine (F) also suggest that CV ash was subjected to a more fluorine‐rich eruption plume and underwent more acidic processing. Particle size does not appear to be a primary control on Fe release. We estimate that the CV eruption had a much larger impact on dissolved Fe (DFe) concentration in the surface ocean than the La Soufrière eruption because of differences in soluble Fe content and particle deposition velocity. These differences may help explain why some eruptions elicit a biological response in the ocean while others do not.

Elliott, Hope E. [Univ. of Miami, FL (United State↗

Sum-of-Fractions Methodology for Actinides in Water- and Polyethylene-Moderated and -Reflected Systems

Sum-of-fractions is a method intended to make sure a subcritical margin for aqueous solutions and slurries of fissionable isotopes exists. The method indicates that a system is subcritical if the sum of the ratios of the mass of each isotope (in a mixture) to its individual minimum subcritical mass limit is less than or equal to one. Historically, the basis of the sum-of-fractions has been derived from allowances given in the American National Standards Institute (ANSI)/ American Nuclear Society (ANS)-8.15-1981. However, the allowance was removed in ANSI/ANS-8.15-2014 due to a lack of technical basis. A methodology was developed to assess the validity of using the sum-of-fractions for water- or polyethylene-moderated systems for the following nuclides: 232 U, 233 U, 234 U, 235 U, 237 Np, 236 Pu, 238 Pu, 239 Pu, 240 Pu, 241 Pu, 242 Pu, 241 Am, 242 m Am, 243 Am, 242 Cm, 243 Cm, 244 Cm, 245 Cm, 246 Cm, 247 Cm, 249 Cf, and 251 Cf. The methodology uses available benchmark data for mixtures of 233 U, 235 U, and 239 Pu to establish the calculational margin, and a mass limit reduction to establish the margin of subcriticality. Water- or polyethylene-moderated and -reflected mixtures containing the nuclides are evaluated with the code system, SCALE 6.2.4. Including the calculational margin, subcritical mass limits for each nuclide were computed for optimally water- or polyethylene-moderated and fully reflected systems. These masses were used to create nuclide mixtures in which the sum of the mass to subcritical mass limit ratios is one. The various nuclide mixtures were modeled over a range of moderation and demonstrate the keff does not exceed the calculational margin. For additional assurance of subcriticality, a significant mass reduction is applied to each computed minimum critical mass of the nuclides without adequate benchmark data consistent with the method in ANSI/ANS-8.15-2014.

07 ISOTOPE AND RADIATION SOURCES↗

Sum-of-Fractions Method

Sum-of-fractions is a method intended to make sure a subcritical margin for aqueous solutions and slurries of fissionable isotopes exists. The method indicates that a system is subcritical if the sum of the ratios of the mass of each isotope (in a mixture) to its individual minimum subcritical mass limit is less than or equal to one. Historically, the basis of the sum-of-fractions has been derived from allowances given in the American National Standards Institute (ANSI)/ American Nuclear Society (ANS)-8.15-1981. However, the allowance was removed in ANSI/ANS-8.15-2014 due to a lack of technical basis. A methodology was developed to assess the validity of using the sum-of-fractions for water- or polyethylene-moderated systems for the following nuclides: 232U, 233U, 234U, 235U, 237Np, 236Pu, 238Pu, 239Pu, 240Pu, 241Pu, 242Pu, 241Am, 242mAm, 243Am, 242Cm, 243Cm, 244Cm, 245Cm, 246Cm, 247Cm, 249Cf, and 251Cf. The methodology uses available benchmark data for mixtures of 233U, 235U, and 239Pu to establish the calculational margin, and a mass limit reduction to establish the margin of subcriticality. Water- or polyethylene-moderated and -reflected mixtures containing the nuclides are evaluated with the code system, SCALE 6.2.4. Including the calculational margin, subcritical mass limits for each nuclide were computed for optimally water- or polyethylene-moderated and fully reflected systems. These masses were used to create nuclide mixtures in which the sum of the mass to subcritical mass limit ratios is one. The various nuclide mixtures were modeled over a range of moderation and demonstrate the keff does not exceed the calculational margin. For additional assurance of subcriticality, a significant mass reduction is applied to each computed minimum critical mass of the nuclides without adequate benchmark data consistent with the method in ANSI/ANS-8.15-2014.

criticality safety, Actinide↗

Development of genetic tools for heterologous protein expression in a pentose‐utilizing environmental isolate of Pseudomonas putida

Abstract Pseudomonas putida has emerged as a promising host for the conversion of biomass‐derived sugars and aromatic intermediates into commercially relevant biofuels and bioproducts. Most of the strain development studies previously published have focused on P. putida KT2440, which has been engineered to produce a variety of non‐native bioproducts. However, P. putida is not capable of metabolizing pentose sugars, which can constitute up to 25% of biomass hydrolysates. Related P. putida isolates that metabolize a larger fraction of biomass‐derived carbon may be attractive as complementary hosts to P. putida KT2440. Here we describe genetic tool development for P. putida M2, a soil isolate that can metabolize pentose sugars. The functionality of five inducible promoter systems and 12 ribosome binding sites was assessed to regulate gene expression. The utility of these expression systems was confirmed by the production of indigoidine from C6 and C5 sugars. Chromosomal integration and expression of non‐native genes was achieved by using chassis‐independent recombinase‐assisted genome engineering (CRAGE) for single‐step gene integration of biosynthetic pathways directly into the genome of P. putida M2. These genetic tools provide a foundation to develop hosts complementary to P. putida KT2440 and expand the ability of this versatile microbial group to convert biomass to bioproducts.

59 BASIC BIOLOGICAL SCIENCES↗

Formation and fate of freshwater on an ice floe in the Central Arctic

The melt of snow and sea ice during the Arctic summer is a significant source of relatively fresh meltwater. The fate of this freshwater, whether in surface melt ponds or thin layers underneath the ice and in leads, impacts atmosphere–ice–ocean interactions and their subsequent coupled evolution. Here, we combine analyses of datasets from the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition (June–July 2020) for a process study on the formation and fate of sea ice freshwater on ice floes in the Central Arctic. Our freshwater budget analyses suggest that a relatively high fraction (58 %) is derived from surface melt. Additionally, the contribution from stored precipitation (snowmelt) outweighs by 5 times the input from in situ summer precipitation (rain). The magnitude and rate of local meltwater production are remarkably similar to those observed on the prior Surface Heat Budget of the Arctic Ocean (SHEBA) campaign, where the cumulative summer freshwater production totaled around 1 m during both. A relatively small fraction (10 %) of freshwater from melt remains in ponds, which is higher on more deformed second-year ice (SYI) compared to first-year ice (FYI) later in the summer. Most meltwater drains laterally and vertically, with vertical drainage enabling storage of freshwater internally in the ice by freshening brine channels. In the upper ocean, freshwater can accumulate in transient meltwater layers on the order of 0.1 to 1 m thick in leads and under the ice. The presence of such layers substantially impacts the coupled system by reducing bottom melt and allowing false bottom growth; reducing heat, nutrient, and gas exchange; and influencing ecosystem productivity. Regardless, the majority fraction of freshwater from melt is inferred to be ultimately incorporated into the upper ocean (75 %) or stored internally in the ice (14 %). Terms such as the annual sea ice freshwater production and meltwater storage in ponds could be used in future work as diagnostics for global climate and process models. For example, the range of values from the CESM2 climate model roughly encapsulate the observed total freshwater production, while storage in melt ponds is underestimated by about 50 %, suggesting pond drainage terms as a key process for investigation.

54 ENVIRONMENTAL SCIENCES↗