Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation 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.

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

At least 1,153 records · Page 64

Sea surface warming and ocean-to-atmosphere feedback driven by large-scale offshore wind farms under seasonally stratified conditions

Offshore wind farms may induce changes in the upper ocean and near-surface atmosphere through coupled ocean-atmosphere feedbacks. Yet, the role of air-sea interactions mediated by offshore wind farms remains poorly understood. Using fully coupled ocean-atmosphere-wave model simulations for seasonally stratified conditions along the US East Coast, we show that simulated cumulative reductions in wind stress due to large-scale wind farm clusters lead to sea surface warming of 0.3° to 0.4°C and a shallower mixed layer. This warming drives upward heat fluxes, destabilizing the atmospheric boundary layer and enhancing wind stress, which partially offsets wake-induced wind deficits. These wake-ocean interactions influence near-surface meteorology and air-sea fluxes, suggesting that a coupled modeling approach may be necessary for assessing potential oceanographic impacts of offshore wind developments. However, ocean coupling exerts limited influence on winds at turbine-relevant heights or within downstream wakes, resulting in minimal impact on long-term energy. These findings suggest that models without ocean coupling may be adequate for wind energy applications.

17 WIND ENERGY↗

Properties and Controlling Processes of Aerosol and Cloud Condensation Nuclei in Marine Boundary Layer Over Eastern North Atlantic (Final Report)

Remote marine low cloud systems have large spatial and temporal coverages. Because of their relatively low optical thickness and background aerosol concentrations, marine low clouds are particularly susceptible to perturbations in aerosols associated with anthropogenic emissions. Recent studies show that a large fraction of the global aerosol indirect forcing can be attributed to changes in marine low clouds, despite their relatively long distance from most anthropogenic sources. Currently, there remain large uncertainties in the magnitude of the global aerosol indirect radiative forcing. One major contribution to this large uncertainty derives from poor understanding of the aerosols under natural conditions, the perturbation by anthropogenic emissions, and the processes that drive them. In this project, we took advantage of the comprehensive, multi-dimensional characterization of meteorological parameters, trace gases, aerosol, cloud, and drizzle fields during the Aerosol and Cloud Experiments in Eastern North Atlantic (ACE-ENA) and the long-term measurements at the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic (ENA) site on Graciosa Island. By synergistically combining in situ measurements, remote sensing retrievals, and model simulations, we quantified aerosol properties in the ENA, their vertical profiles, and variations with season, synoptic state, and airmass. We investigated key aerosol processes, including entrainment of free troposphere aerosol into marine boundary layer, new particle formation, growth of Aitken mode particles, and the influences of these processes on the cloud condensation nuclei (CCN) population in marine boundary layer. The representation of key aerosol processes in the Energy Exascale Earth System Model (E3SM) was evaluated using the observations. This project has led to 26 publications, which are listed at the end of this report.

54 ENVIRONMENTAL SCIENCES↗

Robustness of energy performance of Zero-Net-Energy (ZNE) homes

Zero-net-energy (ZNE) homes produce an adequate amount of energy on-site to meet their energy demand based on source energy for an entire year. California building energy efficiency standards require new residential buildings started in 2020 to be ZNE. For various reasons, a home designed as ZNE may not achieve ZNE performance in real operation. Here we aimed to quantify the robustness of the energy performance of ZNE homes due to weather variability, climate change, and the uncertainty of occupant behavior. A single-family ZNE house, based on the optimal cost-effective design in three California climate zones, was used to develop the EnergyPlus simulation models. Weather variations were considered from a combination of the historical 30 years’ actual meteorological year (AMY) weather data, typical year weather data in TMY3, and future weather data based on Intergovernmental Panel on Climate Change scenarios. Three scenarios of occupant behavior from the energy perspective were defined to represent the uncertainty about occupants’ activities, comfort requirements, and their adaptive interactions with buildings and systems. In terms of annual source energy, the simulation results of the ZNE homes showed: (1) a decrease of 23–38 percent for occupants with energy austerity behavior and an increase of 120–130 percent for occupants with energy wasteful behavior, compared with the baseline assumption of normal occupants; (2) a variation range of –15 percent to +14 percent for the results using 30-year AMY weather data compared with the baseline results using TMY3 weather data; (3) an increase of 10–13 percent with future weather in Fresno and Riverside and a decrease of 15 percent with San Francisco; and (4) climate change can reduce the gap between the austerity and wasteful consumption. These findings provide insights into how ZNE homes may perform in reality and inform architects, engineers, occupants, and policymakers to pay more attention to occupant behavior on design, operation, and regulations of ZNE homes to ensure energy performance robustness.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

