Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sub scale”

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 379 records · Page 21

The three-point function as a probe of models for large-scale structure

We analyze the consequences of models of structure formation for higher order (n-point) galaxy correlation functions in the mildly nonlinear regime. Several variations of the standard Omega = 1 cold dark matter model with scale-invariant primordial perturbations have recently been introduced to obtain more power on large scales, R(sub p) is approximately 20/h Mpc, e.g., low matter-density (nonzero cosmological constant) models, 'tilted' primordial spectra, and scenarios with a mixture of cold and hot dark matter. They also include models with an effective scale-dependent bias, such as the cooperative galaxy formation scenario of Bower et al. We show that higher-order (n-point) galaxy correlation functions can provide a useful test of such models and can discriminate between models with true large-scale power in the density field and those where the galaxy power arises from scale-dependent bias: a bias with rapid scale dependence leads to a dramatic decrease of the the hierarchical amplitudes Q(sub J) at large scales, r is greater than or approximately R(sub p). Current observational constraints on the three-point amplitudes Q(sub 3) and S(sub 3) can place limits on the bias parameter(s) and appear to disfavor, but not yet rule out, the hypothesis that scale-dependent bias is responsible for the extra power observed on large scales.

Frieman, Joshua A.↗

Angstrom level profilometry for sub-millimeter to meter scale surface errors

A versatile noncontacting profilimetry approach is defined which involves the measurement of the local curvature of a test piece by simultaneously measuring its slope at two slightly displaced locations. Several extensions of the measurement technique are emphasized including beam expansion for long scans, measurement of circularity and cone angle of near cylindrical optics, and measurement of absolute flatness.

Glenn, Paul↗

Orbital Processing of Eutectic Rod-Like Arrays

The eutectic is one of only three solidification classes that exist. The others are isostructural and peritectic-class reactions, respectively. Simplistically, in a binaryeutectic phase diagram, a single liquid phase isothermally decomposes to two solid phases in a cooperative manner. The melting point minimum at the eutectic composition, isothermal solidification temperature, near-isocompositional solidification and refined solidification microstructure lend themselves naturally to such applications as brazing and soldering; industries that eutectic alloys dominate. Interest in direct process control of microstructures has led, more recently, to in-situ eutectic directional solidification with applications in electro-magnetics and electro-optics. In these cases, controlled structural refinement and the high aspect ratio and regularity of the distributed eutectic phases is highly significant to the fabrication and application of these in-situ natural composites. The natural pattern formation and scaling of the dispersed phase on a sub-micron scale has enormous potential application, since fabricating bulk materials on this scale mechanically has proven to be particularly difficult. It is thus of obvious importance to understand the solidification of eutectic materials since they are of great commercial significance. The dominant theory that describes eutectic solidification was derived for diffusion-controlled growth of alloys where both solid eutectic phases solidify metallically, i.e. without faceting at the solidification interface. Both high volume fraction (lamellar) and low volume fraction (rod-like) regular metallic arrays are treated by this theory. Many of the useful solders and brazements, however, and most of the regular in-situ composites are characterized by solidification reactions that are faceted/non-faceted in nature, rather than doubly non-faceted (metallic). Further, diffusion-controlled growth conditions are atypical terrestrially since gravitationally-driven convection is pervasive. As a consequence, it is important to determine whether these faceted/non-faceted composites behave in the same manner as their doubly non-faceted counterparts, particularly in the presence of convection. Prior analytical convective sensitivity testing of this theory predicted insensitivity. Prior experimental testing of this theory offered broad-based agreement between theory and experiment, though most results were for high volume fraction lamellar eutectics that solidified without faceting at the solidification interface. Directional solidification experiments of low volume fraction rod eutectics under damped (microgravity or magnetic field) conditions, however, have demonstrated significant sensitivity, challenging this fundamental theory. More recent theories have been proposed which introduce kinetic undercooling, faceting, fluid shear of the solute redistribution zone and the possibility that the interface composition is not the same as the bulk liquid composition. This program tests the established and proposed analytical theories and addresses the origins of discrepancies between the experimental and analytical results.

Larson, David J., Jr.↗

A first-principles derivation of Doppler noise expected from solar wind density fluctuations

