Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “complexes”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

First-Row Transition Metals Complexes with Fused Oxazolidine (FOX) Ligands

In this comprehensive study, we report the syntheses and crystal structures of 6 fused oxazolidine (FOX) bicycles and 22 complexes with FOX ligands. The structures are compared to the few FOX complexes found in literature. For our investigations, we developed synthetic strategies to selectively synthesize chiral or achiral FOX bicycles with pyridyl substituents. The coordination chemistry of the chiral and achiral FOX compounds with first-row transition metals is discussed in an indepth analysis of the crystal structures. We found that the binding modes of the ligands are dependent on the configuration of the FOX backbone and on the steric demand of the substituents which were introduced to the pyridyl groups. Here, the configuration of the FOX ligands is retained during the coordination of metal(II) ions such as Mn 2+ , Fe 2+ , Co 2+ , Ni 2+ , Cu 2+ . Reactions of these metal ions with chiral FOX ligands lead to the formation of chiral complexes. In contrast, stronger Lewis acidic metal(III) ions, such as Fe 3+ or Al 3+ , induce an isomerization of the chiral FOX ligands to their achiral diastereomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

General bounds on holographic complexity

We prove a positive volume theorem for asymptotically AdS spacetimes: the maximal volume slice has nonnegative vacuum-subtracted volume, and the vacuum-subtracted volume vanishes if and only if the spacetime is identically pure AdS. Under the Complexity=Volume proposal, this constitutes a positive holographic complexity theorem. The result features a number of parallels with the positive energy theorem, including the assumption of an energy condition that excludes false vacuum decay (the AdS weak energy condition). Our proof is rigorously established in broad generality in four bulk dimensions, and we provide strong evidence in favor of a generalization to arbitrary dimensions. Our techniques also yield a holographic proof of Lloyd’s bound for a class of bulk spacetimes. We further establish a partial rigidity result for wormholes: wormholes with a given throat size are more complex than AdS-Schwarzschild with the same throat size.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Linear growth of circuit complexity from Brownian dynamics

How rapidly can a many-body quantum system generate randomness? Using path integral methods, we demonstrate that Brownian quantum systems have circuit complexity that grows linearly with time. In particular, we study Brownian clusters of N spins or fermions with time-dependent all-to-all interactions, and calculate the Frame Potential to characterize complexity growth in these models. In both cases the problem can be mapped to an effective statistical mechanics problem which we study using path integral methods. Within this framework it is straightforward to show that the kth Frame Potential comes within ϵ of the Haar value after a time of order t ~ kN + k log k + log ϵ –1 . Using a bound on the diamond norm, this implies that such circuits are capable of coming very close to a unitary k-design after a time of order t ~ kN. We also consider the same question for systems with a time-independent Hamiltonian and argue that a small amount of time-dependent randomness is sufficient to generate a k-design in linear time provided the underlying Hamiltonian is quantum chaotic. These models provide explicit examples of linear complexity growth that are analytically tractable and are directly applicable to practical applications calling for unitary k-designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Chapter One - Selectivity in the activation of C–H bonds by rhodium and iridium complexes

Trispyrazolylborate complexes of rhodium and iridium have been extensively investigated over the past 3 decades with special attention to their ability to activate C–H bonds. The rhodium complexes of tris-(3,5-dimethylpyrazolyl)borate have been the subject of numerous thermodynamic investigations that provide information about rhodium-metal carbon bond strengths. This insight arises as a result of the reversibility of C–H activation with rhodium. In the case of iridium, C–H bond activation reactions are also widespread, and some of these also show reversibility. Access to some key trispyrazolylborate iridium(I) and iridium(III) starting materials has given way to a multitude of studies of reactions with small molecules in which C–H bonds are made and broken reversibly. The stability of Fischer carbenes plays a role in the observed products. Here, the reactivity of trispyrazolylborate complexes of rhodium and iridium compounds over the past decade (since 2010) are summarized here. Related reactions of X–H bonds with these trispyrazolylborate compounds are also included for completeness (X = O, N, S, B, Si).

02 PETROLEUM↗
DOE OSTI2024

Ligand design in lanthanide complexes for luminescence, therapy, and sensing

