Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modern grid”

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 523 records · Page 29

High Resolution Nature Runs and the Big Data Challenge

NASA's Global Modeling and Assimilation Office at Goddard Space Flight Center is undertaking a series of very computationally intensive Nature Runs and a downscaled reanalysis. The nature runs use the GEOS-5 as an Atmospheric General Circulation Model (AGCM) while the reanalysis uses the GEOS-5 in Data Assimilation mode. This paper will present computational challenges from three runs, two of which are AGCM and one is downscaled reanalysis using the full DAS. The nature runs will be completed at two surface grid resolutions, 7 and 3 kilometers and 72 vertical levels. The 7 km run spanned 2 years (2005-2006) and produced 4 PB of data while the 3 km run will span one year and generate 4 BP of data. The downscaled reanalysis (MERRA-II Modern-Era Reanalysis for Research and Applications) will cover 15 years and generate 1 PB of data. Our efforts to address the big data challenges of climate science, we are moving toward a notion of Climate Analytics-as-a-Service (CAaaS), a specialization of the concept of business process-as-a-service that is an evolving extension of IaaS, PaaS, and SaaS enabled by cloud computing. In this presentation, we will describe two projects that demonstrate this shift. MERRA Analytic Services (MERRA/AS) is an example of cloud-enabled CAaaS. MERRA/AS enables MapReduce analytics over MERRA reanalysis data collection by bringing together the high-performance computing, scalable data management, and a domain-specific climate data services API. NASA's High-Performance Science Cloud (HPSC) is an example of the type of compute-storage fabric required to support CAaaS. The HPSC comprises a high speed Infinib and network, high performance file systems and object storage, and a virtual system environments specific for data intensive, science applications. These technologies are providing a new tier in the data and analytic services stack that helps connect earthbound, enterprise-level data and computational resources to new customers and new mobility-driven applications and modes of work. In our experience, CAaaS lowers the barriers and risk to organizational change, fosters innovation and experimentation, and provides the agility required to meet our customers' increasing and changing needs

big data analysis↗

Agricultural Breadbaskets Shift Poleward Given Adaptive Farmer Behavior Under Climate Change

Modern food production is spatially concentrated in global “breadbaskets”. A major unresolved question is whether these peak production regions will shift poleward as the climate warms, allowing some recovery of potential climaterelated losses. While agricultural impacts studies to date have focused on currently cultivated land, the Global Gridded Crop Model Intercomparison Project (GGCMI) Phase 2 experiment allows us to assess changes in both yields and the location of peak productivity regions under warming. We examine crop responses under projected end-of-century warming using 7 process-based models simulating 5 major crops (maize, rice, soybeans, and spring and winter wheat) with a variety of adaptation strategies. We find that in no-adaptation cases, when planting date and cultivar choices are held fixed, regions of peak production remain stationary and yield losses can be severe, since growing seasons contract strongly with warming. When adaptations in management practices are allowed (cultivars that retain growing season length under warming and modified planting dates), peak productivity zones shift poleward and yield losses are largely recovered. While most growing-zone shifts are ultimately limited by geography, breadbaskets studied here move poleward over 600 km on average by end of the century under RCP8.5. These results suggest that agricultural impacts assessments can be strongly biased if restricted in spatial area or in the scope of adaptive behavior considered. Accurate evaluation of food security under climate change requires global modeling and careful treatment of adaptation strategies.

adaptation↗

Spatiotemporal mapping of mesoscopic liquid dynamics

The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Therefore, compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Deep Reinforcement Learning for Distribution System Operations: A Tutorial and Survey

