Search NASA⌕ Search

SEARCH · Search NASA

Results for “radiosonde”

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 91 records · Page 5

Enhanced Boundary Layer Height Detection Using Ceilometer, Surface Meteorology, and Radiation Products With a Random Forest Ensemble Method

This study develops and evaluates a Random Forest (RF) model for estimating planetary boundary layer height (PBLH) using 9 years of data from the Atmospheric Radiation Measurement Southern Great Plains (ARM SGP) user facility, with potential application in the NOAA Surface Radiation (SURFRAD) Network. The model integrates ceilometer, surface meteorology, and radiation measurements, and is trained using thermodynamic PBLH estimates derived from radiosondes. This approach aims to bridge gaps between aerosol-based and thermodynamic-based PBLH estimates. The RF model outperformed traditional methods during daytime and better captured transition periods, demonstrating improved accuracy and robustness. At ARM SGP, it showed a substantial reduction in both bias and RMSE, with a bias near zero (−4.9 m) compared with traditional Haar Wavelet (HW) (70.9 m) and Vaisala BL-View software (124.1 m), and an RMSE of 303.2 m, lower than both BL-View (566.9 m) and HW (404.6 m). During daytime hours, RF consistently outperformed both alternatives, maintaining lower bias and RMSE across all periods. At a second evaluation site, RF achieved the lowest overall RMSE (323.7 m), similar to HW (326.4 m) and significantly better than BL-View (738.3 m). However, all models showed reduced accuracy under stable nighttime conditions, limiting the reliability of PBLH estimates. Key predictors for the model included the lifting condensation level height (LCLH), aerosol gradients, and month for seasonal variability. The study underscores the potential of integrating machine learning with multiple data sets such as surface energy and thermodynamic data to advance PBLH estimation.

boundary layer height↗

Atmospheric temperature, water vapour and liquid water path from two microwave radiometers during MOSAiC

The microwave radiometers HATPRO (Humidity and Temperature Profiler) and MiRAC-P (Microwave Radiometer for Arctic Clouds - Passive) continuously measured radiation emitted from the atmosphere throughout the Multidisciplinary drifting Observatory for the Study of the Arctic Climate (MOSAiC) expedition on board the research vessel Polarstern. From the measured brightness temperatures, we have retrieved atmospheric variables using statistical methods in a temporal resolution of 1 s covering October 2019 to October 2020. The integrated water vapour (IWV) is derived individually from both radiometers. In addition, we present the liquid water path (LWP), temperature and absolute humidity profiles from HATPRO. To prove the quality and to estimate uncertainty, the data sets are compared to radiosonde measurements from Polarstern. The comparison shows an extremely good agreement for IWV, with standard deviations of 0.08–0.19 kg m -2 (0.39–1.47 kg m -2 ) in dry (moist) situations. The derived profiles of temperature and humidity denote uncertainties of 0.7–1.8 K and 0.6–0.45 gm -3 in 0–2 km altitude.

54 ENVIRONMENTAL SCIENCES↗

Quantification and assessment of the atmospheric boundary layer height measured during the AWAKEN experiment by a scanning LiDAR

The atmospheric boundary layer (ABL) height plays a key role in many atmospheric processes as one of the dominant flow length scales. However, a systematic quantification of the ABL height over the entire range of scales (i.e., with periods ranging from one minute to one year) is still lacking in literature. In this work, the ABL height is quantified based on high-resolution measurements collected by a scanning pulsed Doppler LiDAR during the recent American WAKE experimeNt (AWAKEN) campaign. The high availability of ABL height estimates (≈2200 collected over one year and each of them based on 10-min averaged statistics) allows to robustly assess five different ABL height models, i.e., one for convective thermal conditions and four for stable conditions. Thermal condition is quantified by a stability parameter spanning three orders of magnitude and probed by near-ground 3D sonic anemometry. The free-atmosphere stability, quantified by the Brunt–Väisälä frequency, is both calculated from simultaneous radiosonde measurements and obtained from the best fit of two of the chosen ABL height models. Good agreement is found between the data and three of the chosen models, quantified by mean absolute errors on the ABL height between 281 and 585 m. Furthermore, the seasonal variability of the convective ABL height model parameters (−15% to +23% with respect to the year baseline) agrees with the variability of buoyancy-generated turbulence caused by the variation in solar radiation throughout the year.

