Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial statistics”

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 1,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Registration and Fusion of Multiple Source Remotely Sensed Image Data

Earth and Space Science often involve the comparison, fusion, and integration of multiple types of remotely sensed data at various temporal, radiometric, and spatial resolutions. Results of this integration may be utilized for global change analysis, global coverage of an area at multiple resolutions, map updating or validation of new instruments, as well as integration of data provided by multiple instruments carried on multiple platforms, e.g. in spacecraft constellations or fleets of planetary rovers. Our focus is on developing methods to perform fast, accurate and automatic image registration and fusion. General methods for automatic image registration are being reviewed and evaluated. Various choices for feature extraction, feature matching and similarity measurements are being compared, including wavelet-based algorithms, mutual information and statistically robust techniques. Our work also involves studies related to image fusion and investigates dimension reduction and co-kriging for application-dependent fusion. All methods are being tested using several multi-sensor datasets, acquired at EOS Core Sites, and including multiple sensors such as IKONOS, Landsat-7/ETM+, EO1/ALI and Hyperion, MODIS, and SeaWIFS instruments. Issues related to the coregistration of data from the same platform (i.e., AIRS and MODIS from Aqua) or from several platforms of the A-train (i.e., MLS, HIRDLS, OMI from Aura with AIRS and MODIS from Terra and Aqua) will also be considered.

LeMoigne, Jacqueline↗
NASA NTRS2017

Finding the Needles in the Haystacks: High-Fidelity Models of the Modern and Archean Solar System for Simulating Exoplanet Observations

We present two state-of-the-art models of the solar system, one corresponding to the present day and one to the Archean Eon 3.5 billion years ago. Each model contains spatial and spectral information for the star, the planets, and the interplanetary dust, extending to 50 au from the Sun and covering the wavelength range 0.3-2.5 micron. In addition, we created a spectral image cube representative of the astronomical backgrounds that will be seen behind deep observations of extrasolar planetary systems, including galaxies and Milky Way stars. These models are intended as inputs to high-fidelity simulations of direct observations of exoplanetary systems using telescopes equipped with high-contrast capability. They will help improve the realism of observation and instrument parameters that are required inputs to statistical observatory yield calculations, as well as guide development of post-processing algorithms for telescopes capable of directly imaging Earth-like planets.

extrasolar planetary systems↗
DOE OSTI2020

Three-dimensional super line-localization in low signal-to-noise microscope images via prior-apprised unsupervised learning (PAUL)

Biological processes such as processive enzyme turnover and intracellular cargo tracking involve the dynamic motion of a small "article" along a curvilinear biopolymer track. To understand these processes that occur across multiple length and time scales, one must acquire both the trajectory of the particle and the position of the track along which it moves, possibly by combining high-resolution single-particle tracking with conventional microscopy. Yet, usually there is a significant resolution mismatch between these modalities: while the tracked particle is localized with a precision of 10 nm, the image of the surroundings is limited by optical difraction, with 200 nm lateral and 500 nm axial resolutions. Compared to the particle's trajectory, the surrounding curvilinear structure appears as a blurred and noisy image. This disparity in the spatial resolutions of the particle trajectory and the surrounding curvilinear structure image makes data reconstruction, as well as interpretation, particularly challenging. Analysis is further complicated when the curvilinear structures are oriented arbitrarily in 3D space. Here, we present a prior-apprised unsupervised learning (PAUL) approach to extract information from 3D images where the underlying features resemble a curved line such as a filament or microtubule. This three-stage framework starts with a Hessian-based feature enhancement, which is followed by feature registration, where local line segments are detected on repetitively sampled subimage tiles. In the final stage, statistical learning, segments are clustered based on their geometric relationships. Principal curves are then approximated from each segment group via statistical tools including principal component analysis, bootstrap and kernel transformation. This procedure is characterized on simulated images, where sub-voxel medium deviations from true curves have been achieved. The 3D PAUL approach has also been implemented for successful line localization in experimental 3D images of gold nanowires obtained using a multifocal microscope. Lastly, this work not only bridges the resolution gap between two microscopy modalities, but also allows us to conduct 3D super line-localization imaging experiments, without using super-resolution techniques.

42 ENGINEERING↗
DOE OSTI2024

Fire–precipitation interactions amplify the quasi-biennial variability in fires over southern Mexico and Central America