Here, the rapid evolution of modern electric power distribution systems into complex networks of interconnected active devices, distributed generation (DG), and storage poses increasing difficulties for system operators. The large-scale integration of distributed energy resources (DERs) and the rapid exchange of measurement data via communication networks present major opportunities for advancing grid operations but also introduce greater uncertainty, higher data dimensionality, more complex network and device models, and challenging control and optimization problems. Deep reinforcement learning (DRL) algorithms are promising in addressing these challenges. However, they have not been effectively adapted for power systems applications, requiring extensive customization for implementation and evaluation. This has resulted in reproducibility challenges and a steep learning curve for researchers new to applying DRL algorithms to the power systems domain. To bridge these gaps, this tutorial aims to serve as a valuable resource for researchers interested in exploring learning-based algorithms to operate active power distribution networks. Specifically, this work presents a generalized process for translating sequential decision-making problems in power distribution systems into Markov decision process (MDP) formulations, illustrated through concrete grid service examples. Additionally, we introduce a simple environment design strategy to develop and evaluate example DRL algorithms for distribution system applications, complete with an included code repository to guide users through environment construction.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Development of an unsteady aerodynamic analysis for finite-deflection subsonic cascades

An unsteady potential flow analysis, which accounts for the effects of blade geometry and steady turning, was developed to predict aerodynamic forces and moments associated with free vibration or flutter phenomena in the fan, compressor, or turbine stages of modern jet engines. Based on the assumption of small amplitude blade motions, the unsteady flow is governed by linear equations with variable coefficients which depend on the underlying steady low. These equations were approximated using difference expressions determined from an implicit least squares development and applicable on arbitrary grids. The resulting linear system of algebraic equations is block tridiagonal, which permits an efficient, direct (i.e., noniterative) solution. The solution procedure was extended to treat blades with rounded or blunt edges at incidence relative to the inlet flow.

Verdon, J. M.↗

Equilibrium gas flow computations. II - An analysis of numerical formulations of conservation laws

Modern numerical techniques employing properties of flux Jacobian matrices are extended to general, equilibrium gas laws. Generalizations of the Beam-Warming scheme, Steger-Warming and van Leer flux-vector splittings, and Roe's approximate Riemann solver are presented for three-dimensional, time-varying grids. The approximations inherent in previous generalizations are discussed.

Vinokur, Marcel↗

Alaska Meteorology, Energy, and Transmission (MET) Toolkit

The Alaska MET (Meteorology, Energy, and Transmission) Toolkit is the National Laboratory of the Rockies' (NLR) new flagship atmospheric dataset, designed to support comprehensive long-term planning and operations across the entire power sector. Serving as the regional counterpart to CONUS-wide HRRR MET Toolkit, this dataset provides a comprehensive, high-fidelity meteorological record covering Alaska.The Alaska MET Toolkit is delivered at an hourly resolution on a standardized 2-km horizontal grid. This dataset is repackaged from the National Oceanic and Atmospheric Administration's (NOAA) operational High-Resolution Rapid Refresh for Alaska (HRRR-AK) forecasts. Spanning from 2019 to 2025, it overcomes the technical barriers of native weather models by providing spatial regridding from the native 3-km HRRR-AK horizontal resolution to a 2-km grid, temporal gap-filling, and vertical interpolation at key energy-relevant heights. By delivering highly accurate, validation-backed data across a comprehensive suite of atmospheric variables - including temperature, pressure, humidity, and wind characteristics - the Alaska MET Toolkit provides a highly accessible and strictly standardized foundation for modern power system modeling.

17 WIND ENERGY↗

Combining model-based and model-free methods for stochastic control of distributed energy resources

Modern distribution systems are experiencing a fast transformation with the growing penetration of distributed energy resources (DERs). Along with the economic and environmental benefits of DERs, challenges arise to address the uncertainties caused by their inherent volatility. If properly coordinated, however, DERs have the potential to provide the controllability that grid operators need. Here, we propose a hierarchical control framework that combines the model-based and model-free methods for stochastic DER control in distribution systems. The upper-level scheduler considers a chance-constrained optimal power flow problem (model-based) that schedules DER setpoints to minimize the operational cost and maintain the operating reserve. The lower-level distributed DER controllers absorb real-time disturbances and uncertainties using the extremum seeking control (model-free) to achieve grid objectives. The combination of model-based and model-free methods allows us to take the advantages of both methods to effectively manage the uncertainty in distribution systems. The proposed work is demonstrated on the IEEE 13-node feeder.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Small-Signal Stability Constrained Optimal Power Flow of Inverter-Dominated Power Systems with Flexible Operation Mode Selection