17 WIND ENERGY↗

ARM-DOE/adi-vap-pblhtsonde

The Atmospheric Radiation Measurement Program developed a value added product that could run autonomously to determine planetary boundary height using radiosondes. This product has produced data for 117 unique sites and facilities for more than 10 years. The initial software was released around 2013 and data has been produced since 2001 for some sites.

Sivaraman, Chitra↗

AWARE: The Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment

The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) West Antarctic Radiation Experiment (AWARE) performed comprehensive meteorological and aerosol measurements and ground-based atmospheric remote sensing at two Antarctic stations using the most advanced instrumentation available. A suite of cloud research radars, lidars, spectral and broadband radiometers, aerosol chemical and microphysical sampling equipment, and meteorological instrumentation was deployed at McMurdo Station on Ross Island from December 2015 through December 2016. A smaller suite of radiometers and meteorological equipment, including radiosondes optimized for surface energy budget measurement, was deployed on the West Antarctic Ice Sheet between 4 December 2015 and 17 January 2016. AWARE provided Antarctic atmospheric data comparable to several well-instrumented high Arctic sites that have operated for many years and that reveal numerous contrasts with the Arctic in aerosol and cloud microphysical properties. These include persistent differences in liquid cloud occurrence, cloud height, and cloud thickness. Antarctic aerosol properties are also quite different from the Arctic in both seasonal cycle and composition, due to the continent’s isolation from lower latitudes by Southern Ocean storm tracks. Antarctic aerosol number and mass concentrations are not only non-negligible but perhaps play a more important role than previously recognized because of the higher sensitivities of clouds at the very low concentrations caused by the large-scale dynamical isolation. Antarctic aerosol chemical composition, particularly organic components, has implications for local cloud microphysics. The AWARE dataset, fully available online in the ARM Program data archive, offers numerous case studies for unique and rigorous evaluation of mixed-phase cloud parameterization in climate models.

54 ENVIRONMENTAL SCIENCES↗

Global Climate

The 2019 global lower tropospheric temperature (LTT), which encompasses the atmosphere from the surface to ~10 km, ranked second warmest in seven datasets and first or third in the remaining two (Fig. 2.7). These records extend back to 1958 using radiosonde (balloon-borne instrumentation) data and one reanalysis dataset (JRA55), which demonstrate reasonable agreement with the 40+ year satellite record (since late 1978) and two other reanalysis datasets (since 1979 and 1980, ERA5 and MERRA2, respectively). Aweak El Niño contributed to increased global temperatures as 2019 values were +0.44° to +0.68°C higher than the 1981–2010 average (depending on the dataset), being just slightly cooler (~0.07°C on average) than the record warm year of 2016. At least fourof the five globally complete datasets (ERA5, MERRA2, JRA55, RSS, UAH) recorded each of the four months—June, September, November, and December—as experiencing their warmest monthly global LTT.

58 GEOSCIENCES↗

TAMU TRACER: Targeted Mobile Measurements to Isolate the Impacts of Aerosols and Meteorology on Deep Convection

