SEARCH · Engineering Papers

Results for “data gap analysis”

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

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗

DOE OSTI2022

REGAL International Program: Analysis of experimental data for depletion code validation

The Rod-Extremity and Gadolinia AnaLysis (REGAL) Program is a joint international effort to expand the nuclide inventory experimental data for irradiated nuclear fuel, with a specific focus on addressing two challenging needs associated with the characterization of modern, high duty, nuclear fuel. The first challenge is filling the gaps in experimental nuclide inventory data for gadolinia (UO 2 –Gd 2 O 3 ) fuel rods. The huge absorption cross sections of Gd-155 and Gd-157 in the Gd dopant in these rods lead to atypical spatial self-shielding patterns and have an impact on the neutronic environment within the fuel assembly compared to regular UO 2 fuel rods. Additionally, the second challenge is investigating the impact of burnup gradients at rod extremities on fuel composition and neutron leakage, to provide relevant experimental data for assessing computational capabilities to model such impact. A benchmark has been defined as a first step in the development of best-estimate models in the preliminary phase of the experimental data evaluation. Comparison of experimental results obtained in Phase I of the program for two measured pressurized water reactor (PWR) samples, one UO 2 and one UO 2 –Gd 2 O 3 sample, with calculated results obtained with different computational tools based on the defined benchmark are presented and discussed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Chemical control of competing electron transfer pathways in iron tetracyano-polypyridyl photosensitizers

Photoinduced intramolecular electron transfer dynamics following metal-to-ligand charge-transfer (MLCT) excitation of [Fe(CN) 4 (2,2'-bipyridine)] 2- (1), [Fe(CN) 4 (2,3-bis(2-pyridyl)pyrazine)] 2- (2) and [Fe(CN) 4 (2,2'-bipyrimidine)] 2- (3) were investigated in various solvents with static and time-resolved UV-Visible absorption spectroscopy and Fe 2p3d resonant inelastic X-ray scattering (RIXS). This series of polypyridyl ligands, combined with the strong solvatochromism of the complexes, enables the 1 MLCT vertical energy to be varied from 1.64 eV to 2.64 eV and the 3 MLCT lifetime to range from 180 fs to 67 ps. The 3 MLCT lifetimes in 1 and 2 decrease exponentially as the MLCT energy increases, consistent with electron transfer to the lowest energy triplet metal-centred ( 3 MC) excited state, as established by the Tanabe–Sugano analysis of the Fe 2p3d RIXS data. In contrast, the 3 MLCT lifetime in 3 changes non-monotonically with MLCT energy, exhibiting a maximum. This qualitatively distinct behaviour results from a competing 3 MLCT → ground state (GS) electron transfer pathway that exhibits energy gap law behaviour. The 3 MLCT → GS pathway involves nuclear tunnelling for the high-frequency polypyridyl breathing mode (hν = 1530 cm -1 ), which is most displaced for complex 3, making this pathway significantly more efficient. Our study demonstrates that the excited state relaxation mechanism of Fe polypyridyl photosensitizers can be readily tuned by ligand and solvent environment. Furthermore, our study reveals that extending charge transfer lifetimes requires control of the relative energies of the 3 MLCT and the 3 MC states and suppression of the intramolecular distortion of the acceptor ligand in the 3 MLCT excited state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Reliability Assessment of Cooling Fans for PV Inverters: Testing, Modeling, and Case Studies

The reliability of photovoltaic (PV) inverters is critical for long-term solar system performance, with cooling fan failures frequently leading to costly downtime. While much research exists on general cooling fan reliability, little attention has been given to fans operating within PV inverters and their unique environmental challenges. Here, this article proposes a comprehensive methodology to address this gap. First, a failure mode and effects analysis is performed on fans to identify the key failure mechanisms in PV applications, their corresponding stressors, and the models necessary for lifetime prediction. Second, an accelerated life test is designed and conducted to collect valuable experimental data for PV inverter fans in a reasonable amount of time. Third, a mathematical conversion of dynamic mission profiles into effective constant stress levels is derived. Fourth, case studies are given, showcasing lifetime estimates that account for geographic variations in mission profile data. The results demonstrate that this integrated approach leads to an accurate reliability assessment for PV inverter cooling fans.

accelerated life testing (ALT)↗

DOE OSTI2023

Probing rutile solid-phase crystallization of atomically mixed Mn-alloyed TiO 2 coatings through XANES analysis

