Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Power Law”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Joint ROSAT-Compton GRO observations of the X-ray bright Seyfert galaxy IC 4329A

We report a simultaneous ROSAT and Gamma Ray Observatory (GRO) observation of the X-ray-bright Seyfert galaxy IC 4329A. For the GRO Oriented Scintillation Spectrometer Experiment (OSSE) detector, we also present the sum of the data for this and earlier observations. The overall spectrum is very well described as a power law with an energy spectral index of approximately 1 absorbed at low energies plus a strong Compton reflection component, typical for Seyfert 1 galaxies. The low energy absorption can be well described by a sum of a neutral column density of approximately 3 x 10(exp 21)sq cm, most of which is associated with the edge-on galactic disk of IC 4329A, plus an edgelike feature at approximately 700 eV; this feature implies either complex absorption (due to additional ionized material, or due to a partial covering), or a soft excess. The data only weakly constrain the presence of a high-energy cutoff in the underlying power law; they are compatible with an exponential cutoff at any energy E(sub c) approximately greater than 100 keV. The relative steepness of the OSSE data, with the power-law energy index of 1.6 +/- 0.2, can be accounted for entirely by the contribution of the high-energy tail of the reflection component when E(sub c) approaches infinity. (We find that the definite cutoff at an energy E(sub c)approximately 130 keV suggested in the recently published analysis of the OSSE data for this subject is due to a data reduction error.) Including nonsimultaneous Ginga observations with 2 keV fluxes matching well that of ROSAT gives us likely broad-band X-ray/gamma-ray spectra of the object from approximately 0.1 keV up to several hundred keV. We also report the ROSAT spectrum of the companion object to the Seyfert galaxy, the elliptical galaxy IC 4329.

Madejski, G. M.↗
DOE OSTI2022

Comparison of temperature and doping dependence of elastoresistivity near a putative nematic quantum critical point

Abstract Strong electronic nematic fluctuations have been discovered near optimal doping for several families of Fe-based superconductors, motivating the search for a possible link between these fluctuations, nematic quantum criticality, and high temperature superconductivity. Here we probe a key prediction of quantum criticality, namely power-law dependence of the associated nematic susceptibility as a function of composition and temperature approaching the compositionally tuned putative quantum critical point. To probe the ‘bare’ quantum critical point requires suppression of the superconducting state, which we achieve by using large magnetic fields, up to 45 T, while performing elastoresistivity measurements to follow the nematic susceptibility. We performed these measurements for the prototypical electron-doped pnictide, Ba(Fe 1− x Co x ) 2 As 2 , over a dense comb of dopings. We find that close to the putative quantum critical point, the elastoresistivity appears to obey power-law behavior as a function of composition over almost a decade of variation in composition. Paradoxically, however, we also find that the temperature dependence for compositions close to the critical value cannot be described by a single power law.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1988

On the spectrum of the secondary auroral electrons

Results are presented from a statistical study of the spectral characteristics of the (power-law-portion) secondary electrons associated with inverted-V electron precipitation events (IVEs). A total of 106 IVEs observed at low altitudes (below 1000 km) by the DE 2 spacecraft are included in data base. Spectral parameters for the electrons between 20 eV and 100 eV are plotted against the different parameters pertaining to the locations of the IVEs and the primary beams to determine their empirical dependences. It is found that the differential spectrum of the secondary electrons can be modeled by a power law, J proportional to E exp -gamma, with an averaged gamma of 1.85. This spectral parameter is relatively insensitive to the variation of locations of the IVEs, but it depends on the primary beam parameters. These results are compared with predictions, and theories on the generation of the low-energy auroral power-law secondary electrons are assessed. The analysis shows that collisional processses are the dominant source of the secondary electrons.

Fung, S. F.↗
NASA NTRS1994

Thin shell formation in radiative shocks. 1: Supernova remnants in low-density media

