Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “inverse scattering”

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 955 records · Page 53

Dependence of Wave-Breaking Statistics on Wind Stress and Wave Development

Incidence of wave breaking for pure wind driven waves has been studied on Lake Washington at wind speeds up to 8 m/s. Video recordings were employed to identify and categorize the breaking events in terms of micro-scale, spilling and plunging breakers. These events were correlated with the magnitude of the wave spectrum measured with a resistance wire wave gauge and band pass filtered between 6 and 10 Hz. An equivalent percentage of breaking crests were found for spilling and plunging events. Wave forcing as measured by wind stress (or friction velocity, u(sub *), squared) and by inverse wave age, u(sub *)/Cp where Cp is the phase velocity of the waves at the peak of the frequency spectrum, were found to be good prerictors of percentage of breaking crests. When combined in a two parameter regression, those two variables gave small standard deviation and had a high correlation coefficient (66 percent). The combination of u(sub *)(exp 2) and u(sub *)/Cp can be understood in physical terms. Furthermore, for the larger values of u(sub *)(exp 2) the dependence of wave braking and wave age was stronger than at the low end of the values u(sub *)(exp 2) and u(sub *)/Cp. Thus, both the level of wave development as determined by inverse wave age, which we may term relative wind effectiveness for wave forcing and the wind forcing on the water surface determine the incidence of wave breaking. Substituting U(sub 10)(sup 3.75) (which is the dependence of whitecap cover found by Monahan and coworkers) an equivalent correlation was found to the prediction by u(sub *)(exp 2). Slightly better standard deviation value and higher correlation coefficient were found by using a Reynolds number as predictor. A two-parameter regression involving u(sub *)(exp 2) and a Reynold's number proposed by Toba and his colleagues which relates u(sub *)(exp 2) and peak wave frequency, improves the correlation even more but is less easy to interpret in physical terms. The equivalent percentage of breaking crests obtained in our previous study was reported at 8.6 percent for a short record obtained at U(sub 10N) of about 6 m/s. Typical values in the current study for similar conditions are 6 percent, which is consistent with the previous study in view of the scatter. In that study we did not have a video recording system, so the observed breaking may include more of the micro-scaic breaking events, and the value, 8.6 percent, is well within the range of highly probable sampling variability.

Katsaros, Kristina B.↗

Optimal Estimation Inversion of Ionospheric Electron Density from GNSS-POD Limb Measurements: Part II-Validation and Comparison Using NmF2 and hmF2

A growing number of SmallSat/CubeSat constellations with high-rate (50–100 Hz) global navigation satellite system radio occultations (GNSS-RO) as well as low-rate (1 Hz) precise orbit determination (GNSS-POD) limb-viewing capabilities provide unprecedented spatial and temporal sampling rates for ionospheric studies. In the F-region electron density (N 𝑒 ) retrieval process, instead of the conventional onion-peeling (OP) inversion, an optimal estimation (OE) inversion technique was recently developed using total electron content measurements acquired by GNSS-POD link. The new technique is applied to data acquired from the COSMIC-1, COSMIC-2, and Spire constellations. Although both OE and OP techniques use the Abel weighting function in N 𝑒 inversion, OE significantly differs in its performance, especially in the lower F- and E-regions. In this work, we evaluate and compare newly derived data sets using F2 peak properties with other space-based and ground-based observations. We determine the F2 peak N 𝑒 (NmF2) and its altitude (hmF2), and compare them with the OP-retrieved values. Good agreement is observed between the two techniques for both NmF2 and hmF2. In addition, we also utilize autoscaled F2 peak measurements from a number of worldwide Digisonde stations (∼30). The diurnal sensitivity and latitudinal variability of the F2 peak between the two techniques are carefully studied at these locations. Good agreement is observed between OE-retrieved NmF2 and Digisonde-measured NmF2. However, significant differences appear between OE-retrieved hmF2 and Digisonde-measured hmF2. During the daytime, Digisonde-measured hmF2 remains ∼25–45 km below the OE-retrieved hmF2, especially at mid and high latitudes. We also incorporate F-region N 𝑒 measurements from two incoherent scatter radar observations at high latitudes, located in the North American (Millstone Hill) and European (EISCAT at Tromso) sectors. The radar measurements show good agreement with OE-retrieved values. Although there are several possible sources of error in the ionogram-derived N 𝑒 profiles, our further analysis on F1 and F2 layers indicates that the low Digisonde hmF2 is caused by the autoscaled method, which tends to detect a height systematically below the F2 peak when the F1 layer is present.

