SEARCH · Engineering Papers

Results for “lifetime model”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,459 records · Page 81

All sources NASA NTRS DOE OSTI

DOE OSTI2022

NIR fluorescence lifetime macroscopic imaging with a time-gated SPAD camera

The performance of SwissSPAD2 (SS2), a large scale, widefield time-gated CMOS SPAD imager developed for fluorescence lifetime imaging, has recently been described in the context of visible range and fluorescence lifetime imaging microscopy (FLIM) of dyes with lifetimes in the 2.5 – 4 ns range. Here, we explore its capabilities in the NIR regime relevant for small animal imaging, where its sensitivity is lower and typical NIR fluorescent dye lifetimes are much shorter (1 ns or less). We carry out this study in a simple macroscopic imaging setup based on a compact NIR picosecond pulsed laser, an engineered diffuser-based illumination optics, and NIR optimized imaging lens suitable for well-plate or small animal imaging. Because laser repetition rates can vary between models, but the synchronization signal frequency accepted by SS2 is fixed to 20 MHz, we first checked that a simple frequency-division scheme enables data recording for different laser repetition rates. Next, we acquired data using different time gate widths, including gates with duration longer than the laser period, and analyzed the resulting data using both standard nonlinear least-square fit (NLSF) and phasor analysis. We show that the fixed synchronization rate and large gate widths characterizing SS2 (10 ns and over) are not an obstacle to accurately extracting lifetime in the 1 ns range and to distinguishing between close lifetimes. As a result, SS2 and similar very large gated SPAD imagers appear as a versatile alternative to other widefield time-resolved detectors for NIR fluorescence lifetime imaging, including preclinical molecular applications.

FLIM↗

DOE OSTI2026

Overview of the ASDEX Upgrade results

After a 26-month vent ASDEX Upgrade (AUG) went back in operation with a newly designed upper W-divertor suitable for alternative divertor configurations (featuring in-vessel coils and cryo-pump). Parameter scans and an extensive set of measurements were obtained and their interpretation is ongoing. Prompted by the ITER wall change, dedicated experiments on non-boronized plasma startup were contrasted to that employing asymmetric and more symmetric boronizations. The asymmetric boronization proved to be as beneficial as the more symmetric one, which is in contrast to previous model calculations assuming perfect sticking of boron (measurements suggest sticking ≈ 0.3). In the startup phase also the impurity influxes at the outboard limiters were investigated contrasting the unboronized case featuring cold edges (low- Z radiation) to the boronized case, in which the lifetime of the boron layers could be estimated. Pedestal stability investigations revealed that the quasi-continuous exhaust (QCE) regime is obtained when ballooning modes are active in the vicinity of the separatrix and the global peeling-ballooning stability is high enough. Thus, at high enough shaping and high gas flux both can be achieved and QCE is a consequence. The closely related enhanced D-alpha (EDA) mode is not clearly distinguishable from QCE, e.g. the quasi-coherent mode characteristic for EDA also shows up in QCE. In QCE the impurity transport is behaving benign as could be measured for Ne with a novel analysis method making use of a comprehensive set of CXRS measurements. For high radiative fractions the regime of the X-point radiator (XPR) is accessible at AUG and the understanding of its access conditions and behaviour is further developed. Due to the localized radiative cooling at the X-point the XPR can be well diagnosed and thus controlled. For negative triangularity shapes, further experiments at increased shaping resulted in strongly heated L-mode plasmas avoiding ELMs. Two integrated modelling approaches towards ITER suggested that core W-accumulation will be no issue for ITER and that the fusion yield in ITER may be Q = 12 (i.e. ITPA20-IL scaling is too pessimistic). Further, investigations of the ITER ramp-down in AUG provide insights into maintaining position control. Various aspects of shattered pellet injection were investigated in AUG and one of the results show that with increasing Ne fraction the radiation during the current quench increases and the current decay becomes faster.

ASDEX Upgrade↗

DOE OSTI2023

The Dynamics of Megafire Smoke Plumes in Climate Models: Why a Converged Solution Matters for Physical Interpretations