The level of Doppler noise (DN) expected from solar wind (SW) density fluctuations (DF) is derived beginning with the expression for refractive index variations. The calculation takes account of up and downlink paths and of the method actually used to produce the DN values. The usual assumptions that the DF are frozen in, that the large-scale radial variation can be separated from the DF, that the DF power spectrum is a power law with outer scale k sub 0 and that the DF are homogeneous on scales less than 2c delta t, delta t = sample time, are made. The result agrees quite well with previous observations of DN. Corrections for the finite number of points used in the actual algorithm are discussed.

Callahan, P. S.↗

Hard Gamma Ray Emission from the Starburst Galaxy NGC 253

We have completed the study to search for hard gamma ray emission from the starburst galaxy NGC 253. Since supernovae are thought to provide the hard gamma ray emission from the Milky Way, starburst galaxies, with their extraordinarily high supernova rates, are prime targets to search for hard gamma ray emission. We conducted a careful search for hard gamma ray emission from NGC 253 using the archival data from the EGRET experiment aboard the CGRO. Because this starburst galaxy happens to lie near the South Galactic Pole, the Galactic gamma ray background is minimal. We found no significant hard gamma ray signal toward NGC 253, although a marginal signal of about 1.5 sigma was found. Because of the low Galactic background, we obtained a very sensitive upper limit to the emission of greater than 100 MeV gamma-rays of 8 x 10(exp -8) photons/sq cm s. Since we expected to detect hard gamma ray emission, we investigated the theory of gamma ray production in a dense molecular medium. We used a leaky-box model to simulate diffusive transport in a starburst region. Since starburst galaxies have high infrared radiation fields, we included the effects of self-Compton scattering, which are usually ignored. By modelling the expected gamma-ray and synchrotron spectra from NGC 253, we find that roughly 5 - 15% of the energy from supernovae is transferred to cosmic rays in the starburst. This result is consistent with supernova acceleration models, and is somewhat larger than the value derived for the Galaxy (3 - 10%). Our calculations match the EGRET and radio data very well with a supernova rate of 0.08/ yr, a magnetic field B approx. greater than 5 x 10(exp -5) G, a density n approx. less than 100/sq cm, a photon density U(sub ph) approx. 200 eV/sq cm, and an escape time scale tau(sub 0) approx. less than 10 Myr. The models also suggest that NGC 253 should be detectable with only a factor of 2 - 3 improvement in sensitivity. Our results are consistent with the standard picture of gamma-ray acceleration by supernovae.

Jackson, James M.↗

Land-Atmosphere Interactions: Successes, Problems and Prospects

After two decades of active research, a much better understanding of the broader role of biospheric processes on the local climate has emerged. A surface-albedo increase, particularly in desert border regions of the subtropics (as well as the deforested tropical regions), leads to a net surface energy deficit, which in turn leads to a relative sinking and reduced rainfall. On the other hand, studies of the influence of altered ratios of evapotranspiration and sensible fluxes, in situations where the net solar income is unchanged, show that evapotranspiration is a more desirable flux for increased precipitation and vitality of the biosphere. Besides providing water vapor and convective available potential energy (CAPE) to the lower troposphere, evapotranspiration helps in building larger CAPE before "turning on" the moist-convection. Larger CAPE in the lower troposphere enables convection to reach into the deeper atmosphere thereby heating the upper troposphere; indeed, moist-convection is also accompanied by the evaporation of falling precipitation that cools and moistens the lower atmosphere. While convective, as opposed to stratiform, precipitation reduces the fractional cloud cover; it also allows more solar radiation to reach the surface thereby invigorating surface fluxes. These, together with moist convection and associated downdrafts help to maintain the characteristic upper temperature limit(s) of the moist-land as well as oceanic regions. Regardless of the above understanding, several important problems continue to hinder the accurate simulation of a realistic land atmosphere interaction in a numerical model (both GCM and/or Meso-scale models). Among the unsolved problems are parameterization of sub-grid scale land processes that include small-scale variability of soil moisture, snow-cover and snow-physics, the biodiversity of the biosphere, orography, local drainage characteristics under natural conditions, and surface flow over the natural terrain. A well-known non-linear response of surface fluxes to these variations makes the problem of parameterizing land-atmosphere interaction processes hard-to-address and simulate, particularly in a GCM. In our presentation, we will discuss how orographic, snow-cover, and water table interactions can be included into a Simple Biosphere Model such as SiB/SSiB. Figure I shows how, in the Russian region, spring snowmelt affects the soil moisture profile. Corresponding figure 2 shows how interaction with the water table decreases the natural evapotranspiration in the Sahel region simulation. While these simulations need better validation with data, the simulations reveal that surface processes are sensitive to these parameterizations. With these developments, we continue to advance our understanding of the interaction of land with the atmosphere aloft, but the intrinsic variability of the newer parameters, e. g., hydraulic properties of the soil, diminish the positive influences of these advances on the improved climate simulation with GCMs.