GNSS radio occultation↗

Comparison of improved TMD and CGC frameworks in forward quark dijet production

For studying small- x gluon saturation in forward dijet production in high-energy dilute-dense collisions, the improved TMD (ITMD) factorization formula was recently proposed. In the Color Glass Condensate (CGC) framework, it represents the leading term of an expansion in inverse powers of the hard scale. It contains the leading-twist TMD factorization formula relevant for small gluon’s transverse momentum k t , but also incorporates an all-order resummation of kinematical twists, resulting in a proper matching to high-energy factorization at large k t . In this paper, we evaluate the accuracy of the ITMD formula quantitatively, for the case of quark dijet production in high-energy proton-proton( p + p ) and proton-nucleus ( p + A ) collisions at LHC energies. We do so by comparing the quark-antiquark azimuthal angle Δ ϕ distribution to that obtained with the CGC formula. For a dijet with each quark momentum p t much larger than the target saturation scale, Q s , the ITMD formula is a good approximation to the CGC formula in a wide range of azimuthal angle. It becomes less accurate as the jet p t ’s are lowered, as expected, due to the presence of genuine higher-twists contributions in the CGC framework, which represent multi-body scattering effects absent in the ITMD formula. We find that, as the hard jet momenta are lowered, the accuracy of ITMD start by deteriorating at small angles, in the high-energy-factorization regime, while in the TMD regime near Δ ϕ = π , very low values of p t are needed to see differences between the CGC and the ITMD formula. In addition, the genuine twists corrections to ITMD become visible for higher values of p t in p + A collisions, compared to p + p collisions, signaling that they are enhanced by the target saturation scale.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Coaction and double-copy properties of configuration-space integrals at genus zero

We investigate configuration-space integrals over punctured Riemann spheres from the viewpoint of the motivic Galois coaction and double-copy structures generalizing the Kawai-Lewellen-Tye (KLT) relations in string theory. For this purpose, explicit bases of twisted cycles and cocycles are worked out whose orthonormality simplifies the coaction. We present methods to efficiently perform and organize the expansions of configuration-space integrals in the inverse string tension α' or the dimensional-regularization parameter ϵ of Feynman integrals. Generating-function techniques open up a new perspective on the coaction of multiple polylogarithms in any number of variables and analytic continuations in the unintegrated punctures. We present a compact recursion for a generalized KLT kernel and discuss its origin from intersection numbers of Stasheff polytopes and its implications for correlation functions of two-dimensional conformal field theories. We find a non-trivial example of correlation functions in (p, 2) minimal models, which can be normalized to become uniformly transcendental in the p → ∞ limit.

Conformal Field Theory↗

A Laser Frequency Transverse Modulation Might Compensate for the Spectral Broadening Due to Large Electron Energy Spread in Thomson Sources

Compact laser plasma accelerators generate high-energy electron beams with increasing quality. When used in inverse Compton backscattering, however, the relatively large electron energy spread jeopardizes potential applications requiring small bandwidths. We present here a novel interaction scheme that allows us to compensate for the negative effects of the electron energy spread on the spectrum, by introducing a transverse spatial frequency modulation in the laser pulse. Such a laser chirp, together with a properly dispersed electron beam, can substantially reduce the broadening of the Compton bandwidth due to the electron energy spread. We show theoretical analysis and numerical simulations for hard X-ray Thomson sources based on laser plasma accelerators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Development and Testing of the New Surface LER Climatology for OMI UV Aerosol Retrievals