Phase transformation via dopant introduction is an emerging strategy for tuning the electronic and catalytic properties of oxides. However, developing a general theory of ternary alloys and their microstructure evolution remains challenging due to the knowledge gaps in how point defects such as dopants lead to the formation of complex or metastable phases. We investigate the influence of Mn concentration and annealing temperature on the microstructure evolution of amorphous Mn-alloyed TiO 2 synthesized via Atomic Layer Deposition. In conclusion, our findings reveal a new crystallization pathway resulting in a direct amorphous TiO 2 to rutile transformation without an intermediate anatase phase and our multi-parameter thermal processing conditions show potential for AI-assisted data acquisition and analysis for studying complex oxide coatings.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Systems Analysis Approach to Polyethylene Terephthalate and Olefin Plastics Supply Chains in the Circular Economy: A Review of Data Sets and Models

The environmental and economic impacts of implementing a circular economy in plastic waste supply chains are not well understood. The proposed systems analysis framework assesses environmental, social, and economic impacts of plastic waste supply chains in a circular economy. The first objective of this article is to identify datasets, models and knowledge gaps associated with waste plastic supply chain processes, mainly in the U.S. Our literature review indicated that the best datasets exist for virgin plastic resin production, mechanical recycling, landfilling, and incineration, with materials recovery facility being intermediate, and with chemical recycling the lowest. The second objective of this article is to develop an illustrative application of the framework by conducting a preliminary systems analysis of PET bottles with closed-loop recycling. Here, the preliminary systems analysis of PET bottles utilized a linear programming optimization method. Our optimization model indicated that both chemical and mechanical recycling processes are needed to achieve a true circular economy of PET bottles with the least greenhouse gas emissions, specifically reductions of 24% when compared with the linear economy. Good quality and standardized life cycle assessment and techno-economic analysis studies are needed to better understand the environmental, economic, and social impacts of advanced sorting and chemical recycling technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Polyolefin plastic waste hydroconversion to fuels, lubricants, and waxes: a comparative study

Hydroconversion technologies have surged to the forefront of deconstructing plastic waste. Recent studies have been performed over several catalysts with varying conditions and plastics that make comparisons difficult. We compile and compare data from the literature by introducing various metrics and perform a simple energy analysis. We draw mechanistic similarities to and differences from the past literature on small alkane hydroconversion and leverage the former to propose standard approaches to tune product selectivity. We exemplify the plastics materials gap and the challenges it creates. Lastly, we discuss the current limitations and suggest future work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

6⁢𝐻-perovskite dimer lattice with antiferromagnetic interactions: Ba 3⁢ Zn 1−𝑥⁢ Ca 𝑥 ⁢Ru 2 ⁢O 9

Here, we investigate the magnetic behavior of the 6⁢𝐻-perovskite dimer lattice Ba 3 ⁢Zn 1−𝑥 ⁢Ca 𝑥 ⁢Ru 2 ⁢O 9 using analytical theory, density functional theory, inelastic neutron scattering, and modeling of historical magnetization and neutron-scattering data. A dimer mean-field theory built upon classical Luttinger–Tisza analysis generates a phase diagram revealing a transition from a nonmagnetic singlet to a finite-moment ground state as interdimer couplings increase. A (generalized) linear spin-wave theory captures multiplet mixing, excitation gap closing, and fluctuation-induced moment suppression. Density-functional-theory calculations on selected compounds, together with neutron spectroscopy of dilute Ba 3 ⁢Zn⁢(Ru 1−𝑥 ⁢Sb 𝑥 ) 2 ⁢O 9 , confirm the exchange hierarchy, enabling quantification of previously published experiments within this framework. Our results identify three mechanisms for magnetic moment suppression—quantum fluctuations, ligand hybridization, and nonmagnetic-singlet/magnetic-multiplet mixing.