More frequent atmospheric rivers slow the seasonal recovery of Arctic sea ice

In recent decades, Arctic sea ice coverage experienced a drastic decline in winter, when sea ice is expected to recover following the melting season. It is unclear to what extent atmospheric processes such as atmospheric rivers (ARs), intense corridors of moisture transport, contribute to this reduced recovery of sea ice. Here, in this study, using observations and climate model simulations we find a robust frequency increase in ARs in early winter over the Barents-Kara Seas and the central Arctic (ABK) between 1979 and 2021. The moisture carried by more frequent ARs has intensified surface downward longwave radiation and rainfall, caused stronger melting of thin, fragile ice cover, and slowed the seasonal recovery of sea ice, accounting for 34% of the sea ice cover decline in ABK. A series of model ensemble experiments suggests that, in addition to a uniform AR increase in response to anthropogenic warming, tropical Pacific variability also contributes to the observed Arctic AR changes.

54 ENVIRONMENTAL SCIENCES↗

Reducing Sample Size While Improving Equity in Vaccine Clinical Trials: A Machine Learning-Based Recruitment Methodology with Application to Improving Trials of Hepatitis C Virus Vaccines in People Who Inject Drugs

Despite the availability of direct-acting antivirals that cure individuals infected with the hepatitis C virus (HCV), developing a vaccine is critically needed in achieving HCV elimination. HCV vaccine trials have been performed in populations with high incidence of new HCV infection such as people who inject drugs (PWID). Developing strategies of optimal recruitment of PWID for HCV vaccine trials could reduce sample size, follow-up costs and disparities in enrollment. We investigate trial recruitment informed by machine learning and evaluate a strategy for HCV vaccine trials termed PREDICTEE—Predictive Recruitment and Enrichment method balancing Demographics and Incidence for Clinical Trial Equity and Efficiency. PREDICTEE utilizes a survival analysis model applied to trial candidates, considering their demographic and injection characteristics to predict the candidate’s probability of HCV infection during the trial. The decision to recruit considers both the candidate’s predicted incidence and demographic characteristics such as age, sex, and race. We evaluated PREDICTEE using in silico methods, in which we first generated a synthetic candidate pool and their respective HCV infection events using HepCEP, a validated agent-based simulation model of HCV transmission among PWID in metropolitan Chicago. We then compared PREDICTEE to conventional recruitment of high-risk PWID who share drugs or injection equipment in terms of sample size and recruitment equity, with the latter measured by participation-to-prevalence ratio (PPR) across age, sex, and race. Comparing conventional recruitment to PREDICTEE found a reduction in sample size from 802 (95%: 642–1010) to 278 (95%: 264–294) with PREDICTEE, while also reducing screening requirements by 30%. Simultaneously, PPR increased from 0.475 (95%: 0.356–0.568) to 0.754 (95%: 0.685–0.834). Even when targeting a dissimilar maximally balanced population in which achieving recruitment equity would be more difficult, PREDICTEE is able to reduce sample size from 802 (95%: 642–1010) to 304 (95%: 288–322) while improving PPR to 0.807 (95%: 0.792–0.821). PREDICTEE presents a promising strategy for HCV clinical trial recruitment, achieving sample size reduction while improving recruitment equity.

60 APPLIED LIFE SCIENCES↗

Design, Optimization, and In-Depth Understanding of Front and Rear Junction Screen-Printed Double-Side Passivated Contacts Solar Cells