Ozone Monitoring Instrument (OMI) onboard Aura satellite retrieved aerosols properties using UV part of solar spectrum. The OMI near UV aerosol algorithm (OMAERUV) is a global inversion scheme which retrieves aerosol properties both over ocean and land. The current version of the algorithm makes use of TOMS derived Lambertian Equivalent Reflectance (LER) climatology. A new monthly climatology of surface LER at 354 and 388 nm have been developed. This will replace TOMS LER (380 nm and 354nm) climatology in OMI near UV aerosol retrieval algorithm. The main objectives of this study is to produce high resolution (quarter degree) surface LER sets as compared to existing one degree TOMS surface LERs, to product instrument and wavelength consistent surface climatology. Nine years of OMI observations have been used to derive monthly climatology of surface LER. MODIS derived aerosol optical depth (AOD) have been used to make aerosol corrections on OMI wavelengths. MODIS derived BRDF adjusted reflectance product has been also used to capture seasonal changes in the surface characteristics. Finally spatial and temporal averaging techniques have been used to fill the gaps around the globes, especially in the regions with consistent cloud cover such as Amazon. After implementation of new surface data in the research version of algorithm, comparisons of AOD and single scattering albedo (SSA) have been performed over global AERONET sites for year 2007. Preliminary results shows improvements in AOD retrievals globally but more significance improvement were observed over desert and bright locations. We will present methodology of deriving surface data sets and will discuss the observed changes in retrieved aerosol properties with respect to reference AERONET measurements.

Aerosol↗

RSPECT: A PROSPECT-based model incorporating the real structure of rice leaves

Radiative transfer models (RTMs) describe how light is absorbed, scattered, and transmitted within leaves by simulating mechanistic light propagation processes. The PROSPECT model is based on measurable parameters (the leaf biochemical content) and a non-measurable parameter (the leaf anatomical structure represented by the leaf structure parameter (N)). The effect of N on the optical properties of leaves has been investigated through a number of local and global sensitivity analyses. Other studies have directly evaluated the effect of the leaf anatomical structure on spectral reflectance, particularly in the near infrared region. However, the relationship between N and the anatomical structure is unclear. Here, in this study, we leveraged eLeaf, a ray tracing-based 3D rice leaf simulator, to establish relationships between leaf anatomical features and spectral properties, enabling us to replace N in the PROSPECT-4 model with measurable leaf anatomical parameters and develop the RSPECT model. The leaf thickness at minor vein, leaf thickness at bulliform cells, mesophyll thickness at minor vein, and distance between two minor veins could be used to predict N effectively. The RSPECT model achieved spectral simulation accuracy comparable to PROSPECT-4 and was more suitable for parameter inversion of the physical and chemical properties of rice leaves, with relative root mean square errors of 7.4% for chlorophyll content, 5.6% for equivalent water thickness, and 7.5% for dry matter content. In conclusion, RSPECT improves radiative transfer modeling by integrating measurable anatomical features and provides a framework for extending this approach to other vegetation types.

hyperspectral↗

Monte Carlo Bayesian Inference on a Statistical Model of Sub-gridcolumn Moisture Variability Using High-resolution Cloud Observations : Sensitivity Tests and Results - Part II

