Search NASA⌕ Search

SEARCH · Search NASA

Results for “index models”

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 109 records · Page 6

Enhanced light absorption for solid-state brown carbon from wildfires due to organic and water coatings

Abstract Wildfires emit solid-state strongly absorptive brown carbon (solid S-BrC, commonly known as tar ball), critical to Earth’s radiation budget and climate, but their highly variable light absorption properties are typically not accounted for in climate models. Here, we show that from a Pacific Northwest wildfire, over 90% of particles are solid S-BrC with a mean refractive index of 1.49 + 0.056 i at 550 nm. Model sensitivity studies show refractive index variation can cause a ~200% difference in regional absorption aerosol optical depth. We show that ~50% of solid S-BrC particles from this sample uptake water above 97% relative humidity. We hypothesize these results from a hygroscopic organic coating, potentially facilitating solid S-BrC as nuclei for cloud droplets. This water uptake doubles absorption at 550 nm and the organic coating on solid S-BrC can lead to even higher absorption enhancements than water. Incorporating solid S-BrC and water interactions should improve Earth’s radiation budget predictions.

54 ENVIRONMENTAL SCIENCES↗

Impact of Vegetation Assimilation on Flash Drought Characteristics across the Continental United States

Predicting and managing the impacts of flash droughts is difficult owing to their rapid onset and intensification. Flash drought monitoring often relies on assessing changes in root-zone soil moisture. However, the lack of widespread soil moisture measurements means that flash drought assessments often use process-based model data like that from the North American Land Data Assimilation System (NLDAS). Such reliance opens flash drought assessment to model biases, particularly from vegetation processes. Here, we examine the influence of vegetation on NLDAS-simulated flash drought characteristics by comparing two experiments covering 1981–2017: open loop (OL), which uses NLDAS surface meteorological forcing to drive a land surface model using prognostic vegetation, and data assimilation (DA), which instead assimilates near-real-time satellite-derived leaf area index (LAI) into the land surface model. The OL simulation consistently underestimates LAI across the United States, causing relatively high soil moisture values. Both experiments produce similar geographic patterns of flash droughts, but OL produces shorter duration events and regional trends in flash drought occurrence that are sometimes opposite to those in DA. Across the Midwest and Southern United States, flash droughts are 4 weeks (about 70%) longer on average in DA than OL. Moreover, across much of the Great Plains, flash drought occurrence has trended upward according to the DA experiment, opposite to the trend in OL. This sensitivity of flash drought to the representation of vegetation suggests that representing plants with greater fidelity could aid in monitoring flash droughts and improve the prediction of flash drought transitions to more persistent and damaging long-term droughts.

Drought↗

Sex comparisons in physiological and cognitive performance during hypoxic challenge

Within the tactical aviation community, human performance research lags in considering potential psychophysiological differences between male and female aviators due to little inclusion of females during the design and development of aircraft systems. A poor understanding of how male and female aviators differ with respect to human performance results in unknown potential sex differences on aeromedically relevant environmental stressors, perchance leading to suboptimal performance, safety, and health guidelines. For example, previous hypoxia studies have excluded female participants or lacked a sizeable sample to examine sex comparisons. As such, progress toward sensor development and improving hypoxia familiarization training are stunted due to limited knowledge of how individual differences, including sex, may or may not underlie hypoxia symptoms and performance impairment. Investigating sex differences bridges the gap between aerospace medicine and operational health, and addressing hypoxia is one of many facets yet to be studied. In the current study, we retrospectively examined N = 6 hypoxia studies with male-female participant samples (total, N = 189; male, n = 118; female, n = 71). We explored sex as a predictor of physiological response, sensory deficits, the severity of cognitive performance declines, and symptom manifestation via linear and binary logistic regression models. We found that the female sex predicted lower peripheral oxygen saturation and the likelihood of headache reporting in response to hypoxic challenge, yet explained little variance when combined with age and body mass index. The sensory and cognitive performance models did not converge, suggesting high intra-individual variability. Together, sex, age, and body mass index were not the most robust predictors in responses to hypoxic challenge; we cannot infer this for sensory deficits and cognitive performance within an experimentally induced hypoxic environment. The findings have implications for improving hypoxia familiarization training, monitoring sensor development, and emergency response and recovery protocols in case of a hypoxia occurrence suitable for all aircrew. We recommend continuing to elucidate the impact of sex and intrapersonal differences in hypoxia and other aeromedically relevant stressors in tactical aviation.