This paper explores the onset of thin-shell formation in interstellar shocks associated with supernova explosions. We outline a simple but useful scheme that indicates the time at which thin shell formation begins for supernova remnants (SNRs) evolving in a range of interstellar environments, extending the previous analytical models to arbitrary power-law density media. The result depends on the gas cooling properties and the shock velocity and radius. This is then applied to the specific case of SNRs in low-density media. The procedure for defining the time for the onset of shell formation, t(sub sf), equates the value of the adiabat, kappa = p/rho(exp gamma), to zero using the known time dependence of the shock radius and velocity. For the case of a power-law density ambient medium of the form rho(r) = Br(exp -omega), it is found that shell formation can be prevented when the ambient density drops faster than a critical rate. For a cooling function of the form Lambda = Lambda(sub 0) tau(exp beta), with beta = -0.5 (appropriate for line cooling), shell formation never occurs for omega greater than or equal to 9/5. The shell formation time is then computed for spherical shocks in a power-law density medium. For omega = 0, the onset of shell formation is found to be at t(sub sf) approx. equal to 2.87 x 10(exp 4) E51(exp 3/14) n(sub 0 exp -4.7) yr, which agrees well with previous estimates derived by other means. We compare the analytical shell formation time with the results of detailed numerical models for omega = 0 and three different ambient densities and find good agreement. The extension of the criterion for the onset of thin shell formation using the ratio of cooling to swept-up column density is also described. This method provides a useful approximation for cases when the exact solution is not known.

Franco, Jose↗
NASA NTRS1994

The cosmological dependence of cluster density profiles

We use N-body simulations to study the shape of mean cluster density and velocity profiles in the nonlinear regime formed via gravitational instability. The dependence of the final structure on both cosmology and initial density field is examined, using a grid of cosmologies and scale-free initial power spectra P(k) varies as k(exp n). Einstein-de Sitter, open (Omega(sub 0) = 0.2 and 0.1) and flat, low density (Omega(sub 0) = 0.2 lambda(sub 0) = 0.8) models are examined, with initial spectral indices n = -2, -1 and 0. For each model, we stack clusters in an appropriately scaled manner to define an average density profile in the nonlinear regime. The profiles are well fit by a power law rho(r) varies as r(exp -alpha) for radii whereat the local density contrast is between 100 and 3000. This covers 99% of the cluster volume. We find a clear trend toward steeper slopes (larger alphas) with both increasing n and decreasing Omega(sub 0). The Omega(sub 0) dependence is partially masked by the n dependence; there is degeneracy in the values of alpha between the Einstein-de Sitter and flat, low-density cosmologies. However, the profile slopes in the open models are consistently higher than the Omega = 1 values for the range of n examined. Cluster density profiles are thus potentially useful cosmological diagnostics. We find no evidence for a constant density core in any of the models, although the density profiles do tend to flatten at small radii. Much of the flattening is due to the force softening required by the simulations. An attempt is made to recover the unsoftened profiles assuming angular momentum invariance. The recovered profiles in Einstein-de Sitter cosmologies are consistent with a pure power law up to the highest density contrasts (10(exp 6)) accessible with our resolution. The low-density models show significant deviation from a power law above density contrasts approximately 10(exp 5). We interpret this curvature as reflecting the non-scale-invariant nature of the background cosmology in these models. These results are at the limit of our resolution and so should be tested in the future using simulations with larger numbers of particles. Such simulations will also provide insight on the broader problem of understanding, in a statistical sense, the full phase space structure of collapsed, cosmological halos.

Crone, Mary M.↗
DOE OSTI2023

Kibble–Zurek mechanism of Ising domains

The formation of topological defects after a symmetry-breaking phase transition is an overarching phenomenon that encodes the underlying dynamics. The Kibble-Zurek mechanism (KZM) describes these non-equilibrium dynamics of second-order phase transitions and predicts a power-law relationship between the cooling rates and the density of topological defects. It has been verified as a successful model in a wide variety of physical systems, including structure formation in the early Universe and condensed-matter materials. However, it is uncertain if the KZM mechanism is also valid for topologically trivial Ising domains, one of the most common and fundamental types of domain in condensed-matter systems. Here we show that the cooling rate dependence of Ising domain density follows the KZM power law in two different three-dimensional structural Ising domains: ferro-rotation domains in NiTiO 3 and polar domains in BiTeI. However, although the KZM slope of NiTiO 3 agrees with the prediction of the 3D Ising model, the KZM slope of BiTeI exceeds the theoretical limit, providing an example of steepening KZM slope with long-range dipolar interactions. Finally, our results demonstrate the validity of KZM for Ising domains and reveal an enhancement of the power-law exponent for transitions of non-topological quantities with long-range interactions. The Kibble-Zurek mechanism is shown to apply to structural Ising domains in three-dimensional materials. Long-range interactions modify the critical exponents away from theoretical predictions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Theory of Trotter Error with Commutator Scaling

