Variable-resolution Simulations Research Articles

Fast EM-driven nature-inspired optimization of antenna input characteristics using response features and variable-resolution simulation models

Utilization of optimization technique is a must in the design of contemporary antenna systems. Often, global search methods are necessary, which are associated with high computational costs when conducted at the level of full-wave electromagnetic (EM) models. In this study, we introduce an innovative method for globally optimizing reflection responses of multi-band antennas. Our approach uses surrogates constructed based on response features, smoothing the objective function landscape processed by the algorithm. We begin with initial parameter space screening and surrogate model construction using coarse-discretization EM analysis. Subsequently, the surrogate evolves iteratively into a co-kriging model, refining itself using accumulated high-fidelity EM simulation results, with the infill criterion focusing on minimizing the predicted objective function. Employing a particle swarm optimizer (PSO) as the underlying search routine, extensive verification case studies showcase the efficiency and superiority of our procedure over benchmarks. The average optimization cost translates to just around ninety high-fidelity EM antenna analyses, showcasing excellent solution repeatability. Leveraging variable-resolution simulations achieves up to a seventy percent speedup compared to the single-fidelity algorithm.

The 1996 Mid-Atlantic Winter Flood: Exploring Climate Risk through a Storyline Approach

Abstract This article explores the application of thermodynamic perturbations to a historical midlatitude, wintertime, rain-on-snow flood event to evaluate how similar events may evolve under different climate forcings. In particular, we generate a hindcast of the 1996 Mid-Atlantic flood using an ensemble of 14-km variable-resolution simulations completed with the U.S. Department of Energy’s global Energy Exascale Earth System Model (E3SM). We show that the event is skillfully reproduced over the Susquehanna River Basin (SRB) by E3SM when benchmarked against in situ observational data and high-resolution reanalyses. In addition, we perform five counterfactual experiments to simulate the flood under preindustrial conditions and four different levels of warming as projected by the Community Earth System Model Large Ensemble. We find a nonlinear response in simulated surface runoff and streamflow as a function of atmospheric warming. This is attributed to changing contributions of liquid water input from a shallower initial snowpack (decreased snowmelt), increased surface temperatures and rainfall rates, and increased soil water storage. Flooding associated with this event peaks from around +1 to +2 K of global average surface warming and decreases with additional warming beyond this. There are noticeable timing shifts in peak runoff and streamflow associated with changes in the flashiness of the event. This work highlights the utility of using storyline approaches for communicating climate risk and demonstrates the potential nonlinearities associated with hydrologic extremes in areas that experience ephemeral snowpack, such as the SRB.

Technical descriptions of the experimental dynamical downscaling simulations over North America by the CAM–MPAS variable-resolution model

Abstract. Comprehensive assessment of climate datasets is important for communicating model projections and associated uncertainties to stakeholders. Uncertainties can arise not only from assumptions and biases within the model but also from external factors such as computational constraint and data processing. To understand sources of uncertainties in global variable-resolution (VR) dynamical downscaling, we produced a regional climate dataset using the Model for Prediction Across Scales (MPAS; dynamical core version 4.0) coupled to the Community Atmosphere Model (CAM; version 5.4), which we refer to as CAM–MPAS hereafter. This document provides technical details of the model configuration, simulations, computational requirements, post-processing, and data archive of the experimental CAM–MPAS downscaling data. The CAM–MPAS model is configured with VR meshes featuring higher resolutions over North America as well as quasi-uniform-resolution meshes across the globe. The dataset includes multiple uniform- (240 and 120 km) and variable-resolution (50–200, 25–100, and 12–46 km) simulations for both the present-day (1990–2010) and future (2080–2100) periods, closely following the protocol of the North American Coordinated Regional Climate Downscaling Experiment. A deviation from the protocol is the pseudo-warming experiment for the future period, using the ocean boundary conditions produced by adding the sea surface temperature and sea-ice changes from the low-resolution version of the Max Planck Institute Earth System Model (MPI-ESM-LR) in the Coupled Model Intercomparison Project Phase 5 to the present-day ocean state from a reanalysis product. Some unique aspects of global VR models are evaluated to provide background knowledge to data users and to explore good practices for modelers who use VR models for regional downscaling. In the coarse-resolution domain, strong resolution sensitivity of the hydrological cycles exists over the tropics but does not appear to affect the midlatitude circulations in the Northern Hemisphere, including the downscaling target of North America. The pseudo-warming experiment leads to similar responses of large-scale circulations to the imposed radiative and boundary forcings in the CAM–MPAS and MPI-ESM-LR models, but their climatological states in the historical period differ over various regions, including North America. Such differences are carried to the future period, suggesting the importance of the base state climatology. Within the refined domain, precipitation statistics improve with higher resolutions, and such statistical inference is verified to be negligibly influenced by horizontal remapping during post-processing. Limited (≈50 % slower) throughput of the current code is found on a recent many-core/wide-vector high-performance computing system, which limits the lengths of the 12–46 km simulations and indirectly affects sampling uncertainty. Our experience shows that global and technical aspects of the VR downscaling framework require further investigations to reduce uncertainties for regional climate projection.