As the climate system warms, megafires have become more frequent with devastating effects. A byproduct of these events is the creation of smoke plumes that can rise into the stratosphere and spread across the globe where they reside for many months. To gain a deeper understanding of the plume dynamics, global climate simulations of a megafire were performed at a wide range of grid spacings from 2.0° down to 7 km, including a 7 km nonhydrostatic experiment. The analysis focuses on how the resolved dynamics affects the specification of the plume characteristics such as injection height and black carbon (BC) mass. Prior studies initialize the smoke plume at one or a few grid points and this is shown here to produce severely dissipative dynamics. In order to validate such simulations with observations, enhancements of the plume characteristics to offset the dissipation is necessary. Using a numerically converged simulation, sensitivity tests show that to approximate the observed stratospheric lifetime, a reduction in BC fraction by 50% is necessary for external mixtures. The vorticity dynamics of the plume is also analyzed with a Lagrangian budget to understand the mechanisms responsible for the evolution of a collocated anticyclonic vortex. The results can be distilled down into a simple conceptual model. As the plume rises, the air diverges at the top of the updraft where the largest concentrations of smoke are found. This divergence induces a dilution of the background cyclonic absolute vorticity producing an anticyclonic vortex. Vortex decay occurs from opposite arguments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Modeling and Understanding of Rear Junction Double-Side Passivated Contact Solar Cells with Selective Area TOPCon on Front

Device modeling is performed to propose > 25% efficient industry-compatible rear junction double-side passivated contacts solar cell structure with full area p-TOPCon on the rear and selective area n-TOPCon under the front grid pattern (selective TOPCon). Here, this design enables the use of thicker TOPCon (>100nm) on the front for traditional screen-printed contacts without incurring metal-induced damage, high parasitic absorption loss, and compromise in lateral transport or carrier collection on the front. Rear junction design with appropriate bulk lifetime and resistivity combination eliminate the need for heavy doping in the front field region because carriers can flow through the bulk Si without appreciable FF loss. High VOC is maintained because high-quality Si surface passivation in the field region by Al 2 O 3 /SiN gives J 0 comparable to the TOPCon. Our device modeling specifies the practically achievable properties and parameters for each region, including full area rear p-TOPCon, selective area front n-TOPCon, bulk and contacts, to achieve 25.4% efficiency screen-printed bifacial rear junction selective TOPCon cells.

14 SOLAR ENERGY↗

NASA NTRS2004

Starting Conditions for Hydrothermal Systems Underneath Martian Craters: Hydrocode Modeling

Mars is the most Earth-like of the Solar System s planets, and the first place to look for any sign of present or past extraterrestrial life. Its surface shows many features indicative of the presence of surface and sub-surface water, while impact cratering and volcanism have provided temporary and local surface heat sources throughout Mars geologic history. Impact craters are widely used ubiquitous indicators for the presence of sub-surface water or ice on Mars. In particular, the presence of significant amounts of ground ice or water would cause impact-induced hydrothermal alteration at Martian impact sites. The realization that hydrothermal systems are possible sites for the origin and early evolution of life on Earth has given rise to the hypothesis that hydrothermal systems may have had the same role on Mars. Rough estimates of the heat generated in impact events have been based on scaling relations, or thermal data based on terrestrial impacts on crystalline basements. Preliminary studies also suggest that melt sheets and target uplift are equally important heat sources for the development of a hydrothermal system, while its lifetime depends on the volume and cooling rate of the heat source, as well as the permeability of the host rocks. We present initial results of two-dimensional (2D) and three-dimensional (3D) simulations of impacts on Mars aimed at constraining the initial conditions for modeling the onset and evolution of a hydrothermal system on the red planet. Simulations of the early stages of impact cratering provide an estimate of the amount of shock melting and the pressure-temperature distribution in the target caused by various impacts on the Martian surface. Modeling of the late stage of crater collapse is necessary to characterize the final thermal state of the target, including crater uplift, and distribution of the heated target material (including the melt pool) and hot ejecta around the crater.

Pierazzo, E.↗

NASA NTRS2008

Assessment and Consequences of the Delayed Breakup of the Antarctic Polar Vortex in Two Versions of the GEOS Chemistry-Climate Model

