Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surface moisture flux”

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 433 records · Page 24

The environmental influence on tropical cyclone precipitation

The intensity, spatial, and temporal changes in precipitation were examined in three North Atlantic hurricanes during 1989 (Dean, Gabrielle, and Hugo) using precipitation estimates made from Special Sensor Microwave/Imager (SSM/I) measurements. In addition, analyses from a barotropic hurricane forecast model and the European Centre for Medium-Range Weather Forecast model were used to examine the relationship between the evolution of the precipitation in these tropical cyclones and external forcing. The external forcing parameters examined were (1) mean climatological sea surface temperatures, (2) vertical wind shear, (3) environmental tropospheric water vapor flux, and (4) upper-tropospheric eddy relative angular momentum flux convergence. The analyses revealed that (1) the SSM/I precipitation estimates were able to delineate and monitor convective ring cycles similar to those observed with land-based and aircraft radar and in situ measurements; (2) tropical cyclone intensification was observed to occur when these convective rings propagated into the inner core of these systems (within 111 km of the center) and when the precipitation rates increased; (3) tropical cyclone weakening was observed to occur when these inner-core convective rings dissipated; (4) the inward propagation of the outer convective rings coincided with the dissipation of the inner convective rings when they came within 55 km of each other; (5) in regions with the combined warm sea surface temperatures (above 26 C) and low vertical wind shear (less than 5 m/s), convective rings outside the region of strong lower-tropospheric inertial stability could be initiated by strong surges of tropospheric moisture, while convective rings inside the region of strong lower-tropospheric inertial stability could be enhanced by upper-tropospheric eddy relative angular momentum flux convergence.

Rodgers, Edward B.↗

Prototype Hanford Barrier Performance Monitoring Report: Fiscal Year 2024

This report provides an annual update on performance monitoring of the Prototype Hanford Barrier (PHB) at the Hanford Site. The PHB is part of a long-term study that serves as the scientific basis for many of the evapotranspiration-capillary barrier designs currently used globally and for future engineered barrier designs planned for remedial actions over waste and demolition sites at the Hanford Site. The PHB allows us to identify potential future issues and develop better monitoring techniques before these engineered barriers are constructed. Surface barriers like the PHB are essential for preventing water infiltration and curbing the spread of contaminants to groundwater. Effective monitoring of soil moisture levels above and below these barriers is crucial given that performance metrics could span up to 1,000 years. From July 2023 to June 2024, the water flux (as measured by tipping buckets) through the 2-m-thick silt loam layer of the PHB – a fine silty material that stores water under high capillary tension – remained well below the 0.5-mm-per-year performance threshold, demonstrating its effectiveness in preventing water penetration. In 2024, degraded tipping bucket gauges were replaced to ensure the accuracy of future data. Additionally, neutron probes revealed that the wetting front from the rainy season only penetrated to a maximum depth of 1.2 m into the silt loam. This further demonstrated the barrier's efficiency, as the wetting front did not fully penetrate the 2-m-thick capillary barrier. The western gravel slope of the PHB, composed of Hanford Site pit gravel of varying sizes, demonstrated low flux rates. In contrast, the eastern riprap slope exhibited higher flux rates. Zhang (2017) found that the riprap side slope of the PHB had the highest drainage rates in January and the lowest in late summer or early fall, indicating significant summer evaporation. Drainage rates increased significantly under enhanced precipitation conditions, far exceeding the design criterion, which could lead to water infiltration into the waste zone. The study introduced the “edge effect,” where elevated drainage rates from the riprap may migrate laterally beneath the barrier, compromising its ability to isolate waste, and recommended expanding the barrier and conducting further research to mitigate this issue. This report provides a review on surface barrier edge effects and their potential impact to subsurface contaminant migration.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

How similar are temperature and humidity profiles in the unstable boundary layer?