Fires have great ecological, social, and economic impact. However, fire prediction and management remain challenges due to a limited understanding of their roles in the Earth system. Fires over southern Mexico and Central America (SMCA) are a good example of this, greatly impacting local air quality and regional climate. Here we report that the spring peak (April–May) of fire activities in this region has a distinct quasi-biennial signal based on multiple satellite datasets measuring different fire characteristics. The variability is initially driven by quasi-biennial variations in precipitation. Composite analysis indicates that strong fire years correspond to suppressed ascending motion and weakened precipitation over the SMCA. The anomalous precipitation over the SMCA is further found to be mostly related to the East Pacific–North Pacific (EP-NP) pattern 2 months prior to the fire season. The positive phase of the EP-NP leads to enhanced precipitation over the eastern US but suppressed precipitation over the SMCA, similar to the spatial pattern of precipitation differences between strong and weak fire years. Meanwhile, the quasi-biennial signals in precipitation and fires appear to be amplified by their interactions through a positive feedback loop at short timescales. Model simulations show that in strong fire years, more aerosol particles are released and transported downstream over the Gulf of Mexico and the eastern US, where suspended light-absorbing aerosols warm the atmosphere and cause the ascending motion of the air aloft. Subsequently, a compensating downward motion is formed over the region of the fire source and ultimately suppresses precipitation and intensifies fires. Statistical analysis shows the different durations of the two-way interaction, where the fire suppression effect of precipitation lasts for more than 20 d, while fire leads to a decrease in precipitation at shorter timescales (3–5 d). This study demonstrates the importance of fire–climate interactions in shaping the fire activities on an interannual scale and highlights how precipitation–fire interactions at short timescales contribute to the interannual variability in both fire and precipitation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1989

Context dependent anti-aliasing image reconstruction

Image Reconstruction has been mostly confined to context free linear processes; the traditional continuum interpretation of digital array data uses a linear interpolator with or without an enhancement filter. Here, anti-aliasing context dependent interpretation techniques are investigated for image reconstruction. Pattern classification is applied to each neighborhood to assign it a context class; a different interpolation/filter is applied to neighborhoods of differing context. It is shown how the context dependent interpolation is computed through ensemble average statistics using high resolution training imagery from which the lower resolution image array data is obtained (simulation). A quadratic least squares (LS) context-free image quality model is described from which the context dependent interpolation coefficients are derived. It is shown how ensembles of high-resolution images can be used to capture the a priori special character of different context classes. As a consequence, a priori information such as the translational invariance of edges along the edge direction, edge discontinuity, and the character of corners is captured and can be used to interpret image array data with greater spatial resolution than would be expected by the Nyquist limit. A Gibb-like artifact associated with this super-resolution is discussed. More realistic context dependent image quality models are needed and a suggestion is made for using a quality model which now is finding application in data compression.

Beaudet, Paul R.↗
NASA NTRS2007

Single-Shot Scalar-Triplet Measurements in High-Pressure Swirl-Stabilized Flames for Combustion Code Validation

In support of NASA ARMD's code validation project, we have made significant progress by providing the first quantitative single-shot multi-scalar data from a turbulent elevated-pressure (5 atm), swirl-stabilized, lean direct injection (LDI) type research burner operating on CH4-air using a spatially-resolved pulsed-laser spontaneous Raman diagnostic technique. The Raman diagnostics apparatus and data analysis that we present here were developed over the past 6 years at Glenn Research Center. From the Raman scattering data, we produce spatially-mapped probability density functions (PDFs) of the instantaneous temperature, determined using a newly developed low-resolution effective rotational bandwidth (ERB) technique. The measured 3-scalar (triplet) correlations, between temperature, CH4, and O2 concentrations, as well as their PDF s, also provide a high-level of detail into the nature and extent of the turbulent mixing process and its impact on chemical reactions in a realistic gas turbine injector flame at elevated pressures. The multi-scalar triplet data set presented here provides a good validation case for CFD combustion codes to simulate by providing both average and statistical values for the 3 measured scalars.

Kojima, Jun↗
NASA NTRS2017

Hypervolatiles in a Jupiter-Family Comet: Observations of 45P/Honda-Mrkos-Pajdusakova Using iSHELL at the NASA-IRTF