Given the intermittence and low inertia nature of inverter-based resources (IBRs), modern power systems with high penetration of IBRs challenge the conventional optimal power flow (OPF) analysis and the system may experience unexpected failures if stability constraints are not incorporated. This study proposes a small-signal stability-constrained OPF (SSSC-OPF) with flexible operation mode selection between grid-forming (GFM) and grid-following (GFL) modes for IBRs to address these challenges. The approach aims to maintain system stability with a sufficient stability margin while minimizing operation costs. The effectiveness of the proposed method is validated through extensive case studies on the IEEE 14-bus system. The results demonstrate that the proposed method is able to support system-level power flow analysis, reduce generation costs, and ensure stability under various disturbances.

grid-following↗

Field Measurement and Analysis of Frequency and RoCoF for Low-Inertia Power Systems

The frequency measurement and rate of change of frequency (RoCoF) calculations are important for active power, frequency control, and frequency-related protection in power system operations. With the rapid decrease in inertia under modern power systems, however, multiple power variations introduced by variable renewable energy bring large RoCoF and introduce oscillations in frequency measurements. To shed light on this issue, this article proposes a frequency measurement and RoCoF analysis technology, as well as the sensitivity analysis for the real low-inertia power grid. First, a real-world event in an island power grid is analyzed as an example, and an event numerical model is estimated considering four uncertainties, including event RoCoF, oscillation magnitude, oscillation frequency, and underfrequency load-shedding threshold. Then, a sensitivity analysis is performed using the Sobol indices method and quasi-Monte Carlo simulations to identify the dominant uncertainties based on the event model. Here, the dominant uncertainties are compared to provide a guideline for the frequency estimation algorithm selection under different event cases. This algorithm is compared with an enhanced zero-crossing algorithm and gives a 3.41% better RoCoF error rate. Finally, an experiment is conducted using universal grid analyzers to validate the simulation results and the theoretical analysis.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Downscaling atmospheric chemistry simulations with physically consistent deep learning

Abstract. Recent advances in deep convolutional neural network (CNN)-based super resolution can be used to downscale atmospheric chemistry simulations with substantially higher accuracy than conventional downscaling methods. This work both demonstrates the downscaling capabilities of modern CNN-based single image super resolution and video super-resolution schemes and develops modifications to these schemes to ensure they are appropriate for use with physical science data. The CNN-based video super-resolution schemes in particular incur only 39 % to 54 % of the grid-cell-level error of interpolation schemes and generate outputs with extremely realistic small-scale variability based on multiple perceptual quality metrics while performing a large (8×10) increase in resolution in the spatial dimensions. Methods are introduced to strictly enforce physical conservation laws within CNNs, perform large and asymmetric resolution changes between common model grid resolutions, account for non-uniform grid-cell areas, super-resolve lognormally distributed datasets, and leverage additional inputs such as high-resolution climatologies and model state variables. High-resolution chemistry simulations are critical for modeling regional air quality and for understanding future climate, and CNN-based downscaling has the potential to generate these high-resolution simulations and ensembles at a fraction of the computational cost.

54 ENVIRONMENTAL SCIENCES↗

Experimental Investigation of Separated and Transitional Boundary Layers Under Low-Pressure Turbine Airfoil Conditions