In this chapter, we explore the research done by the de Bettencourt-Dias group in luminescent lanthanide complexes. The group has isolated ligands that sensitize the luminescence of lanthanide ions, and that bear additional functional groups that impart unique properties to the resulting complexes, such as the ability to generate cytotoxic singlet oxygen, or report on, through changes in emission intensity, viscosity or temperature at the nanoscale. The energy transfer processes are accessed through one-photon excitation or two-photon excitation through the ligand, as well as upconversion. Moreover, the group explored unique coordination chemistry of lanthanides, taking advantage of secondary interactions to isolate outer-sphere complexes of these ions.

Barber, Patrick S.↗
DOE OSTI2026

Emergence of complex-regular eutectic patterns in Al–Ge: Observations from correlative nano-imaging

Eutectic solidification exemplifies nonequilibrium pattern formation, making it a well-studied moving boundary problem. Yet the mechanisms behind the formation of complex-regular microstructures – particularly in highly anisotropic systems with a significant volume fraction of a faceted phase – remain poorly understood. Our understanding of such systems is made complicated by the nonlinear interface kinetics and unique growth dynamics characteristic of faceted phases. To address these challenges, we investigate a model Al–Ge eutectic system, where the faceted Ge phase constitutes a substantial volume fraction (~0.35) and where the two solid phases arrange into so-called “fishbone” or “feather” complex-regular patterns. Using synchrotron-based x-ray nano-imaging and nanotomography with high spatial resolution (22 nm per pixel), we capture in real-time the evolution of the solid–liquid interfaces and the resulting three-dimensional microstructures in this faceted/non-faceted eutectic system. By integrating these observations with electron backscattered diffraction, we elucidate the crystallographic biases on the solidification process and the mechanisms driving the formation of such complex-regular microstructures. These findings inform a new growth model for irregular eutectics in (near-)symmetrical phase diagrams, offering insight on advanced microstructural design and processing strategies. More broadly, we demonstrate how interfacial curvature is generated in irregular eutectic alloys and how it depends on the volume fraction of the faceted phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Process optimization of complex geometries using feed forward control for laser powder bed fusion additive manufacturing

Additive manufacturing (AM) enables the fabrication of complex designs that are difficult to create by other means. Metal parts manufactured by laser powder bed fusion (LPBF) can incorporate intricate design features and demonstrate desirable mechanical properties. However, printing a part that is qualified for its intended application often involves reprinting and discarding many parts to eliminate defects, improve dimensional accuracy, and increase repeatability. The process of iteratively converging on the appropriate build parameters increases the time and cost of creating functional LPBF manufactured parts. This work describes a fast, scalable method for part-scale process optimization of arbitrary geometries. Additionally, the computational approach uses feature extraction to identify scan vectors in need of parameter adaptation and applies results from simulation-based feed forward control models. This method provides a framework to quickly optimize complex parts through the targeted application of models with a range of fidelity and by automating the transfer of optimization strategies to new part designs. The computational approach and algorithmic framework are described, a software package is implemented, the method is applied to parts with complex features, and parts are printed on a customized open architecture LPBF machine.

3D printing↗
DOE OSTI2025

Data-driven equation-free dynamics applied to many-protein complexes: The microtubule tip relaxation

Microtubules (MTs) constitute the largest components of the eukaryotic cytoskeleton and play crucial roles in various cellular processes, including mitosis and intracellular transport. The property allowing MTs to cater to such diverse roles is attributed to dynamic instability, which is coupled to the hydrolysis of GTP (guanosine-5'-triphosphate) to GDP (guanosine-5'-diphosphate) within the β-tubulin monomers. Understanding the equilibrium dynamics and the structural features of both GDP- and GTP-complexed MT tips, especially at an all-atom level, remains challenging for both experimental and computational methods because of their dynamic nature and the prohibitive computational demands of simulating large, many-protein systems. This study employs the “equation-free” multiscale computational method to accelerate the relaxation of all-atom simulations of MT tips toward their putative equilibrium conformation. Using large MT lattice systems (14 protofilaments × 8 heterodimers) comprising ~21-38 million atoms, we applied this multiscale approach to leapfrog through time and nearly double the computational efficiency in realizing relaxed all-atom conformations of GDP- and GTP-complexed MT tips. Commencing from an initial 4 μs unbiased all-atom simulation, we interleave coarse projective “equation-free” jumps with short bursts of all-atom molecular dynamics simulation to realize an additional effective simulation time of 1.875 μs. Our 5.875 μs of effective simulation trajectories for each system expose the subtle yet essential differences in the structures of MT tips as a function of whether β-tubulin monomer is complexed with GDP or GTP, as well as the lateral interactions within the MT tip, offering a refined understanding of features underlying MT dynamic instability. Furthermore, the approach presents a robust and generalizable framework for future explorations of large biomolecular systems at atomic resolution.