In this work, detailed numerical modeling is performed for front junction (FJ) and rear junction (RJ) n-type Si solar cells with screen-printed double-side poly-Si based tunnel oxide passivated contacts (TOPCon). Fundamental understanding of both structures is crucial to explore factors that limit the efficiency potential of the two designs. By investigating several key parameters such as front poly sheet resistance and thickness, bulk material properties, and carrier transport in our simulation model, we determine and explain why RJ cells outperform FJ cells. Our findings reveal that FJ suffers from present technological limitations of p-poly based passivated contacts, namely, i) large recombination observed in textured p-TOPCon emitter and ii) low solid solubility and hole mobility in boron-doped poly-Si resulting in very high sheet resistance of front p-poly layers which must be made thin to reduce absorption losses. In contrast, RJ desensitized the cell efficiency to front sheet resistance allowing application of ultra-thin front n-poly layers and is therefore ideally suited for double-side TOPCon cells. Lastly, we have developed a roadmap for efficiency projections of commercial-type screen-printed RJ and FJ topologies reaching ~24% and >22.5% efficiencies respectively, quantifying and explaining the impact of technological innovations on the performance of each design.

14 SOLAR ENERGY↗

Hydrodynamic description of D meson production in high-energy heavy-ion collisions *

The large values and constituent-quark-number scaling of the elliptic flow of low-p T D mesons imply that charm quarks, initially produced through hard processes, might be partially thermalized through strong interactions with quark-gluon plasma (QGP) in high-energy heavy-ion collisions. To quantify the degree of thermalization of low-p T charm quarks, we compare the D 0 meson spectra and elliptic flow from a hydrodynamic model to experimental data as well as transport model simulations. We use an effective charm chemical potential at the freeze-out temperature to account for the initial charm quark production from hard processes and assume that they are thermalized in the local comoving frame of the medium before freeze-out. D 0 mesons are sampled statistically from the freeze-out hyper-surface of the expanding QGP as described by the event-by-event (3+1)D viscous hydrodynamic model CLVisc. Both the hydrodynamic and transport models can describe the elliptic flow of D 0 mesons at P T < 3 GeV/c as measured in Au+Au collisions at $sqrts_{NN} = 200$ GeV. Though the experimental data on spectra are consistent with the hydrodynamic result at small p T ~ 1 GeV/c, they deviate from the hydrodynamic model at high transverse momentum, GeV/c. The diffusion and parton energy loss mechanisms in the transport model can describe the measured spectra reasonably well within the theoretical uncertainty. Finally, our comparative study indicates that charm quarks only approach local thermal equilibrium at small p T , even though they acquire sizable elliptic flow comparable to light-quark hadrons at both small and intermediate p T .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamic P-Q Capability and Abnormal Operation Analysis of a Wind Turbine with Doubly-Fed Induction Generator

After several accidents occurred in wind farms in the United States and around the world, the subsynchronous oscillation (SSO) issues in grid-connected wind farms have gained serious attention. Particularly, these issues have caused significant challenges to wind turbines with a doubly fed induction generator (DFIG) because it is connected to the grid via both its stator and rotor paths. Traditionally, a P-Q capability chart is utilized to assure the safe operation boundary for a synchronous generator. But the energy conversion characteristics of a DFIG wind turbine are completely different. A critical factor to affect the reliable operation of a DFIG wind turbine is the rated current and pulse-width modulation (PWM) saturation constraints of its power converters. These constraint conditions can be affected by the wind turbine rotating speed and grid conditions. However, a detailed study of DFIG P-Q capability from these perspectives has not been conducted, which has hindered adequate understanding of many abnormal wind turbine operations reported in the literature and the development of advanced control technologies to overcome the challenges. The proposed study in this paper considers vector control implementation to DFIG power electronic converters in the dq reference frame, and the models and algorithms developed for the P-Q capability study have addressed specific DFIG power converter constraints that are different from those of a traditional synchronous generator. The paper especially focuses on exploring the dynamic natures of DFIG P-Q capability under uncertain and variable conditions to explore the root causes of many abnormal operations of DFIG wind turbines reported in the literature. The proposed study is validated through an electromagnetic transient simulation model of a grid-connected DFIG wind turbine. The proposed study has the potential to lead to the development of new DFIG control technologies that can help overcome the challenges for many abnormal operations of DFIG wind turbines.

17 WIND ENERGY↗

Separating the Influences of Low-Latitude Warming and Sea Ice Loss on Northern Hemisphere Climate Change