We used the new high spectral resolution cross-dispersed facility spectrograph, iSHELL, at the NASA Infrared Telescope Facility on Maunakea, HI, to observe Jupiter-family comet (JFC) 45P/Honda-Mrkos-Pajdusakova. We report water production rates, as well as production rates and abundance ratios relative to H2O, for eight trace parent molecules (native ices), CO, CH4, H2CO, CH3OH, HCN, NH3, C2H2, and C2H6, on 2 days spanning UT 2017 January 6/7 and 7/8, shortly following perihelion. Trace species were measured simultaneously with H2O and/or OH prompt emission, a proxy for H2O production, thereby providing a robust and consistent means of establishing the native ice composition of 45P. Its favorable geocentric radial velocity (approximately minus 35 kilometers per second) permitted sensitive measures of the "hypervolatiles" CO and CH4, which are substantially undercharacterized in JFCs. Our results represent the most precise ground-based measures of CO and CH4 to date in a JFC, providing a foundation for building meaningful statistics regarding their abundances. The abundance ratio for CH4 in 45P (0.79 percent plus or minus 0.06 percent relative to H2O) was consistent with its median value as measured among Oort Cloud comets, whereas CO (0.60 percent plus or minus 0.04 percent) was strongly depleted. Compared with all measured comets, HCN (0.049 percent plus or minus 0.012 percent) was strongly depleted, CH3OH (3.6 percent plus or minus 0.3 percent) was enriched, and the remaining species were consistent with their respective median abundances. The volatile composition measured for 45P could indicate processing of ices prior to their incorporation into its nucleus. Spatial analysis of emissions suggests enhanced release of more volatile species into the sunward-facing hemisphere of the coma.

DiSanti, Michael A.↗
NASA NTRS2019

Can We Use the QA4ECV Black-sky Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) using AVHRR Surface Reflectance to Assess Terrestrial Global Change?

NOAA platforms provide the longest period of terrestrial observation since the 1980s. The progress in calibration, atmospheric corrections and physically based land retrieval offers the opportunity to reprocess these data for extending terrestrial product time series. Within the Quality Assurance for Essential Climate Variables (QA4ECV) project, the black-sky Joint Research Centre (JRC)-fraction of absorbed photosynthetically active radiation (FAPAR) algorithm was developed for the AVHRR sensors on-board NOAA-07 to -16 using the Land Surface Reflectance Climate Data Record. The retrieval algorithm was based on the radiative transfer theory, and uncertainties were included in the products. We proposed a time and spatial composite for providing both 10-day and monthly products at 0.05º × 0.05º. Quality control and validation were achieved through benchmarking against third-party products, including Sea-Viewing Wide Field-of-View Sensor (SeaWiFS) datasets produced with the same retrieval algorithm. Past ground-based measurements, providing a proxy of FAPAR, showed good agreement of seasonality values over short homogeneous canopies and mixed vegetation. The average difference between SeaWiFS and QA4ECV monthly products over 2002–2005 is about 0.075 with a standard deviation of 0.091. We proposed a monthly linear bias correction that reduced these statistics to 0.02 and 0.001. The complete harmonized long-term time series was then used to address its fitness for the purpose of analysis of global terrestrial change.

AVHRR↗
DOE OSTI2021

Identifying geological structures through microseismic cluster and burst analyses complementing active seismic interpretation

At the Decatur carbon capture and storage site (IL, USA) CO 2 has been injected from 2011–2014 and from 2017 to present near the base of the Lower Mt. Simon Sandstone saline reservoir, resulting in microseismicity. Microseismicity is mainly located in the basement and distributed in distinct spatial clusters. The lack of significant impedance contrasts within the basement makes the interpretation of active-source seismic reflection data challenging, however, recent reprocessing allowed to resolve faults above and at the top of the basement. These faults generally do not coincide with the location of microseismic events and their continuation to the general depth of the seismic events cannot be assumed. This paper shows how the interpretation of the microseismicity can complement structural interpretations of active-source seismic reflection data. In particular, we analyze clusters and bursts (abrupt increases) of microseismicity, identify unresolved, smaller-scale weaknesses and extract statistical parameters. These parameters allow comparisons with the interpreted faults, and with fracture sets intercepted by boreholes. During injection at the Decatur site, the injection pressure was kept far below fracture pressure, nevertheless, seismic events were induced and spread far beyond the expected extent of the CO 2 plume. We argue that local stress transfers related to the CO 2 injection reactivated pre-existing fractures within the critically stressed basement. Finally, we conducted a slip tendency analysis for faults interpreted from active seismic, selected cluster, bursts and nodal planes from focal mechanisms to determine if the interpreted structures are optimally oriented with respect to the stress regime. Our results suggest that the orientation of fractures close to the injection well, generally shows slight deviations from the optimal orientation for slip. This might indicate either slight local deviations of the maximum horizontal stress azimuth from the average direction used in the analysis, or the lack of optimally oriented fractures at this location.