The Lie-Trotter formula, together with its higher-order generalizations, provides a simple approach to decomposing the exponential of a sum of operators. Despite significant effort, the error scaling of such product formulas remains poorly understood. We develop a theory of Trotter error that overcomes the limitations of truncating the Baker-Campbell-Hausdorff expansion. Our analysis directly exploits the commutativity of operator summands, producing tighter error bounds for both real- and imaginary-time evolutions. Whereas previous work achieves similar goals for systems with geometric locality or Lie-algebraic structure, our approach holds in general. We give a host of improved algorithms for digital quantum simulation and quantum Monte Carlo methods, nearly matching or even outperforming the best previous results. Our applications include: (i) a simulation of second-quantized plane-wave electronic structure, nearly matching the interaction-picture algorithm of Low and Wiebe; (ii) a simulation of $k$-local Hamiltonians almost with induced one-norm scaling, faster than the qubitization algorithm of Low and Chuang; (iii) a simulation of rapidly decaying power-law interactions, outperforming the Lieb-Robinson-based approach of Tran et al.; (iv) a hybrid simulation of clustered Hamiltonians, dramatically improving the result of Peng, Harrow, Ozols, and Wu; and (v) quantum Monte Carlo simulations of the transverse field Ising model and quantum ferromagnets, tightening previous analyses of Bravyi and Gosset. We obtain further speedups using the fact that product formulas can preserve the locality of the simulated system. Specifically, we show that local observables can be simulated with complexity independent of the system size for power-law interacting systems, which implies a Lieb-Robinson bound nearly matching a recent result of Tran et al. Our analysis reproduces known tight bounds for first- and second-order formulas. We further investigate the tightness of our bounds for higher-order formulas. For quantum simulation of a one-dimensional Heisenberg model with an even-odd ordering of terms, our result overestimates the complexity by only a factor of $5$. Our bound is also close to tight for power-law interactions and other orderings of terms. This suggests that our theory can accurately characterize Trotter error in terms of both the asymptotic scaling and the constant prefactor.

quantum computing, numerical analysis↗
DOE OSTI2022

Interpreting radial correlation Doppler reflectometry using gyrokinetic simulations

A linear response, local model for the DBS amplitude applied to gyrokinetic simulations shows that radial correlation Doppler reflectometry measurements (RCDR, Schirmer et al 2007 Plasma Phys. Control. Fusion 49 1019) are not sensitive to the average turbulence radial correlation length, but to a correlation length that depends on the binormal wavenumber $k_{\perp}$ selected by the Doppler backscattering (DBS) signal. Nonlinear gyrokinetic simulations show that the turbulence naturally exhibits a nonseparable power law spectrum in wavenumber space, leading to a power law dependence of the radial correlation length with binormal wavenumber $l_r \sim C k_{\perp}^{-\alpha} (\alpha \approx 1)$ which agrees with the inverse proportionality relationship between the measured $l_r$ and $k_{\perp} $ observed in experiments (Fernández-Marina et al 2014 Nucl. Fusion 54 072001). This new insight indicates that RCDR characterizes the eddy aspect ratio in the perpendicular plane to the magnetic field. It also motivates future use of a nonseparable turbulent spectrum to quantitatively interpret RCDR and potentially other turbulence diagnostics. The radial correlation length is only measurable when the radial resolution at the cutoff location $W_\text n$ satisfies $W_\text n \ll l_r$, while the measurement becomes dominated by $W_\text n$ for $W_\text n \gg l_r$. This suggests that $l_r$ is likely to be inaccessible for electron-scale DBS measurements (${k_{\perp}}{\rho_{\text s}} \gt 1$). The effect of $W_\text n$ on ion-scale radial correlation lengths could be nonnegligible.

radial correlation Doppler reflectometry↗
DOE OSTI2022

Crustal Strain Rates in the Western United States and Their Relationship with Earthquake Rates