Wu, Jiangbo [University of Chicago, IL (United Sta↗
DOE OSTI2024

Familial Alzheimer mutations stabilize synaptotoxic γ-secretase-substrate complexes

Mutations that cause familial Alzheimer’s disease (FAD) are found in amyloid precursor protein (APP) and presenilin, the catalytic component of γ-secretase, that together produce amyloid β-peptide (Aβ). Nevertheless, whether Aβ is the primary disease driver remains controversial. We report here that FAD mutations disrupt initial proteolytic events in the multistep processing of APP substrate C99 by γ-secretase. Cryoelectron microscopy reveals that a substrate mimetic traps γ-secretase during the transition state, and this structure aligns with activated enzyme-substrate complex captured by molecular dynamics simulations. In silico simulations and in cellulo fluorescence microscopy support stabilization of enzyme-substrate complexes by FAD mutations. Neuronal expression of C99 and/or presenilin-1 in Caenorhabditis elegans leads to synaptic loss only with FAD-mutant transgenes. Designed mutations that stabilize the enzyme-substrate complex and block Aβ production likewise led to synaptic loss. Collectively, these findings implicate the stalled process—not the products—of γ-secretase cleavage of substrates in FAD pathogenesis.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Hydrothermal solubility of Dy hydroxide as a function of pH and stability of Dy hydroxyl aqueous complexes from 25 to 250 °C

The rare earth elements (REE) have important applications in green energy technologies. The formation of mineral deposits in geologic systems commonly involves hydrothermal fluids which can mobilize the REE. However, the REE speciation is not well known as a function of pH. The thermodynamic properties of REE hydroxyl complexes used in geochemical models are based on the Helgeson-Kirkham-Flowers (HKF) equation of state parameters which were derived by extrapolation of low temperature experimental and estimated data. In this study, Dy hydroxide solubility experiments are combined with available literature data to improve these models from 25 to 250 °C and optimize the thermodynamic properties of Dy 3+ and Dy hydroxyl complexes using GEMSFITS. Batch-type solubility experiments were conducted from 150 to 250 °C and at saturated water vapor pressure in perchloric acid solutions with initial pH values of 2 to 5 in 0.5 pH unit increments. The measured solubility of Dy hydroxide is retrograde with temperature and decreases with pH. The logarithm of total dissolved Dy molality ranges from –2.3 to –5.3 at 150 °C (pH 4.7–5.5), from –2.4 to –5.6 at 200 °C (pH 3.9–5.1), and from –3.7 to –6.9 at 250 °C (pH of 3.4 and 5.0). The optimized standard partial molal Gibbs energies of formation (Δ f G° T ) derived for Dy 3+ and DyOH 2+ display a close to linear relationship with temperature, fitting with previous optimizations based on DyPO 4 solubility data in the literature. A comparison of the optimized ΔfG°T values for aqueous Dy species with predictions from available HKF parameters indicates significant differences ranging from +11 to –26 kJ/mol between 25 and 250 °C. The experimental fits are used to derive the Dy hydroxide solubility products (K s0 ) and formation constants for the hydrolysis of Dy (β n with n = 1 to 3; Dy 3+ + nOH – = DyOH n 3-n ) as a function of temperature. The optimization method presented yields accurate thermodynamic properties for the Dy 3+ aqua ions and the DyOH 2+ species at the acidic to mildly acidic pH studied whereas more experimental work is needed at near-neutral and alkaline conditions to better constrain the other hydroxyl complexes. Furthermore, the optimized thermodynamic data have a significant impact on geochemical modeling of the mobility and solubility of REE minerals in acidic hydrothermal fluids.

58 GEOSCIENCES↗
DOE OSTI2025

Spectrophotometric determination of the stability of La hydroxyl complexes at near neutral to alkaline pH from 25 to 75 °C

The hydrolysis of rare earth elements (REE) potentially controls their mobility during fluid-rock interaction in a broad range of pH and temperature conditions. However, there is still a lack of thermodynamic data for modeling accurately the stability of REE hydroxyl complexes in hydrothermal aqueous fluids. Here, in this study, UV–Vis spectrophotometric experiments were conducted from 25 to 75 °C in near-neutral to alkaline NaOH-bearing aqueous solutions with varying lanthanum (La) concentrations (0 to ∼0.23 mmol/kg). The color indicator m-cresol purple was used to determine in situ pH and derive the average OH− ligand number ($\overrightarrow{n}$) and formation constants for the La hydroxyl complexes (LaOH 2+ , La(OH) 2 + , and La(OH) 3 0 ). From 25 to 50 °C, $\overrightarrow{n}$ ranges between ∼1 and 2 at pH from 7.0 to 9.3. At 75 °C, $\overrightarrow{n}$ ranges between ∼1.5 and 3 at pH from 6.3 to 8.8. These results suggest the predominance of LaOH 2+ and La(OH) 2 + complexes from 25 to 50 °C, and an increased predominance of La(OH) 3 0 at 75 °C. The cumulative formation constants (β n °, n = 1 to 3) are derived for the reaction La 3+ + nOH − = La(OH) n 3-n , and fitted between 25 and 250 °C by combining the UV–Vis and literature solubility data. The resulting logβ n ° are expressed as function of temperature (T in Kelvin): logβ 1 ° = −1.786 + 0.0133 T + 1.049·10 3 /T; logβ 2 ° = −5.797 + 0.0267 T + 2.713·10 3 /T; logβ 3 ° = 6.435 + 0.0223 T + 512.7/T. A comparison between these new fits and existing extrapolations using the Helgeson-Kirkham-Flowers equation of state indicates significant differences in the predicted hydrolysis of La. The latter extrapolations should therefore be updated for the hydrolysis of REE.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Peridynamic neural operators: A data-driven nonlocal constitutive model for complex material responses

Neural operators, which can act as implicit solution operators of hidden governing equations, have recently become popular tools for learning the responses of complex real-world physical systems. Nevertheless, most neural operator applications have thus far been data-driven and neglect the intrinsic preservation of fundamental physical laws in data. Here, in this work, we introduce a novel integral neural operator architecture called the Peridynamic Neural Operator (PNO) that learns a nonlocal constitutive law from data. This neural operator provides a forward model in the form of state-based peridynamics, with objectivity and momentum balance laws automatically guaranteed. As applications, we demonstrate the expressivity and efficacy of our model in learning complex material behaviors from both synthetic and experimental data sets. We also compare the performances with baseline models that use predefined constitutive laws. We show that, owing to its ability to capture complex responses, our learned neural operator achieves improved accuracy and efficiency. Moreover, by preserving the essential physical laws within the neural network architecture, the PNO is robust in treating noisy data. The method shows generalizability to different domain configurations, external loadings, and discretizations.

42 ENGINEERING↗
DOE OSTI2025

CASM Monte Carlo: Calculations of the thermodynamic and kinetic properties of complex multicomponent crystals

Monte Carlo techniques play a central role in statistical mechanics approaches that connect macroscopic thermodynamic and kinetic properties to the electronic structure of a material. This paper describes the implementation of Monte Carlo techniques for the study of multicomponent crystalline materials within the Clusters Approach to Statistical Mechanics (CASM) software suite, and demonstrates their use in model systems to calculate free energies and kinetic coefficients, study phase transitions, and construct phase diagrams from first principles. Many crystal structures are complex, with multiple sublattices occupied by differing sets of chemical species, along with the presence of vacancies or interstitial species. This imposes constraints on concentration variables, the form of thermodynamic potentials, and the values of kinetic transport coefficients. The framework used by CASM to formulate thermodynamic potentials and kinetic transport coefficients accounting for arbitrarily complex crystal structures is presented and demonstrated with examples of increasing complexity. Additionally, an overview of the capabilities of the CASM software specific to Monte Carlo methods is given, and a new CASM software package is introduced, casm-flow, which helps automate the setup, submission, management, and analysis of Monte Carlo simulations.

Cluster expansion↗
DOE OSTI2021

Azo–hydrazone tautomerism in organometallic complexes triggered by a -Re(CO) 3 (L) core: A spectroscopic and theoretical study

The spectroscopic properties and tautomeric behavior of five novel Rhenium (I) tricarbonyl complexes bearing an azo ligand are presented. The organic ligand is stable in solution as the hydrazone tautomer. It remains as hydrazone in different medium conditions (solvent, concentration, pH, etc.) without the formation of detectable amounts of the azo tautomer. However, the complexation of this ligand to the strong electron-withdrawing fragment -Re(CO) 3 X (X = Cl – , Br – ) causes tautomerism to appear in the organic moiety. Two well-defined regions in the electronic spectra for both tautomers were observed, allowing the azo/enol-keto/hydrazone equilibrium to be followed. TD-DFT calculations indicate that for the keto/hydrazone form, the main absorption band is attributed to an IL transition. In contrast, the azo/enol species shows a major contribution of the 1 MLCT (dπ(Re) → dπ*(NN)) transition. The tautomeric equilibrium is easily shifted using solvents with different dielectric constant and hydrogen bond donor/acceptor (HBD/HBA) abilities. Also, the tautomerization process is deeply influenced by the electronic properties of the axial ligand. Therefore, this behavior represents a different strategy for the design of novel materials with optical properties. Furthermore, tautomerism is affected by the concentration of complexes and the presence of water in solvents. Finally, the acid-base behavior of these compounds in a 50% v/v ethanol-buffer system showed that the azo/enol tautomer is stabilized in acidic media. At the same time, an increase of pH promotes tautomerization toward the keto/hydrazone, followed by forming the anionic form at pH > 8.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Spectroscopy and dynamics of the v = 1 levels derived from the OH(D) stretching modes in the I‾∙HDO complex using CW and time-resolved, resonant two-photon infrared excitation of the cryogenically cooled ions

The vibrational energy levels of the two isotopomers adopted by the I‾∙HDO ion-molecule complex occur such that the OH(D) stretching fundamentals span its dissociation energy, thus enabling a spectroscopic investigation of the dynamics displayed by a system prepared in the vicinity of the dissociation threshold. This regime is explored using infrared photoexcitation of the mass-selected complexes cooled in a cryogenic radiofrequency (Paul) ion trap. Survey spectra are obtained at modest resolution using two-color, IR-IR photodissociation with nanosecond lasers to establish the level structure and unimolecular decay dynamics of the v = 1 and 2 levels of the bound OH(D) oscillator. The v = 1 levels are prepared by fixed frequency excitation in the trap and the absorption spectra arising from this excited state are probed by photofragmentation of the complex with a second pulsed IR laser after a variable delay time (0 to 20 ms). The transitions to levels above the dissociation threshold for I‾ + HDO formation are observed to be sharp (~5 cm -1 FWHM). At low pressure, the bound OD (v = 1) population relaxes very slowly (~3 ms), consistent with resonant fluorescence in the IR. The collisional quenching rate constants of this level by the He buffer gas were estimated to be on the order 3 x 10 -12 cm 3 /s based on a crude Stern-Volmer analysis. Here, the rotational fine structure and linewidths (≲0.01 cm -1 FWHM) arising from transitions of the non-bonded OH stretch fundamental of the OD-bound isotopomer that lies just above the dissociation energy are determined using single photon photodissociation by excitation of the 10 K ions with a single-frequency, CW IR laser in the ion trap.

Cryogenic ion infrared spectroscopy↗
DOE OSTI2022

Structure of a dimeric photosystem II complex from a cyanobacterium acclimated to far-red light

Photosystem II (PSII) is the water-splitting enzyme central to oxygenic photosynthesis. To drive water oxidation, light is harvested by accessory pigments, mostly chlorophyll (Chl) a molecules, which absorb visible light (400–700 nm). Some cyanobacteria facultatively acclimate to shaded environments by altering their photosynthetic machinery to additionally absorb far-red light (FRL, 700–800 nm), a process termed far-red light photoacclimation or FaRLiP. During far-red light photoacclimation, FRL-PSII is assembled with FRL-specific isoforms of the subunits PsbA, PsbB, PsbC, PsbD, and PsbH, and some Chl-binding sites contain Chls d or f instead of the usual Chl a. The structure of an apo-FRL-PSII monomer lacking the FRL-specific PsbH subunit has previously been determined, but visualization of the dimeric complex has remained elusive. Here, we report the cryo-EM structure of a dimeric FRL–PSII complex. The site assignments for Chls d and f are consistent with those assigned in the previous apo-FRL-PSII monomeric structure. All sites that bind Chl d or Chl f at high occupancy exhibit a FRL-specific interaction of the formyl moiety of the Chl d or Chl f with the protein environment, which in some cases involves a phenylalanine sidechain. The structure retains the FRL-specific PsbH2 subunit, which appears to alter the energetic landscape of FRL-PSII, redirecting energy transfer from the phycobiliprotein complex to a Chl f molecule bound by PsbB2 that acts as a bridge for energy transfer to the electron transfer chain. Collectively, these observations extend our previous understanding of the structure-function relationship that allows PSII to function using lower energy FRL.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Review of additive manufacturing and densification techniques for the net- and near net-shaping of geometrically complex silicon nitride components

Silicon nitride (Si 3 N 4 ) is an advantageous material due its unique combination of mechanical, thermal, chemical, and electrical properties both at ambient and elevated temperatures. Because of these properties there are a wide range of applications for Si 3 N 4 components. Applications include heat exchangers, environmental barrier coatings, osteointegration scaffolds, radomes, and integrated circuitry. Furthermore, such applications often require geometric complexity for efficient and/or effective operation. However, traditional ceramics processing methods such as hot-pressing or die extrusion are typically limited to simple axis-symmetric shapes. With the advent of additive manufacturing, there has been significant advancement into the forming of geometrically complex Si 3 N 4 components. This review documents additive manufacturing advancements that have demonstrated, or are capable of, fabricating Si 3 N 4 components with complex geometry.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Accounting for uncertainty in complex alluvial aquifer modeling by Bayesian multi-model approach

Alluvial aquifers by nature are complex caused by varied depositional environments. Developing a reliable groundwater model to represent an alluvial aquifer is non-trivial. Also, relying on a single best calibrated model may not be sufficient because of an inadequate choice of model parameter values. To better understand groundwater dynamics and improve model prediction reliability, this study presents a Bayesian multi-model uncertainty quantification (BMMUQ) framework to account for model parameter uncertainty in complex alluvial groundwater modeling. The methodology was applied to the agriculturally intensive Mississippi River alluvial aquifer (MRAA), Northeast Louisiana. An aquifer architecture was first constructed using 7,259 well logs in the MRAA area which covers three fluvial deposits (alluvium, braided-stream terrace, and braided-stream terrace-loess). A 12-layer MODFLOW model was then developed to address the alluvial aquifer complexity and well calibrated through a genetic algorithm. This study quantified model parameter uncertainty in hydraulic conductivity and specific storage of sand facies. Bayesian model averaging (BMA) with the Expectation Maximization (EM) algorithm was adopted to derive posterior model weights and head variances of 50 alternative conceptual groundwater flow models, and thereby obtains BMA ensemble model predictions instead of only relying on the best calibrated conceptual model. Overall, results show that an estimated around 950 million m3 of groundwater storage loss occurs in 2015 with respect to the beginning of 2004, due to high groundwater demand for irrigation in the MRAA area. Explicitly quantifying model uncertainty can produce more reliable groundwater level predictions from BMA ensemble model. The presented groundwater modeling framework improves our understanding of the MRAA and provides a valuable tool to assist agricultural water management.

54 ENVIRONMENTAL SCIENCES↗