Sud, Y. C.↗

Revisiting Scaling Relations for Giant Radio Halos in Galaxy Clusters

Many galaxy clusters host megaparsec-scale radio halos, generated by ultrarelativistic electrons in the magnetized intracluster medium. Correlations between the synchrotron power of radio halos and the thermal properties of the hosting clusters were established in the last decade, including the connection between the presence of a halo and cluster mergers. The X-ray luminosity and redshift-limited Extended GMRT Radio Halo Survey provides a rich and unique dataset for statistical studies of the halos. We uniformly analyze the radio and X-ray data for the GMRT cluster sample, and use the new Planck Sunyaev-Zel'dovich (SZ) catalog to revisit the correlations between the power of radio halos and the thermal properties of galaxy clusters. We find that the radio power at 1.4 GHz scales with the cluster X-ray (0.1-2.4 keV) luminosity computed within R(sub 500) as P(sub 1.4) approx. L(2.1+/-0.2) - 500). Our bigger and more homogenous sample confirms that the X-ray luminous (L(sub 500) > 5 × 10(exp 44) erg/s)) clusters branch into two populations-radio halos lie on the correlation, while clusters without radio halos have their radio upper limits well below that correlation. This bimodality remains if we excise cool cores from the X-ray luminosities. We also find that P(sub 1.4) scales with the cluster integrated SZ signal within R(sub 500), measured by Planck, as P(sub 1.4) approx. Y(2.05+/-0.28) - 500), in line with previous findings. However, contrary to previous studies that were limited by incompleteness and small sample size, we find that "SZ-luminous" Y(sub 500) > 6×10(exp −5) Mpc(exp 2) clusters show a bimodal behavior for the presence of radio halos, similar to that in the radio-X-ray diagram. Bimodality of both correlations can be traced to clusters dynamics, with radio halos found exclusively in merging clusters. These results confirm the key role of mergers for the origin of giant radio halos, suggesting that they trigger the relativistic particle acceleration.

radiation mechanisms↗

Hubble Space Telescope Medium Deep Survey. 2: Deconvolution of Wide Field Camera field galaxy images in the 13 hour + 43 deg field

