Search NASA⌕ Search

SEARCH · Search NASA

Results for “dominant balance”

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 253 records · Page 14

Shear Stress Dependence of Force Networks in 3D Dense Suspensions

The geometric organization and force networks of 3D dense suspensions that exhibit both shear thinning and thickening have been examined as a function of varying strength of interparticle attractive interactions using lubrication flow discrete element simulations. Significant rearrangement of the geometric topology does not occur at either the local or global scale as these systems transition across the shear thinning and shear thickening regimes. In contrast massive rearrangements in the balance of attractive, lubrication, and contact forces are observed with interesting behavior of network growth and competition. In agreement with prior work, in shear thinning regions the attractive force is dominant, however as the shear thickening region is approached there is growth of lubrication forces. Lubrication forces oppose the attraction forces, but as viscosity continues to increase under increasing shear stress, the lubrication forces are dominated by contact forces that also resist attraction. Contact forces are the dominant interactions during shear thickening and are an order of magnitude higher than their values in the shear-thinning regime. At high attractive interaction strength, contact networks can form even under shear thinning conditions, however high shear stress is still required before contact networks become the driving mechanism of shear thickening. Analysis of the contact force network during shear thickening generally indicate a uniformly spreading network that rapidly forms across empty domains, however the growth patterns have structure that are significantly dependent upon the strength of interparticle interactions, indicating subtle variations in the mechanism of shear thickening.

Edens, Lance E.↗

Hysteresis between gas breakdown and plasma discharge

In direct-current (DC) discharge, it is well known that hysteresis is observed between the Townsend (gas breakdown) and glow regimes. Forward and backward voltage sweep is performed using a one-dimensional particle-in-cell Monte Carlo collision (PIC-MCC) model considering a ballast resistor. When increasing the applied voltage after reaching the breakdown voltage (V b ), transition from Townsend to glow discharges is observed. When decreasing the applied voltage from the glow regime, the discharge voltage (V d ) between the anode–cathode gap can be smaller than the breakdown voltage, resulting in a hysteresis, which is consistent with experimental observations. Next, the PIC-MCC model is used to investigate the self-sustaining voltage (V s ) in the presence of finite initial plasma densities between the anode and cathode gap. It is observed that the self-sustaining voltage coincides with the discharge voltage obtained from the backward voltage sweep. In addition, the self-sustaining voltage decreases with increased initial plasma density and saturates above a certain initial plasma density, which indicates a change in plasma resistivity. The decrease in self-sustaining voltage is associated with the electron heat loss at the anode for the low pd (rarefied) regime. In the high pd (collisional) regime, the ion energy loss toward the cathode due to the cathode fall and the inelastic collision loss of electrons in the bulk discharge balance out. Finally, it is demonstrated that the self-sustaining voltage collapses to a singular value, despite the presence of a initial plasma, for microgaps when field emission is dominant, which is also consistent with experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Resolving Salt-Induced Agglomeration of Laponite Suspensions Using X-ray Photon Correlation Spectroscopy and Molecular Dynamics Simulations

Linking the physics of the relaxation behavior of viscoelastic fluids as they form arrested gel states to the underlying chemical changes is essential for developing predictive controls on the properties of the suspensions. In this study, 3 wt.% laponite suspensions are studied as model systems to probe the influence of salt-induced relaxation behavior arising from the assembly of laponite nanodisks. X-ray Photon Correlation Spectroscopy (XPCS) measurements show that laponite suspensions prepared in the presence of 5 mM concentrations of CaCl2, MgCl2 and CsCl salts accelerate the formation of arrested gel states, with CaCl2 having a significant impact followed by CsCl and MgCl2 salts. The competing effects of ion size and charge on relaxation behavior are noted. For example, the relaxation times of laponite suspensions in the presence of Mg2+ ions are slower compared to Cs+ ions despite the higher charge, suggesting that cation size dominates in this scenario. The faster relaxation behavior of laponite suspensions in the presence of Ca2+ ions compared to Cs+ ions shows that a higher charge dominates the size of the ion. The trends in relaxation behavior are consistent with the cluster formation behavior of laponite suspensions and the electrostatic interactions predicted from MD simulations. Charge balance is achieved by the intercalation of the cations at the negatively charged surfaces of laponite suspensions. These studies show that the arrested gel state of laponite suspensions is accelerated in the presence of salts, with ion sizes and charges having a competing effect on relaxation behavior.