58 GEOSCIENCES↗
DOE OSTI2022

Optically Induced Anisotropy in Time-Resolved Scattering: Imaging Molecular-Scale Structure and Dynamics in Disordered Media with Experiment and Theory

Time-resolved scattering experiments enable imaging of materials at the molecular scale with femtosecond time resolution. However, in disordered media they provide access to just one radial dimension thus limiting the study of orientational structure and dynamics. Here we introduce a rigorous and practical theoretical framework for predicting and interpreting experiments combining optically induced anisotropy and time-resolved scattering. Using impulsive nuclear Raman and ultrafast x-ray scattering experiments of chloroform and simulations, we demonstrate that this framework can accurately predict and elucidate both the spatial and temporal features of these experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Evaluation of Turbulence and Dispersion in Multiscale Atmospheric Simulations over Complex Urban Terrain during the Joint Urban 2003 Field Campaign

Abstract This paper evaluates the representation of turbulence and its effect on transport and dispersion within multiscale and microscale-only simulations in an urban environment. These simulations, run using the Weather Research and Forecasting Model with the addition of an immersed boundary method, predict transport and mixing during a controlled tracer release from the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma. This work extends the results of a recent study through analysis of turbulence kinetic energy and turbulence spectra and their role in accurately simulating wind speed, direction, and tracer concentration. The significance and role of surface heat fluxes and use of the cell perturbation method in the numerical simulation setup are also examined. Our previous study detailed the model development necessary for our multiscale simulations, examined model skill at predicting wind speeds and tracer concentrations, and demonstrated that dynamic downscaling from mesoscale to microscale through a sequence of nested simulations can improve predictions of transport and dispersion relative to a microscale-only simulation forced by idealized meteorology. Here, predictions are compared with observations to assess qualitative agreement and statistical model skill at predicting wind speed, wind direction, tracer concentration, and turbulent kinetic energy at locations throughout the city. We also investigate the scale distribution of turbulence and the associated impact on model skill, particularly for predictions of transport and dispersion. Our results show that downscaled large-scale turbulence, which is unique to the multiscale simulations, significantly improves predictions of tracer concentrations in this complex urban environment. Significance Statement Simulations of atmospheric transport and mixing in urban environments have many applications, including pollution modeling for urban planning or informing emergency response following a hazardous release. These applications include phenomena with spatial scales spanning from millimeters to kilometers. Most simulations resolve flow only within the urban area of interest, omitting larger scales of turbulence and regional influences. This study examines a method that resolves both the small and large-scale flow features. We evaluate simulation accuracy by comparing predictions with observations from an experiment involving the release of a tracer gas in Oklahoma City, Oklahoma, with emphasis on correctly modeling turbulent fluctuations. Our results demonstrate the importance of resolving large-scale flow features when predicting transport and dispersion in urban environments.

42 ENGINEERING↗
NASA NTRS1980

An error analysis of basic kinematic quantities

A description is presented of the computational procedures used to obtain spatial fields of relative vorticity, advection of vorticity, horizontal divergence, and vertical motion. An investigation is conducted of the effectiveness of these procedures to reduce the effects of deliberately introduced random errors, giving attention to the confidence which can be placed on values of the derived parameters. It is concluded that fields of 500 mb vorticity and 700 mb vertical velocity, which are frequently used in synoptic scale forecasting schemes, show minimal effects of the deliberately introduced errors. Vorticity advection is found to be the most sensitive term to input data errors. Large fluctuations in fields of this parameter are evident both statistically and qualitatively.

Berecek, E. M.↗
NASA NTRS2020

Impact of Convectively Detrained Ice Crystals on the Tropical Upper Troposphere and Lower Stratosphere