Pajerowski, Daniel M. [Oak Ridge National Laborato↗

DOE OSTI2025

A Modelica Implementation of an Organic Rankine Cycle

Organic Rankine cycle (ORC) systems generate power from low-grade heat sources, such as geothermal sources and industrial waste heat. A key feature is that a working fluid is selected to match the temperature of the source. With the vast pool of candidate working fluids comes the challenge of developing a large number of robust thermodynamic media models. We implemented a subcritical ORC model in Modelica that uses working fluid data records and interpolation schemes in lieu of thermodynamic medium evaluation for energy recovery estimation. This is a component model that can be integrated into a larger energy system model. It does not require detailed thermodynamic, heat transfer, or machine analysis. Our ORC model fills a gap where working fluids are ready to choose or easy to add, and at the same time can be integrated into an energy system.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2026

Phase change material-to-refrigerant heat exchangers: Experimental validation and uneven melting analysis

Thermal energy storage (TES) using phase change materials (PCMs) enables load shifting and reduces mismatches between the building's thermal demand and the heat pump (HP) system's thermal capacity. While most TES studies have focused on single-phase heat transfer fluids, PCM-to-refrigerant heat exchangers (PRHX), in which the refrigerant undergoes liquid-vapor phase change, remain largely underexplored. Here, to address this gap, this study systematically investigates a shell-and-tube PRHX operating as the condenser in an HP through experiments and simulations. A two-dimensional enthalpy-based finite-volume PRHX model was developed and validated against experimental data, including refrigerant and PCM temperature profiles, PRHX capacity, and system-level performance. Results show three characteristic stages: Stage I, the condenser temperature rapidly increases with PCM in the solid state absorbing sensible heat; Stage II, stable operation during PCM phase change; and Stage III, performance decline once part of the PCM becomes fully melted. Heat transfer analysis revealed that uneven PCM melting along the condenser length was driven primarily by variations in the refrigerant-side heat transfer coefficient, rather than local approach temperature differences. Further study demonstrated that minimizing refrigerant outlet subcooling led to more uniform PCM melting and extended Stages I and II by up to 89%. This paper provides a validated PRHX model, clarifies the mechanisms of uneven PCM melting, and highlights subcooling control as an effective strategy to improve TES performance in HP systems without requiring secondary loops.

Graphite matrix↗

DOE OSTI2021

A simplified integrated framework for predicting the economic impacts of feedstock variations in a catalytic fast pyrolysis conversion process

Feedstock attributes of lignocellulosic biomass, such as particle size, compositional makeup, and moisture content, can vary substantially even within pre-processed materials and have a significant effect on conversion in fast pyrolysis-based processes. However, the economic impacts of these attributes are not well understood. To address this, biomass deconstruction phenomena captured with a versatile particle-scale simulation were linked to techno-economic impacts via reduced-order models. Parametric analysis of the particle-scale model, which was validated using literature data, was used in combination with multiple linear regression models to develop correlations between feedstock attributes and yields of pyrolysis oil, gas, and char. Yields were then correlated with the minimum fuel selling price (MFSP) using a techno-economic model, bridging the gap between physics-based biomass conversion simulations and predictions of MFSP for a catalytic fast-pyrolysis process. Empirical correlations derived from the literature regarding the impact of mineral matter (ash) on oil yield were also considered. The model correlations deployed in the integrated framework capture the impacts of variation in feedstock attributes on the MFSP. Variations in ash were shown to have the biggest impact, varying MFSP by -13%/+22% due to catalytic effects and lower relative amounts of convertible lignocellulosic material. It was also found that, if ash can be controlled to low levels, the increased extractives in forest residues can help compensate for some yield losses associated with increased ash. As a result, other inputs considered (particle size, moisture content, and reactor temperature) had relatively negligible effects on process economics within the ranges analyzed considering particle-scale effects alone.

BIOMASS FUELS↗

DOE OSTI2023

Analysis of Fire-Induced Circulations during the FireFlux2 Experiment

Despite recent advances in both coupled fire modeling and measurement techniques to sample the fire environment, the fire–atmosphere coupling mechanisms that lead to fast propagating wildfires remain poorly understood. This knowledge gap adversely affects fire management when wildland fires propagate unexpectedly rapidly and shift direction due to the fire impacts on local wind conditions. In this work, we utilized observational data from the FireFlux2 prescribed burn and numerical simulations performed with a coupled fire–atmosphere model WRF-SFIRE to assess the small-scale impacts of fire on local micrometeorology under moderate wind conditions (10–12 m/s). The FireFlux2 prescribed burn provided a comprehensive observational dataset with in situ meteorological observations as well as IR measurements of fire progression. To directly quantify the effects of fire–atmosphere interactions, two WRF-SFIRE simulations were executed. One simulation was run in a two-way coupled mode in which the heat and moisture fluxes emitted from the fire were injected into the atmosphere, and the other simulation was performed in a one-way coupled mode for which the atmosphere was not affected by the fire. The difference between these two simulations was used to analyze and quantify the fire impacts on the atmospheric circulation at different sections of the fire front. The fire-released heat fluxes resulted in vertical velocities as high as 10.8 m/s at the highest measurement level (20 m above ground level) gradually diminishing with height and dropping to 7.9 m/s at 5.77 m. The fire-induced horizontal winds indicated the strongest fire-induced flow at the lowest measurement levels (as high as 3.3 m/s) gradually decreasing to less than 1 m/s at 20 m above ground level. The analysis of the simulated flow indicates significant differences between the fire-induced circulation at the fire head and on the flanks. The fire-induced circulation was much stronger near the fire head than at the flanks, where the fire did not produce particularly strong cross-fire flow and did not significantly change the lateral fire progression. However, at the head of the fire the fire-induced winds blowing across the front were the strongest and significantly accelerated fire progression. The two-way coupled simulation including the fire-induced winds produced 36.2% faster fire propagation than the one-way coupled run, and more realistically represented the fire progression.

Benik, Jeremy T.↗

DOE OSTI2021

Techno-Economic Analysis of decentralized preprocessing systems for fast pyrolysis biorefineries with blended feedstocks in the southeastern United States

This study evaluated the economic feasibility of fast pyrolysis biorefineries fed with blended pine residues and switchgrass in the Southeastern U.S. with different supply chain design. Previous techno-economic analyses (TEA) have focused on either blended biomass or decentralized preprocessing without investigating the impacts of varied process parameters, technology options, and real-world biomass distribution. This study fills the literature gap by modeling scenarios for different biomass blending ratios, biorefinery and preprocessing site (so-called depot) capacities, and alternative preprocessing technologies. High-resolution, real-world geospatial data were analyzed using Geographic Information Systems to facilitate supply chain design and TEA. For a decentralized system, the minimum fuel selling price (MFSP) of biofuel was 3.92–4.33 per gallon gasoline equivalent (GGE), while the MFSP for the centralized biorefinery at the same capacities ranged between 3.75–4.02/GGE. Implementing a high moisture pelleting process depot rather than a conventional pelleting process lowered the MFSP by 0.03–0.17/GGE. Scenario analysis indicated decreased MFSP with increasing biorefinery capacities but not necessarily with increasing depot size. Medium-size depots (500 OMDT/day) achieved the lowest MFSP. Here, this analysis identified the optimal blending ratios for two preprocessing technologies at varied depot sizes. Counterintuitively, increasing the proportion of higher cost switchgrass reduced the MFSP for large biorefineries (>5000 ODMT/day), but increased the MFSP for small biorefineries (1000–2500 ODMT/day). Although the decentralized systems have a higher MFSP based on current analysis, it has other potential benefits such as mitigated supply chain risks and improved feedstock quality that are difficult to be quantified in this TEA.

09 BIOMASS FUELS↗

DOE OSTI2022

Passive and active peer effects in the spatial diffusion of residential solar panels: A case study of the Las Vegas Valley

This research analyzes the role of peer influences on the adoption of residential rooftop solar photovoltaic panels (PV) within the context of the Diffusion of Innovation Theory. PV literature indicates that adopters are influenced by word of mouth (WOM) information exchange with peers, i.e., active peer effects, while other studies suggest that living near households with visible rooftop PV installations influences adoption, i.e., passive peer effects. We bridge the gap in this literature by conducting a mixed methods analysis. We develop and administer a survey to Las Vegas Valley (LVV) residents to identify current and potential PV adopters' perceptions of peer-effects and consumer intention variables. We conduct a spatial analysis of Google's Project Sunroof data to identify LVV neighborhoods in the later stages of the PV diffusion process, i.e., those with the highest PV adoption rates. Key results show that current PV adopters living in early diffusion areas report significantly higher active and passive peer effects compared to adopters in later diffusion areas. Potential adopters in later diffusion areas report higher passive peer effects than those in early diffusion areas. Overall, because LVV has a low PV adoption rate (<3%), short term strategies aimed at increasing PV adoption should emphasize WOM active peer effects. Here, we caution against long-term green marketing strategies focusing solely on peer-effects as the PV market matures.

14 SOLAR ENERGY↗

DOE OSTI2025

Portable Distortion Free Large Solid Angle Coverage Detector

Neutron single crystal diffractometers require large solid angle coverage for optimum performance. This can be achieved by tiling flat detectors in a cylindrical or spherical geometry, but results in large gaps in detector coverage and parallax distortion. The detector edges exhibit degraded resolution, distortion, and gamma rejection. Since the detector edge regions are a significant fraction of the detector active area, they must be removed from the experimental data set, requiring extra beam time to collect enough analysis data. A spherical detector with a continuous surface would effectively address this issue while eliminating most boundary ‘dead’ areas. Here we report on the development of a novel hemispherical shaped neutron detector using seamlessly tiled readout modules to form the desired shape. The heart of the detector is a specially developed curved neutron scintillator coupled to high resolution silicon photomultiplier (SiPM) Anger cameras via custom made fiber optic tapers (FOTs). The detector has been assembled and initial tests have been conducted at the High Flux Isotope Reactor (HFIR) beamlines at Oak Ridge National Laboratory (ORNL). Here, in this work, we describe details of the scintillator design, fabrication and characterization, evaluation of individual detector modules, the details of the detector design implementation, and evaluation of the assembled detector at ORNL beamlines.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

Tidal Resource Gaps Analysis Technical Report

The tidal resource gaps project was created to address a growing body of evidence that models underpredict tidal current speeds compared to measurements at a number of the top-ranking tidal energy sites. In response, this project compared opportunistic tidal power density measurements from 16 tidal energy hot spot sites with estimates from resource assessment data to identify discrepancies. To improve the accuracy of resource estimates from model data, updated data from eight improved model simulations were obtained. Model improvements included grid refinement, domain coupling, and the use of unstructured or nested grids. New resource estimates were calculated from the updated models, and these data were used to update the tidal hot spots-a list of promising tidal energy sites around the United States.

16 TIDAL AND WAVE POWER↗

DOE OSTI2020

A Multiyear Gridded Data Ensemble of Surface Biogenic Carbon Fluxes for North America: Evaluation and Analysis of Results

Accurate and fine-scale estimates of biogenic carbon fluxes are critical for measuring and monitoring the biosphere's responses and feedback to the climate system. Currently available data products from flux towers and model-intercomparison projects struggle to adequately represent spatiotemporal dynamics of surface biogenic carbon fluxes, and to quantify their uncertainties, which also are crucial to atmospheric inversion systems. To address these gaps, we introduce a new perturbed-parameter model ensemble with the CASA model to estimate surface biogenic carbon fluxes at monthly and 3-hourly scales for North America at ~500-m and 5-km resolutions. We first use the Extended Fourier Amplitude Sensitivity Testing to choose the three most sensitive parameters to be perturbed, maximum light use efficiency (E max ), optimal temperature of photosynthesis (T opt ), and temperature response of respiration (Q 10 ). The initial range for each parameter is broadly sampled for the L1 ensemble, but then we pruned E max with site-level primary productivity to derive an L2 ensemble with narrower uncertainty ranges. Ensembles are strongly correlated with site-level results at both monthly and 3-hourly scales, and the spread across L1/L2 ensemble members encompasses the range of AmeriFlux observations. Monthly variability in the L2 ensemble mean is 85% of the observed variability. The L2 ensemble outperforms diverse data products with the highest Taylor skill scores at diurnal to annual scales. Finally, the ensemble's seasonality agrees well with other models for most biome types and in high and middle latitudes, but inconsistencies are found in subtropical and tropical ecoregions and for annual totals over North America.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Floods and Heavy Precipitation at the Global Scale: 100‐Year Analysis and 180‐Year Reconstruction

Floods and heavy precipitation have disruptive impacts worldwide, but their historical variability remains only partially understood at the global scale. This article aims at reducing this knowledge gap by jointly analyzing seasonal maxima of streamflow and precipitation at more than 3,000 stations over a 100‐year period. The analysis is based on Hidden Climate Indices (HCIs). Like standard climate indices (e.g., Nino 3.4, NAO), HCIs are used as covariates explaining the temporal variability of data, but unlike them, HCIs are estimated from the data. In this work, a distinction is made between common HCIs, that affect both heavy precipitation and floods, and specific HCIs, that exclusively affect one or the other. Overall, HCIs do not show noticeable autocorrelation, but some are affected by noticeable trends. In particular, strong and wide‐ranging trends are identified in precipitation‐specific HCIs, while trends affecting flood‐specific HCIs are weaker and have more localized effects. A probabilistic model is then derived to link HCIs and large‐scale atmospheric variables (pressure, wind, temperature) and to reconstruct HCIs since 1836 using the 20CRv3 reanalysis. In turn this allows estimating the probability of occurrence of floods and heavy precipitation at the global scale. This 180‐year reconstruction highlights flood hot‐spots and hot‐moments in the distant past, well before the establishment of perennial monitoring networks. Finally, the approach presented in this study is generic and paves the way for an improved characterization of historical variability by making a better use of long but highly irregular station data sets.

54 ENVIRONMENTAL SCIENCES↗