In mid-winter, winds circle the globe at speeds greater than 200 km/hr (approximately 130mph) in the middle atmosphere. This strong jet bounds the region known as the polar vortex. The presence of the Antarctic polar vortex is a key ingredient in the formation of the 'ozone hole', because the air inside the vortex is cold and isolated from lower latitudes, creating ideal conditions for large-scale chemical ozone depletion. Many atmospheric models are not able to reproduce observed winds in the middle atmosphere. Specifically, the polar vortices tend to break down too late and peak wind speeds are higher than observed. Hurwitz et al. find that the delayed break-up of the Antarctic polar vortex is due to weaker-than-observed wave driving from the lower atmosphere during the October-November period. The delayed break-up of the Antarctic polar vortex changes the temperature structure of the middle atmosphere, which biases the amount of chemical ozone depletion that can occur in late winter and spring. Also, the extended lifetime of the polar vortex strengthens the 'overturning' circulation cell in the middle atmosphere, changing the amount of ozone, methane and other chemical species that is transported from low to high latitudes. As greenhouse gas concentrations continue to rise, the atmospheric temperature structure and resulting wind structure are expected to change. Clearly, if models cannot duplicate the observed late 20th century high-latitude winds, their ability to simulate the polar vortices in future must be poor. Understanding model weaknesses and improving the modeled polar vortices will be necessary for accurate predictions of ozone recovery in the coming century.

Hurwitz, M.M.↗

DOE OSTI2026

Prognostic simulations of mixed-phase clouds with model AC-1D v1.0: the impact of aerosol types and freezing parameterizations on ice crystal budgets

Mixed-phase clouds at high latitudes contribute to the uncertainty in predicting cloud feedbacks and climate sensitivity, mainly due to the complexity of microphysical processes that influence the partitioning between the supercooled liquid and ice phases, and hence, cloud radiative effects on regional scales. Particularly in Arctic mixed-phase clouds, the activation of ice-nucleating particles (INPs) from various aerosol populations remains a leading source of uncertainty. We developed an aerosol-cloud one-dimensional (AC-1D) model, which provides a novel framework to prognostically treat INP and ice crystal budgets while explicitly accounting for polydisperse and multicomponent aerosol that activate INPs following different freezing parameterizations. The AC-1D model is informed by large-eddy simulations to probe the impact of INP representation on predicted ice crystal number concentrations (N i ) and ice crystal budgets in mixed-phase Arctic stratus. We apply three immersion freezing (IMF) parameterizations, two time-independent (singular) and one time-dependent (classical nucleation theory), to predict the evolution of the INP reservoir and resulting ice crystal budget from polydisperse mineral dust, organic (humic-like substances), and sea spray aerosol particle size distributions. Our analysis focuses on how variations in aerosol number concentration and cloud system parameters such as cloud cooling rate, cloud-top entrainment rate, and ice crystal fall speed influence the INP reservoir and ice crystal budgets. Furthermore, this study investigates the competitive ice nucleation dynamics in mixed aerosol environments and provides a process-level quantification of the INP budget terms, which directly controls ice crystal budgets. For all studied case scenarios, the aerosol types and associated particle size distributions significantly impact INP and N i , and the choice between a time-dependent and a singular freezing description yields orders-of-magnitude differences in the predicted INP and N i over the 10 h simulation time, reflecting typical cloud lifetimes. Our results show that the influence of cloud cooling, INP entrainment, and sedimentation varies significantly depending on the chosen freezing parameterization. These findings underscore the critical need for robust IMF parameterizations and precise cloud system observations to enhance the accuracy of models in predicting mixed-phase cloud structure and evolution.

Arctic clouds↗

DOE OSTI2023

Impact of Global Warming on U.S. Summertime Mesoscale Convective Systems: A Simple Lagrangian Parcel Model Perspective