The role of convectively detrained ice crystals on the humidity of the tropical upper troposphere and lower stratosphere (UTLS) is investigated in simulations of cirrus clouds along trajectories launched from the 378K potential temperature level in the tropics. The one-dimensional (vertical) cloud model tracks individual ice crystals through their lifecycle beginning with detrainment from convection, followed by deposition growth, sedimentation and sublimation. Convective influence of the parcels is diagnosed by tracing the trajectories through time-dependent fields of convective cloud-top height adjusted to match the CloudSAT and CALIPSO statistics. Model simulations of UTLS water vapor and cloud fields are evaluated and constrained by comparison with Aura MLS and CALIPSO measurements. Preliminary results indicate sensitivity of the detrained ice crystal lifecycle to atmospheric conditions downstream of convection. Specifically, cooling (high relative humidity and supersaturation) downstream of convection leads to deposition growth and sedimentation of detrained ice crystals, resulting in net dehydration of the UTLS. In contrast, warming (low relative humidity and subsaturation) downstream of convection leads to sublimation of detrained ice crystals and subsequent hydration. As such, the impact of detrained ice crystals on the humidity of the UTLS exhibits distinct spatial variability. Detrained ice crystals predominantly dehydrate the UTLS in the tropical mean. Sensitivities to the convectively detrained ice crystal size and concentration are also examined using measurements from the StatoClim aircraft campaign. The importance of convectively detrained ice crystals will be discussed within the context of the overall contribution of convection to the lower stratospheric humidity.

Ueyama, Rei↗
NASA NTRS2005

Global Survey and Statistics of Radio-Frequency Interference in AMSR-E Land Observations

Radio-frequency interference (RFI) is an increasingly serious problem for passive and active microwave sensing of the Earth. To satisfy their measurement objectives, many spaceborne passive sensors must operate in unprotected bands, and future sensors may also need to operate in unprotected bands. Data from these sensors are likely to be increasingly contaminated by RFI as the spectrum becomes more crowded. In a previous paper we reported on a preliminary investigation of RFI observed over the United States in the 6.9-GHz channels of the Advanced Microwave Scanning Radiometer (AMSR-E) on the Earth Observing System Aqua satellite. Here, we extend the analysis to an investigation of RFI in the 6.9- and 10.7-GHz AMSR-E channels over the global land domain and for a one-year observation period. The spatial and temporal characteristics of the RFI are examined by the use of spectral indices. The observed RFI at 6.9 GHz is most densely concentrated in the United States, Japan, and the Middle East, and is sparser in Europe, while at 10.7 GHz the RFI is concentrated mostly in England, Italy, and Japan. Classification of RFI using means and standard deviations of the spectral indices is effective in identifying strong RFI. In many cases, however, it is difficult, using these indices, to distinguish weak RFI from natural geophysical variability. Geophysical retrievals using RFI-filtered data may therefore contain residual errors due to weak RFI. More robust radiometer designs and continued efforts to protect spectrum allocations will be needed in future to ensure the viability of spaceborne passive microwave sensing.

microwave remote sensing↗
DOE OSTI2023

Transboundary effects from idealized regional geoengineering