Abstract We present a suite of strain rate models for the western United States based on geologic and geodetic data. The geologic data consist of Quaternary fault-slip rates and the geodetic data consists of a new compilation of Global Positioning System (GPS) velocities derived from continuous, semicontinuous, and campaign measurements. We remove postseismic deformation from the GPS time series in order for our geodetic strain rate model to best capture the interseismic strain accumulation rate. We present models based on either geologic or geodetic data, but also create a hybrid model. Although there are some differences between the models, the large-scale features are the same, with the noticeable exception for the Pacific Northwest where interseismic strain is naturally more distributed than the long-term strain release. We also present a map of earthquake rate densities based on mainshocks, and the result has similar spatial features similar to the strain rate models (at least in the southwestern United States). We perform a general correlation analysis between strain rate and seismicity rate (south of Cascadia) and find a change in linearity between seismicity and strain rates from slow to faster deforming areas with seismicity rates relatively lower for the latter. The extent of that change depends a bit on assumptions made on the declustering and completeness of the catalog, but the finding of a change in slope is robust across the different strain rate models. Linearity for all areas is only expected when Gutenberg–Richter parameters and parameters involved in the conversion from strain to moment rate are uniform across the study area. We discuss these qualifications, but find no single satisfactory explanation for our observation. Moreover, when considering a rather short time and space, theoretical considerations of sampling from a power-law distribution actually predict there to be a power law instead of a linear relationship, generally consistent with our observation.

Geochemistry & Geophysics↗
DOE OSTI2025

X-Ray Polarization of the Magnetar 1E 1841−045

We report on IXPE and NuSTAR observations beginning 40 days after the 2024 outburst onset of magnetar 1E 1841−045, marking the first IXPE observation of a magnetar in an enhanced state. Our spectropolarimetric analysis indicates that both a blackbody (BB) plus double power-law (PL) and a double blackbody plus power-law spectral model fit the phase-averaged intensity data well, with a hard PL tail (Γ = 1.19 and 1.35, respectively) dominating above ≈5 keV. For the former model, we find the soft PL (the dominant component at soft energies) exhibits a polarization degree (PD) of ≈30% while the hard PL displays a PD of ≈40%. Similarly, the cool BB of the 2BB+PL model possesses a PD of ≈15% and a hard PL PD of ≈57%. For both models, each component has a polarization angle (PA) compatible with celestial north. Model-independent polarization analysis supports these results, wherein the PD increases from ≈15% to ≈70% in the 2–3 keV and 6–8 keV ranges, respectively, while the PA remains nearly constant. We find marginal evidence for phase-dependent variability of the polarization properties, namely a higher PD at phases coinciding with the hard X-ray pulse peak. We compare the hard X-ray PL to the expectation from resonant inverse Compton scattering (RICS) and secondary pair cascade synchrotron radiation from primary high-energy RICS photons; both present reasonable spectropolarimetric agreement with the data, albeit the latter does so more naturally. We suggest that the soft PL X-ray component may originate from a Comptonized corona in the inner magnetosphere.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The contribution of Humboldt Glacier, northern Greenland, to sea-level rise through 2100 constrained by recent observations of speedup and retreat

Abstract. Humboldt Glacier, northern Greenland, has retreated and accelerated through the 21st century, raising concerns that it could be a significant contributor to future sea-level rise. We use a data-constrained ensemble of three-dimensional higher-order ice sheet model simulations to estimate the likely range of sea-level rise from the continued retreat of Humboldt Glacier. We first solve for basal traction using observed ice thickness, bed topography, and ice surface velocity from the year 2007 in a PDE-constrained (partial differential equation) optimization. Next, we impose calving rates to match mean observed retreat rates from winter 2007–2008 to winter 2017–2018 in a transient calibration of the exponent in the power-law basal friction relationship. We find that power-law exponents in the range of 1/7–1/5 – rather than the commonly used 1/3–1 – are necessary to reproduce the observed speedup over this period. We then tune an iceberg calving parameterization based on the von Mises stress yield criterion in another transient-calibration step to approximate both observed ice velocities and terminus position in 2017–2018. Finally, we use the range of basal friction relationship exponents and calving parameter values to generate the ensemble of model simulations from 2007–2100 under three climate forcing scenarios from CMIP5 (two RCP8.5 forcings, Representative Concentration Pathway) and CMIP6 (one SSP5-8.5 forcing, Shared Socioeconomic Pathway). Our simulations predict 5.2–8.7 mm of sea-level rise from Humboldt Glacier, significantly higher than a previous estimate (∼ 3.5 mm) and equivalent to a substantial fraction of the 40–140 mm predicted by ISMIP6 from the whole Greenland Ice Sheet. Our larger future sea-level rise prediction results from the transient calibration of our basal friction law to match the observed speedup, which requires a semi-plastic bed rheology. In many simulations, our model predicts the growth of a sizable ice shelf in the middle of the 21st century. Thus, atmospheric warming could lead to more retreat than predicted here if increased surface melt promotes hydrofracture of the ice shelf. Our data-constrained simulations of Humboldt Glacier underscore the sensitivity of model predictions of Greenland outlet glacier response to warming to choices of basal shear stress and iceberg calving parameterizations. Further, transient calibration of these parameterizations, which has not typically been performed, is necessary to reproduce observed behavior. Current estimates of future sea-level rise from the Greenland Ice Sheet could, therefore, contain significant biases.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1973