36 MATERIALS SCIENCE↗

Observation and simulation of the radial evolution and stream structure of solar wind turbulence

This paper presents a brief overview of the observed evolution in a variety of quantities describing the turbulent evolution of the interplanetary plasma and describes simulation results consistent with many features of the evolution. The turbulence is manifested through a dissipation at small scales in the inner heliosphere with a corresponding evolution in the breakpoint between a relatively flat and a Kolmogoroff spectrum; an evolution from kinetically to (slightly) magnetically dominated energy of the plasma fluctuations; a general decrease in the cross helicity or 'Alfvenicity'; changes in the anisotropy of the fluctuations; and the increasing predominance of quasi-pressure-balanced structures in the compressive component of the fluctuations. MHD simulations with shear layers either side of a central current sheet show that even in the absence of compressibility the lack of a mean field along the direction of the main flow in the current sheet leads to rapid nonlinear evolution and the observed characteristics of 'Elsasser spectra' of the fields in the inner heliosphere. Adding compressibility to the simulations does not greatly change the 'incompressive' quantities but leads in addition to observed correlations between a measure of compression and other quantities.

Roberts, D. A.↗

An Approach to In-Situ Observations of Volcanic Plumes

Volcanoes have long been recognized as playing a dominant role in the birth, and possibly the death, of biological populations. They are possible sources of primordial gases, provide conditions sufficient for creating amino acids, strongly affect the heat balance in the atmosphere, and have been shown to sustain life (in oceanic vents.) Eruptions can have profound effects on local flora and fauna, and for very large eruptions, may alter global weather patterns and cause entire species to fail. Measurements of particulates, gases, and dynamics within a volcanic plume are critical to understanding both how volcanoes work and how plumes affect populations, environment, and aviation. Volcanic plumes and associated eruption columns are a miasma of toxic gases, corrosive condensates, and abrasive particulates that makes them hazardous to nearby populations and poses a significant risk to all forms of aviation. Plumes also provide a mechanism for sampling the volcanic interior, which, for hydrothermal environments, may host unique biological populations.

Smythe, W. D.↗

Structural Properties of [N1888][TFSI] Ionic Liquid: A Small Angle Neutron Scattering and Polarizable Molecular Dynamics Study

In this study, we investigate the quaternary ammonium-based ionic liquid (QAIL), methyltrioctylammonium bis(trifluoromethylsulfonyl)imide, [N 1888 ][TFSI], utilizing small angle neutron scattering (SANS) measurements and polarizable molecular dynamics (MD) simulations to characterize the shortand long-range liquid structure. Scattering structure factors show signatures of three length scales in reciprocal space indicative of alternating polarity (k ~ 0.44 Å –1 ), charge (k ~ 0.75 Å –1 ), and neighboring or adjacent (k ~ 1.46 Å –1 ) domains. Excellent agreement between simulation and experimental scattering structure factors validates various simulation analyses that provide detailed atomistic characterization of the different length scale correlations. The first solvation shell structure is illustrated by obtaining radial, angular, dihedral, and combined distribution functions, where two dominant spatial motifs, N + ···N – and N + ···O – , compete for optimal packing around the polar head of the [N 1888 ] + cation. Intermediate and long-range structures are governed by the balance between local electroneutrality and octyl chain networking, respectively. By computing the charge-correlation structure factor, S ZZ , and the spatial extent of the octyl chain network using graph theory, the bulk-phase structure of [N 1888 ][TFSI] is characterized in terms of electrostatic screening and apolar domain formation length scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Arbuscular mycorrhizal fungi equalize differences in plant fitness and facilitate plant species coexistence through niche differentiation

