Consequences of the Large-Scale Subsidence Rate on the Stably Stratified Atmospheric Boundary Layer Over the Arctic Ocean, as seen in Large-Eddy Simulations

The analysis of surface heat fluxes and sounding profiles from SHEBA indicated possible significant effects of subsidence on the structure of stably-stratified ABLs (Mirocha et al. 2005). In this study the influence of the large-scale subsidence rate on the stably stratified atmospheric boundary layer (ABL) over the Arctic Ocean during clear sky, winter conditions is investigated using a large-eddy simulation model. Simulations are conducted while varying the subsidence rate between 0, 0.001 and 0.002 ms{sup -1}, and the resulting quasi-equilibrium ABL structure and evolution are examined. Simulations conducted without subsidence yield ABLs that are deeper, more strongly mixed, and cool much more rapidly than were observed. The addition of a small subsidence rate significantly improves agreement between the simulations and observations regarding the ABL height, potential temperature profiles and bulk heating rates. Subsidence likewise alters the shapes of the surface-layer flux, stress and shear profiles, resulting in increased vertical transport of heat while decreasing vertical momentum transport. A brief discussion of the relevance of these results to parameterization of the stable ABL under subsiding conditions in large-scale numerical weather and climate prediction models is presented

An improved WRF for urban-scale and complex-terrain applications

Simulations of atmospheric flow through urban areas must account for a wide range of physical phenomena including both mesoscale and urban processes. Numerical weather prediction models, such as the Weather and Research Forecasting model (WRF), excel at predicting synoptic and mesoscale phenomena. With grid spacings of less than 1 km (as is required for complex heterogeneous urban areas), however, the limits of WRF's terrain capabilities and subfilter scale (SFS) turbulence parameterizations are exposed. Observations of turbulence in urban areas frequently illustrate a local imbalance of turbulent kinetic energy (TKE), which cannot be captured by current turbulence models. Furthermore, WRF's terrain-following coordinate system is inappropriate for high-resolution simulations that include buildings. To address these issues, we are implementing significant modifications to the ARW core of the Weather Research and Forecasting model. First, we are implementing an improved turbulence model, the Dynamic Reconstruction Model (DRM), following Chow et al. (2005). Second, we are modifying WRF's terrain-following coordinate system by implementing an immersed boundary method (IBM) approach to account for the effects of urban geometries and complex terrain. Companion papers detailing the improvements enabled by the DRM and the IBM approaches are also presented (by Mirocha et al., paper 13.1, and K.A. Lundquist et al., paper 11.1, respectively). This overview of the LLNL-UC Berkeley collaboration presents the motivation for this work and some highlights of our progress to date. After implementing both DRM and an IBM for buildings in WRF, we will be able to seamlessly integrate mesoscale synoptic boundary conditions with building-scale urban simulations using grid nesting and lateral boundary forcing. This multi-scale integration will enable high-resolution simulations of flow and dispersion in complex geometries such as urban areas, as well as new simulation capabilities in regions of complex terrain

Development of an Immersed Boundary Method to Resolve Complex Terrain in the Weather Research and Forecasting Model