Modern low-pressure turbine airfoils are subject to increasingly stronger pressure gradients as designers impose higher loading in an effort to improve efficiency and to reduce part count. The adverse pressure gradients on the suction side of these airfoils can lead to boundary-layer separation, particularly under cruise conditions. Separation bubbles, notably those which fail to reattach, can result in a significant degradation of engine efficiency. Accurate prediction of separation and reattachment is hence crucial to improved turbine design. This requires an improved understanding of the transition flow physics. Transition may begin before or after separation, depending on the Reynolds number and other flow conditions, has a strong influence on subsequent reattachment, and may even eliminate separation. Further complicating the problem are the high free-stream turbulence levels in a real engine environment, the strong pressure gradients along the airfoils, the curvature of the airfoils, and the unsteadiness associated with wake passing from upstream stages. Because of the complicated flow situation, transition in these devices can take many paths that can coexist, vary in importance, and possibly also interact, at different locations and instances in time. The present work was carried out in an attempt to systematically sort out some of these issues. Detailed velocity measurements were made along a flat plate subject to the same nominal dimensionless pressure gradient as the suction side of a modern low-pressure turbine airfoil ('Pak-B'). The Reynolds number based on wetted plate length and nominal exit velocity, Re, was varied from 50;000 to 300; 000, covering cruise to takeoff conditions. Low, 0.2%, and high, 7%, inlet free-stream turbulence intensities were set using passive grids. These turbulence levels correspond to about 0.2% and 2.5% turbulence intensity in the test section when normalized with the exit velocity. The Reynolds number and free-stream turbulence level do not have a significant effect on the location of boundary-layer separation unless they are high enough to induce transition upstream of separation. The location and extent of the transition zone, in contrast, depend strongly on Re and TI. The beginning of reattachment closely follows the onset of transition. Under low free-stream turbulence conditions the boundary layer is laminar at separation and then begins to exhibit fluctuations in a finite frequency band in the shear layer over the separation bubble. These fluctuations are due to instability waves. The fluctuations grow in magnitude, higher harmonics are generated, and finally lead to a breakdown to turbulence. Transition begins in the shear layer, but quickly spreads to the near wall region and causes the boundary layer to reattach. The transition is rapid and the resulting turbulence contains a full range of high and low frequencies. Under high free-stream turbulence conditions, slowly growing low-frequency fluctuations are induced in the pretransitional boundary layer by the free-stream. The separation bubbles are considerably thinner than in the low TI cases, resulting in thinner boundary layers at the end of the test wall. At Re=50,000 and 100,000, the pre-transitional boundary layer separates at about the same location as in the low TI cases. Transition occurs through a bypass mode, begins upstream of the corresponding low-TI location, and proceeds in a manner similar to that of an attached boundary layer. Under high TI at Re=200,000 and 300,000, transition begins before separation. The boundary layer may separate, but if it does the separation bubble is very short and does not significantly affect the downstream development of the boundary layer. A comparison is made to previous work in a simulated cascade.

Hultgren, Lennart S.↗

A critical review of cyber-physical security for building automation systems

Modern Building Automation Systems (BASs), as the brain that enable the smartness of a smart building, often require increased connectivity both among system components as well as with outside entities, such as the cloud, to enable low-cost remote management, optimized automation via outsourced cloud analytics, and increased building-grid integrations. As smart buildings move towards open communication technologies, providing access to BASs through the building's intranet, or even remotely through the Internet, has become a common practice. However, increased connectivity and accessibility come with increased cyber security threats. BASs were historically developed as closed environments with limited cyber-security considerations. As a result, BASs in many buildings are vulnerable to cyber-attacks that may cause adverse consequences, such as occupant discomfort, excessive energy usage, and unexpected equipment downtime. Therefore, there is a strong need to advance the state-of-the-art in cyber-physical security for BASs and provide practical solutions for attack mitigation in buildings. However, an inclusive and systematic review of BAS vulnerabilities, potential cyber-attacks with impact assessment, detection & defense approaches, and cyber resilient control strategies is currently lacking in the literature. This review paper fills the gap by providing a comprehensive up-to-date review of cyber-physical security for BASs at three levels in commercial buildings: management level, automation level, and field level. The general BASs vulnerabilities and protocol-specific vulnerabilities for the four dominant BAS protocols (i.e., BACnet, KNX, LonWorks, and Modbus) are reviewed, followed by a discussion on four attack targets and seven potential attack scenarios. Furthermore, the impact of cyber-attacks on BASs is summarized as signal corruption, signal delaying, and signal blocking. The typical cyber-attack detection and defense approaches are identified at the three levels. Cyber resilient control strategies for BASs under attack are categorized into passive and active resilient control schemes. Open challenges and future opportunities are finally discussed.