Mycorrhizal fungi are essential to the establishment of the vast majority of plant species but are often conceptualized with contradictory roles in plant community assembly. On the one hand, host-specific mycorrhizal fungi may allow a plant to be competitively dominant by enhancing growth. On the other hand, host-specific mycorrhizal fungi with different functional capabilities may increase nutrient niche partitioning, allowing plant species to coexist. Here, to resolve the balance of these two contradictory forces, we used a controlled greenhouse study to manipulate the presence of two main types of mycorrhizal fungus, ectomycorrhizal fungi and arbuscular mycorrhizal fungi, and used a range of conspecific and heterospecific competitor densities to investigate the role of mycorrhizal fungi in plant competition and coexistence. We find that the presence of arbuscular mycorrhizal fungi equalizes fitness differences between plants and stabilizes competition to create conditions for host species coexistence. Furthermore, our results show how below-ground mutualisms can shift outcomes of plant competition and that a holistic view of plant communities that incorporates their mycorrhizal partners is important in predicting plant community dynamics.

09 BIOMASS FUELS↗

A Scalable Multigrid Reduction Framework for Multiphase Poromechanics of Heterogeneous Media

Simulation of multiphase poromechanics involves solving a multiphysics problem in which multiphase flow and transport are tightly coupled with the porous medium deformation. To capture this dynamic interplay, fully implicit methods, also known as monolithic approaches, are usually preferred. The main bottleneck of a monolithic approach is that it requires solution of large linear systems that result from the discretization and linearization of the governing balance equations. Because such systems are nonsymmetric, indefinite, and highly ill-conditioned, preconditioning is critical for fast convergence. Recently, most efforts in designing efficient preconditioners for multiphase poromechanics have been dominated by physics-based strategies. Current state-of-the-art “black-box” solvers such as algebraic multigrid (AMG) are ineffective because they cannot effectively capture the strong coupling between the mechanics and the flow subproblems, as well as the coupling inherent in the multiphase flow and transport process. In this work, we develop an algebraic framework based on multigrid reduction (MGR) that is suited for tightly coupled systems of PDEs. Using this framework, the decoupling between the equations is done algebraically through defining appropriate interpolation and restriction operators. One can then employ existing solvers for each of the decoupled blocks or design a new solver based on knowledge of the physics. We demonstrate the applicability of our framework when used as a “black-box” solver for multiphase poromechanics. Here, we show that the framework is flexible to accommodate a wide range of scenarios, as well as efficient and scalable for large problems.

97 MATHEMATICS AND COMPUTING↗

A GPU-based compressible combustion solver for applications exhibiting disparate space and time scales

High-speed chemically active flows pose significant computational challenges due to their disparate space and time scales, with stiff chemistry often dominating simulation time. While modern scientific computing programs achieve exascale performance by leveraging graphics processing units (GPUs), existing GPU-based compressible combustion solvers face critical limitations in memory management, load balancing, and handling the highly localized nature of chemical reactions. To this end, we present a high-performance compressible reacting flow solver built on the AMReX framework and optimized for multi-GPU settings. Here, our approach addresses three GPU performance bottlenecks: memory access patterns through column-major storage optimization, computational workload variability via a bulk-sparse integration strategy for chemical kinetics, and multi-GPU load distribution for adaptive mesh refinement applications. The solver adapts existing matrix-based chemical kinetics formulations to multi-grid contexts. Using representative combustion applications, including 2D and 3D detonations and a 3D jet-in-crossflow configuration, we demonstrate 1.4–5× performance improvements over initial implementations on an in-house cluster of NVIDIA H100 GPUs, and near-ideal weak scaling on the Frontier supercomputer (Oak Ridge Leadership Computing Facility) with up to 1024 AMD Instinct MI250X GPUs. Roofline analysis reveals substantial improvements in arithmetic intensity for both convection (∼ 10 ×) and chemistry (∼ 4 ×) routines, confirming efficient utilization of GPU memory bandwidth and computational resources.