Rapid yield optimization of miniaturized microwave passives by response features and variable-fidelity EM simulations

The operation of high-frequency devices, including microwave passive components, can be impaired by fabrication tolerances but also incomplete knowledge concerning operating conditions (temperature, input power levels) and material parameters (e.g., substrate permittivity). Although the accuracy of manufacturing processes is always limited, the effects of parameter deviations can be accounted for in advance at the design phase through optimization of suitably selected statistical performance figures. Perhaps the most popular one is the yield, which provides a straightforward assessment of the likelihood of fulfilling performance conditions imposed upon the system given the assumed deviations of designable parameters. The latter are typically quantified by means of probability distributions pertinent to the fabrication process. The fundamental obstacle of the yield-driven design is its high computational cost. The primary mitigation approach nowadays is the employment of surrogate modeling methods. Yet, a construction of reliable metamodels becomes problematic for systems featuring a large number of degrees of freedom. Our work proposes a technique for fast yield optimization of microwave passives, which relies on response feature technology as well as variable-fidelity simulation models. Utilization of response features enables efficient handling of issues related to the system response nonlinearities. Meanwhile, the incorporation of variable-resolution simulations allows for accelerating the yield estimation process, which translates into remarkably low overall cost of the optimizing the yield. Our approach is verified with the use of three microstrip couplers. Comprehensive benchmarking demonstrates its superiority in terms of computational efficiency over the state-of-the-art algorithms, whereas reliability is corroborated by electromagnetic-driven Monte Carlo simulations.

Expedited Gradient-Based Design Closure of Antennas Using Variable-Resolution Simulations and Sparse Sensitivity Updates

Numerical optimization has been playing an increasingly important role in the design of contemporary antenna systems. Due to the shortage of design-ready theoretical models, optimization is mainly based on electromagnetic (EM) analysis, which tends to be costly. Numerous techniques have evolved to abate this cost, including surrogate-assisted frameworks for global optimization, or sparse sensitivity updates for speeding up local search. In the latter, CPU-heavy updates of the system response sensitivity through finite differentiation are suppressed based on, e.g., the magnitude of design variability during the optimization run. Another approach is to incorporate variable-resolution simulations. Recently, a technique exploiting a continuous spectrum of admissible model fidelity levels has been reported, thereby allowing for a considerable reduction of computational expenditures. Seeking further savings, this work introduces an accelerated gradient-based algorithm with sparse sensitivity updates and variable-resolution EM simulations. Our technique is validated using four broadband antennas, and demonstrated to offer substantial (around 80%) savings over the benchmark while maintaining acceptable design quality.

Impact of Dust-Cloud-Radiation-Precipitation Dynamical Feedback on Subseasonal-to-Seasonal Variability of the Asian Summer Monsoon in Global Variable-Resolution Simulations With MPAS-CAM5

In this study, we investigate the effects of increased dust emission from the Middle East deserts on subseasonal-to-seasonal (S2S) variability of the Asian summer monsoon (ASM). Numerical experiments are performed using the Model for Prediction Across Scales (MPAS) coupled with the Community Atmosphere Model (CAM5) physics, with regional refinement at 30 km grid resolution over South Asia and the surrounding regions. Result shows that increased dust emission and transport from the Middle East/West Asia region induces a strong dust-cloud-radiation-precipitation-circulation feedback, resulting in a colder surface over the desert regions and the western Tibetan Plateau, but warmer and moister troposphere with enhanced cloudiness and precipitation over the Pakistan/Northwest India (PNWI) region. The latter changes are amplified by the dust aerosol induced Elevated Heat Pump (EHP) mechanism along the West Himalayas/Iranian Plateau foothill regions, most pronounced during May-June. During July-August, cloud radiation feedback further enhances the warming of the upper troposphere, and cooling of the land surface over the PNWI and adjacent regions over West Asia. The upper tropospheric heating and increased precipitation over PNWI spur a large-scale anomalous Rossby wavetrain and a northward displacement of the subtropical jetstream, manifesting in a contraction of the South Asian High and westward shift of the Western Pacific Subtropical High. As a result, the entire ASM precipitation-cloud system is displaced westward. Precipitation and cloudiness are intensified over northwest and western India and west Asia, but suppressed over southern and central East Asia. Analyses of the S2S variability of the upper level vorticity balance suggests that heating by Middle East dust plays an important role in exciting, and anchoring a teleconnection pattern through interactions among dust-cloud radiation, precipitation heating over the PNWI, the development of an upper level Rossby wavetrain and the northward shift of the boreal summer jetstream over Eurasia.