Flow and dispersion processes in urban areas are profoundly influenced by the presence of buildings which divert mean flow, affect surface heating and cooling, and alter the structure of turbulence in the lower atmosphere. Accurate prediction of velocity, temperature, and turbulent kinetic energy fields are necessary for determining the transport and dispersion of scalars. Correct predictions of scalar concentrations are vital in densely populated urban areas where they are used to aid in emergency response planning for accidental or intentional releases of hazardous substances. Traditionally, urban flow simulations have been performed by computational fluid dynamics (CFD) codes which can accommodate the geometric complexity inherent to urban landscapes. In these types of models the grid is aligned with the solid boundaries, and the boundary conditions are applied to the computational nodes coincident with the surface. If the CFD code uses a structured curvilinear mesh, then time-consuming manual manipulation is needed to ensure that the mesh conforms to the solid boundaries while minimizing skewness. If the CFD code uses an unstructured grid, then the solver cannot be optimized for the underlying data structure which takes an irregular form. Unstructured solvers are therefore often slower and more memory intensive than their structured counterparts. Additionally, urban-scale CFD models are often forced at lateral boundaries with idealized flow, neglecting dynamic forcing due to synoptic scale weather patterns. These CFD codes solve the incompressible Navier-Stokes equations and include limited options for representing atmospheric processes such as surface fluxes and moisture. Traditional CFD codes therefore posses several drawbacks, due to the expense of either creating the grid or solving the resulting algebraic system of equations, and due to the idealized boundary conditions and the lack of full atmospheric physics. Meso-scale atmospheric boundary layer simulations, on the other hand, are performed by numerical weather prediction (NWP) codes, which cannot handle the geometry of the urban landscape, but do provide a more complete representation of atmospheric physics. NWP codes typically use structured grids with terrain-following vertical coordinates, include a full suite of atmospheric physics parameterizations, and allow for dynamic synoptic scale lateral forcing through grid nesting. Terrain following grids are unsuitable for urban terrain, as steep terrain gradients cause extreme distortion of the computational cells. In this work, we introduce and develop an immersed boundary method (IBM) to allow the favorable properties of a numerical weather prediction code to be combined with the ability to handle complex terrain. IBM uses a non-conforming structured grid, and allows solid boundaries to pass through the computational cells. As the terrain passes through the mesh in an arbitrary manner, the main goal of the IBM is to apply the boundary condition on the interior of the domain as accurately as possible. With the implementation of the IBM, numerical weather prediction codes can be used to explicitly resolve urban terrain. Heterogeneous urban domains using the IBM can be nested into larger mesoscale domains using a terrain-following coordinate. The larger mesoscale domain provides lateral boundary conditions to the urban domain with the correct forcing, allowing seamless integration between mesoscale and urban scale models. Further discussion of the scope of this project is given by Lundquist et al. [2007]. The current paper describes the implementation of an IBM into the Weather Research and Forecasting (WRF) model, which is an open source numerical weather prediction code. The WRF model solves the non-hydrostatic compressible Navier-Stokes equations, and employs an isobaric terrain-following vertical coordinate. Many types of IB methods have been developed by researchers; a comprehensive review can be found in Mittal and Iaccarino [2005]. To the authors knowledge, this is the first IBM approach that is able to use a pressure-based coordinate. The immersed boundary method presented here uses direct forcing, first suggested by Mohd-Yusof [1997], to impose a no-slip boundary condition. Additionally, the WRF model has been modified to include a no-slip bottom boundary condition enabling direct comparisons with the IBM solution for problems with gently sloping terrain. The accuracy and efficiency of the immersed boundary solver is examined within the context of a two-dimensional Witch of Agnesi hill. Results are also presented for two-dimensional flow over several blocks of New York City, which demonstrate the IB method's ability to handle extremely complex terrain with sharp corners and steep terrain gradients

Using Satellite-Derived Fire Arrival Times for Coupled Wildfire-Air Quality Simulations at Regional Scales of the 2020 California Wildfire Season

Wildfire frequency has increased in the Western US over recent decades, driven by climate change and a legacy of forest management practices. Consequently, human structures, health, and life are increasingly at risk due to wildfires. Furthermore, wildfire smoke presents a growing hazard for regional and national air quality. In response, many scientific tools have been developed to study and forecast wildfire behavior, or test interventions that may mitigate risk. In this study, we present a retrospective analysis of 1 month of the 2020 Northern California wildfire season, when many wildfires with varying environments and behavior impacted regional air quality. We simulated this period using a coupled numerical weather prediction model with online atmospheric chemistry, and compare two approaches to representing smoke emissions: an online fire spread model driven by remotely sensed fire arrival times and a biomass burning emissions inventory. First, we quantify the differences in smoke emissions and timing of fire activity, and characterize the subsequent impact on estimates of smoke emissions. Next, we compare the simulated smoke to surface observations and remotely sensed smoke; we find that despite differences in the simulated smoke surface concentrations, the two models achieve similar levels of accuracy. We present a detailed comparison between the performance and relative strengths of both approaches, and discuss potential refinements that could further improve future simulations of wildfire smoke. Finally, we characterize the interactions between smoke and meteorology during this event, and discuss the implications that increases in regional smoke may have on future meteorological conditions

Simulating atmosphere flow for wind energy applications with WRF-LES

Forecasts of available wind energy resources at high spatial resolution enable users to site wind turbines in optimal locations, to forecast available resources for integration into power grids, to schedule maintenance on wind energy facilities, and to define design criteria for next-generation turbines. This array of research needs implies that an appropriate forecasting tool must be able to account for mesoscale processes like frontal passages, surface-atmosphere interactions inducing local-scale circulations, and the microscale effects of atmospheric stability such as breaking Kelvin-Helmholtz billows. This range of scales and processes demands a mesoscale model with large-eddy simulation (LES) capabilities which can also account for varying atmospheric stability. Numerical weather prediction models, such as the Weather and Research Forecasting model (WRF), excel at predicting synoptic and mesoscale phenomena. With grid spacings of less than 1 km (as is often required for wind energy applications), however, the limits of WRF's subfilter scale (SFS) turbulence parameterizations are exposed, and fundamental problems arise, associated with modeling the scales of motion between those which LES can represent and those for which large-scale PBL parameterizations apply. To address these issues, we have implemented significant modifications to the ARW core of the Weather Research and Forecasting model, including the Nonlinear Backscatter model with Anisotropy (NBA) SFS model following Kosovic (1997) and an explicit filtering and reconstruction technique to compute the Resolvable Subfilter-Scale (RSFS) stresses (following Chow et al, 2005).We are also modifying WRF's terrain-following coordinate system by implementing an immersed boundary method (IBM) approach to account for the effects of complex terrain. Companion papers presenting idealized simulations with NBA-RSFS-WRF (Mirocha et al.) and IBM-WRF (K. A. Lundquist et al.) are also presented. Observations of flow through the Altamont Pass (Northern California) wind farm are available for validation of the WRF modeling tool for wind energy applications. In this presentation, we use these data to evaluate simulations using the NBA-RSFS-WRF tool in multiple configurations. We vary nesting capabilities, multiple levels of RSFS reconstruction, SFS turbulence models (the new NBA turbulence model versus existing WRF SFS turbulence models) to illustrate the capabilities of the modeling tool and to prioritize recommendations for operational uses. Nested simulations which capture both significant mesoscale processes as well as local-scale stable boundary layer effects are required to effectively predict available wind resources at turbine height