Abstract Mesoscale convective systems (MCSs) are the dominant rainfall producer over the United States during the warm season, causing natural disasters and severe weather every year. Global climate models have large uncertainty in projecting precipitation changes in the future climate. Here, a simple Lagrangian parcel model is used to investigate the impact of global warming on MCS initiation and growth. The single-column parcel model projects a mean precipitation decrease over the central United States and an increase to its east, in agreement with the CMIP5 model projection. It also highlights the crucial role of current climate mean-state model bias in exaggerating the change in future mean precipitation projection by 25%. As for convective population, the model captures the decreased occurrence frequency of weak to moderate convection and increased frequency of strong convection due to the increased CAPE and CIN, in agreement with convection-permitting regional simulations. Novel parameterizations of gust-front propagation speed and subsidence strength are developed as guided by cloud-resolving simulations. The multicolumn parcel model employing those parameterizations captures readily the cold pool–induced upscale growth feature. It simulates smaller mesoscale clusters over the central United States under global warming due to gust-front slowdown and subsidence strength enhancement, which are further attributed to land aridity–induced weakening of initial accumulated precipitation and strengthening of updraft speed, respectively. That said, mesoscale clusters could become bigger under more favorable conditions in future climate, including boundary layer moistening, convection lifetime lengthening, and cold pool mechanical-lifting enhancement, which require further investigations to improve mechanistic understanding of future MCS changes.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2017

Development of a 13 kW Hall Thruster Propulsion System Performance Model for AEPS

The Advanced Electric Propulsion System (AEPS) program will develop a flight 13kW Hall thruster propulsion system based on NASA's HERMeS thruster. The AEPS system includes the Hall Thruster, the Power Processing Unit (PPU) and the Xenon Flow Controller (XFC). These three primary components must operate together to ensure that the system generates the required combinations of thrust and specific impulse at the required system efficiencies for the desired system lifetime. At the highest level, the AEPS system will be integrated into the spacecraft and will receive power, propellant, and commands from the spacecraft. Power and propellant flow rates will be determined by the throttle set points commanded by the spacecraft. Within the system, the major control loop is between the mass flow rate and thruster current, with time-dependencies required to handle all expected transients, and additional, much slower interactions between the thruster and cathode temperatures, flow controller and PPU. The internal system interactions generally occur on shorter timescales than the spacecraft interactions, though certain failure modes may require rapid responses from the spacecraft. The AEPS system performance model is designed to account for all these interactions in a way that allows evaluation of the sensitivity of the system to expected changes over the planned mission as well as to assess the impacts of normal component and assembly variability during the production phase of the program. This effort describes the plan for the system performance model development, correlation to NASA test data, and how the model will be used to evaluate the critical internal and external interactions. The results will ensure the component requirements do not unnecessarily drive the system cost or overly constrain the development program. Finally, the model will be available to quickly troubleshoot any future unforeseen development challenges.

Stanley, Steven↗

DOE OSTI2022

Spark plasma sintered, MoNbTi-based multi-principal element alloys with Cr, V, and Zr

MoNbTi, MoNbTiZr, CrMoNbTiZr, and MoNbTiVZr multi-principal element alloys (MPEAs) were fabricated via spark plasma sintering (SPS) and investigated for use in high-strength applications. The fabrication method by SPS and powder metallurgy differs from those presented in the prior literature, where most MPEAs are fabricated using arc melting (AM) methods. Here cryogenic milling was used to maximize potential defect sinks (grain boundaries) for radiation resistance and to increase ultimate tensile strength through the Hall-Petch effect. SPS was chosen for consolidation in order to maintain a fine-grained structure during densification. Each alloy was characterized using x-ray diffraction and scanning electron microscopy for phase identification and compositional homogeneity. The base ternary alloy MoNbTi presented a predominantly single BCC system with minor cubic phases. The introduction of additional alloying elements—Zr, V, and Cr—heightened the phase complexity and increased the fractions of a secondary BCC phase and an HCP phase from Zr. The addition of Cr induced a larger fraction of the Laves phase to form. The addition of V caused the precipitation of small Mo inclusions. Thermodynamic analysis was performed to understand the separation of phases in each alloy. Discrepancies among the phase predictions generated by thermodynamic models, phases previously presented in the literature, and the characterization results suggest that MPEA fabrication methods, especially solid-state methods, require significant investigation to ensure that alloys can remain stable throughout their anticipated service lifetimes.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Topological protection of coherence in disordered open quantum systems