Abstract Analyzing a multimodel ensemble of coupled climate model simulations forced with Arctic sea ice loss using a two-parameter pattern-scaling technique to remove the cross-coupling between low- and high-latitude responses, the sensitivity to high-latitude sea ice loss is isolated and contrasted to the sensitivity to low-latitude warming. Despite some differences in experimental design, the Northern Hemisphere near-surface atmospheric sensitivity to sea ice loss is found to be robust across models in the cold season; however, a larger intermodel spread is found at the surface in boreal summer, and in the free tropospheric circulation. In contrast, the sensitivity to low-latitude warming is most robust in the free troposphere and in the warm season, with more intermodel spread in the surface ocean and surface heat flux over the Northern Hemisphere. The robust signals associated with sea ice loss include upward turbulent and longwave heat fluxes where sea ice is lost, warming and freshening of the Arctic Ocean, warming of the eastern North Pacific Ocean relative to the western North Pacific with upward turbulent heat fluxes in the Kuroshio Extension, and salinification of the shallow shelf seas of the Arctic Ocean alongside freshening in the subpolar North Atlantic Ocean. In contrast, the robust signals associated with low-latitude warming include intensified ocean warming and upward latent heat fluxes near the western boundary currents, freshening of the Pacific Ocean, salinification of the North Atlantic, and downward sensible and longwave fluxes over the ocean.

54 ENVIRONMENTAL SCIENCES↗

Local limits of detection for anthropogenic aerosol-cloud interactions

Ship tracks are quasi-linear cloud patterns produced from the interaction of ship emissions with low boundary layer clouds. They are visible throughout the diurnal cycle in satellite images from space-borne assets like the Advanced Baseline Imagers (ABI) aboard the National Oceanic and Atmospheric Administration Geostationary Operational Environmental Satellites (GOES-R). However, complex atmospheric dynamics often make it difficult to identify and characterize the formation and evolution of tracks. Ship tracks have the potential to increase a cloud's albedo and reduce the impact of global warming. Thus, it is important to study these patterns to better understand the complex atmospheric interactions between aerosols and clouds to improve our climate models, and examine the efficacy of climate interventions, such as marine cloud brightening. Over the course of this 3-year project, we have developed novel data-driven techniques that advance our ability to assess the effects of ship emissions on marine environments and the risks of future marine cloud brightening efforts. The three main innovative technical contributions we will document here are a method to track aerosol injections using optical flow, a stochastic simulation model for track formations and an automated detection algorithm for efficient identification of ship tracks in large datasets.

54 ENVIRONMENTAL SCIENCES↗

On the relationship between precipitation extreme and local temperature over eastern China based on convection permitting simulations: roles of different moisture processes and precipitation types

The Clausius–Clapeyron (CC) scaling, which indicates a roughly 7% increase in saturated water vapor per 1 °C increase in temperature, can serve as a strong constraint linking the intensity of precipitation extremes and local temperature. However, the relationship between precipitation extreme and local temperature (referred to as the PE-T relationship) does not always follow the CC scaling and is highly dependent on climate regimes. In this study, we investigated the impacts of different moisture processes and precipitation types on the PE-T relationship over eastern China during the summertime based on convection-permitting model simulations. Consistent with observations, the simulated intensity of precipitation extremes increases with temperature at a rate close to CC (double-CC) scaling below (above) 20 °C. When the temperature exceeds 25 °C, precipitation intensity starts to drop. Precipitation extremes are mainly contributed by the stratiform, MCS (i.e., mesoscale convective system) convective, and non-MCS convective precipitation at low (< 20 °C), medium (20–25 °C), and high (> 25 °C) temperatures, respectively, suggesting that the double-CC scaling occurs when convective types become dominant, while the negative scaling at high temperatures is attributed to the reduced horizontal scale of convection. Corresponding to the reduced intensity of precipitation at high temperatures, there are stronger divergence and subsidence in the low-level atmosphere, which is probably caused by the net cooling associated with the enhanced melting and evaporation of falling hydrometeors due to the lower relative humidity in the low-level atmosphere. Overall, our findings contribute to a deeper understanding of the temperature dependence of precipitation extremes in eastern China.

54 ENVIRONMENTAL SCIENCES↗

Stochastic reconstruction of fracture network pattern using spatial point processes