Modeling extreme precipitation over East China with a global variable-resolution modeling framework (MPASv5.2): impacts of resolution and physics

Abstract. The non-hydrostatic atmospheric Model for Prediction Across Scales (MPAS-A), a global variable-resolution modeling framework, is applied at a range of resolutions from hydrostatic (60, 30, 16 km) to non-hydrostatic (4 km) scales using regional refinement over East Asia to simulate an extreme precipitation event. The event is triggered by a typical wind shear in the lower layer of the Meiyu front in East China on 25–27 June 2012 during the East Asian summer monsoon season. The simulations are evaluated using ground observations and reanalysis data. The simulated distribution and intensity of precipitation are analyzed to investigate the sensitivity to model configuration, resolution, and physics parameterizations. In general, simulations using global uniform-resolution and variable-resolution meshes share similar characteristics of precipitation and wind in the refined region with comparable horizontal resolution. Further experiments at multiple resolutions reveal the significant impacts of horizontal resolution on simulating the distribution and intensity of precipitation and updrafts. More specifically, simulations at coarser resolutions shift the zonal distribution of the rain belt and produce weaker heavy precipitation centers that are misplaced relative to the observed locations. In comparison, simulations employing 4 km cell spacing produce more realistic features of precipitation and wind. The difference among experiments in modeling rain belt features is mainly due to the difference in simulated wind shear formation and evolution during this event. Sensitivity experiments show that cloud microphysics have significant effects on modeling precipitation at non-hydrostatic scales, but their impacts are relatively small compared to that of convective parameterizations for simulations at hydrostatic scales. This study provides the first evidence supporting the use of convection-permitting global variable-resolution simulations for studying and improving forecasting of extreme precipitation over East China and motivates the need for a more systematic study of heavy precipitation events and the impacts of physics parameterizations and topography in the future. The key points are as follows. Model for Prediction Across Scales (MPAS) simulations at global uniform and variable resolutions share similar characteristics of precipitation and wind in the refined region. Numerical experiments reveal significant impacts of resolution on simulating the distribution and intensity of precipitation and updrafts. This study provides evidence supporting the use of convection-permitting global variable-resolution simulation to study extreme precipitation.

Effects of Localized Grid Refinement on the General Circulation and Climatology in the Community Atmosphere Model

Abstract Using the spectral element (SE) dynamical core within the National Center for Atmospheric Research–Department of Energy Community Atmosphere Model (CAM), a regionally refined nest at 0.25° (~28 km) horizontal resolution located over the North Atlantic is embedded within a global 1° (~111 km) grid. A 23-yr simulation using Atmospheric Model Intercomparison Project (AMIP) protocols and default CAM, version 5, physics is compared to an identically forced run using the global 1° (~111 km) grid without refinement. The addition of a refined patch over the Atlantic basin does not noticeably affect the global circulation. In the area where the refinement is located, large-scale precipitation increases with the higher resolution. This increase is partly offset by a decrease in precipitation resulting from convective parameterizations, although total precipitation is also slightly higher at finer resolutions. Equatorial waves are not significantly impacted when traversing multiple grid spacings. Despite the grid transition region bisecting northern Africa, local zonal jets and African easterly wave activity are highly similar in both simulations. The frequency of extreme precipitation events increases with resolution, although this increase is restricted to the refined patch. Topography is better resolved in the nest as a result of finer grid spacing. The spatial patterns of variables with strong orographic forcing (such as precipitation, cloud, and precipitable water) are improved with local refinement. Additionally, dynamical features, such as wind patterns, associated with steep terrain are improved in the variable-resolution simulation when compared to the uniform coarser run.

Aquaplanet Experiments Using CAM’s Variable-Resolution Dynamical Core