Here, we consider topological protection mechanisms in dissipative quantum systems in the presence of quenched disorder, with the intent to prolong the coherence time of a fiducial qubit. The qubit is part of a network of other qubits and dissipative cavities whose coupling parameters are tunable, such that topological edge states can be stabilized. The evolution of the fiducial qubit is entirely determined by a non-Hermitian Hamiltonian which thus emerges from a bona fide physical process. Even in the presence of disorder, a winding number W can be defined and evaluated in real space, as long as certain symmetries are preserved. Hence we can construct the topological phase diagrams of noisy open quantum models, such as the non-Hermitian disordered Su-Schrieffer-Heeger dimer model and a trimer model that includes longer-range couplings. For finite-size systems we find that there are precisely W modes localized at one end of the chain. In such topological phases the qubit's coherence lifetime is exponentially large in the system size. In the presence of competing disorder parameters, interesting reentrance phenomena of topologically nontrivial sectors are observed. This means that in certain parameter regions, increasing disorder drastically increases the coherence time of the fiducial qubit.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Unraveling higher-order contributions to spin excitations probed using resonant inelastic x-ray scattering

Resonant inelastic x-ray scattering (RIXS) is an evolving tool for investigating the spin dynamics of strongly correlated materials, which complements inelastic neutron scattering. In isotropic spin-$\frac{1}{2}$ Heisenberg antiferromagnetic (HAFM) spin chains, both techniques have observed non-spin-conserving (NSC) excitations confined to the two-spinon phase space. However, a recent O K-edge RIXS study of the one-dimensional HAFM Sr 2 CuO 3 observed spin-conserving (SC) four-spinon excitations outside the two-spinon phase space. Here, we demonstrate that analogous four-spinon excitations can also be accessed at the Cu L 3 edge in the related material SrCuO 2 . Through detailed modeling, we establish that these excitations appear in both the SC and NSC channels of the Cu L 3 edge, and are only captured by higher-order terms in the ultrashort core-hole lifetime expansion. Finally, since these terms encode information about spin-spin correlations extending beyond nearest neighbors, our results offer different possibilities for studying nonlocal spin correlations in quantum magnets.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Attosecond Time-Domain Measurement of Core-Level-Exciton Decay in Magnesium Oxide

Excitation of ionic solids with extreme ultraviolet pulses creates localized core-level excitons, which in some cases couple strongly to the lattice. In this paper, core-level-exciton states of magnesium oxide are studied in the time domain at the Mg L 2 , 3 edge with attosecond transient reflectivity spectroscopy. Attosecond pulses trigger the excitation of these short-lived quasiparticles, whose decay is perturbed by time-delayed near-infrared pulses. Combined with a few-state theoretical model, this reveals that the infrared pulse shifts the energy of bright (dipole-allowed) core-level-exciton states as well as induces features arising from dark core-level excitons. We report coherence lifetimes for the two lowest core-level excitons of 2.3 ± 0.2 and 1.6 ± 0.5 fs and show that these are primarily a consequence of strong exciton-phonon coupling, disclosing the drastic influence of structural effects in this ultrafast relaxation process.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

3D Frequency Domain Reflectometry Digital Twin of an Electrical Cable: A First Glance

Electrical cables within nuclear power plants (NPPs) are critical components required for power, control, and instrumentation systems which may be exposed to stressors, such as elevated temperatures and gamma radiation. Such stressors can lead to a reduction in the remaining useful life of electrical cables, jeopardizing the safety of NPP systems. To evaluate the effect of stressors on the degradation of electrical cables, electrical reflectometry methods are commonly employed. Frequency domain reflectometry (FDR) is a non-destructive electrical reflectometry method that uses transmission line theory to detect degradation or impedance changes within electrical cables. However, in most cases FDR is only applied to de-energized cables, limiting the application in NPPs as the cable system must be taken offline. In this work, we explore the development of an FDR digital twin to predict the degradation of an electrical cable exposed to elevated temperature, which is expected to reduce the need for offline FDR. A 3-conductor low-voltage electrical cable was selected for evaluation of the digital twin. The fully three-dimensional digital twin was developed in COMSOL using the RF module. A cable length of 30-m and frequency bandwidth of 400 MHz was selected to mimic real-world application of FDR. Over a 1-m region, the permittivity of the insulation was varied by up to 20% to model thermal degradation. The results demonstrate accurate detection of the insulation damage region, supporting further investigation of the FDR digital twin using real-world data and machine learning for predictive damage estimation or remaining lifetime.