Difficulty in using observations to isolate the impacts of aerosols from meteorology on deep convection often stems from inability to resolve the spatiotemporal variations in the environment serving as the storm’s inflow region. During the DOE TRacking Aerosol Convection interactions ExpeRiment (TRACER) in June-September 2022, a Texas A&M University (TAMU) team conducted a mobile field campaign to characterize the meteorological and aerosol variability in airmasses that serve as inflow to convection across the ubiquitous mesoscale boundaries associated with the sea- and bay-breezes in the Houston, Texas, region. These boundaries propagate inland over the fixed DOE Atmospheric Radiation Measurement (ARM) sites. However, convection occurs on either or both the continental or maritime sides or along the boundary. The maritime and continental airmasses serving as convection inflow may be quite distinct, with different meteorological and aerosol characteristics that fixed-site measurements cannot simultaneously sample. Thus, a primary objective of TAMU TRACER was to provide mobile measurements similar to those at the fixed sites, but in the opposite airmass across these moving mesoscale boundaries. TAMU TRACER collected radiosonde, lidar, aerosol, cloud condensation nuclei (CCN), and ice nucleating particle (INP) measurements on 29 enhanced operations days covering a variety of maritime, continental, outflow, and pre-frontal airmasses. This paper summarizes the TAMU TRACER deployment and measurement strategy, instruments, available datasets, and provides sample cases highlighting differences between these mobile measurements and those made at the ARM sites. We also highlight the exceptional TAMU TRACER undergraduate student participation in high impact learning activities through forecasting and field deployment opportunities.

54 ENVIRONMENTAL SCIENCES↗

Improving Process Level Understanding of Boundary Layer Winds over the Northeast U.S. Shelf: The Third Wind Forecast Improvement Project (WFIP3)

The third Wind Forecast Improvement Project (WFIP3), a U.S. Department of Energy and National Oceanic and Atmospheric Administration sponsored investigation, sought to improve understanding of the physical phenomena in the atmosphere and ocean that dictate the structure and variability of wind and thermodynamic fields within the marine atmospheric boundary layer. WFIP3 focused on mesoscale and submesoscale flows -- including sea breezes, low-level jets, low-level clouds, and coastal storms -- and the ability of advanced numerical model parameterizations to represent them within fully coupled oceanic and atmospheric modeling systems and foundational weather forecast models. WFIP3 conducted a comprehensive 18-month observational study over the Northeast U.S. outer continental shelf, a high-use coastal zone, using a 3D multiscale sensor array to highly resolve the temporal, vertical, and horizontal structure of the coupled atmospheric and oceanic boundary layers. Multiple land-based study sites adjacent to the coastal ocean observed surface meteorology and vertical profiles of atmospheric properties via passive infrared and microwave radiometers, active lidars and radars, and radiosondes. At sea, an array of surface flux buoys and two vertical profiling lidar buoys observed both atmospheric and oceanic properties, augmented by land-based oceanographic radar systems and routine ship-based surveys. Intensive observations of the marine atmospheric boundary layer over the ocean was done from an air-sea interaction flux tower and extended deployments of a large autonomous barge platform. Numerous critical forecasting phenomena were observed that are being evaluated within regional coupled and uncoupled modeling systems, including the National Oceanic and Atmospheric Administration's foundational High-Resolution Rapid Refresh forecast model.

Kirincich, Anthony↗

Climatology of Severe Local Storm Environments and Synoptic-Scale Features over North America in ERA5 Reanalysis and CAM6 Simulation

Severe local storm (SLS) activity is known to occur within specific thermodynamic and kinematic environments. These environments are commonly associated with key synoptic-scale features—including southerly Great Plains low-level jets, drylines, elevated mixed layers, and extratropical cyclones—that link the large-scale climate to SLS environments. This work analyzes spatiotemporal distributions of both extreme values of SLS environmental parameters and synoptic-scale features in the ERA5 reanalysis and in the Community Atmosphere Model, version 6 (CAM6), historical simulation during 1980–2014 over North America. Compared to radiosondes, ERA5 successfully reproduces SLS environments, with strong spatiotemporal correlations and low biases, especially over the Great Plains. Both ERA5 and CAM6 reproduce the climatology of SLS environments over the central United States as well as its strong seasonal and diurnal cycles. ERA5 and CAM6 also reproduce the climatological occurrence of the synoptic-scale features, with the distribution pattern similar to that of SLS environments. Compared to ERA5, CAM6 exhibits a high bias in convective available potential energy over the eastern United States primarily due to a high bias in surface moisture and, to a lesser extent, storm-relative helicity due to enhanced low-level winds. Additionally, composite analysis indicates consistent synoptic anomaly patterns favorable for significant SLS environments over much of the eastern half of the United States in both ERA5 and CAM6, though the pattern differs for the southeastern United States. Overall, our results indicate that both ERA5 and CAM6 are capable of reproducing SLS environments as well as the synoptic-scale features and transient events that generate them.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Information Content of a Synergy of Ground-Based and Space-Based Infrared Sounders. Part I: Clear-Sky Environments