60 APPLIED LIFE SCIENCES↗

Modeling inter‐reader variability in clinical target volume delineation for soft tissue sarcomas using diffusion model

Abstract Background Accurate delineation of the clinical target volume (CTV) is essential in the radiotherapy treatment of soft tissue sarcomas. However, this process is subject to inter‐reader variability due to the need for clinical assessment of risk and extent of potential microscopic spread. This can lead to inconsistencies in treatment planning, potentially impacting treatment outcomes. Most existing automatic CTV delineation methods do not account for this variability and can only generate a single CTV for each case. Purpose This study aims to develop a deep learning‐based technique to generate multiple CTV contours for each case, simulating the inter‐reader variability in the clinical practice. Methods We employed a publicly available dataset consisting of fluorodeoxyglucose positron emission tomography (FDG‐PET), x‐ray computed tomography (CT), and pre‐contrast T1‐weighted magnetic resonance imaging (MRI) scans from 51 patients with soft tissue sarcoma, along with an independent validation set containing five additional patients. An experienced reader drew a contour of the gross tumor volume (GTV) for each patient based on multi‐modality images. Subsequently, two additional readers, together with the first one, were responsible for contouring three CTVs in total based on the GTV. We developed a diffusion model‐based deep learning method that is capable of generating arbitrary number of different and plausible CTVs to mimic the inter‐reader variability in CTV delineation. The proposed model incorporates a separate encoder to extract features from the GTV masks, leveraging the critical role of GTV information in accurate CTV delineation. Results The proposed diffusion model demonstrated superior performance with the highest Dice Index (0.902 compared to values below 0.881 for state‐of‐the‐art models) and the best generalized energy distance (GED) (0.209 compared to values exceeding 0.221 for state‐of‐the‐art models). It also achieved the second‐highest recall and precision metrics among the compared ambiguous image segmentation models. Results from both datasets exhibited consistent trends, reinforcing the reliability of our findings. Additionally, ablation studies exploring different model structures and input configurations highlighted the significance of incorporating prior GTV information for accurate CTV delineation. Conclusions The proposed diffusion model successfully generates multiple plausible CTV contours for soft tissue sarcomas, effectively capturing inter‐reader variability in CTV delineation.