42 ENGINEERING↗

The structure and energy balance of cool star atmospheres

The atmospheric structure and energy balance phenomena associated with magnetic fields in the Sun are reviewed and it is shown that similar phenomena occur in cool stars. The evidence for the weakening or disappearance of transition regions and coronae is discussed together with the appearance of extended cool chromospheres with large mass loss, near V-R = 0.80 in the H-R diagram. Like the solar atmosphere, these atmospheres are not homogeneous and there is considerable evidence for plage regions with bright TR emission lines that overlie dark (presumably magnetic) star spots. The IUE observations are providing important information on the energy balance in these atmospheres that should guide theoretical calculations of the nonradiative heating rate. Recent high dispersion spectra are providing unique information concerning which components of close binary systems are the dominant contributors to the observed emission. A recent unanticipated discovery is that the transition lines are redshifted (an antiwind) in DRa (G2 Ib) and perhaps other stars. Finally, the G and K giants and supergiants are classified into three groups depending on whether their atmospheres are dominated by closed magnetic flux tubes, open field geometries, or a predominately open geometry with a few closed flux tubes embedded.

Linsky, J. L.↗

Mechanics of a gaseous film barrier to lubricant wetting of elastohydrodynamically lubricated conjunctions

Two analytical models, one based on simple hydrodynamic lubrication and the other on soft elastohydrodynamic lubrication, are presented and compared to delineate the dominant physical parameters that govern the mechanics of a gaseous film between a small droplet of lubricant and the outer race of a ball bearing. Both models are based on the balance of gravity forces, air drag forces, and air film lubrication forces and incorporate a drag coefficient C sub D and a lubrication coefficient C sub L to be determined from experiment. The soft elastohydrodynamic lubrication (EHL) model considers the effects of droplet deformation and solid-surface geometry; the simpler hydrodynamic lubrication (HL) model assumes that the droplet remains essentially spherical. The droplet's angular position depended primarily on the ratio of gas inertia to droplet gravity forces and on the gas Reynolds number and weakly on the ratio of droplet gravity forces to surface tension forces (Bond number) and geometric ratios for the soft EHL. An experimental configuration in which an oil droplet is supported by an air film on the rotating outer race of a ball bearing within a pressure-controlled chamber produced measurements of droplet angular position as a function of outer-race velocity droplet size and type, and chamber pressure.

Prahl, J. M.↗

Coordination and control – limits in standard representations of multi-reservoir operations in hydrological modeling

Abstract. Major multi-reservoir cascades represent a primary mechanism for dealing with hydrologic variability and extremes within institutionally complex river basins worldwide. These coordinated management processes fundamentally reshape water balance dynamics. Yet, multi-reservoir coordination processes have been largely ignored in the increasingly sophisticated representations of reservoir operations within large-scale hydrological models. The aim of this paper is twofold, namely (i) to provide evidence that the common modeling practice of parameterizing each reservoir in a cascade independently from the others is a significant approximation and (ii) to demonstrate potential unintended consequences of this independence approximation when simulating the dynamics of hydrological extremes in complex reservoir cascades. We explore these questions using the Water Balance Model, which features detailed representations of the human infrastructure coupled to the natural processes that shape water balance dynamics. It is applied to the Upper Snake River basin in the western US and its heavily regulated multi-reservoir cascade. We employ a time-varying sensitivity analysis that utilizes the method of Morris factor screening to explicitly track how the dominant release rule parameters evolve both along the cascade and in time according to seasonal high- and low-flow events. This enables us to address aim (i) by demonstrating how the progressive and cumulative dominance of upstream releases significantly dampens the ability of downstream reservoir rules' parameters to influence flow conditions. We address aim (ii) by comparing simulation results with observed reservoir operations during critical low-flow and high-flow events in the basin. Our time-varying parameter sensitivity analysis with the method of Morris clarifies how independent single-reservoir parameterizations and their tacit assumption of independence leads to reservoir release behaviors that generate artificial water shortages and flooding, whereas the observed coordinated cascade operations avoided these outcomes for the same events. To further explore the role of (non-)coordination in the large deviations from the observed operations, we use an offline multi-reservoir water balance model in which adding basic coordination mechanisms drawn from the observed emergency operations is sufficient to correct the deficiencies of the independently parameterized reservoir rules from the hydrological model. These results demonstrate the importance of understanding the state–space context in which reservoir releases occur and where operational coordination plays a crucial role in avoiding or mitigating water-related extremes. Understanding how major infrastructure is coordinated and controlled in major river basins is essential for properly assessing future flood and drought hazards in a changing world.