A comparison was made between temperature and humidity profiles obtained from radiosonde measurements under unstable conditions over hilly terrain in northeastern Kansas during the summer and fall of 1987. The summer profiles were found to be similar in a layer over the height range given by z = 65(plus or minus 45) to 124(plus or minus 54)m, where z is the height above the ground surface reference level. Evidence for this similarity was the closeness of the scaled ratio of the temperature and humidity gradients at any given time to the mean Bowen ratio determined at ground-based surface flux stations. Above the layer, in the outer region of the atmospheric boundary layer, it was found that the profiles of humidity and potential temperature are completely dissimilar in terms of the surface Bowen ratio. The profiles of the fall campaign when soil moisture content was low and vegetation senescent displayed a less pronounced similarity.

Sugita, Michiaki↗

Estimation of bare soil evaporation from a multi-frequency airborne SAR

An experiment was performed with the NASA/JPL airborne radar polarimeter in Sept. 1989 in an agricultural area near Fresno, CA. Based on 1st-order surface backscattering models, a physically based algorithm for retrieval of soil moisture and surface roughness was developed. It is shown that copolarization ratio is sensitive to soil moisture but not to soil roughness at high incidence angles (between 40 and 60 deg). The derived soil moisture was employed to derive a two-layer heat and energy flux model to estimate evaporation from bare soils.

Soares, Joao V.↗

Implications of summertime marine stratocumulus on the North American climate

This study focuses on the effects of summertime stratocumulus over the eastern Pacific. This cloud is linked to the semi-permanent sub-tropical highs that dominate the low-level circulation over the Pacific and Atlantic. Subsidence on the eastern flank of these highs creates an inversion based about 800 m above sea level that caps moist air near the surface. This air overlies cool waters driven by upwelling along the coastal regions of North America. Strong surface north-westerlies mix the boundary layer enough to saturate the air just below the capping inversion. Widespread stratocumulus is thus formed. All calculations were carried out using the GENESIS general circulation model that was run at MSFC. Among the more important properties of the model is that it includes radiative forcing due to absorption of solar radiation and the emission of infrared radiation, interactive clouds (both stratocumulus and cumulus types), exchanges of heat and moisture with the lower boundary. Clouds are interactive in the sense that they impact the circulation by modifying the fields of radiative heating and turbulent fluxes of heat and moisture in the boundary layer. In turn, clouds are modified by the winds through the advection of moisture. In order to isolate the effects of mid- and high-latitude stratocumulus, two runs were made with the model: one with and the other without stratocumulus. The runs were made for a year, but with perpetual July conditions, i.e., solar forcing was fixed. The diurnal solar cycle, however, was allowed for. The sea surface temperature distribution was fixed in both runs to represent climatological July conditions. All dependent variables were represented at 12 surfaces of constant sigma = p/p(sub O), where p is pressure and p(sub O) is surface pressure. To facilitate analysis, model output was transformed to constant pressure surfaces. Structures no smaller in size than 7.5 degrees longitude and 4.5 degrees in latitude were resolved. Smaller features of the circulation were parameterized. The model thus captures synoptic- and planetary-scale circulation features.

Clark, John H. E.↗

Size-Dependent Characteristics of Surface-Rooted Three-Dimensional Convective Objects in Continental Shallow Cumulus Simulations

A segmentation algorithm is applied to high resolution simulations of shallow continental convection to identify individual convective 3D objects within the convective boundary layer, in order to investigate which properties of the objects vary with the object width. The study analyses the geometry of the objects, along with their profiles of vertical velocity and total water, to assess various assumptions often used in spectral mass-flux convection schemes. The methodology of this paper is unique in that we use (a) a novel application of the watershed algorithm to detect individual objects in the constantly evolving continental boundary layer efficiently, and (b) an unprecedentedly large number of simulations being analyzed. In total, 26 days of LASSO simulations at the Atmospheric Radiation Measurement Southern Great Plains site are analyzed, yielding roughly one million objects. Plume-like surface-rooted objects are found to be omnipresent, the vertical extent of which is strongly dependent on the object width. The vertical velocity and moisture anomaly profiles of the widest objects are roughly consistent with the classic buoyancy-driven rising plume model. The kinematic and thermodynamic properties of the objects vary with object width. This width dependence is strongest above cloud base, but much weaker below. Finally the impact of neglecting the contribution of covariances to the vertical moisture flux, which is commonly used in mass-flux parameterizations, is investigated. The average effect of neglecting covariances increases linearly with object width, leading to a 20% flux underestimation for 2 km wide objects. Implications of the results for spectral convection scheme development are briefly discussed.