A Space-based Observational Strategy for Characterizing the First Stars and Galaxies Using the Redshifted 21 cm Global Spectrum

© 2017. The American Astronomical Society. All rights reserved. The redshifted 21 cm monopole is expected to be a powerful probe of the epoch of the first stars and galaxies (10 < z < 35). The global 21 cm signal is sensitive to the thermal and ionization state of hydrogen gas and thus provides a tracer of sources of energetic photons-primarily hot stars and accreting black holes-which ionize and heat the high redshift intergalactic medium (IGM). This paper presents a strategy for observations of the global spectrum with a realizable instrument placed in a low-altitude lunar orbit, performing night-time 40-120 MHz spectral observations, while on the farside to avoid terrestrial radio frequency interference, ionospheric corruption, and solar radio emissions. The frequency structure, uniformity over large scales, and unpolarized state of the redshifted 21 cm spectrum are distinct from the spectrally featureless, spatially varying, and polarized emission from the bright foregrounds. This allows a clean separation between the primordial signal and foregrounds. For signal extraction, we model the foreground, instrument, and 21 cm spectrum with eigenmodes calculated via Singular Value Decomposition analyses. Using a Markov Chain Monte Carlo algorithm to explore the parameter space defined by the coefficients associated with these modes, we illustrate how the spectrum can be measured and how astrophysical parameters (e.g., IGM properties, first star characteristics) can be constrained in the presence of foregrounds using the Dark Ages Radio Explorer (DARE)

HERA Phase i Limits on the Cosmic 21 cm Signal: Constraints on Astrophysics and Cosmology during the Epoch of Reionization

Recently, the Hydrogen Epoch of Reionization Array (HERA) has produced the experiment's first upper limits on the power spectrum of 21 cm fluctuations at z ∼ 8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization from these limits. We find that the IGM must have been heated above the adiabatic-cooling threshold by z ∼ 8, independent of uncertainties about IGM ionization and the radio background. Combining HERA limits with complementary observations constrains the spin temperature of the z ∼ 8 neutral IGM to 27 K 630 K (2.3 K 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the cosmic microwave background dominates the z ∼ 8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones. The z ∼ 10 limits require even earlier heating if dark-matter interactions cool the hydrogen gas. If an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities of L r,ν /SFR &gt; 4 × 1024 W Hz-1 yr and L X /SFR &lt; 7.6 × 1039 erg s-1 yr. The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent Experiment to Detect the Global EOR Signature (EDGES) measurement. The framework described here provides a foundation for the interpretation of future HERA results

Improved Constraints on the 21 cm EoR Power Spectrum and the X-Ray Heating of the IGM with HERA Phase I Observations

We report the most sensitive upper limits to date on the 21 cm epoch of reionization power spectrum using 94 nights of observing with Phase I of the Hydrogen Epoch of Reionization Array (HERA). Using similar analysis techniques as in previously reported limits, we find at 95% confidence that Δ2(k = 0.34 h Mpc−1) ≤ 457 mK2 at z = 7.9 and that Δ2(k = 0.36 h Mpc−1) ≤ 3496 mK2 at z = 10.4, an improvement by a factor of 2.1 and 2.6, respectively. These limits are mostly consistent with thermal noise over a wide range of k after our data quality cuts, despite performing a relatively conservative analysis designed to minimize signal loss. Our results are validated with both statistical tests on the data and end-to-end pipeline simulations. We also report updated constraints on the astrophysics of reionization and the cosmic dawn. Using multiple independent modeling and inference techniques previously employed by HERA Collaboration, we find that the intergalactic medium must have been heated above the adiabatic cooling limit at least as early as z = 10.4, ruling out a broad set of so-called “cold reionization” scenarios. If this heating is due to high-mass X-ray binaries during the cosmic dawn, as is generally believed, our result’s 99% credible interval excludes the local relationship between soft X-ray luminosity and star formation and thus requires heating driven by evolved low-metallicity stars