Fracture spatial patterns can strongly affect fluid flow in the subsurface. Proximity and distribution of fractures control reservoir flow behavior over various length scales. In many studies, however, simplified geometrical patterns are generated for fractures which may lead to unrealistic subsurface models. Here we introduce a new method for characterization and modeling of fracture spatial patterns based on outcrop observations. We use Ripley's K-function to characterize the arrangement of fracture barycenters and intersection points over various length scales. In addition, we use semivariograms to quantify spatial correlation in fracture intensity maps. Using this information, we develop a stochastic algorithm that generates two-dimensional fracture network realizations with spatial properties similar to those of a real fracture network measured in the field. Numerical simulation models indicate that the generated fracture realizations exhibit similar flow behaviors as that of the original fracture network. Furthermore, such modeling tools expand and improve our capability in building representative fracture models and in quantification of uncertainty in naturally and hydraulically fractured reservoirs.

58 GEOSCIENCES↗

Computing the Paschen curve for argon with speed-limited particle-in-cell simulation

Upon inclusion of collisions, the speed-limited particle-in-cell (SLPIC) simulation method successfully computed the Paschen curve for argon. The 1D3V simulations modeled an electron cascade across an argon-filled capacitor, including electron-neutral ionization, electron-neutral elastic collisions, electron-neutral excitation, and ion-induced secondary electron emission. In electrical breakdown, the timescale difference between ion and electron motion makes traditional PIC methods computationally slow. To decrease this timescale difference and speed up computation, here we used SLPIC, a time-domain algorithm that limits the speed of the fastest electrons in the simulation. The SLPIC algorithm facilitates a straightforward, fully kinetic treatment of dynamics and collisions. SLPIC was as accurate as PIC, but ran up to 200 times faster. SLPIC accurately computed the Paschen curve for argon over three orders of magnitude in pressure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Arctic Mixed-Phase Cloud Base Ice Precipitation Properties Over the NSA Site