54 ENVIRONMENTAL SCIENCES↗

A new convective parameterization scheme for global models

A convective adjustment scheme is presented for parameterization of the vertical temperature and moisture structure for large-scale numerical weather models. The temperature and moisture structure is adjusted toward observed quasi-equilibrium structures to maintain a realistic vertical temperature and moisture structure in grid-scale global models during convection. Two reference thermodynamic structures are modeled analytically to account for shallow and deep convection. The thermodynamic structure obtained is permitted to change as a function of large-scale advective terms, radiation and surface fluxes. The reference profiles are applied when an energy integral provides a description of the shallow or deep convective processes occurring. Comparisons with GATE, BOMEX and ATEX data and a 10 day forecast for weather over the tropical ocean illustrate the climate and precipitation forecast accuracy improvements available with the adjustment scheme relative to standard parameterization schemes.

Miller, M. J.↗

Water cycle research associated with the CaPE hydrometeorology project (CHymP

One outgrowth of the Convection and Precipitation/Electrification (CaPE) experiment that took place in central Florida during July and August 1991 was the creation of the CaPE Hydrometeorology Project (CHymP). The principal goal of this project is to investigate the daily water cycle of the CaPE experimental area by analyzing the numerous land and atmosphere in situ and remotely sensed data sets that were generated during the 40-days of observations. The water cycle comprises the atmospheric branch. In turn, the atmospheric branch comprises precipitation leaving the base of the atmospheric volume under study, evaporation and transpiration entering the base, the net horizontal fluxes of water vapor and cloud water through the volume and the conversion of water vapor to cloud water and vice-versa. The sum of these components results in a time rate of change in the water and liquid water (or ice) content of the atmospheric volume. The components of the land branch are precipitation input to and evaporation and transpiration output from the surface, net horizontal fluxes of surface and subsurface water, the sum of which results in a time rate of change in surface and subsurface water mass. The objective of CHymP is to estimate these components in order to determine the daily water budget for a selected area within the CaPE domain. This work began in earnest in the summer of 1992 and continues. Even estimating all the budget components for one day is a complex and time consuming task. The discussions below provides a short summary of the rainfall quality assessment procedures followed by a plan for estimating the horizontal moisture flux.

Duchon, Claude E.↗

Relationships between evaprorative fraction and remotely sensed vegetation index and microwave brightness temperature for semiarid rangelands

Measurements of the microwave brightness temperature (TB) with the Pushbroom Microwave Radiometer (PBMR) over the Walnut Gulch Experiment Watershed were made on selected days during the MONSOON 90 field campaign. The PBMR is an L-band instrument (21-cm wavelength) that can provide estimates of near-surface soil moisture over a variety of surfaces. Aircraft observations in the visible and near-infrared wavelengths collected on selected days also were used to compute a vegetation index. Continuous micrometeorological measurements and daily soil moisture samples were obtained at eight locations during experimental period. Two sites were instrumented with time domain reflectometry probes to monitor the soil moisture profile. The fraction of available energy used for evapotranspiration was computed by taking the ratio of latent heat flux (LE) to the sum of net radiation (Rn) and soil heat flux (G). This ratio is commonly called the evaporative fraction (EF) and normally varies between 0 and 1 under daytime convective conditions with minimal advection. A wide range of environmental conditions existed during the field campaign, resulting in average EF values for the study area varying from 0.4 to 0.8 and values of TB ranging from 220 to 280 K. Comparison between measured TB and EF for the eight locations showed an inverse relationship. Other days were included in the analysis by estimating TB with the soil moisture data. Because transpiration from the vegetation is more strongly coupled to root zone soil moisture, significant scatter in this relationship existed at high values of TB or dry near-surface soil moisture conditions. The variation in EF under dry near-surface soil moisture conditions was correlated to the amount of vegetation cover estimated with a remotely sensed vegetation index. These findings indicate that information obtained from optical and microwave data can be used for quantifying the energy balance of semiarid areas. The microwave data can indicate when soil evaporation is significantly contributing to EF, while the optical data is helpful for quantifying the spatial variation in EF due to the distribution of vegetation cover.

Kustas, W. P.↗

Interannual variability of spring and summer monsoon growing season carbon exchange at a semiarid savanna over nearly two decades

Eddy covariance measurements of land-atmosphere energy, carbon, and water exchange now span multiple decades at some sites, supporting an improved understanding of flux interannual variability (IAV) and its ecophysiological and physical controls. Most eddy covariance IAV studies have focused on temperate forest ecosystems, where carbon fluxes are large and flux records are longest – but also where IAV is much lower than in dryland regions, which have been identified as an essential driver of the trend and variability in the global terrestrial carbon sink. In this study, we leveraged 19 years of continuous micrometeorological measurements at the AmeriFlux US-SRM mesquite savanna site in southern Arizona, USA to quantify the IAV, trends, and drivers of carbon fluxes during the distinct spring and summer growing seasons. We also assessed the ability of modern satellite and land surface models to capture the IAV of seasonal water and carbon fluxes. Annual net ecosystem production (NEP) was small and highly variable (23 +/- 64 gC m –2 yr –1 ). Precipitation and associated measures of water availability determined most of the variability in NEP, largely through their influence on annual and seasonal gross ecosystem productivity (GEP) as opposed to ecosystem respiration (ER). Root-zone soil moisture captured between 73% (spring) and 85% (summer) of GEP variability and between 73% (spring) and 58% (summer) of ER variability. Throughout the study period, soil moisture and greenness increased with associated increases in GEP, ER and NEP. These trends were strongly influenced by very productive and wet summer growing seasons during the last two years, which were characterized by abundant understory grass cover. Typically, less than half of the variability in growing season GEP and evapotranspiration was captured by satellite-based estimates and land surface model simulations with local site forcing and calibration, highlighting the ongoing utility of long-term datasets to support careful model testing and improvement.

54 ENVIRONMENTAL SCIENCES↗

Simulation of mesoscale convective response

Results are reported from a transfer of thermodynamic and dynamic data downscale from a two-dimensional mesoscale model to a two-dimensional cloud model and a subsequent examination of the differential convective response of the cloud model to a mesoscale structure. The mesoscale model included a high resolution PBL formulation, with convective forces expressed in a profile of exchange coefficients over the height of the PBL and the magnitude of the surface heat flux. Account was also taken of forcing by long- and short-wave radiation, surface forcing by the surface energy budget, the terrain shape, and possible wave reflection at the top boundary. The slab-symmetric cloud model possessed subgrid-scale features, five types of moisture terms, and a subroutine for accretion processes. The mesoscale environment was observed to have a significant impact on convective response, i.e., the vertical velocity and the cloud water. Various factors which were not included in the study and which must be considered are discussed.

Mcnider, R. T.↗

Measuring surface fluxes in CAPE

Two stations (site 1612 and site 2008) were operated by the University of Georgia group from 6 July 1991 to 18 August 1991. The following data were collected continuously: surface energy fluxes (i.e., net radiation, soil heat fluxes, sensible heat flux and latent heat flux), air temperature, vapor pressure, soil temperature (at 1 cm depth), and precipitation. Canopy reflectance and light interception data were taken three times at each site between 6 July and 18 August. Soil moisture content was measured twice at each site.

Kanemasu, E. T.↗

Cloud Modeling Using Field Project Data for the Study of Precipitation Processes

The use of cloud-resolving models (CRMs) in the study of precipitation process and their relation to the large-scale environment can be generally categorized into two approaches. The first approach is so called "cloud ensemble modeling". In this approach, many clouds of different size in various stages of their lifecycles can be present at any model simulation time. Large-scale effects are derived from observations and imposed into the model as the main forcing. The advantage of this approach is that the modeled convection will be forced to have the same intensity, thermodynamic budget and organization as the obserations.This approach will also allow CRMs to perform multi-day or multi-week time integrations. The second approach usually requires initial temperature and water vapor profiles that have a medium to large CAPE, and open lateral boundary conditions are used. The modeled clouds could be termed "self-forced convection". Model improvements, such as in the microphysics, are achieved using the second approach. In cloud ensemble modeling, accurate large-scale advective tendencies for temperature and water vapor are the main forcing for the CRMs. We found that the large-scale advective terms for temperature and water vapor are not always consistent, For example, large-scale forcing could indicate strong drying which would produce cooling in the model through evaporation but not contain large-scale advective heating to compensate. This discrepancy in forcing would cause differences between the observed and modeled latent heating profiles. Good measurements of other quantities (i.e., surface fluxes and radiation) are also required to perform variational objective analysis that computes and minimizes a "cost function" that constrains the difference between the large-scale advective forcing in temperature and water vapor. With self-forced convection, accurate vertical distributions of temperature, moisture (water vapor), and horizontal winds are required. The timing of the measurements relative to cloud development is crucial (i.e., prior to cloud triggering). Microphysical measurements (i.e., the cloud number concentration and size distribution) can also be used in this second approach but are of secondary importance with cloud ensemble modeling. In this paper, data collected during TRMM field campaigns (FCs; i.e., SCSMEX, LBA and KWAJEX) which were aimed at validating TRMM products (i.e., rainfall and the vertical distribution of latent heating) will be used to examine the impact of errors in the initial conditions (e.g.., soundings an large-scale forcing) on simulated rainfall distributions and brightness. Rainfall and precipitations simulated from a CRM will also be compared with those estimated by a schocastic model.

Tao, W.-K.↗

EOS Aqua AMSR-E Arctic Sea Ice Validation Program: Arctic2003 Aircraft Campaign Flight Report

In March 2003 a coordinated Arctic sea ice validation field campaign using the NASA Wallops P-3B aircraft was successfully completed. This campaign was part of the program for validating the Earth Observing System (EOS) Aqua Advanced Microwave Scanning Radiometer (AMSR-E) sea ice products. The AMSR-E, designed and built by the Japanese National Space Development Agency for NASA, was launched May 4, 2002 on the EOS Aqua spacecraft. The AMSR-E sea ice products to be validated include sea ice concentration, sea ice temperature, and snow depth on sea ice. This flight report describes the suite of instruments flown on the P-3, the objectives of each of the seven flights, the Arctic regions overflown, and the coordination among satellite, aircraft, and surface-based measurements. Two of the seven aircraft flights were coordinated with scientists making surface measurements of snow and ice properties including sea ice temperature and snow depth on sea ice at a study area near Barrow, AK and at a Navy ice camp located in the Beaufort Sea. Two additional flights were dedicated to making heat and moisture flux measurements over the St. Lawrence Island polynya to support ongoing air-sea-ice processes studies of Arctic coastal polynyas. The remaining flights covered portions of the Bering Sea ice edge, the Chukchi Sea, and Norton Sound.

Cavalieri, D. J.↗

Characterizing the Effects of Irrigation in the Middle East and North Africa Using Remotely Sensed Vegetation and Water Cycle Observations

A majority of the countries in the Middle East and North Africa (MENA) region suffer from water scarcity due in part to widespread rainfall deficits, unprecedented levels of water demand, and the inefficient use of renewable freshwater resources. Since a majority of the water withdrawal in the MENA is used for irrigation, there is a desperate need for improved understanding of irrigation practices and agricultural water use in the region. Here, satellite-derived irrigation maps and crop-type agricultural data are applied to the Land Data Assimilation System for the MENA region (MENA LDAS), designed to provide regional, gridded fields of hydrological states and fluxes relevant for water resources assessments. Within MENA-LDAS, the Catchment Land Surface Model (CLSM) simulates the location, timing, and amount of water applied through agricultural irrigation practices over the region from 2002-2012. In addition to simulating the irrigation impact on evapotranspiration, soil moisture, and runoff, we also investigate regional changes in terrestrial water storage (TWS) observed from the Gravity Recovery and Climate Experiment (GRACE) and simulated by CLSM.

Bolten, John↗

NASA/GEWEX SRB Surface Radiation Budget Fluxes Through the Prism of Weather States

The NASA/GEWEX Surface Radiation Budget Version 4 (SRB Rel4) product is now publicly available, covering the period from July 1983 through June 2017. This supersedes the SRB Rel3 (Stackhouse et al., 2011) product which has been used in the community for a wide variety of applications, including climate model validation, agriculture, solar energy, and architecture. SRB Rel4 uses the newly recalibrated and processed ISCCP HXS product as its primary input for cloud and radiance data, replacing ISCCP DX with a ninefold increase in pixel count (10km instead of 30km). This version retains a 1°x1° resolution but benefits from a much larger number of samples per grid box than the SRB Rel3. ISCCP also provides an atmospheric temperature and moisture dataset known as nnHIRS which we use here, along with Seaflux and Landflux surface and near-surface meteorological parameters. Rel4 incorporates several important algorithm improvements. These include recalculated shortwave (SW) atmospheric transmissivities and reflectivities yielding a somewhat less transmissive atmosphere. Both shortwave and longwave (LW) now also include variable aerosol composition and radiative properties, allowing for the use of a detailed aerosol history from the Max Planck Institute Aerosol Climatology (MACv1). LW and SW algorithms now produce pristine sky fluxes, allowing the aerosol flux effects to be quantified. For SW, ocean albedo and snow/ice albedo are improved from Release 3. Total solar irradiance is now variable, and reduced to an average of 1361 Wm-2. The radiative treatment of ice cloud is improved. For the LW, a climatological monthly varying spectral surface emissivity is added. Here we evaluate the SRB Rel4 top of atmosphere (TOA) and surface fluxes in the context of weather states and oceanic regions. Weather states (Tselioudis et al., 2013, Tselioudis et al., 2021) are defined by K-means clustering 2-D histograms of satellite-retrieved cloud optical depths and cloud top pressures. The most recent (Tselioudis et al., 2021) weather states derived from the International Satellite Cloud Climatology Project (ISCCP) H-series data results in eight cloud weather states and one clear sky weather state, for a total of nine weather states. We examine the statistics of SRB Rel4 fluxes for each weather state and a variety of relevant land and ocean regions. We compare to similar analyses performed on CERES EBAF and SYN1Deg flux products, and validate against Baseline Surface Radiometer Network (BSRN) land stations and Pacific Marine Environmental Laboratory (PMEL) ocean buoy data. It is demonstrated that the weather state concept is valuable when examining the comparative strengths and weaknesses of satellite radiative flux algorithms.

Atmospheric Radiation↗

Definition of a Technology Validation Mission for P-band Reflectometry using Signals of Opportunity

Root-Zone Soil Moisture (RZSM) (moisture profile in the top meter of soil) and Snow Water Equivalent (SWE) (total snow pack water content) are identified as priority target variables in the ESAS 2017 decadal survey [1] with critical roles in hydrology and water management. RZSM estimates are vital for understanding multiple Earth system processes and forecasting (for example, droughts [2]). Simultaneous knowledge of surface and RZSM could enable a breakthrough in estimating key unobserved hydrologic fluxes and reduce uncertainty in net ecosystem exchange (NEE), carbon balance [3] discharge estimates, and crop yield forecasts [4] .With the high albedo and insulating properties of snow, monitoring, SWE accumulation would provide a key constraint on the potential runoff during spring ablation while monitoring SWE disappearance rates would provide a key constraint on SWE partition into runoff vs. infiltration/recharge. [5] demonstrated that knowledge of early-spring SWE generally contributes most to streamflow forecast skill in the Western U.S. SWE is also a source of water storage that provides the water resources during spring snowmelt. Despite such potentially transformative contributions, accurate RZSM and SWE measurements are unattainable with current technology. While active/passive L-band methods (e.g. SMAP, SMOS) can reliably retrieve surface soil moisture in the top 5 cm of soil [6], [7]. RZSM estimates are only available through model assimilation of brightness temperatures with a radiative transfer and land surface models [8]. SWE estimation uses multi-frequency passive microwave techniques (e.g. [9]-[11]), which have significant problems with deeper snow and in forested and mountainous environments [12]. Signals of opportunity (SoOp) in P-band (200-400 MHz) is a new remote sensing technique with the capability of estimating both essential hydrologic variables, RZSM and SWE, circumventing many of the aforementioned limitations under all weather conditions day and night. SoOp is the re-utilization of existing powerful satellite transmissions within bands allocated for communications or navigation. P-band SoOp sensitivity to soil moisture has been demonstrated in an airborne experiment over Oklahoma in 2016 [13]. Recent theory [14] and experiments [15] have also confirmed that the reflection coefficient phase is proportional to SWE.

Garrison, J. L.↗

Urban weather modeling using WRF: linking physical assumptions, code implementation, and observational needs

The Weather Research and Forecasting (WRF) model includes urban schemes that simulate the influence of urban surfaces on the atmosphere using parameterizations for flux, and radiative exchanges. Three core schemes – the Bulk urban parameterization, Single-Layer Urban Canopy Model (SLUCM), and Multi-Layer Urban Canopy Model (MLUCM) – represent increasing levels of complexity. Although the parameterizations within these urban schemes are described in the literature, their specific implementation remains poorly documented, thus slowing down model development efforts. This manuscript provides a roadmap to the three urban schemes in WRF version 4.5.2, presenting equations using the same symbols as in the model code, along with references to code lines, and including graphics and explanations that connect the code to its physical foundations. Our thorough review of the urban parameterizations implemented in WRF version 4.5.2 highlighted a handful of parameters that may introduce discontinuities in simulations: (i) in the SLUCM, a 1 mm h −1 rain rate threshold is employed to switch between two minimum moisture availability parameterizations, thus impacting latent heat flux calculations; (ii) in the SLUCM a threshold is used to partition shortwave radiation into direct and diffuse components; (iii) in all three urban schemes, the bulk Richardson number is employed to select the similarity function, which influences the vertical distribution of heat and momentum. We also identified a highly simplified treatment of the radiative balance on roof surfaces. The implications of these simplifications can be assessed through targeted observations across relevant conditions, including varying precipitation rates, cloud cover, and transitions between stability regimes. Furthermore, the widespread application of the Monin-Obukhov similarity theory in these urban schemes warrants model evaluation under highly stable and unstable conditions and in heterogeneous urban settings with variable land cover and building heights on scales finer than model resolution. To address these challenges, we offer guidance on observational strategies, emphasizing the need for multi-parametric measurements to capture potential compensating biases and multi-height measurements that align with the levels where quantities are diagnostic and prognosed in the model (i.e., the lowest atmospheric level of the WRF model). Finally, our inspection of the code revealed implementation bugs that have now been corrected in WRF versions 4.6.0 and 4.6.1. Sensitivity tests over the Atlanta urban area show that these corrections affect surface temperatures, underscoring the importance of performing rigorous documentation and verification of the implementation of parameterizations in model code.

54 ENVIRONMENTAL SCIENCES↗