54 ENVIRONMENTAL SCIENCES↗

Influence of interfacial responses of Berea Sandstone in low-salinity waterflooding environments

Wetting conditions and fluid-fluid interfacial properties are essential to understanding oil recovery in the context of low-salinity and smart waterflooding. Many have deemed wettability alteration to be the dominant interaction, while some highlight the importance of fluid-fluid interactions. In some rock-fluid systems, benefits from the use of additives can arise to favor hydrocarbon production. Understanding the balance between the fluid-fluid and rock-fluid interactions when utilizing additives is crucial to evaluate benefits in improved efficiency. Ultimately, this understanding will shed light on recovery mechanisms. Regarding oil recovery, previous studies indicate that the injection of a naphthenic acid (NA) blend dissolved in an injection brine improves oil recovery in Berea Sandstone, particularly under low-salinity conditions. Interfacial rheology results for the brine-oil interface correlated well with the recovery response. This observation was speculated to serve as a potential improved recovery mechanism through oil-trapping suppression. However, the said studies did not evaluate potentially concurrent rock wettability alteration during recovery or the initial wetting condition. The purpose of this work is to determine the production behavior in systems with NAs, namely NA-10 (3-Cyclohexanepropinoic acid), NA-14 (Cyclopentane carboxylic acid), NA-18 (Cyclohexane carboxylic acid), and NA-mix (Naphthenic acid mixture), at different initial wetting conditions under spontaneous imbibition. In addition, we tested the NAs ability to alter wettability in Berea Sandstone and determine each “additive’s” capacity to improve recovery. Here, time-domain nuclear magnetic resonance (TD-NMR), contact angle, and Amott cell tests, i.e. spontaneous imbibition, were used to analyze the influence of the selected naphthenic acids on oil displacement under different initial wetting conditions. Our results indicate that individual naphthenic acids indeed have potential as additives in low-salinity waterflooding for oil-wet systems while negatively impacting recovery in water-wet systems. However, results suggest that the evolution of the interfacial properties and wettability conditions are equally important on the recovery behavior.

54 ENVIRONMENTAL SCIENCES↗

An accelerated approach to inline non-LTE modeling

A goal across many physics and engineering disciplines is to use radiation-hydrodynamics codes to accurately simulate systems that are expected to deviate from local thermodynamic equilibrium (LTE), such as laser-driven hohlraums, stellar coronae, and supernova ejecta. The computational expense of traditional non-LTE (NLTE) opacity calculations can dominate the cost of simulating these systems and render NLTE physics untenable in most cases. In this paper, we present a new parallel implementation for NLTE opacity calculations that employs advanced load balancing, vectorization, and memory management to significantly reduce the computational costs of these calculations. Using the radiation-hydrodynamics code xRAGE , we demonstrate this new implementation on three test problems and present order-of-magnitude speed improvements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Doped GaInAs/GaP Quantum Well Superlattice Solar Cells With 27.5% Efficiency