Spencer, Mychal P.↗

DOE OSTI2021

Connect the Dots: In Situ 4-D Seismic Monitoring of CO 2 Storage With Spatio-Temporal CNNs

4-D seismic imaging has been widely used in CO 2 sequestration projects to monitor the fluid flow in the volumetric subsurface region that is not sampled by wells. Ideally, real-time monitoring and near-future forecasting would provide site operators with great insights to understand the dynamics of the subsurface reservoir and assess any potential risks. However, due to obstacles such as high deployment cost, availability of acquisition equipment, exclusion zones around surface structures, only very sparse seismic imaging data can be obtained during monitoring. That leads to an unavoidable and growing knowledge gap over time. The operator needs to understand the fluid flow throughout the project lifetime and the seismic data are only available at a limited number of times. This is insufficient for understanding reservoir behavior. To overcome those challenges, we have developed spatio-temporal neural-network-based models that can produce high-fidelity interpolated or extrapolated images effectively and efficiently. Specifically, our models are built on an autoencoder, and incorporate the long short-term memory (LSTM) structure with a new loss function regularized by optical flow. We validate the performance of our models using real 4-D post-stack seismic imaging data acquired at the Sleipner CO 2 sequestration field. We employ two different strategies in evaluating our models. Numerically, we compare our models with different baseline approaches using classic pixel-based metrics. We also conduct a blind survey and collect a total of 20 responses from domain experts to evaluate the quality of data generated by our models. Finally, via both numerical and expert evaluation, we conclude that our models can produce high-quality 2-D/3-D seismic imaging data at a reasonable cost, offering the possibility of real-time monitoring or even near-future forecasting of the CO 2 storage reservoir.

4-D seismic imaging↗

NASA NTRS1972

Analysis of stress-rupture data from S-glass composites.

Both composite and free strands have been tested to study matrix effectiveness. Glass from different lots were used, and processing parameters were varied to demonstrate the effects of material variability. The statistical patterns that emerge from these experiments show that lifetime distributions are characteristically skewed. As a class, the lifetime distributions are well approximated by the reduced form of the Weibull distribution. The strength retention data display little early degradation. The retained-strength distributions can also be classed as Weibull distributions. A statistical model applied to the data gives predictions of progressive failure, strength retention path, and expected life for two basic types of internal fracture processes.

Robinson, E. Y.↗

NASA NTRS1973

An engineering model 30 cm ion thruster

Thruster development at Hughes Research Laboratories and NASA Lewis Research Center has brought the 30-cm mercury bombardment ion thruster to the state of an engineering model. This thruster has been designed to have sufficient internal strength for direct mounting on gimbals, to weigh 7.3 kg, to operate with a corrected overall efficiency of 71%, and to have 10,000 hours lifetime. Subassemblies, such as the ion optical system, isolators, etc., have been upgraded to meet launch qualification standards. This paper presents a summary of the design specifications and performance characteristics which define the interface between the thruster module and the remainder of the propulsion system.

Poeschel, R. L.↗

NASA NTRS1974

Isotopic composition of low energy beryllium in the galactic cosmic rays

The multielement charged-particle telescope on the IMP-5 satellite obtained data essentially continuously for more than three years in the interplanetary medium under stable instrument conditions which have made it possible to separate the H, He, and Be isotopes in the galactic cosmic radiation using the technique of double dE/dx vs residual E and particle range measurements. Special emphasis is placed on demonstrating that in the energy range from 50 to 150 MeV/nucleon the instrument clearly resolved Be-7 vs Be-9. A detailed analysis of approximately 100 Be events collected over the mission lifetime yields the ratios (Be-7/Be) = 0.50 plus or minus 0.07, (Be-9/Be) = 0.41 plus or minus 0.10, and (Be-10/Be) = 0.09 plus or minus 0.10. These results are consistent with current models of galactic cosmic ray propagation after taking account of the effects of solar modulation.

Garcia-Munoz, M.↗