97 MATHEMATICS AND COMPUTING↗

The Global Modeling and Assimilation Office (GMAO) 4d-Var and its Adjoint-based Tools

The fifth generation of the Goddard Earth Observing System (GEOS-5) Data Assimilation System (DAS) is a 3d-var system that uses the Grid-point Statistical Interpolation (GSI) system developed in collaboration with NCEP, and a general circulation model developed at Goddard, that includes the finite-volume hydrodynamics of GEOS-4 wrapped in the Earth System Modeling Framework and physical packages tuned to provide a reliable hydrological cycle for the integration of the Modern Era Retrospective-analysis for Research and Applications (MERRA). This MERRA system is essentially complete and the next generation GEOS is under intense development. A prototype next generation system is now complete and has been producing preliminary results. This prototype system replaces the GSI-based Incremental Analysis Update procedure with a GSI-based 4d-var which uses the adjoint of the finite-volume hydrodynamics of GEOS-4 together with a vertical diffusing scheme for simplified physics. As part of this development we have kept the GEOS-5 IAU procedure as an option and have added the capability to experiment with a First Guess at the Appropriate Time (FGAT) procedure, thus allowing for at least three modes of running the data assimilation experiments. The prototype system is a large extension of GEOS-5 as it also includes various adjoint-based tools, namely, a forecast sensitivity tool, a singular vector tool, and an observation impact tool, that combines the model sensitivity tool with a GSI-based adjoint tool. These features bring the global data assimilation effort at Goddard up to date with technologies used in data assimilation systems at major meteorological centers elsewhere. Various aspects of the next generation GEOS will be discussed during the presentation at the Workshop, and preliminary results will illustrate the discussion.

Todling, Ricardo↗

Tools And Methods to Analyze Plant Outage Schedule and Assist Schedulers in Improving Outage Resilience

Refueling outages of nuclear power plants (NPPs) are considered one of the most critical phases throughout the plant lifetime. In such instances, tens of thousands of activities (e.g., maintenance, surveillance) are performed in a short amount of time (typically 2-3 weeks unless major backfitting or modernization projects are carried out) by a large number of crews (e.g., electricians, mechanics) that are hired as contractors. As a consequence, a plant outage can be expensive not only in terms of costs (e.g., contractor labor, material), but also in terms of loss generation since the plant is taken off the grid during the full outage duration (an indicative metric is about 1.2M$/day of loss of revenue). Thus, there is a continuous need to decrease the economic impact of outages on plant finances. This can be done by: decreasing the frequency of plant outages (e.g., from 18 to 24 months), reducing the time to complete the outage, and reducing the risk of outage delays. The Optimization of Outage Activities project under the Risk Informed Systems Analysis Pathway (RISA) sponsored by Department of Energy (DOE) Light Water Reactor Sustainability (LWRS) Program focuses on developing tools and methods to support NPPs with outage schedule optimization. The developed tools and methods are designed to analyze plant outage schedule with the goal of identify critical elements in the schedule that might pose a high risk of delays. These methods and tools can be considered resource-centric in the sense that they address outage challenges as a resource optimization problem. In this context, resources are either time and crews; outage delays occurs when either (or both) resources are insufficient to complete the set of tasks assigned at a specific time instant of the outage. This report provides details on how plant resources (time and crews) can be allocated in such a way that delays are minimized. In this respect, two classes of methods have been developed: the first one focuses on the time resource and how variability of the time to complete outage tasks may impact outage delays. The second one integrates available resources to assess when dailies activities should be performed such that the risk of outage delays are minimized.