Celestial diffuse gamma radiation above 30 MeV observed by SAS-2

The Small Astronomy Satellite (SAS)-2, launched on November 15, 1972, carried into orbit a 32-deck magnetic-core digitized spark chamber gamma ray telescope to study celestial gamma radiation in the energy range above 30 MeV. In the study of several regions with b sub 2 15 deg, a finite, diffuse flux of gamma rays with a steep energy spectrum in the energy region from 35 to 200 MeV is observed. Representing the energy spectrum by a power law of the form dJ/dE = AE to - alpha power over this energy range, alpha is found along with the integral flux above 100 MeV. Combining this result with existing low energy gamma ray data yields an energy spectrum which is not a simple power law in energy, as in the X-ray region, but which demonstrates first an increase and then a decrease in slope, consistent within uncertainties with that predicted by cosmological theories, including the continuous production of high energy gamma rays primarily from neutral pi mesons throughout the history of the universe.

Fichtel, C. E.↗
NASA NTRS1981

Temperature-dependent creep buckling of plates

Time-dependent lateral deflection of flat rectangular plates is predicted by the Norton-Bailey (Norton 1929, Bailey 1935) power law for material creep. The plates have a through-thickness steady-state temperature distribution, and the effects are considered by using Maxwell's law to modify the power creep law. Equations are derived for creep exponents of 3 and 5, using the sandwich plate element to predict creep buckling of plates. Predictions of creep buckling with a temperature variation between the inner and outer plate surfaces are found to be somewhat dependent on the creep buckling relationship assumed. When significant scatter justifies a variation in the creep constants up to an order of magnitude, discrepancies in predictions using the two exponents are reasonable, and for one engineering material, the predictions have the same degree of agreement with experimental data as have the respective creep laws.

Ross, D. A.↗
NASA NTRS1986

The magnetosphere of Uranus - Hot plasma and radiation environment

Inferences are drawn on the morphology and composition of the Uranus magnetosphere based on low-energy charged particle data collected by Voyager 2. Proton and electron energies in the magnetosphere attained energies of 4 and 1.2 MeV, respectively, although electron intensities surpassed the proton intensities at most energy levels. Protons dominated in the ion energy regime 0.6-1.0 MeV. The ion and electron spectra were Maxwellian below about 200 keV and had a power law distribution at energies over 590 keV. The power law was reduced by a factor of nearly three inside the orbit of Miranda. The proton population is dense enough to polymerize CO and CH4 ice surfaces within 10,000-100,000 yr. The data indicated that the particles are swept out at least to the orbit of Titania by the satellites. The morphology of the magnetosphere closely resembles that around Jupiter, except that plasma sheet distorsion from particle loading is negligible in regions within 15 Uranus radii.

Krimigis, S. M.↗
NASA NTRS1997

RXTE Observation of Cygnus X-1: Spectra and Timing

We present first results from the analysis of an RXTE observation of Cyg X-1 in its low state, taken about two months after the end of the high state. With Gamma approx. equal to 1.45 the spectrum is considerably harder than previous low-state measurements. The observed spectrum can be explained by a Comptonization spectrum as that emitted from a spherical corona surrounded by a cold accretion disk. The optical depth of the corona is between 2 and 2.5 and the temperature is between 60 and 80 keV. Temporal analysis shows a typical Root Mean Square (RMS) noise of approximately 25%. The Pulse Shape Discrimination (PSD) can be described as consisting of a flat component followed by an 1/f power-law, followed by an f(sup -1.6) power-law. The lag of the hard photons with respect to the soft photons is consistent with prior observations. The coherence function is remarkably close to unity from 0.01 Hz to 10 Hz.

Wilms, J.↗
NASA NTRS2001

Models for X-Ray Emission from Isolated Pulsars