Part I presented a Monte Carlo Bayesian method for constraining a complex statistical model of GCM sub-gridcolumn moisture variability using high-resolution MODIS cloud data, thereby permitting large-scale model parameter estimation and cloud data assimilation. This part performs some basic testing of this new approach, verifying that it does indeed significantly reduce mean and standard deviation biases with respect to the assimilated MODIS cloud optical depth, brightness temperature and cloud top pressure, and that it also improves the simulated rotational-Ramman scattering cloud optical centroid pressure (OCP) against independent (non-assimilated) retrievals from the OMI instrument. Of particular interest, the Monte Carlo method does show skill in the especially difficult case where the background state is clear but cloudy observations exist. In traditional linearized data assimilation methods, a subsaturated background cannot produce clouds via any infinitesimal equilibrium perturbation, but the Monte Carlo approach allows finite jumps into regions of non-zero cloud probability. In the example provided, the method is able to restore marine stratocumulus near the Californian coast where the background state has a clear swath. This paper also examines a number of algorithmic and physical sensitivities of the new method and provides guidance for its cost-effective implementation. One obvious difficulty for the method, and other cloud data assimilation methods as well, is the lack of information content in the cloud observables on cloud vertical structure, beyond cloud top pressure and optical thickness, thus necessitating strong dependence on the background vertical moisture structure. It is found that a simple flow-dependent correlation modification due to Riishojgaard (1998) provides some help in this respect, by better honoring inversion structures in the background state.

Monte Carlo Bayesian inference↗

The nakhlite meteorites: Augite-rich igneous rocks from Mars