Quantum wells can extend the absorption range of a solar cell and are typically placed in the intrinsic region of the device to enable efficient carrier collection via drift. Thick intrinsic regions are needed for significant absorption in the low bandgap quantum wells, resulting in a large depletion region recombination and ultimately a solar cell with a low fill factor (FF). However, in quantum well superlattice solar cells where tunneling plays a dominant role in carrier transport, collection by carrier diffusion may be possible, relieving the requirement for quantum wells to be placed in the intrinsic region. Here, we investigate doping in stress-balanced quantum well superlattice solar cells using thin 2 nm GaP barriers. Doping reduces J02 depletion region recombination and improves the solar cell FF up to 86.7%, but very high doping eventually reduces the carrier collection, leading to a tradeoff in efficiency. Finally, we show that the barrier thickness also plays an important role in carrier collection, and we demonstrate high efficiency 27.5% single-junction devices with a doped superlattice.

13 HYDRO ENERGY↗

The need for carbon-emissions-driven climate projections in CMIP7

Abstract. Previous phases of the Coupled Model Intercomparison Project (CMIP) have primarily focused on simulations driven by atmospheric concentrations of greenhouse gases (GHGs), for both idealized model experiments and climate projections of different emissions scenarios. We argue that although this approach was practical to allow parallel development of Earth system model simulations and detailed socioeconomic futures, carbon cycle uncertainty as represented by diverse, process-resolving Earth system models (ESMs) is not manifested in the scenario outcomes, thus omitting a dominant source of uncertainty in meeting the Paris Agreement. Mitigation policy is defined in terms of human activity (including emissions), with strategies varying in their timing of net-zero emissions, the balance of mitigation effort between short-lived and long-lived climate forcers, their reliance on land use strategy, and the extent and timing of carbon removals. To explore the response to these drivers, ESMs need to explicitly represent complete cycles of major GHGs, including natural processes and anthropogenic influences. Carbon removal and sequestration strategies, which rely on proposed human management of natural systems, are currently calculated in integrated assessment models (IAMs) during scenario development with only the net carbon emissions passed to the ESM. However, proper accounting of the coupled system impacts of and feedback on such interventions requires explicit process representation in ESMs to build self-consistent physical representations of their potential effectiveness and risks under climate change. We propose that CMIP7 efforts prioritize simulations driven by CO2 emissions from fossil fuel use and projected deployment of carbon dioxide removal technologies, as well as land use and management, using the process resolution allowed by state-of-the-art ESMs to resolve carbon–climate feedbacks. Post-CMIP7 ambitions should aim to incorporate modeling of non-CO2 GHGs (in particular, sources and sinks of methane and nitrous oxide) and process-based representation of carbon removal options. These developments will allow three primary benefits: (1) resources to be allocated to policy-relevant climate projections and better real-time information related to the detectability and verification of emissions reductions and their relationship to expected near-term climate impacts, (2) scenario modeling of the range of possible future climate states including Earth system processes and feedbacks that are increasingly well-represented in ESMs, and (3) optimal utilization of the strengths of ESMs in the wider context of climate modeling infrastructure (which includes simple climate models, machine learning approaches and kilometer-scale climate models).

54 ENVIRONMENTAL SCIENCES↗

Infrared coronal emission lines and the possibility of their maser emission in Seyfert nuclei