Abstract A variable-resolution option has been added within the spectral element (SE) dynamical core of the U.S. Department of Energy (DOE)–NCAR Community Atmosphere Model (CAM). CAM-SE allows for static refinement via conforming quadrilateral meshes on the cubed sphere. This paper investigates the effect of mesh refinement in a climate model by running variable-resolution (var-res) simulations on an aquaplanet. The variable-resolution grid is a 2° (~222 km) grid with a refined patch of 0.25° (~28 km) resolution centered at the equator. Climatology statistics from these simulations are compared to globally uniform runs of 2° and 0.25°. A significant resolution dependence exists when using the CAM version 4 (CAM4) subgrid physical parameterization package across scales. Global cloud fraction decreases and equatorial precipitation increases with finer horizontal resolution, resulting in drastically different climates between the uniform grid runs and a physics-induced grid imprinting in the var-res simulation. Using CAM version 5 (CAM5) physics significantly improves cloud scaling at different grid resolutions. Additional precipitation at the equator in the high-resolution mesh results in collocated zonally anomalous divergence in both var-res simulations, although this feature is much weaker in CAM5 than CAM4. The equilibrium solution at each grid spacing within the var-res simulations captures the majority of the resolution signal of the corresponding globally uniform grids. The var-res simulation exhibits good performance with respect to wave propagation, including equatorial regions where waves pass through grid transitions. In addition, the increased frequency of high-precipitation events in the refined 0.25° area within the var-res simulations matches that observed in the global 0.25° simulations.

Exploring a Global Multiresolution Modeling Approach Using Aquaplanet Simulations*

Abstract Results from aquaplanet experiments performed using the Model for Prediction across Scales (MPAS) hydrostatic dynamical core implemented within the Department of Energy (DOE)–NCAR Community Atmosphere Model (CAM) are presented. MPAS is an unstructured-grid approach to climate system modeling that supports both quasi-uniform and variable-resolution meshing of the sphere based on conforming grids. Using quasi-uniform simulations at resolutions of 30, 60, 120, and 240 km, the authors evaluate the performance of CAM-MPAS via its kinetic energy spectra, general circulation, and precipitation characteristics. By analyzing an additional variable-resolution simulation with grid spacing that varies from 30 km in a spherical, continental-sized equatorial region to 240 km elsewhere, the CAM-MPAS’s potential for use as a regional climate simulation tool is explored. Similar to other quasi-uniform aquaplanet simulations, tropical precipitation increases with resolution, indicating the resolution sensitivity of the physical parameterizations. Comparison with the finite volume (FV) dynamical core suggests a weaker tropical circulation in the CAM-MPAS simulations, which is evident in reduced tropical precipitation and a weaker Hadley circulation. In the variable-resolution simulation, the kinetic energy spectrum within the high-resolution region closely resembles the quasi-uniform 30-km simulation, indicating a robust simulation of the fluid dynamics. As suggested by the quasi-uniform simulations, the CAM4 physics behave differently in the high and low resolution regions. A positive precipitation anomaly occurs on the western edge of the high-resolution region, exciting a Gill-type response; this zonal asymmetry represents the errors incurred in a variable resolution setting. When paired with a multiresolution mesh, the aquaplanet test case offers an exceptional opportunity to examine the response of physical parameterizations to grid resolution.

Simulating Global and North American Climate Using the Global Environmental Multiscale Model with a Variable-Resolution Modeling Approach

Abstract Results from two simulations using the Global Environmental Multiscale (GEM) model in a variable-resolution modeling approach are evaluated. Simulations with a highly resolved domain positioned over North America and over the tropical Pacific–eastern Indian Ocean are assessed against the GEM uniform grid control run, 40-yr ECMWF Re-Analysis (ERA-40), and available observations in terms of regional and global climate and interannual variability. It is found that the variable-resolution configurations realistically simulate global and regional climate over North America with seasonal means and variability generally closer to ERA-40 or observations than the control run. Systematic errors of the control run are still present within the variable-resolution simulations but alleviated to some extent over their respective highly resolved domains. Additionally, there is some evidence of performance deterioration due to the increased resolution. There is little evidence that an increased resolution over the tropical Pacific–eastern Indian Ocean, with better-resolved local processes (e.g., convection and equatorial waves), has a significant impact on the extratropical time mean fields. However, in terms of simulating the Northern Hemisphere atmospheric flow anomaly associated with the dominant mode of sea surface temperature interannual variability in the equatorial eastern Pacific (i.e., El Niño), both stretched configurations have more realistic teleconnection patterns than the control run.

Global variable resolution simulations of Mississippi River basin rains of summer 1993

A variable resolution global model is used to simulate a 13 day period of heavy rainfall during the 1993 Mississippi River basin floods. The model allows interaction between the atmosphere of the inner, fine resolution region that covers most of the United States and the outer global domain. We describe the impact of flow generated by latent heating over the central United States and its interaction with the global atmosphere. For this case, global constraints on the divergent flow impose substantial modifications upon the rainfall evolution within the central region. The most accurate simulations of rainfall occur when both local divergence and local thermal modifications interact with their global counterparts. The moisture budget is strongly influenced by the imposition of both divergence and thermal fields, and the temporal and spatial distribution of horizontal moisture flux and rainfall change substantially with the degree of outer state prespecification.