Dong, Yafei [Yale Biomedical Imaging Institute Yal↗

Geodyn Material Library: Pseudocap models for dry porous tocks

This report describes the second edition of the Pseudocap Strength models for porous rocks implemented in GEODYN material library. The first model was developed in 2007 and calibrated for concrete. Then, the model parameters were calibrated based on triaxial tests reported for limestones and sandstones of various porosities. In these models some key parameters were chosen as functions of the reference porosity,Φ. Since then multiple modifications were implemented in the model, therefore, it has been recalibrated for some common porous materials such as limestones, sandstones, alluvium, tuffs and granite. Two types of models are described in this repot. The first type (called Pseudocap Model or PM) is for rocks from a specific location. Parameters were calibrated for several specific geologic materials. The second type (called Generic Pseudocap Model or GPM) is useful for the sites where only basic information (rock type, porosity) is available. Generic models include built-in correlations between porosities and other mechanical properties observed for certain rock types. The models of both types were validated by comparing not only to quasi-static triaxial tests for these materials but also to shock Hugoniot data and spherical explosion data for some materials. All models were derived in the frame of isotropic plasticity. They are designed to be used in explicit finite element/difference codes. Tangent stiffness tensor is not provided but can be calculated numerically for the model to be used in implicit finite element codes. For an isotropic material the stress can be decomposed into volumetric and deviatoric parts. The volumetric part is modeled using an Equation of state (EOS) which calculates the pressure and the bulk sound speed as functions of the internal specific energy and density. Here a simple, Mie-Gruneisen EOS is presented, but tabulated EOS (LEOS) provided by the library can be used as well. The stress is limited by the yield surface which depends on three invariants of the stress tensor and specific internal energy. In addition, to capture the strain-rate dependence a simple multiplier is used for the yield surface which depends on the equivalent plastic strain rate. The failure surface (the ultimate yield, Y f , defined later) is chosen in the Hoek-Brown form, commonly used in rock mechanics. It includes measurable parameters such as Unconfined Compressive Strength (UCS) as well as scale parameters characterizing the quality of the rock such as GSI (Geologic Strength Index). Thus, even though the model is calibrated for small samples it offers a way to extrapolate the strength to the field scale using geological characterization of the rock mass. The model captures effects of brittle-ductile transition in rocks by introducing a cap multiplier to the yield function. The rate of dilatancy (bulking) is proportional to the slope of the yield surface affected by the cap. Therefore, it takes place only at low confinements when the pressure is less than the brittle-ductile transition pressure, P BD . On the contrary, the porous compaction takes place at pressures higher than P BD . The cap moves as the porosity is compacted or new porosity is generated due to dilatancy. The porous compaction is modeled using an evolution equation which includes deviatoric stress so that the onset of compaction corresponds to the cap surface. The model captures effects of shear-enhanced compaction which is an important for porous rocks. Section 2 describes the modeling framework and Section 3 presents the model calibration procedure. Section 4 compares experimental data for various rocks versus model predictions. The model parameters used for this comparison are given in Appendix. The files with material constants are available with the latest GEODYN material library distribution.

58 GEOSCIENCES↗

Experimental and modeling investigation of binary liquid mixtures

Binary mixtures of liquids may be encountered in industrial or remote sensing scenarios and present challenges to positive identification compared to neat single-component liquids. Our investigation examines whether one can predict the optical properties of the mixture, i.e. its complex index of refraction, by assuming a linear superposition of the real and imaginary components of the index of refraction in proportion to the ratio of each constituent. To investigate this hypothesis various liquid mixtures were created using mass ratios. The mixtures were then characterized as to their complex index of refraction and used in numerical modeling calculations of thin liquid mixture films on surfaces and compared with composed mixtures using linear n and k synthetic mixtures where the n and k components of the complex index of refraction were combined in similar ratios. The comparison of modeling and experimental results is presented with recommendations for further investigation.

infrared (IR) spectroscopy, Liquid mixtures, compl↗

Using synthetic infrared spectra derived from n/k optical constants for standoff detection of chemical deposits

We report results from a recent field experiment to test the validity of using physics-based synthetic infrared spectra to serve as endmembers in a spectral database targeted at chemical deposits. Specifically, the optical constants n and k, (the real and imaginary part of the refractive index) were used to first model infrared reflectance spectra for different thicknesses of chemical layers (e.g. acetaminophen, methylphosphonic acid – MPA, etc.) on various conducting and insulating substrates such as aluminum, wood, and glass. In the experimental portion of the research, thin films of the solid and liquid analytes were deposited onto such substrates to form micron-thick layers of the analytes at different thicknesses: Standoff data from an imaging instrument were then recorded and analyzed to not only identify the different analytes, but also quantify the layer/deposit thickness. To gauge success, the detection results using the synthetic data were compared to the results from hemispherical reflectance (HRF) spectra that were collected for the same sample planchets measured in the field via standoff methods. Preliminary results indicate good agreement between the synthetic reference data as compared to the lab-measured HRF data in terms of their ability to quantitatively reduce longwave infrared data. Specifically, modeled IR spectra for acetaminophen on an aluminum planchet at various thicknesses (1, 2, 5, 10, 15, and 20 µm) were synthesized and compared with standoff field reflectance data as well as HRF laboratory reflectance spectra for two samples: a 5.2 µm- and 12.8 µm-thick layer of acetaminophen on aluminum. Using a first-order approximation, analysis of the field data estimates the thicknesses of the samples to be 2 and 10 µm for the two samples, respectively, while the HRF laboratory data yields thickness estimates of between 5-10 µm and 10 µm, respectively. Both yield reasonable estimates, with the uncertainty most likely due to factors yet to be accounted for in the synthetic spectra such as light scattering.

complex index of refraction, optical constants, hy↗

Validation of the High-Resolution Salish Sea Tidal Hydrodynamic Model

In this study, a tidal hydrodynamic model was developed and validated to simulate tidal currents in Puget Sound, Washington, to support tidal energy resource characterization using the unstructured-grid, Finite Volume Community Ocean Model (FVCOM). The Salish Sea tidal hydrodynamic model was driven by tides along two open boundaries at the entrance of the Strait of Juan de Fuca and north end of Georgia Strait, and river flows from 19 major rivers in the Salish Sea. To simulate the tidal current in Puget Sound, a high-resolution model grid is required to accurately represent the complex coastlines and bathymetry. The spatial resolution of the model grid varies from ~10 m near river boundaries and ~30 m in small tidal channels and estuaries to near 1000 m inside Georgia Strait and at the open boundaries. Model validation was carried out by comparing simulated and observed water levels at 12 tidal stations and currents at 135 Acoustic Doppler Current Profiler stations in the model domain. A set of model performance metrics, including root mean square error, scatter index, bias, and linear correlation coefficient, were used to quantify the model skills in simulating the tidal hydrodynamics in Puget Sound. Error statistics showed an overall good agreement between simulated and observed tidal elevations and currents, which demonstrated that the Puget Sound tidal model can be used to accurately characterize the tidal stream energy resource in Puget Sound.

16 TIDAL AND WAVE POWER↗

An Interpretable Index of Social Vulnerability to Environmental Hazards

Index-based measures of social vulnerability to environmental hazards are commonly modeled from composites of population-level risk factors. These models overlook individual context in communities' experiences of environmental hazards, producing metrics that may hinder spatial decision support for mitigating and responding to hazards. This paper introduces an interpretable, high-resolution model for generating an individual-oriented social vulnerability index (IOSVI) for the United States built on synthetic populations that couples individual and social determinants of vulnerability. The IOSVI combines an individual vulnerability index (IVI) that ranks individuals in an area’s synthetic population based on intersecting risk factors, with a social vulnerability index (SVI) based on the population’s cumulative distribution of IVI scores. Interpretability of the IOSVI procedure is demonstrated through examples of national, metropolitan, and neighborhood (census tract) level spatial variation in index scores and IVI themes, as well as an exploratory analysis examining risk factors affecting a specific sub-population (military veterans) in areas of high social and environmental vulnerability.

Tuccillo, Joe↗

Seasonal grassland productivity forecast for the U.S. Great Plains using Grass–Cast

Every spring, ranchers in the drought-prone U.S. Great Plains face the same difficult challenge —trying to estimate how much forage will be available for livestock to graze during the upcoming summer grazing season. To reduce this uncertainty in predicting forage availability, we developed an innovative new grassland productivity forecast system, named Grass-Cast, to provide science-informed estimates of growing season above ground net primary production (ANPP). Grass-Cast uses over 30 yr of historical data including weather and the satellite-derived normalized vegetation difference index (NDVI)—combined with ecosystem modeling and seasonal precipitation forecasts—to predict if rangelands in individual counties are likely to produce below-normal, near-normal, or above-normal amounts of grass biomass (lbs/ac). Grass-Cast also provides a view of rangeland productivity in the broader region, to assist in larger scale decision-making—such as where forage resources for grazing might be more plentiful if a rancher’s own region is at risk of drought. Grass-Cast is updated approximately every two weeks from April through July. Each Grass-Cast forecast provides three scenarios of ANPP for the upcoming growing season based on different precipitation outlooks. Near real-time 8-d NDVI can be used to supplement Grass-Cast in predicting cumulative growing season NDVI and ANPP starting in mid-April for the Southern Great Plains and mid-May to early June for the Central and Northern Great Plains. Here, we present the scientific basis and methods for Grass-Cast along with the county-level production forecasts from 2017 and 2018 for ten states in the U.S. Great Plains. The correlation between early growing season forecasts and the end-of growing season ANPP estimate is >50% by late May or early June. In a retrospective evaluation, we compared Grass-Cast end-of-growing season ANPP results to an independent dataset and found that the two agreed 69% of the time over a 20-yr period. Although some predictive tools exist for forecasting upcoming growing season conditions, none predict actual productivity for the entire Great Plains. The Grass-Cast system could be adapted to predict grassland ANPP outside of the Great Plains or to predict perennial biofuel grass production.

54 ENVIRONMENTAL SCIENCES↗

Electrocoagulation Combined with Ultrafiltration Membranes as Pretreatment for RO Desalination of Synthetic Cooling Tower Blowdown Water

Electrocoagulation (EC), an electrochemical water treatment process, is commonly used to remove particulate and colloidal matter from water. Here, in this study, we demonstrate that EC, when coupled with membrane filtration, is also capable of removing dissolved species, such as Ca + , Mg 2+ , and SiO x . The removal of such species is important for downstream membrane-based desalination treatment of the water that can suffer from reduced performance due to membrane scaling. Here, we describe how EC can be combined with a low-pressure membrane (LPM) system to offer efficient (and potentially universal) pretreatment for downstream membrane desalination. Synthetic water, simulating cooling tower blowdown (CTBD) with elevated concentrations of hardness and silicates (Ca: 418 ppm, Mg: 63 ppm, SiO 2 : 50 ppm) is treated using EC coupled to ultrafiltration (UF) to remove Mg up to 30 ± 1%, Ca up to 29 ± 1%, and silica up to 99 ± 1%. We evaluated the effectiveness of the EC-UF pretreatment system in reducing downstream RO scaling using thermodynamic modeling to predict the saturation index (SI) at the RO membrane/water interface. An SI value below zero (SI < 0) indicates under-saturated conditions (with respect to a particular mineral) where mineral scaling does not take place, which correlates with improved water recovery. Our findings suggest that the EC-UF pretreatment system was able to increase water recovery by up to 30%, compared to 0% recovery without pretreatment, under optimal conditions (feed solution pH of 7 and an EC charge loading of 1800 C/L). Finally, we conducted an economic analysis showing that implementing an EC-UF system for CTBD water could yield a cost benefit of up to $\$$14.13 per m 3 compared to direct brine disposal.

cooling tower blowdown↗

Single-cell, whole-embryo phenotyping of mammalian developmental disorders

Mouse models are a critical tool for studying human diseases, particularly developmental disorders. However, conventional approaches for phenotyping may fail to detect subtle defects throughout the developing mouse. Here we set out to establish single-cell RNA sequencing of the whole embryo as a scalable platform for the systematic phenotyping of mouse genetic models. We applied combinatorial indexing-based single-cell RNA sequencing to profile 101 embryos of 22 mutant and 4 wild-type genotypes at embryonic day 13.5, altogether profiling more than 1.6 million nuclei. The 22 mutants represent a range of anticipated phenotypic severities, from established multisystem disorders to deletions of individual regulatory regions. We developed and applied several analytical frameworks for detecting differences in composition and/or gene expression across 52 cell types or trajectories. Some mutants exhibit changes in dozens of trajectories whereas others exhibit changes in only a few cell types. We also identify differences between widely used wild-type strains, compare phenotyping of gain- versus loss-of-function mutants and characterize deletions of topological associating domain boundaries. Notably, some changes are shared among mutants, suggesting that developmental pleiotropy might be ‘decomposable’ through further scaling of this approach. Overall, our findings show how single-cell profiling of whole embryos can enable the systematic molecular and cellular phenotypic characterization of mouse mutants with unprecedented breadth and resolution.

60 APPLIED LIFE SCIENCES↗

xsearch

xsearch is a utility to index and search for Coupled model Intercomparison (CMIP) datasets from data stores that are in a standard directory structure.

Po-Chedley, StephenD [Lawrence Livermore National ↗

Evaluation of satellite Leaf Area Index in California vineyards for improving water use estimation

Remote sensing estimation of evapotranspiration (ET) directly quantifies plant water consumption and provides essential information for irrigation scheduling, which is a pressing need for California vineyards as extreme droughts become more frequent. Many ET models take satellite-derived Leaf Area Index (LAI) as a major input, but how uncertainties of LAI estimations propagate to ET and the partitioning between evaporation and transpiration is poorly understood. Here we assessed six satellite-based LAI estimation approaches using Landsat and Sentinel-2 images against ground measurements from four vineyards in California and evaluated ET sensitivity to LAI in the thermal-based two-source energy balance (TSEB) model. We found that radiative transfer modeling-based approaches predicted low to medium LAI well, but they significantly underestimated high LAI in highly clumped vine canopies (RMSE ~ 0.97 to 1.27). Cubist regression models trained with ground LAI measurements from all vineyards achieved high accuracy (RMSE ~ 0.3 to 0.48), but these empirical models did not generalize well between sites. Red edge bands and the related vegetation index (VI) from the Sentinel-2 satellite contain complementary information of LAI to VIs based on near-infrared and red bands. TSEB ET was more sensitive to positive LAI biases than negative ones. Positive LAI errors of 50% resulted in up to 50% changes in ET, while negative biases of 50% in LAI caused less than 10% deviations in ET. However, even when ET changes were minimal, negative LAI errors of 50% led to up to a 40% reduction in modeled transpiration, as soil evaporation and plant transpiration responded to LAI change divergently. These findings call for careful consideration of satellite LAI uncertainties for ET modeling, especially for the partitioning of water loss between vine and soil or cover crop for effective vineyard irrigation management.

54 ENVIRONMENTAL SCIENCES↗

Automated analysis of unlabeled PV data with Solar Data Tools software: Overview and feature updates

Distributed rooftop PV systems: ubiquitous, yet commonly have unlabeled data Difficult or impossible to form a performance index We developed Solar Data Tools (SDT), an open-source Python library for analyzing PV power (and irradiance) time-series data SDT enables analysis of unlabeled PV data—no model, no meteorological data, no performance index required Takes a statistical signal processing approach Data processing steps are largely pre-defined and automatic regardless of system type—from utility tracking systems to multi-pitch rooftop systems

Meyers-Im, Bennet E↗

Phases of $ \mathcal{N} $ = 2 Sachdev-Ye-Kitaev models

We study N = 2 supersymmetric Sachdev-Ye-Kitaev (SYK) models with com- plex fermions at non-zero background charge. Motivated by multi-charge supersymmetric black holes, we propose a new N = 2 SYK model with multiple U(1) symmetries, integer charges, and a non-vanishing supersymmetric index, realizing features not present in known SYK models. In both models, a conformal solution with a super-Schwarzian mode emerges at low temperatures, signaling the appearance of nearly AdS2/BPS physics. However, in contrast to complex SYK, the fermion scaling dimension depends on the background charge in the conformal limit. For a critical charge, we find a high to low entropy phase transition in which the conformal solution ceases to be valid. This transition has a simple interpretation– the fermion scaling dimension violates the unitarity bound. We offer some comments on a holographic interpretation for supersymmetric black holes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Estimating biotic integrity to capture existence value of freshwater ecosystems

The US Environmental Protection Agency (EPA) uses a water quality index (WQI) to estimate benefits of proposed Clean Water Act regulations. The WQI is relevant to human use value, such as recreation, but may not fully capture aspects of nonuse value, such as existence value. Here, we identify an index of biological integrity to supplement the WQI in a forthcoming national stated preference survey that seeks to capture existence value of streams and lakes more accurately within the conterminous United States (CONUS). We used literature and focus group research to evaluate aquatic indices regularly reported by the EPA’s National Aquatic Resource Surveys. We chose an index that quantifies loss in biodiversity as the observed-to-expected (O/E) ratio of taxonomic composition because focus group participants easily understood its meaning and the environmental changes that would result in incremental improvements. However, available datasets of this index do not provide the spatial coverage to account for how conditions near survey respondents affect their willingness to pay for its improvement. Therefore, we modeled and interpolated the values of this index from sampled sites to 1.1 million stream segments and 297,071 lakes across the CONUS to provide the required coverage. The models explained 13 to 36% of the variation in O/E scores and demonstrate how modeling can provide data at the required density for benefits estimation. We close by discussing future work to improve performance of the models and to link biological condition with water quality and habitat models that will allow us to forecast changes resulting from regulatory options.

59 BASIC BIOLOGICAL SCIENCES↗

Annual Mean Arctic Amplification 1970–2020: Observed and Simulated by CMIP6 Climate Models

While the annual mean Arctic Amplification (AA) index varied between two and three during the 1970–2000 period, it reached values exceeding four during the first two decades of the 21st century. The AA did not change in a continuous fashion but rather in two sharp increases around 1986 and 1999. During those steps the mean global surface air temperature trend remained almost constant, while the Arctic trend increased. Although the “best” CMIP6 models reproduce the increasing trend of the AA in 1980s they do not capture the sharply increasing trend of the AA after 1999 including its rapid step-like increase. We propose that the first sharp AA increase around 1986 is due to external forcing, while the second step close to 1999 is due to internal climate variability, which models cannot reproduce in the observed time.

54 ENVIRONMENTAL SCIENCES↗