Regional geoengineering, by reflecting sunlight over a very limited spatial domain, might be considered as a means to target specific regional impacts of climate change. One of the obvious concerns raised by such approaches is the extent to which the resulting effects would be detectable well beyond the targeted region (e.g. in neighbouring countries). A few studies have explored this question for targeted regions that are still comparatively large. We consider idealized simulations with increased ocean albedo over relatively small domains; the Gulf of Mexico (0.23% of Earth's surface) and over the Australian Great Barrier Reef (0.07%), both with negligible global radiative forcing. Applied over these very small domains, the only statistically significant non-local changes we find are some limited reduction on summer precipitation in Florida in the Gulf of Mexico case (adjacent to the targeted region). The lack of transboundary effects suggests that governance needs for such targeted interventions are quite distinct from those for more global sunlight reflection.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Heterogeneous Multi-Domain Dataset Synthesis to Facilitate Privacy and Risk Assessments in Smart City IoT

The emergence of the Smart Cities paradigm and the rapid expansion and integration of Internet of Things (IoT) technologies within this context have created unprecedented opportunities for high-resolution behavioral analytics, urban optimization, and context-aware services. However, this same proliferation intensifies privacy risks, particularly those arising from cross-modal data linkage across heterogeneous sensing platforms. To address these challenges, this paper introduces a comprehensive, statistically grounded framework for generating synthetic, multimodal IoT datasets tailored to Smart City research. The framework produces behaviorally plausible synthetic data suitable for preliminary privacy risk assessment and as a benchmark for future re-identification studies, as well as for evaluating algorithms in mobility modeling, urban informatics, and privacy-enhancing technologies. As part of our approach, we formalize probabilistic methods for synthesizing three heterogeneous and operationally relevant data streams—cellular mobility traces, payment terminal transaction logs, and Smart Retail nutrition records—capturing the behaviors of a large number of synthetically generated urban residents over a 12-week period. The framework integrates spatially explicit merchant selection using K-Dimensional (KD)-tree nearest-neighbor algorithms, temporally correlated anchor-based mobility simulation reflective of daily urban rhythms, and dietary-constraint filtering to preserve ecological validity in consumption patterns. In total, the system generates approximately 116 million mobility pings, 5.4 million transactions, and 1.9 million itemized purchases, yielding a reproducible benchmark for evaluating multimodal analytics, privacy-preserving computation, and secure IoT data-sharing protocols. To show the validity of this dataset, the underlying distributions of these residents were successfully validated against reported distributions in published research. We present preliminary uniqueness and cross-modal linkage indicators; comprehensive re-identification benchmarking against specific attack algorithms is planned as future work. This framework can be easily adapted to various scenarios of interest in Smart Cities and other IoT applications. By aligning methodological rigor with the operational needs of Smart City ecosystems, this work fills critical gaps in synthetic data generation for privacy-sensitive domains, including intelligent transportation systems, urban health informatics, and next-generation digital commerce infrastructures.

IoT↗
DOE OSTI2023

Investigating the impact of environmental factors on West Nile virus human case prediction in Ontario, Canada

West Nile virus is the most common mosquito borne disease in North America and the leading cause of viral encephalitis. West Nile virus is primarily transmitted between birds and mosquitoes while humans are incidental, dead-end hosts. Climate change may increase the risk of human infections as climatic variables have been shown to affect the mosquito life cycle, biting rate, incubation period of the disease in mosquitoes, and bird migration patterns. We develop a zero-inflated Poisson model to investigate how human West Nile virus case counts vary with respect to mosquito abundance and infection rates, bird abundance, and other environmental covariates. We use a Bayesian paradigm to fit our model to data from 2010–2019 in Ontario, Canada. Our results show mosquito infection rate, temperature, precipitation, and crow abundance are positively correlated with human cases while NDVI and robin abundance are negatively correlated with human cases. We find the inclusion of spatial random effects allows for more accurate predictions, particularly in years where cases are higher. Our model is able to accurately predict the magnitude and timing of yearly West Nile virus outbreaks and could be a valuable tool for public health officials to implement prevention strategies to mitigate these outbreaks.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Limitations of Bin and Bulk Microphysics in Reproducing the Observed Spatial Structure of Light Precipitation

Abstract Coarse-gridded atmospheric models often account for subgrid-scale variability by specifying probability distribution functions (PDFs) of process rate inputs such as cloud and rainwater mixing ratios ( q c and q r , respectively). PDF parameters can be obtained from numerous sources: in situ observations, ground- or space-based remote sensing, or fine-scale modeling such as large-eddy simulation (LES). LES is appealing to constrain PDFs because it generates large sample sizes, can simulate a variety of cloud regimes/case studies, and is not subject to the ambiguities of observations. However, despite the appeal of using model output for parameterization development, it has not been demonstrated that LES satisfactorily reproduces the observed spatial structure of microphysical fields. In this study, the structure of observed and modeled microphysical fields are compared by applying bifractal analysis, an approach that quantifies variability across spatial scales, to simulations of a drizzling stratocumulus field that span a range of domain sizes, drop concentrations (a proxy for mesoscale organization), and microphysics schemes (bulk and bin). Simulated q c closely matches observed estimates of bifractal parameters that measure smoothness and intermittency. There are major discrepancies between observed and simulated q r properties, though, with bulk simulated q r consistently displaying the bifractal properties of observed clouds (smooth, minimally intermittent) rather than rain while bin simulations produce q r that is appropriately intermittent but too smooth. These results suggest fundamental limitations of bulk and bin schemes to realistically represent higher-order statistics of the observed rain structure.

54 ENVIRONMENTAL SCIENCES↗