A model is proposed for the observed combination of power-law and thermal X-rays from rotationally powered pulsars. For gamma-ray pulsars with accelerators very many stellar radii above the neutron star surface, 100 MeV curvature gamma-rays from e(-) or e(+) flowing starward out of such accelerators are converted to e1 pairs on closed field lines all around the star. These pairs strongly affect X-ray emission from near the star in two ways. (1) The pairs are a source of synchrotron emission immediately following their creation in regions where B approx. 10(exp 10) G. This emission, in the photon energy range 0.1 keV less than E(sub X) less than 5 MeV, has a power-law spectrum with energy index 0.5 and X-ray luminosity that depends on the back-flow current, and is typically approx. 10(exp 33) ergs/ s. (2) The pairs ultimately a cyclotron resonance "blanket" surrounding the star except for two holes along the open field line bundles which pass through it. In such a blanket the gravitational pull on e(+,-) pairs toward the star is balanced by the hugely amplified push of outflowing surface emitted X-rays wherever cyclotron resonance occurs. Because of it the neutron star is surrounded by a leaky "hohlraum" of hot blackbody radiation with two small holes, which prevents direct X-ray observation of a heated polar cap of a gamma-ray pulsar. Weakly spin modulated radiation from the blanket together with more strongly spin-modulated radiation from the holes through it would then dominate observed low energy (0.1-10 keV) emission. For non-y-ray pulsars, in which no such accelerators with their accompanying extreme relativistic back-flow toward the star are expected, optically thick e1 resonance blankets should not form (except in special cases very close to the open field line bundle). From such pulsars blackbody radiation from both the warm stellar surface and the heated polar caps should be directly observable. In these pulsars, details of the surface magnetic field evolution, especially of polar cap areas, become relevant to observations. The models are compared to X-ray data from Geminga, PSR 1055-52, PSR 0656+14, PSR 1929+10, and PSR 0950+08.

Wang, F. Y.-H.↗
NASA NTRS1976

Detection of nuclear gamma rays from Centaurus A

Results are reported for an observation of nuclear gamma rays in the energy range between 0.033 and 12.25 MeV from Centaurus A using a balloon-borne actively collimated NaI(Tl) crystal scintillation counter. The observing procedure is outlined, no systematic errors are found in the data, and power-law fits to the source's energy spectrum are attempted. A power law of approximately 0.86E to the -1.9 power photon/sq cm/sec per keV is shown to give an acceptable fit to the continuum, and the detection of two gamma-ray lines at 1.6 and 4.5 MeV, respectively, is discussed. It is found that the low-energy gamma-ray luminosity of Cen A is 9.4 by 10 to the 43rd power erg/sec for a distance of 5 Mpc and that Cen A is apparently variable in low-energy gamma radiation. It is suggested that the broad feature detected at 1.6 MeV may be due to three blended lines (possibly excited Ne-20, Mg-24, and Si-28), the 4.5-MeV line is most likely due to deexcitation of excited C-12, and the nuclear excitation results from either cosmic-ray bombardment of Cen A's interstellar medium or nucleosynthesis within the source.

Hall, R. D.↗
NASA NTRS2014

Use of 2d-video Disdrometer to Derive Mean Density-size and Ze-SR Relations: Four Snow Cases from the Light Precipitation Validation Experiment

The application of the 2D-video disdrometer to measure fall speed and snow size distribution and to derive liquid equivalent snow rate, mean density-size and reflectivity-snow rate power law is described. Inversion of the methodology proposed by Böhm provides the pathway to use measured fall speed, area ratio and '3D' size measurement to estimate the mass of each particle. Four snow cases from the Light Precipitation Validation Experiment are analyzed with supporting data from other instruments such as Precipitation Occurrence Sensor System (POSS), Snow Video Imager (SVI), a network of seven snow gauges and three scanning C9 band radars. The radar-based snow accumulations using the 2DVD-derived Ze-SR relation are in good agreement with a network of seven snow gauges and outperform the accumulations derived from a climatological Ze-SR relation used by the Finnish Meteorological Institute (FMI). The normalized bias between radar-derived and gauge accumulation is reduced from 96% when using the fixed FMI relation to 28% when using the Ze-SR relations based on 2DVD data. The normalized standard error is also reduced significantly from 66% to 31%. For two of the days with widely different coefficients of the Ze-SR power law, the reflectivity structure showed significant differences in spatial variability. Liquid water path estimates from radiometric data also showed significant differences between the two cases. Examination of SVI particle images at the measurement site corroborated these differences in terms of unrimed versus rimed snow particles. The findings reported herein support the application of Böhm's methodology for deriving the mean density-size and Ze-SR power laws using data from 2D-video disdrometer.

radar↗