Abstract Profiles of atmospheric temperature and water vapor from remotely sensed platforms provide critical observations within the temporal and spatial gaps of the radiosonde network. The 2017 National Academies of Science Decadal Survey highlighted that observations of the planetary boundary layer (PBL) from the current space-based observing system are not of the necessary accuracy or resolution for monitoring and predicting high-impact weather phenomena. One possible solution to improving observations of the PBL is supplementing the existing space-based observing system with a network of ground-based profilers. A synthetic information content study is developed utilizing profiles from the Atmospheric Radiation Measurement (ARM) program sites at the Southern Great Plains (SGP), east North Atlantic (ENA), and North Slope of Alaska (NSA) to assess the benefits, in terms of degrees of freedom (DOF), vertical resolution, and uncertainties, of a synergy between the ground-based Atmospheric Emitted Radiance Interferometer (AERI) with space-based hyperspectral infrared (IR) sounders. A combination of AERI with any of the three polar-orbiting IR sounders: the Atmospheric Infrared Sounder (AIRS), the Cross-track Infrared Sounder (CrIS), or the Infrared Atmospheric Sounding Interferometer (IASI), results in a DOF increase of 30%–40% in the surface-to-700-hPa layer compared to the space-based instrument alone. Introducing AERI measurements to the observing system also results in significant improvements to vertical resolution and uncertainties in the bottom 1000 m of the atmosphere compared to CrIS measurements alone. A synergy of CrIS and AERI exceeds the 1-km-vertical-resolution goal set by the Decadal Survey in the lowest 1000 m.

54 ENVIRONMENTAL SCIENCES↗

Convective-storm environments in subtropical South America from high-frequency soundings during RELAMPAGO-CACTI

During the RELAMPAGO-CACTI field experiments in 2018-19, an unprecedented number of balloon-borne soundings were collected in Argentina. Radiosondes were launched from both fixed and mobile platforms, yielding 2712 soundings during the period 15 October 2018-30 April 2019. Approximately 20% of these soundings were collected by highly mobile platforms, strategically positioned for each intensive observing period, and launching approximately once per hour. The combination of fixed and mobile soundings capture both the overall conditions characterizing the RELAMPAGO-CACTI campaign, as well as the detailed evolution of environments supporting the initiation and upscale growth of deep convective storms, including some that produced hazardous hail and heavy rainfall. Episodes of frequent convection were characterized by sufficient quantities of moisture and instability for deep convection, along with deep-layer vertical wind shear supportive of organized or rotating storms. Eleven soundings showed most-unstable convective available potential energy (MUCAPE) exceeding 6000 J kg-1, comparable to the extreme instability observed in other parts of the world with intense deep convection. Parameters used to diagnose severe-storm potential showed that conditions were often favorable for supercells and severe hail, but not for tornadoes, primarily owing to insufficient low-level wind shear. High-frequency soundings also revealed the structure and evolution of the boundary layer leading up to convection initiation, convectively generated cold pools, the South American Low-Level Jet (SALLJ), and elevated nocturnal convection. This sounding dataset will enable improved understanding and prediction of convective storms and their surroundings in subtropical South America, as well as comparisons with other heavily studied regions such as the central United States that have not previously been possible.

Schumacher, Russ S.↗

Low-Level Mesoscale and Cloud-Scale Interactions Promoting Deep Convection Initiation