The seven nakhlite meteorites are augite-rich igneous rocks that formed in flows or shallow intrusions of basaltic magma on Mars. They consist of euhedral to subhedral crystals of augite and olivine (to 1 cm long) in fine-grained mesostases. The augite crystals have homogeneous cores of Mg' = 63% and rims that are normally zoned to iron enrichment. The core-rim zoning is cut by iron-enriched zones along fractures and is replaced locally by ferroan low-Ca pyroxene. The core compositions of the olivines vary inversely with the steepness of their rim zoning - sharp rim zoning goes with the most magnesian cores (Mg' = 42%), homogeneous olivines are the most ferroan. The olivine and augite crystals contain multiphase inclusions representing trapped magma. Among the olivine and augite crystals is mesostasis, composed principally of plagioclase and/or glass, with euhedra of titanomagnetite and many minor minerals. Olivine and mesostasis glass are partially replaced by veinlets and patches of iddingsite, a mixture of smectite clays, iron oxy-hydroxides and carbonate minerals. In the mesostasis are rare patches of a salt alteration assemblage: halite, siderite, and anhydrite/ gypsum. The nakhlites are little shocked, but have been affected chemically and biologically by their residence on Earth. Differences among the chemical compositions of the nakhlites can be ascribed mostly to different proportions of augite, olivine, and mesostasis. Compared to common basalts, they are rich in Ca, strongly depleted in Al, and enriched in magmaphile (incompatible) elements, including the LREE. Nakhlites contain little pre-terrestrial organic matter. Oxygen isotope ratios are not terrestrial, and are different in anhydrous silicates and in iddingsite. The alteration assemblages all have heavy oxygen and heavy carbon, while D/H values are extreme and scattered. Igneous sulfur had a solar-system isotopic ratio, but in most minerals was altered to higher and lower values. High precision analyses show mass-independent fractionations of S isotopes. Nitrogen and noble gases are complex and represent three components: two mantle sources (Chas-E and Chas-S), and fractionated Martian atmosphere. The nakhlites are igneous cumulate rocks, formed from basaltic magma at approx.1.3 Ga, containing excess crystals over what would form from pure magma. After accumulation of their augite and olivine crystals, they were affected (to various degrees) by crystallization of the magma, element diffusion among minerals and magma, chemical reactions among minerals and magma, magma movement among the crystals, and post-igneous chemical equilibration. The extent of these modifications varies, from least to greatest, in the order: MIL03346, NWA817, Y000593, Nakhla = Governador Valadares, Lafayette, and NWA998. Chemical, isotopic, and chronologic data confirm that the nakhlites formed on Mars, most likely in thick lava flows or shallow intrusions. Their crystallization ages, referenced to crater count chronologies for Mars, suggest that the nakhlites formed on the large volcanic constructs of Tharsis, Elysium, or Syrtis Major. The nakhlites were suffused with liquid water, probably at approx.620 ma. This water dissolved olivine and mesostasis glass, and deposited iddingsite and salt minerals in their places. The nakhlites were ejected from Mars at approx.10.75Ma by an asteroid impact and fell to Earth within the last 10,000 years. Although the nakhlites are enriched in incompatible elements, their source mantle was strongly depleted. This depletion event was ancient, as the nakhlites source mantle was fractionated while short-lived radionuclides (e.g., t(sub 1/2 = 9 my) were still active. This differentiation event may have been core formation coupled with a magma ocean, as is inferred for the moon.

Treiman, Allan H.↗

A cloudiness transition in a marine boundary layer

Boundary layer cloudiness plays several important roles in the energy budget of the earth. Low level stratocumulus are highly reflective clouds which reduce the net incoming shortwave radiation at the earth's surface. Climatically, the transition to a small area fraction of scattered cumulus clouds occurs as the air flows over warmer water. Although these clouds reflect less sunlight, they still play an important role in the boundary layer equilibrium by transporting water vapor upwards, and enhancing the surface evaporation. The First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment (FIRE) included a marine stratocumulus experiment off the southern California coast from June 29 to July 19, 1987. The objectives of this experiment were to study the controls on fractional cloudiness, and to assess the role of cloud-top entrainment instability (CTEI) and mesoscale structure in determining cloud type. The focus is one research day, July 7, 1987, when coordinated aircraft missions were flown by four research aircraft, centered on a LANDSAT scene at 1830 UTC. The remarkable feature of this LANDSAT scene is the transition from a clear sky in the west through broken cumulus to solid stratocumulus in the east. The dynamic and thermodynamic structure of this transition in cloudiness is analyzed using data from the NCAR Electra. By averaging the aircraft data, the internal structure of the different cloud regimes is documented, and it is shown that the transition between broken cumulus and stratocumulus is associated with a change in structure with respect to the CTEI condition. However, this results not from sea surface temperature changes, but mostly from a transition in the air above the inversion, and the breakup appears to be at a structure on the unstable side of the wet virtual adiabat.

Betts, Alan K.↗

Size distribution and scattering phase function of aerosol particles retrieved from sky brightness measurements

Ground-based measurements of the solar transmission and sky radiance in a horizontal plane through the Sun are taken in several geographical regions and aerosol types: dust in a desert transition zone in Israel, sulfate particles in Eastern and Western Europe, tropical aerosol in Brazil, and mixed continental/maritime aerosol in California. Stratospheric aerosol was introduced after the eruption of Mount Pinatubo in June 1991. Therefore measurements taken before the eruption are used to analyze the properties of tropospheric aerosol; measurements from 1992 are also used to detect the particle size and concentration of stratospheric aerosol. The measurements are used to retrieve the size distribution and the scattering phase function at large scattering angles of the undisturbed aerosol particles. The retrieved properties represent an average on the entire atmospheric column. A comparison between the retrieved phase function for a scattering angle of 120 deg, with phase function predicted from the retrieved size distribution, is used to test the assumption of particle homogeneity and sphericity in radiative transfer models (Mie theory). The effect was found to be small (20% +/- 15%). For the stratospheric aerosol (sulfates), as expected, the phase function was very well predicted using the Mie theory. A model with a power law distribution, based on the spectral dependence of the optical thickness, alpha, cannot estimate accurately the phase function (up to 50% error for lambda = 0.87 microns). Before the Pinatubo eruption the ratio between the volumes of sulfate and coarse particles was very well correlated with alpha. The Pinatubo stratospheric aerosol destroyed this correlation. The aerosol optical properties are compared with analysis of the size, shape, and composition of the individual particles by electron microscopy of in situ samples. The measured volume size distribution before the injection of stratospheric aerosol consistently show two modes, sulfate particles with r(sub m) less than 0.2 microns and coarse paritcles with r(sub m) greater than 0.7 microns. The 'window' in the tropospheric aerosol in this radius range was used to observe a stable stratospheric aerosol in 1992, with r(sub m) approximately 0.5 microns. A combination of such optical thickness and sky measurements can be used to assess the direct forcing and the climatic impact of aerosol. Systematic inversion for the key aerosol types (sulfates, smoke, dust, and maritime aerosol) of the size distribution and phase function can give the relationship between the aerosol physical and optical properties that can be used to compute the radiative forcing. This forcing can be validated in dedicated field experiments.

Kaufman, Y. J.↗

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗

One-loop matrix elements of effective superstring interactions: α' -expanding loop integrands

In the low-energy effective action of string theories, non-abelian gauge interactions and supergravity are augmented by infinite towers of higher-mass-dimension operators. We propose a new method to construct one-loop matrix elements with insertions of operators D 2k F n and D 2k R n in the tree-level effective action of type-I and type-II superstrings. Inspired by ambitwistor string theories, our method is based on forward limits of moduli-space integrals using string tree-level amplitudes with two extra points, expanded in powers of the inverse string tension α'. Similar to one-loop ambitwistor computations, intermediate steps feature non-standard linearized Feynman propagators which eventually recombine to conventional quadratic propagators. With linearized propagators the loop integrand of the matrix elements obey one-loop versions of the monodromy and KLT relations. We express a variety of four- and five-point examples in terms of quadratic propagators and formulate a criterion on the underlying genus-one correlation functions that should make this recombination possible at all orders in α'. The ultraviolet divergences of the one-loop matrix elements are crosschecked against the non-separating degeneration of genus-one integrals in string amplitudes. Conversely, our results can be used as a constructive method to determine degenerations of elliptic multiple zeta values and modular graph forms at arbitrary weight.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Earth System Modeling and Field Experiments in the Arctic-Boreal Zone - Report from a NASA Workshop

Early climate modeling studies predicted that the Arctic Ocean and surrounding circumpolar land masses would heat up earlier and faster than other parts of the planet as a result of greenhouse gas-induced climate change, augmented by the sea-ice albedo feedback effect. These predictions have been largely borne out by observations over the last thirty years. However, despite constant improvement, global climate models have greater difficulty in reproducing the current climate in the Arctic than elsewhere and the scatter between projections from different climate models is much larger in the Arctic than for other regions. Biogeochemical cycle (BGC) models indicate that the warming in the Arctic-Boreal Zone (ABZ) could lead to widespread thawing of the permafrost, along with massive releases of CO2 and CH4, and large-scale changes in the vegetation cover in the ABZ. However, the uncertainties associated with these BGC model predictions are even larger than those associated with the physical climate system models used to describe climate change. These deficiencies in climate and BGC models reflect, at least in part, an incomplete understanding of the Arctic climate system and can be related to inadequate observational data or analyses of existing data. A workshop was held at NASA/GSFC, May 22-24 2012, to assess the predictive capability of the models, prioritize the critical science questions; and make recommendations regarding new field experiments needed to improve model subcomponents. This presentation will summarize the findings and recommendations of the workshop, including the need for aircraft and flux tower measurements and extension of existing in-situ measurements to improve process modeling of both the physical climate and biogeochemical cycle systems. Studies should be directly linked to remote sensing investigations with a view to scaling up the improved process models to the Earth System Model scale. Data assimilation and observing system simulation studies should be used to guide the deployment pattern and schedule for inversion studies as well. Synthesis and integration of previously funded Arctic-Boreal projects (e.g., ABLE, BOREAS, ICESCAPE, ICEBRIDGE, ARCTAS) should also be undertaken. Such an effort would include the integration of multiple remotely sensed products from the EOS satellites and other resources.

Sellers, Piers↗

Imaging surface topography with coherent x-ray reflectivity: Theory, kinematics, and simulations

A theoretical formalism is described for understanding coherent x-ray reflectivity (CXR) from the surface of a semi-infinite crystal having a variable surface topography, described by the height profile ℎ(𝑥,𝑦). The surface topography is imaged as a complex “effective density,” obtained from the phasing and inversion of the coherent x-ray reflectivity data, measured through a rocking scan centered at a vertical momentum transfer 𝑄$^{0}_{𝑧}$ and a vertical range Δ⁢𝑄 𝑧 . The formalism predicts that the effective density has an amplitude with a maximum located at the surface height for each position within the surface plane. The phase of the effective density has a lateral variation that is controlled by the surface height and a vertical variation that reflects a combination of the interfacial structure and specific choice of measurement conditions. This understanding enables direct observation of nanometer-scale interfacial topography, i.e., ℎ⁡(𝑥,𝑦)⁢𝑐 𝑠 (where 𝑐 𝑠 is the vertical substrate lattice parameter) with Å-scale sensitivity to surface height. Numerical simulations illustrate and confirm the theoretical results. These results show how the interpretation of the interfacial density phase obtained by CXR data inversion (i.e., surface topography with respect to a flat surface) is conceptually similar to that previously known for Bragg coherent diffraction imaging (BCDI) measurements of isolated nanoparticles (i.e., lattice displacements with respect to an ideal crystal lattice). This suggests that CXR can be thought of as a form of dark field imaging with respect to the bright field BCDI approach. An implication of these results is that interfacial imaging may bypass some of the significant challenges associated with BCDI imaging of multiple particles having different orientations.

X-ray imaging↗

NETL Plastic Pipes Project (Final Report)

Plastic or composite pipelines have been the bane of the utility locating industry because they are neither conductive nor magnetic which are the properties traditionally used to locate buried utilities. Ground penetrating radar (GPR) is an effective geophysical tool for locating plastic/composite pipelines where resistive cover allows for adequate penetration of radar energy. However, GPR has limited applicability in areas where the soil cover is conductive due to significant clay and/or salt content. This study examines complementary near-surface geophysical methods that are potentially useful for locating buried plastic/composite pipelines, either singly or in combination. Specifically, this modeling study used computational numerical methods to forward model the response of GPR, resistivity, seismic, gravity gradiometry, and photoacoustic/thermoacoustic imaging methods to plastic/composite pipelines for various scenarios including: (1) pipe diameters ranging between 2 in. to 12 in.; (2) burial depths ranging between 3 ft. to 4 ft.; (3) various degrees in contrast in physical properties (i.e., electrical permittivity, elasticity, resistivity, density); and (4) various experimental acquisition choices (e.g., GPR radar and seismic source frequencies, electrode spacing). Numerical modeling performed herein reconfirmed that GPR is the preferred method for detecting/locating plastic pipelines. A caveat for GPR detection is that the material covering the plastic pipe (trench fill material and adjacent soil) must be sufficiently resistive to allow the two-way propagation to the required depth of investigation and back to the surface. GPR was the only method modeled in this study that can be used to directly detect plastic pipelines of 2-in.-diameter and larger when buried 3-ft-deep. GPR data processing and imaging also can determine pipe depth, pipe diameter, trench dimensions, and moisture conditions. Seismic modeling results suggest that direct detection of a 12-in.-diameter plastic pipe at 3-ft.-depth may be possible under favorable conditions; however, the associated signature would be weak (e.g., surface- to S-wave, backscattered surface-waves, and/or forward scattered surface-waves to S-wave). Direct pipe detection under field conditions with noise and strong lateral geologic heterogeneity is doubtful. Numerical modeling also suggests that plastic pipelines can be indirectly located by detecting the trench in which they are buried. GPR, direct current (DC) resistivity, and seismic methods have the potential to locate the pipeline trench if there is sufficient contrast between the trench-wall and trench-fill materials for the physical property being measured by each method (i.e., electrical permittivity for GPR; resistivity for DC resistivity; or density, compressional velocity, or shear velocity for seismic). Modeling also indicated that currently available (commercial) gravity gradiometers would be unable to directly detect/locate plastic pipelines ≤ 8-in.-diameter when buried 3-ft.-deep given the typical instrument noise floor for field surveying as well as the expected density variations due to geologic heterogeneity. The numerical modeling performed in this project did not identify a universal geophysical technology that can locate buried plastic pipelines in all parts of the United States (although GPR is suggested for all areas with resistive cover). However, the project results suggest that a towed land streamer simultaneously acquiring multiple geophysical data types including multi-offset GPR, multi-channel DC resistivity, seismic geophone- and/or distributed acoustic sensing (DAS), and potentially photoacoustic/thermoacoustic data would be an appropriate platform for locating buried plastic pipeline. Moreover, the complementary multiphysics data acquired by a towed land streamer would permit the use of joint and/or cooperative inversion frameworks for a more rigorous and consistent data interpretation.

42 ENGINEERING↗

Sea Surface Slope Statistics for Intermediate and Shore Scale Ocean Waves Measured Using a Low-Altitude Aircraft

Ocean surface remote sensing techniques often rely on scattering or emission linked to shorter- scale gravity-capillary ocean wavelets. However, it is increasingly apparent that slightly longer wavelengths of O(10 to 500 cm) are vital components in the robust sea surface description needed to link varied global remote sensing data sets. This paper describes a sensor suite developed to examine sea surface slope variations in the field using an aircraft flying at very low altitude (below 30 m) and will also provide preliminary measurements detailing changes in slope characteristics versus sea state and friction velocity. Two-dimensional surface slope is measured using simultaneous range measurements from three compact short-range laser altimeters mounted in an equilateral triangle arrangement with spacing of about 1 m. In addition, all three lasers provide independent wave elevation profiles after GPS-aided correction for aircraft altitude. Laser range precision is 1 cm rms while vertical motion correction is 15 cm rms. The measurements are made along-track at approximately 1 m intervals setting the spatial scale of the measurement to cover waves of intermediate to long scale. Products available for this array then include surface elevation, two-dimensional slope distribution, and the cross- and along-track 1-D slope distributions. To complement the laser, a down-looking mm-wave radar scatterometer is centered within the laser array to measure radar backscatter simultaneously with the laser slope. The radar's footprint is nominally 1 m in diameter. Near-vertical radar backscatter is inversely proportional to the small-scale surface slope variance and to the tilt of the underlying (laser-measured) surface facet. Together the laser and radar data provide information on wave roughness from the longest scales down to about 1 cm. These measurements are complemented by aircraft turbulence probe data that provides robust surface flux information.

Vandemack, Douglas↗

Characterizing a New Surface-Based Shortwave Cloud Retrieval Technique, Based on Transmitted Radiance for Soil and Vegetated Surface Types

This paper presents an approach using the GEneralized Nonlinear Retrieval Analysis (GENRA) tool and general inverse theory diagnostics including the maximum likelihood solution and the Shannon information content to investigate the performance of a new spectral technique for the retrieval of cloud optical properties from surface based transmittance measurements. The cumulative retrieval information over broad ranges in cloud optical thickness (tau), droplet effective radius (r(sub e)), and overhead sun angles is quantified under two conditions known to impact transmitted radiation; the variability in land surface albedo and atmospheric water vapor content. Our conclusions are: (1) the retrieved cloud properties are more sensitive to the natural variability in land surface albedo than to water vapor content; (2) the new spectral technique is more accurate (but still imprecise) than a standard approach, in particular for tau between 5 and 60 and r(sub e) less than approximately 20 nm; and (3) the retrieved cloud properties are dependent on sun angle for clouds of tau from 5 to 10 and r(sub e) less than 10 nm, with maximum sensitivity obtained for an overhead sun.

Shannon information content↗