Cloud-climate feedbacks are still the greatest source of uncertainty in current climate projections. Arctic clouds, which are predominantly stratiform and supercooled, often long-lived, and nearly continuously precipitate ice particles, contribute roughly 10% of the uncertainty attributed to the global cloud feedback. This arctic cloud uncertainty is driven by incomplete observational and theoretical knowledge required to estimate and explain the state and active processes occurring in those clouds. A focus on ice precipitation properties at arctic cloud base rather than the surface deconfounds the product of cloud condensate sink processes from the influence of the atmospheric thermodynamic state below cloud base, rendering cloud-base properties a more appealing target for inference and evaluation of model simulations. This data set provides more than 1800 samples of cloud-base ice precipitation properties over Utqiagvik, North Slope of Alaska, all of which were retrieved using the synthesis of ARM radar and lidar measurements. The retrieved ice precipitation variables in this data set include, among others, the ice number concentration, water content, PSD parameters, precipitation rate, mass-weighted fall velocity, vertical air motion, and effective radius, all of which are highly valuable for model evaluation and a general understanding of polar cloud sink processes. Each variable sample includes its mean value and associated uncertainty. Additional variables based on ARM measurements (liquid layer statistics, etc.) are included in this data set. The retrieval algorithm and analysis of this data set are described in Silber (JGR, 2023, https://doi.org/10.1029/2022JD038202).

54 ENVIRONMENTAL SCIENCES↗

Polarimetric Radar Characteristics of Simulated and Observed Intense Convective Cores for a Midlatitude Continental and Tropical Maritime Environment

This study contrasts midlatitude continental and tropical maritime deep convective cores using polarimetric radar observables and retrievals from selected deep convection episodes during the MC3E and TWPICE field campaigns. The continental convective cores produce stronger radar reflectivities throughout the profiles, while maritime convective cores produce more positive differential reflectivity Z dr and larger specific differential phase K dp above the melting level. Hydrometeor identification retrievals revealed the presence of large fractions of rimed ice particles (snow aggregates) in the continental (maritime) convective cores, consistent with the Z dr and K dp observations. The regional cloud-resolving model simulations with bulk and size-resolved bin microphysics are conducted for the selected cases, and the simulation outputs are converted into polarimetric radar signals and retrievals identical to the observational composites. Both the bulk and the bin microphysics reproduce realistic land and ocean (L-O) contrasts in reflectivity, polarimetric variables of rain drops, and hydrometeor profiles, but there are still large uncertainties in describing Z dr and K dp of ice crystals associated with the ice particle shapes/orientation assumptions. Sensitivity experiments are conducted by swapping background aerosols between the continental and maritime environments, revealing that background aerosols play a role in shaping the distinct L-O contrasts in radar reflectivity associated with raindrop sizes, in addition to the dominant role of background thermodynamics.

54 ENVIRONMENTAL SCIENCES↗

Fluctuations and power-law scaling of dry, frictionless granular rheology near the hard-particle limit

Here the flow of frictionless granular particles is studied with stress-controlled discrete element modeling simulations for systems varying in size from 300 to 100,000 particles. The volume fraction and shear-stress ratio μ are relatively insensitive to system size for a wide range of inertial numbers I. Second-order effects in strain rate, such as normal stress differences, require large system sizes to accurately extract meaningful results, notably a nonmonotonic dependence in the first normal stress difference with strain rate. The rheological response represented by the μ(I) scalar model works well at describing the lower-order aspects of the rheology, except near the quasistatic limit of these stress-controlled flows. The pressure is varied over five decades, and a pressure dependence of the coordination number is observed, which is not captured by the inertial number. Large fluctuations observed for small systems N≤1000 near the quasistatic limit can lead to the arrest of flow resulting in challenges to fitting the data to rheological relationships. The inertial number is also insufficient for capturing the pressure-dependent behavior of property fluctuations. Fluctuations in the flow and microstructural properties are measured in both the quasistatic and inertial regimes, including shear stress, pressure, strain rate, normal stress differences, volume fraction, coordination number, and contact fabric anisotropy. The fluctuations in flow properties scale self-similarly with pressure and system size. A transition in the scaling of fluctuations of stress properties and contact fabric anisotropy are measured and proposed as a quantitative identification of the transition from inertial to quasistatic flow.

36 MATERIALS SCIENCE↗

Multi-timescale operations of nuclear-renewable hybrid energy systems for reserve and thermal product provision

In this paper, an optimal operation strategy of a nuclear-renewable hybrid energy system (N-R HES), in conjunction with a district heating network, is developed within a comprehensive multi-timescale electricity market framework. The grid-connected N-R HES is simulated to explore the capabilities and benefits of N-R HES of providing energy products, different reserve products, and thermal products. An N-R HES optimization and control strategy is formulated to exploit the benefits from the hybrid energy system in terms of both energy and ancillary services. A case study is performed on the customized NREL-118 bus test system with high renewable penetrations, based on a multi-timescale (i.e., three-cycle) production cost model. Both day-ahead and real-time market clearing prices are determined from the market model simulation. In conclusion, the results show that the N-R HES can contribute to the reserve requirements and also meet the thermal load, thereby increasing the economic efficiency of N-R HES (with increased revenue ranging from 1.55% to 35.25% at certain cases) compared to the baseline case where reserve and thermal power exports are not optimized.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Gridded Sub-daily Climate Forcings for North America Based on Daymet and GSWP3 (Daymet-GSWP3)

To support high spatial and temporal resolution land surface modeling, this dataset provides 3-hourly time step historic weather forcing at 1-km spatial resolution for the entire North America. The latest Daymet V4 data provides gridded historic daily weather observations at 1-km spatial resolution from 1980 to 2014. Using sub-daily temporal information from the Global Soil Wetness Project Phase 3 (GSWP3), Daymet was further temporally downscaled to 3-hourly time steps and provided in the format required for land surface model simulations. The process of temporal downscaling preserves the relative magnitude in each sub-daily time step from GSWP3 while maintaining the total and average values from Daymet for each day. This results in a blended 1980-2014 Daymet-GSWP3 dataset. Available variables include surface air temperature, precipitation, specific humidity, shortwave and longwave radiation, wind speed, and pressure. These data can be used as a high-resolution meteorological forcing dataset to support high-resolution land surface modeling where accurate meteorological forcing datasets built from historic observations and/or reanalysis datasets are desirable.

54 ENVIRONMENTAL SCIENCES↗