Data from scanning radars, radiosondes, and vertical profilers deployed during three field campaigns are analyzed to study interactions between cloud-scale updrafts associated with initiating deep moist convection and the surrounding environment. Three cases are analyzed in which the radar networks permitted dual-Doppler wind retrievals in clear air preceding and during the onset of surface precipitation. These observations capture the evolution of (i) the mesoscale and boundary layer flow, and (ii) low-level updrafts associated with deep moist convection initiation (CI) events yielding sustained or short-lived precipitating storms. The elimination of convective inhibition did not distinguish between sustained and unsustained CI events, though the vertical distribution of convective available potential energy may have played a role. The clearest signal differentiating the initiation of sustained versus unsustained precipitating deep convection was the depth of the low-level horizontal wind convergence associated with the mesoscale flow feature triggering CI, a sharp surface wind shift boundary, or orographic upslope flow. The depth of the boundary layer relative to the height of the LFC failed to be a consistent indicator of CI potential. Widths of the earliest detectable low-level updrafts associated with sustained precipitating deep convection were ~3–5 km, larger than updrafts associated with surrounding boundary layer turbulence (~1–3 km wide). It is hypothesized that updrafts of this larger size are important for initiating cells to survive the destructive effects of buoyancy dilution via entrainment.

54 ENVIRONMENTAL SCIENCES↗

Near-Cloud Atmospheric Ingredients for Deep Convection Initiation

A lack of routine environmental observations located near deepening cumulus congestus clouds limits verification of important theorized and simulated updraft–environment interaction processes occurring during deep convection initiation (CI). We analyze radiosonde profiles collected during several hundred CI events near a mountain range in central Argentina during the CACTI field campaign. Statistical analyses illustrate environmental conditions supporting radar-observed CI outcomes that span a spectrum of convective cell depths, widths, and durations, as well as events lacking precipitating convection. Tested environmental factors include a large variety of sounding-derived measurements of CAPE, CIN, moisture, terrain-relative winds, vertical shear, and lifted parcel properties, with supplemental model reanalysis of background larger-scale vertical motion. CAPE and CIN metrics do not consistently differentiate CI success from failure. Only a few environmental factors contain consistent monotonic relationships among the spectrum of cloud depths achieved during CI: (i) the depth and strength of background ascent, and (ii) the component of low-level flow oriented parallel to the ridgeline. These metrics suggest that the ability of the surrounding flow to lift parcels to their LFC and terrain-modified flow are consistently relevant processes for CI. Low- to midlevel relative humidity strongly discriminated between CI and non-CI events, likely reflecting entrainment-driven dilution processes. However, we could not confidently conclude that relative humidity similarly discriminated robust from marginal CI events. Circumstantial evidence was found linking cell width, an important cloud property governing the probability of CI, to LCL height, boundary layer depth, depth and magnitude of the CIN layer, and ambient wind shear.

Storm environments↗

The Potential Roles of Preexisting Airmass Boundaries on a Tornadic Supercell Observed by TORUS on 28 May 2019

On 28 May 2019, a tornadic supercell, observed as part of Targeted Observation by UAS and Radars of Supercells (TORUS) produced an EF-2 tornado near Tipton, Kansas. The supercell was observed to interact with multiple preexisting airmass boundaries. These boundaries and attendant air masses were examined using unoccupied aircraft system (UAS), mobile mesonets, radiosondes, and dual-Doppler analyses derived from TORUS mobile radars. The cool-side air mass of one of these boundaries was found to have higher equivalent potential temperature and backed winds relative to the warm-side air mass; features associated with mesoscale air masses with high theta-e (MAHTEs). It is hypothesized that these characteristics may have facilitated tornadogenesis. The two additional boundaries were produced by a nearby supercell and appeared to weaken the tornadic supercell. This work represents the first time that UAS have been used to examine the impact of preexisting airmass boundaries on a supercell, and it provides insights into the influence environmental heterogeneities can have on the evolution of a supercell.

54 ENVIRONMENTAL SCIENCES↗

Observed Variability in Convective Cell Characteristics and Near-Storm Environments across the Sea- and Bay-Breeze Fronts in Southeast Texas