Energetic emitting regions have traditionally been studied via x-ray, UV and optical emission lines of highly ionized intermediate mass elements. Such lines are often referred to as 'coronal lines' since the ions, when produced by collisional ionization, reach maximum abundance at electron temperatures of approx. 10(exp 5) - 10(exp 6) K typical of the sun's upper atmosphere. However, optical and UV coronal lines are also observed in a wide variety of Galactic and extragalactic sources including the Galactic interstellar medium, nova shells, supernova remnants, galaxies and QSOs. Infrared coronal lines are providing a new window for observation of energetic emitting regions in heavily dust obscured sources such as infrared bright merging galaxies and Seyfert nuclei and new opportunities for model constraints on physical conditions in these sources. Unlike their UV and optical counterparts, infrared coronal lines can be primary coolants of collisionally ionized plasmas with 10(exp 4) less than T(sub e)(K) less than 10(exp 6) which produce little or no optical or shorter wavelength coronal line emission. In addition, they provide a means to probe heavily dust obscured emitting regions which are often inaccessible to optical or UV line studies. In this poster, we provide results from new model calculations to support upcoming Infrared Space Observatory (ISO) and current ground-based observing programs involving infrared coronal emission lines in AGN. We present a complete list of infrared (lambda greater than 1 micron) lines due to transitions within the ground configurations 2s(2)2p(k) and 3s(2)3p(k) (k = 1 to 5) or the first excited configurations 2s2p and 3s3p of highly ionized (x greater than or equal to 100 eV) astrophysically abundant (n(X)/n(H) greater than or equal to 10(exp -6)) elements. Included are approximately 74 lines in ions of O, Ne, Na, Mg, Al, Si, S, Ar, Ca, Fe, and Ni spanning a wavelength range of approximately 1 - 280 microns. We present new results from detailed balance calculations, new critical densities for collisional de-excitation, intrinsic photon rates, branching ratios, and excitation temperatures for the majority of the compiled transitions. The temperature and density parameter space for dominant cooling via infrared coronal lines is presented, and the relationship of infrared to optical coronal lines is discussed.

Greenhouse, Matthew A.↗

Global Ocean Evaporation: How Well Can We Estimate Interannual to Decadal Variability?

Evaporation from the world's oceans constitutes the largest component of the global water balance. It is important not only as the ultimate source of moisture that is tied to the radiative processes determining Earth's energy balance but also to freshwater availability over land, governing habitability of the planet. Here we focus on variability of ocean evaporation on scales from interannual to decadal by appealing to three sources of data: the new MERRA-2 (Modern-Era Retrospective analysis for Research and Applications -2); climate models run with historical sea-surface temperatures, ice and atmospheric constituents (so-called AMIP experiments); and state-of-the-art satellite retrievals from the Seaflux and HOAPS (Hamburg Ocean-Atmosphere Parameters and Fluxes from Satellite) projects. Each of these sources has distinct advantages as well as drawbacks. MERRA-2, like other reanalyses, synthesizes evaporation estimates consistent with observationally constrained physical and dynamical models-but data stream discontinuities are a major problem for interpreting multi-decadal records. The climate models used in data assimilation can also be run with lesser constraints such as with SSTs and sea-ice (i.e. AMIPs) or with additional, minimal observations of surface pressure and marine observations that have longer and less fragmentary observational records. We use the new ERA-20C reanalysis produced by ECMWF embodying the latter methodology. Still, the model physics biases in climate models and the lack of a predicted surface energy balance are of concern. Satellite retrievals and comparisons to ship-based measurements offer the most observationally-based estimates, but sensor inter-calibration, algorithm retrieval assumptions, and short records are dominant issues. Our strategy depends on maximizing the advantages of these combined records. The primary diagnostic tool used here is an analysis of bulk aerodynamic computations produced by these sources and uses a first-order Taylor series analysis of wind speed, SST, near-surface stability and relative humidity variations around climatology to gauge the importance of these components. We find that the MERRA-2 evaporation record is strongly influenced by the availability of wind speed and humidity from passive microwave imagers beginning in the late 1980s as well as by the SST record. The trend over the period 1980 to present is nearly 10%. AMIP or the ERA-20C trends are much smaller. We find that ENSO-related signals involving both wind speed and thermodynamic variability remain the primary signal in the latter and are confirmed by satellite retrievals. We present uncertainty estimates based on the various data sources and discuss the implications for GEWEX water and energy budget science challenges.

Robertson, Franklin R.↗