97 MATHEMATICS AND COMPUTING↗

The E3SM Diagnostics Package (E3SM Diags v2.7): a Python-based diagnostics package for Earth system model evaluation

Abstract. The E3SM Diagnostics Package (E3SM Diags) is a modern, Python-based Earth system model (ESM) evaluation tool (with Python module name e3sm_diags), developed to support the Department of Energy (DOE) Energy Exascale Earth System Model (E3SM). E3SM Diags provides a wide suite of tools for evaluating native E3SM output, as well as ESM data on regular latitude–longitude grids, including output from Coupled Model Intercomparison Project (CMIP) class models. E3SM Diags is modeled after the National Center for Atmospheric Research (NCAR) Atmosphere Model Working Group (AMWG, 2022) diagnostics package. In its version 1 release, E3SM Diags included a set of core essential diagnostics to evaluate the mean physical climate from model simulations. As of version 2.7, more process-oriented and phenomenon-based evaluation diagnostics have been implemented, such as analysis of the quasi-biennial oscillation (QBO), the El Niño–Southern Oscillation (ENSO), streamflow, the diurnal cycle of precipitation, tropical cyclones, ozone and aerosol properties. An in situ dataset from DOE's Atmospheric Radiation Measurement (ARM) program has been integrated into the package for evaluating the representation of simulated cloud and precipitation processes. This tool is designed with enough flexibility to allow for the addition of new observational datasets and new diagnostic algorithms. Additional features include customizable figures; streamlined installation, configuration and execution; and multiprocessing for fast computation. The package uses an up-to-date observational data repository maintained by its developers, where recent datasets are added to the repository as they become available. Finally, several applications for the E3SM Diags module were introduced to fit a diverse set of use cases from the scientific community.

54 ENVIRONMENTAL SCIENCES↗

CFD Sensitivity Analysis of a Modern Civil Transport Near Buffet-Onset Conditions

A computational fluid dynamics (CFD) sensitivity analysis is conducted for a modern civil transport at several conditions ranging from mostly attached flow to flow with substantial separation. Two different Navier-Stokes computer codes and four different turbulence models are utilized, and results are compared both to wind tunnel data at flight Reynolds number and flight data. In-depth CFD sensitivities to grid, code, spatial differencing method, aeroelastic shape, and turbulence model are described for conditions near buffet onset (a condition at which significant separation exists). In summary, given a grid of sufficient density for a given aeroelastic wing shape, the combined approximate error band in CFD at conditions near buffet onset due to code, spatial differencing method, and turbulence model is: 6% in lift, 7% in drag, and 16% in moment. The biggest two contributers to this uncertainty are turbulence model and code. Computed results agree well with wind tunnel surface pressure measurements both for an overspeed 'cruise' case as well as a case with small trailing edge separation. At and beyond buffet onset, computed results agree well over the inner half of the wing, but shock location is predicted too far aft at some of the outboard stations. Lift, drag, and moment curves are predicted in good agreement with experimental results from the wind tunnel.

Rumsey, Christopher L.↗

Increasing Prevalence of Hot Drought Across Western North America Since the 16th Century

Across western North America (WNA), 20th-21st century anthropogenic warming has increased the prevalence and severity of concurrent drought and heat events, also termed hot droughts. However, the lack of independent spatial reconstructions of both soil moisture and temperature limits the potential to identify these events in the past and to place them in a long-term context. We develop the Western North American Temperature Atlas (WNATA), a data-independent 0.5° gridded reconstruction of summer maximum temperatures back to the 16th century. Our evaluation of the WNATA with existing hydroclimate reconstructions reveals an increasing association between maximum temperature and drought severity in recent decades, relative to the past five centuries. The synthesis of these paleo-reconstructions indicates that the amplification of the modern WNA megadrought by increased temperatures and the frequency and spatial extent of compound hot and dry conditions in the 21st century are likely unprecedented since at least the 16th century.

Drought↗