We present isophotal profiles of six faint field galaxies from some of the first deep images taken for the Hubble Space Telescope (HST) Medium Deep Survey (MDS). These have redshifts in the range z = 0.126 to 0.402. The images were taken with the Wide Field Camera (WFC) in `parallel mode' and deconvolved with the Lucy method using as the point-spread function nearby stars in the image stack. The WFC deconvolutions have a dynamic range of 16 to 20 dB (4 to 5 mag) and an effective resolution approximately less than 0.2 sec (FWHM). The multiorbit HST images allow us to trace the morphology, light profiles, and color gradients of faint field galaxies down to V approximately equal to 22 to 23 mag at sub-kpc resolution, since the redshift range covered is z = 0.1 to 0.4. The goals of the MDS are to study the sub-kpc scale morphology, light profiles, and color gradients for a large samole of faint field galaxies down to V approximately equal to 23 mag, and to trace the fraction of early to late-type galaxies as function of cosmic time. In this paper we study the brighter MDS galaxies in the 13 hour + 43 deg MDS field in detail, and investigate to what extent model fits with pure exponential disks or a(exp 1/4) bulges are justified at V approximately less than 22 mag. Four of the six field galaxies have light profiles that indicate (small) inner bulges following r(exp 1/4) laws down to 0.2 sec resolution, plus a dominant surrounding exponential disk with little or no color gradients. Two occur in a group at z = 0.401, two are barred spiral galaxies at z = 0.179 and z = 0.302, and two are rather subluminous (and edge-on) disk galaxies at z = 0.126 and z = 0.179. Our deep MDS images can detect galaxies down to V, I approximately less than 25 to 26 mag, and demonstrate the impressive potential of HST--even with its pre-refurbished optics--to resolve morphological details in galaxies at cosmologically significant distances (v approximately less than 23 mag). Since the median redshift of these galaxies is approximately less than 0.4, the HST resolution allows us to study sub kpc size scales at the galaxy, which cannot be done with stable images over wide fields from the best ground-based sites.

Windhorst, R. A.↗

Consistency Between Sun-Induced Chlorophyll Fluorescence and Gross Primary Production of Vegetation in North America

Accurate estimation of the gross primary production (GPP) of terrestrial ecosystems is vital for a better understanding of the spatial-temporal patterns of the global carbon cycle. In this study,we estimate GPP in North America (NA) using the satellite-based Vegetation Photosynthesis Model (VPM), MODIS (Moderate Resolution Imaging Spectrometer) images at 8-day temporal and 500 meter spatial resolutions, and NCEP-NARR (National Center for Environmental Prediction-North America Regional Reanalysis) climate data. The simulated GPP (GPP (sub VPM)) agrees well with the flux tower derived GPP (GPPEC) at 39 AmeriFlux sites (155 site-years). The GPP (sub VPM) in 2010 is spatially aggregated to 0.5 by 0.5-degree grid cells and then compared with sun-induced chlorophyll fluorescence (SIF) data from Global Ozone Monitoring Instrument 2 (GOME-2), which is directly related to vegetation photosynthesis. Spatial distribution and seasonal dynamics of GPP (sub VPM) and GOME-2 SIF show good consistency. At the biome scale, GPP (sub VPM) and SIF shows strong linear relationships (R (sup 2) is greater than 0.95) and small variations in regression slopes ((4.60-5.55 grams Carbon per square meter per day) divided by (milliwatts per square meter per nanometer per square radian)). The total annual GPP (sub VPM) in NA in 2010 is approximately 13.53 petagrams Carbon per year, which accounts for approximately 11.0 percent of the global terrestrial GPP and is within the range of annual GPP estimates from six other process-based and data-driven models (11.35-22.23 petagrams Carbon per year). Among the seven models, some models did not capture the spatial pattern of GOME-2 SIF data at annual scale, especially in Midwest cropland region. The results from this study demonstrate the reliable performance of VPM at the continental scale, and the potential of SIF data being used as a benchmark to compare with GPP models.

photosynthesis model↗

A Viewing Angle-Kinetic Luminosity Scheme for BL Lacertae Objects

We propose a unified classification for BL Lacertae objects (BLs), focusing on the synchrotron peak frequency of the spectral energy distribution. The unification scheme is based on the angle theta that describes the orientation of the relativistic jet and on the electron kinetic luminosity Lamba(sub kin) of the jet. We assume that Lamba(sub kin) scales with the size of the jet (r) in a self-similar fashion (Lamba(sub kin) is proportional to r(exp 2)), as supported by observational data. The jets are self-similar in geometry and have the same pressure and median magnetic field at the inlet, independent of size. The self-similarity is broken for the highest energy electrons, which radiate mainly at high frequencies, since for large sources they suffer more severe radiative energy losses over a given fraction of the jet length. We calculate the optically thin synchrotron spectrum using an accelerating inner jet model based on simple relativistic gas dynamics and show that it can fit the observed infrared-to-X-ray spectrum of PKS 2155-304. We couple the accelerating jet model to the unification scheme and compare the results to complete samples of BLs. The negative apparent evolution of X-ray-selected BLs is explained as a result of positive evolution of the jet electron kinetic luminosity Lamba(sub kin). We review observational arguments in favor of the existence of scaled-down accretion disks and broad emission-line regions in BLs. The proposed unification scheme can explain the lack of observed broad emission lines in X-ray-selected BLs as well as the existence of those lines preferentially in luminous radio-selected BLs. Finally, we review observational arguments that suggest the extension of this unification scheme to all blazars.

Georganopoulos, Markos↗

The effect of small-scale vertical mixing of horizontal momentum in a general circulation model

Several experiments are described in which the sub-grid-scale vertical eddy viscosity in the GISS global general circulation model was varied. The results show that large viscosities suppress large-scale eddies in middle and high latitudes, but enhance the circulation in the tropical Hadley cell and increase the extent of the tropical easterlies. Comparison with observations shows that the GISS model requires eddy viscosities about 1 sq m per sec or less to give realistic results for middle and high latitudes, and eddy viscosities about 100 sq m per sec to give realistic results for low latitudes. A plausible mechanism for the implied increase in small-scale mixing in low latitudes is cumulus convection.

Stone, P. H.↗

Evaluation of NASA's High-Resolution Global Composition Simulations: Understanding a Pollution Event in the Chesapeake Bay During the Summer 2017 OWLETS Campaign

Recirculation of pollutants due to a bay breeze effect is a key meteorological mechanism impacting air quality near urban coastal areas, but regional and global chemical transport models have historically struggled to capture this phenomenon. We present a case study of a high ozone (O3) episode observed over the Chesapeake Bay during the NASA Ozone Water-Land Environmental Transition Study (OWLETS) in summer 2017. OWLETS included a complementary suite of ground-based and airborne observations, with which we characterize the meteorological and chemical context of this event and develop a framework to evaluate model performance. Two publicly-available NASA global high-resolution coupled chemistry-meteorology models (CCMMs) are investigated: GEOS-CF and MERRA2-GMI. The GEOS-CF R squared value for comparisons between the NASA Sherpa C-23 aircraft measurements to the GEOS-CF resulted in good agreement (R squared: 0.67) on July 19th and fair agreement (R squared: 0.55) for July 20th. Compared to surface observations, we find the GEOS-CF product with a 25 x 25 km squared grid box, at an hourly (R squared: 0.62 to 0.87) and 15-minute (R squared: 0.64 to 0.87) interval for six regional sites outperforms the hourly nominally 50 x 50 km squared gridded MERRA2-GMI (R squared: 0.53 to 0.76) for four of the six sites, suggesting it is better capable of simulating complex chemical and meteorological features associated with ozone transport within the Chesapeake Bay airshed. When the GEOS-CF product was compared to the TOLNet LiDAR observations at both NASA Langley Research Center (LaRC) and the Chesapeake Bay Bridge Tunnel (CBBT), the median differences at LaRC were -6 to 8% and at CBBT were ± 7% between 400 to 2000 m ASL. This indicates that, for this case study, the GEOS-CF is able to simulate surface level ozone diurnal cycles and vertical ozone profiles at small scales between the surface level and 2000 m ASL. Evaluating global chemical model simulations at sub-regional scales will help air quality scientists understand the complex processes occurring at small spatial and temporal scales within complex surface terrain changes, simulating nighttime chemistry and deposition, and the potential to use global chemical transport simulations in support of regional and sub-regional field campaigns.

NASA Ozone Water-Land Environmental Transition Stu↗

Renormalization group estimates of transport coefficients in the advection of a passive scalar by incompressible turbulence

The advection of a passive scalar by incompressible turbulence is considered using recursive renormalization group procedures in the differential sub grid shell thickness limit. It is shown explicitly that the higher order nonlinearities induced by the recursive renormalization group procedure preserve Galilean invariance. Differential equations, valid for the entire resolvable wave number k range, are determined for the eddy viscosity and eddy diffusivity coefficients, and it is shown that higher order nonlinearities do not contribute as k goes to 0, but have an essential role as k goes to k(sub c) the cutoff wave number separating the resolvable scales from the sub grid scales. The recursive renormalization transport coefficients and the associated eddy Prandtl number are in good agreement with the k-dependent transport coefficients derived from closure theories and experiments.

Zhou, YE↗

An Investigation of Differential Deposition for Figure Corrections in Full-Shell Grazing-Incidents X-Ray Optics

We are investigating differential deposition as a way of correcting small figure errors inside full-shell grazing-incidence x-ray optics. The optics in our study are fabricated using the electroformed-nickel-replication technique, and the figure errors arise from fabrication errors in the mandrel, from which the shells are replicated, as well as errors induced during the electroforming process. Combined, these give sub-micron-scale figure deviations which limit the angular resolution of the optics to approx. 10 arcsec. Sub-micron figure errors can be corrected by selectively depositing (physical vapor deposition) material inside the shell. The requirements for this filler material are that it must not degrade the ultra-smooth surface finish necessary for efficient x-ray reflection (approx. 5 A rms), and must not be highly stressed. In addition, a technique must be found to produce well controlled and defined beams within highly constrained geometries, as some of our mirror shells are less than 3 cm in diameter.

Gubarev, Mikhail V.↗

High Tc superconductors: The scaling of Tc with the number of bound holes associated with charge transfer neutralizing the multivalence cations

It is observed that for the known high-T(sub c) Cu-, Tl-, and Bi-based superconductors, T(sub c) scales consistently with the number of bound holes per unit cell which arise from charge transfer excitations of frequency approximately = 3 x 10(exp 13) that neutralized the multivalence cations into diamagnetic states. The resulting holes are established on the oxygens. Extrapolation of this empirical fit in the up-temperature direction suggests a T(sub c) of about 220-230 K at a value of 25 holes/unit cell (approximately the maximum that can be materials-engineered into a high-T(sub c) K2MnF4 or triple Perovskite structure). In the down-temperature direction, the extrapolation gives a T(sub c) in the vicinity of 235 K for the Y-Ba-Cu-O system as well as the known maximum temperature of 23 K for low-T(sub c) materials shown by Nb3Ge. The approach is also consistent with the experimental findings that only multivalence ions which are diamagnetic in their atomic state (Cu, Tl, Bi, Pb, and Sb) associate with high-T(sub c) compounds.

Vezzoli, G. C.↗

Electron Currents and Heating in the Ion Diffusion Region of Asymmetric Reconnection

In this letter the structure of the ion diffusion region of magnetic reconnection at Earths magnetopause is investigated using the Magnetospheric Multiscale (MMS) spacecraft. The ion diffusion region is characterized by a strong DC electric field, approximately equal to the Hall electric field, intense currents, and electron heating parallel to the background magnetic field. Current structures well below ion spatial scales are resolved, and the electron motion associated with lower hybrid drift waves is shown to contribute significantly to the total current density. The electron heating is shown to be consistent with large-scale parallel electric fields trapping and accelerating electrons, rather than wave-particle interactions. These results show that sub-ion scale processes occur in the ion diffusion region and are important for understanding electron heating and acceleration.

Graham, D. B.↗

3D Woven Composite End Ring Structure: Design Development and Verification Testing

The National Aeronautics and Space Administration (NASA) Space Technology and Mission Directorate (STMD) is advancing composite technologies for exploration missions. The Composite Technology for Exploration (CTE) project goal is to advance lightweight joint concepts for Space Launch System (SLS) class structures. To this end CTE is developing and demonstrating a 3D woven carbon /epoxy composite joint design co-bonded in a circumferential end ring configuration, applied to a SLS sized conceptual payload attachment structure. This joint is a lightweight alternative to heavy and costly metallic joints. The 3D woven joint development advances lightweight joint technologies through material and structural design, material and joint analyses, and fabrication, as well as demonstration through building block testing. Prior limited building block tests are mentioned. In this paper, new sub-element-scale test results are described that demonstrate the structural capability well beyond design limit load. Then, results for a sub-component test article with a 889 mm joint subjected to combined loads are presented.

Kenneth N Segal↗

Full-Coverage High-Resolution Daily PM(sub 2.5) Estimation using MAIAC AOD in the Yangtze River Delta of China

Satellite aerosol optical depth (AOD) has been used to assess population exposure to fine particulate matter (PM (sub 2.5)). The emerging high-resolution satellite aerosol product, Multi-Angle Implementation of Atmospheric Correction(MAIAC), provides a valuable opportunity to characterize local-scale PM(sub 2.5) at 1-km resolution. However, non-random missing AOD due to cloud snow cover or high surface reflectance makes this task challenging. Previous studies filled the data gap by spatially interpolating neighboring PM(sub 2.5) measurements or predictions. This strategy ignored the effect of cloud cover on aerosol loadings and has been shown to exhibit poor performance when monitoring stations are sparse or when there is seasonal large-scale missngness. Using the Yangtze River Delta of China as an example, we present a Multiple Imputation (MI) method that combines the MAIAC high-resolution satellite retrievals with chemical transport model (CTM) simulations to fill missing AOD. A two-stage statistical model driven by gap-filled AOD, meteorology and land use information was then fitted to estimate daily ground PM(sub 2.5) concentrations in 2013 and 2014 at 1 km resolution with complete coverage in space and time. The daily MI models have an average R(exp 2) of 0.77, with an inter-quartile range of 0.71 to 0.82 across days. The overall Ml model 10-fold cross-validation R(exp 2) (root mean square error) were 0.81 (25 gm(exp 3)) and 0.73 (18 gm(exp 3)) for year 2013 and 2014, respectively. Predictions with only observational AOD or only imputed AOD showed similar accuracy.Comparing with previous gap-filling methods, our MI method presented in this study performed bette rwith higher coverage, higher accuracy, and the ability to fill missing PM(sub 2.5) predictions without ground PM(sub 2.5) measurements. This method can provide reliable PM(sub 2.5)predictions with complete coverage that can reduce biasin exposure assessment in air pollution and health studies.

PM2.5↗