Abstract During the DOE Atmospheric Radiation Measurement (ARM) Tracking Aerosol Convection Interactions Experiment (TRACER) IOP spanning June–September 2022, two fixed ARM sites and a mobile team concurrently sampled the airmass heterogeneity across sea- and bay-breeze fronts around the greater Houston metropolitan region. Here, we quantify the spatiotemporal variability between maritime (coastal/bay side of breeze fronts) and continental (inland side of breeze fronts) air masses over 15 IOP days characterized by strong sea-breeze forcing. We analyze environmental profile data from 177 radiosondes and use S- and C-band radar data to track and quantify the variability in attributes of more than 2300 shallow and transitioning cells across different air masses. The composite analysis of environmental profiles indicates that during the early afternoon, the sea-breeze maritime air mass exhibits lower convective available potential energy (CAPE) than the bay-breeze maritime air mass. As the sea breeze advances inland with time, CAPE within the maritime air mass exceeds that of the continental air mass to the north of the breeze fronts. In general, maritime cells have a larger mean composite reflectivity and cell widths than continental cells; however, the response varies between shallow and transitioning cells. Mean composite 20-dB Z echo-top heights, however, are similar across air masses for both shallow and transitioning cells. The continental and maritime inflow air mass for transitioning cells has significantly different mean values for mixed-layer entrainment CAPE, lifted condensation level, level of free condensation, boundary layer depth, and diluted equilibrium level. For shallow cells, only total precipitable water shows a significant difference. Significance Statement The greater Houston metropolitan area is a natural laboratory for understanding the individual impacts of background meteorology and aerosols on convective clouds. Due to its proximity to the Gulf Coast and Galveston Bay, the Houston region experiences a diurnal precipitation cycle in the summer, driven by convection triggered from sea- and bay-breeze fronts. These fronts act as a boundary between air masses with distinct thermodynamic and environmental characteristics. Convergence along these fronts and interactions between storm outflow and the fronts facilitate convection initiation in different mesoscale air masses. This study quantifies the heterogeneity among these air masses while investigating their influence on cloud microphysics. We find that the effect of airmass heterogeneity is more pronounced for the bulk microphysical properties in shallow clouds.

Sharma, Milind↗

Evaluation of the Rapid Refresh Numerical Weather Prediction Model over Arctic Alaska

Abstract Despite a need for accurate weather forecasts for societal and economic interests in the U.S. Arctic, thorough evaluations of operational numerical weather prediction in the region have been limited. In particular, the Rapid Refresh Model (RAP), which plays a key role in short-term forecasting and decision-making, has seen very limited assessment in northern Alaska, with most evaluation efforts focused on lower latitudes. In the present study, we verify forecasts from version 4 of the RAP against radiosonde, surface meteorological, and radiative flux observations from two Arctic sites on the northern Alaskan coastline, with a focus on boundary layer thermodynamic and dynamic biases, model representation of surface inversions, and cloud characteristics. We find persistent seasonal thermodynamic biases near the surface that vary with wind direction, and may be related to the RAP’s handling of sea ice and ocean interactions. These biases seem to have diminished in the latest version of the RAP (version 5), which includes refined handling of sea ice, among other improvements. In addition, we find that despite capturing boundary layer temperature profiles well overall, the RAP struggles to consistently represent strong, shallow surface inversions. Further, while the RAP seems to forecast the presence of clouds accurately in most cases, there are errors in the simulated characteristics of these clouds, which we hypothesize may be related to the RAP’s treatment of mixed-phase clouds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterizing warm atmospheric boundary layer over land by combining Raman and Doppler lidar measurements

PBL plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, for PBL characterization and process study, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary. Here, new approaches are developed to provide reliable PBLH and MLH to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) user facility. Compared with widely used lidar aerosol measurements for PBLH, WVMR is a better trace for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to re-construct 2-D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm seasons’ RL and ML measurements at the SGP site for five years to study warm season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable PBLH and MLH dataset will be valuable for PBL process study, model evolution, and PBL parameterization improvement.

54 ENVIRONMENTAL SCIENCES↗