Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical error”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-P-36: Evaluation of MLC Positioning Errors in Dynamic IMRT Treatments by Analyzing Dynalog Files

Purpose: To assess the accuracy of MLC positioning in Varian linear accelerator, in dynamic IMRT technique, from the analysis of dynalog files generated by the MLC controller. Methods: In Clinac accelerators (pre-TrueBeam technology), control system has an approximately 50ms delay (one control cycle time). Then, the system compares the measured position to the planned position corresponding to the next control cycle. As it has been confirmed by Varian technical support, this effect causes that measured positions appear in dynalogs one cycle out of phase with respect to the planned positions. Around 9000 dynalogs have been analyzed, coming from the three linear accelerators of our center (one Trilogy and two Clinac 21EX) equipped with a Millennium 120 MLC. In order to compare our results to recent publications, leaf positioning errors (RMS and 95th percentile) are calculated with and without delay effect. Dynalogs have been analyzed using a in-house Matlab software. Results: The RMS errors were 0.341, 0.339 and 0.348mm for each Linac; being the average error 0.343 mm. The 95th percentiles of the error were 0.617, 0.607 and 0.625; with an average of 0.617mm. A recent multi-institution study carried out by Kerns et al. found a mean leaf RMS error of 0.32mm and a 95th percentile error value of 0.64mm.Without delay effect, mean leaf RMS errors obtained were 0.040, 0.042 and 0.038mm for each treatment machine; being the average 0.040mm. The 95th percentile error values obtained were 0.057, 0.058 and 0.054 mm, with an average of 0.056mm. Conclusion: Results obtained for the mean leaf RMS error and the mean 95th percentile were consistent with the multi-institution study. Calculated error statistics with delay effect are significantly larger due to the speed proportional and systematic leaf offset. Consequently it is proposed to correct this effect in dynalogs analysis to determine the MLC performance.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Stochastic embedding DFT: Theory and application to p -nitroaniline in water

Over this past decade, we combined the idea of stochastic resolution of identity with a variety of electronic structure methods. In our stochastic Kohn-Sham density functional theory (DFT) method, the density is an average over multiple stochastic samples, with stochastic errors that decrease as the inverse square root of the number of sampling orbitals. Here in this paper, we develop a stochastic embedding density functional theory method (se-DFT) that selectively reduces the stochastic error (specifically on the forces) for a selected subsystem(s). The motivation, similar to that of other quantum embedding methods, is that for many systems of practical interest, the properties are often determined by only a small subsystem. In stochastic embedding DFT, two sets of orbitals are used: a deterministic one associated with the embedded subspace and the rest, which is described by a stochastic set. The method agrees exactly with deterministic calculations in the limit of a large number of stochastic samples. We apply se-DFT to study a p-nitroaniline molecule in water, where the statistical errors in the forces on the system (the p-nitroaniline molecule) are reduced by an order of magnitude compared with nonembedding stochastic DFT.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

An Update on Monitoring Stellar Orbits in the Galactic Center

Using 25 years of data from uninterrupted monitoring of stellar orbits in the Galactic Center, we present an update of the main results from this unique data set: a measurement of mass and distance to Sgr A*. Our progress is not only due to the eight-year increase in time base, but also to the improved definition of the coordinate system. The star S2 continues to yield the best constraints on the mass of and distance to Sgr A*; the statistical errors of 0.13×10{sup 6} M{sub ⊙} and 0.12 kpc have halved compared to the previous study. The S2 orbit fit is robust and does not need any prior information. Using coordinate system priors, the star S1 also yields tight constraints on mass and distance. For a combined orbit fit, we use 17 stars, which yields our current best estimates for mass and distance: M=4.28±0.10|{sub stat.}±0.21|{sub sys}×10{sup 6} M{sub ⊙} and R{sub 0}=8.32±0.07|{sub stat.}±0.14|{sub sys} kpc. These numbers are in agreement with the recent determination of R {sub 0} from the statistical cluster parallax. The positions of the mass, of the near-infrared flares from Sgr A*, and of the radio source Sgr A* agree to within 1 mas. In total, we have determined orbits for 40 stars so far, a sample which consists of 32 stars with randomly oriented orbits and a thermal eccentricity distribution, plus eight stars that we can explicitly show are members of the clockwise disk of young stars, and which have lower-eccentricity orbits.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Range-separated stochastic resolution of identity: Formulation and application to second-order Green’s function theory

We develop a range-separated stochastic resolution of identity (RS-SRI) approach for the four-index electron repulsion integrals, where the larger terms (above a predefined threshold) are treated using a deterministic RI and the remaining terms are treated using a SRI. The approach is implemented within a second-order Green’s function formalism with an improved O(N 3 ) scaling with the size of the basis set, N. Moreover, the RS approach greatly reduces the statistical error compared to the full stochastic version [T. Y. Takeshita et al., J. Chem. Phys. 151, 044114 (2019)], resulting in computational speedups of ground and excited state energies of nearly two orders of magnitude, as demonstrated for hydrogen dimer chains and water clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Impact of transport model errors on the global and regional methane emissions estimated by inverse modelling

A modelling experiment has been conceived to assess the impact of transport model errors on methane emissions estimated in an atmospheric inversion system. Synthetic methane observations, obtained from 10 different model outputs from the international TransCom-CH4 model inter-comparison exercise, are combined with a prior scenario of methane emissions and sinks, and integrated into the three-component PYVAR-LMDZ-SACS (PYthon VARiational-Laboratoire de Météorologie Dynamique model with Zooming capability-Simplified Atmospheric Chemistry System) inversion system to produce 10 different methane emission estimates at the global scale for the year 2005. The same methane sinks, emissions and initial conditions have been applied to produce the 10 synthetic observation datasets. The same inversion set-up (statistical errors, prior emissions, inverse procedure) is then applied to derive flux estimates by inverse modelling. Consequently, only differences in the modelling of atmospheric transport may cause differences in the estimated fluxes. Here in our framework, we show that transport model errors lead to a discrepancy of 27 Tg yr -1 at the global scale, representing 5% of total methane emissions. At continental and annual scales, transport model errors are proportionally larger than at the global scale, with errors ranging from 36 Tg yr -1 in North America to 7 Tg yr -1 in Boreal Eurasia (from 23 to 48%, respectively). At the model grid-scale, the spread of inverse estimates can reach 150% of the prior flux. Therefore, transport model errors contribute significantly to overall uncertainties in emission estimates by inverse modelling, especially when small spatial scales are examined. Sensitivity tests have been carried out to estimate the impact of the measurement network and the advantage of higher horizontal resolution in transport models. The large differences found between methane flux estimates inferred in these different configurations highly question the consistency of transport model errors in current inverse systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Measured Energy-Dependent Neutron Attenuation Through the Stacked Printed Circuit Boards

Engineers often stack printed circuit boards (PCBs), populated with semiconductor devices, one behind the other in neutron beams to improve the error statistics in radiation-effects tests. As the neutron beam traverses through the boards, the neutron flux is both attenuated and the energy spectrum changes, which may have a significant impact on the results of such tests. In this work, we have developed a technique to measure the energy-dependent neutron attenuation due to the PCBs and to determine the energy-dependent effective cross section of neutron interactions with a single PCB. With this effective cross section, engineers can calculate the change in the neutron energy spectrum and the total neutron attenuation after any number of identical PCBs in the stack. This information will allow test engineers to partially correct inaccuracies in soft error tests.

42 ENGINEERING↗
DOE OSTI2016

The Sloan Digital Sky Survey Reverberation Mapping Project: An investigation of biases in C iv emission line properties

We investigate the dependence on data quality of quasar properties measured from the C iv emission line region at high redshifts. Our measurements come from 32 epochs of Sloan Digital Sky Survey Reverberation Mapping Project spectroscopic observations of 482 z>1.46 quasars. We compare the differences between measurements made from the single-epoch (SE) and coadded spectra, focusing on the C iv λ1549 emission line because of its importance for studies of high-redshift quasar demographics and physical properties, including black hole masses. In addition to statistical errors increasing (by factors of ∼2–4), we find increasing systematic offsets with decreasing signal-to-noise ratio (S/N). The systematic difference (measurement uncertainty) in our lowest-S/N (<5) subsample between the SE and coadded spectrum (i) C iv equivalent width is 17 Å (31 Å), (ii) centroid wavelength is <1 Å (2 Å), and fractional velocity widths, ΔV/V, characterized by (iii) the line dispersion, σ {sub l}, is 0.104 (0.12), and (iv) the mean absolute deviation (MAD) is 0.072 (0.11). These remain smaller than the 1σ measurement uncertainties for all subsamples considered. The MAD is found to be the most robust line-width characterization. Offsets in the C iv FWHM velocity width and the C iv profile characterized by FWHM/σ {sub l} are only smaller than the statistical uncertainties when S/N > 10, although offsets in lower-S/N spectra exceed the statistical uncertainties by only a factor of ∼1.5 and may depend on the type of functional fit to the line. Characterizing the C iv line profile by the kurtosis is the least robust property investigated, as the median systematic coadded–SE measurement differences are larger than the statistical uncertainties for all S/N subsamples.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Optimizing time integration for accurate recovery of shockwave interface location in radiography

We present simulations and experiments of time integrated radiographic imaging of a moving 1D shock wave front and a quantitative method for determining the statistical error in locating the shock front as a function of integration time and noise in the radiograph. We discuss the trade-off between increasing motion blur, which leads to decreased shock front location certainty, and increasing signal-to-noise, which leads to improved image quality with increasing integration time. We find an optimum integration time between a short integration time, where noise limits the error, and a long integration time, where motion blurring limits the error. This methodology can be used to tune experimental configurations to obtain the highest quality radiograph for a given experimental configuration.

Bremsstrahlung↗
DOE OSTI2020

The completed SDSS-IV extended baryon oscillation spectroscopic survey: pairwise-inverse probability and angular correction for fibre collisions in clustering measurements

ABSTRACT The completed extended Baryon Oscillation Spectroscopic Survey (eBOSS) catalogues contain redshifts of 344 080 quasars at 0.8 < z < 2.2, 174 816 luminous red galaxies between 0.6 < z < 1.0, and 173 736 emission-line galaxies over 0.6 < z < 1.1 in order to constrain the expansion history of the Universe and the growth rate of structure through clustering measurements. Mechanical limitations of the fibre-fed spectrograph on the Sloan telescope prevent two fibres being placed closer than 62 arcsec in a single pass of the instrument. These ‘fibre collisions’ strongly correlate with the intrinsic clustering of targets and can bias measurements of the two-point correlation function resulting in a systematic error on the inferred values of the cosmological parameters. We combine the new techniques of pairwise-inverse probability and the angular upweighting (PIP+ANG) to correct the clustering measurements for the effect of fibre collisions. Using mock catalogues, we show that our corrections provide unbiased measurements, within data precision, of both the projected $\rm {\mathit{ w}_p}\left(\mathit{ r}_p\right)$ and the redshift-space multipole ξ(ℓ = 0, 2, 4)(s) correlation functions down to $0.1\, h^{-1}{\rm Mpc}$, regardless of the tracer type. We apply the corrections to the eBOSS DR16 catalogues. We find that, on scales $s\gtrsim 20\, h^{-1}{\rm Mpc}$ for ξℓ, as used to make baryon acoustic oscillation and large-scale redshift-space distortion measurements, approximate methods such as nearest-neighbour upweighting are sufficiently accurate given the statistical errors of the data. Using the PIP method, for the first time for a spectroscopic program of the Sloan Digital Sky Survey, we are able to successfully access the one-halo term in the clustering measurements down to $\sim 0.1\, h^{-1}{\rm Mpc}$ scales. Our results will therefore allow studies that use the small-scale clustering to strengthen the constraints on both cosmological parameters and the halo occupation distribution models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

What is the physical origin of strong Lyα emission? I. Demographics of Lyα emitter structures

We present the results of structure analyses for a large sample of 426 Lyα emitters (LAEs) at z ∼ 2.2 that are observed with the Hubble Space Telescope/Advanced Camera for Surveys and WFC3-IR during deep extra-galactic legacy surveys. We confirm that the merger fraction and the average ellipticity of LAE's stellar component are 10%-30% and 0.4-0.6, respectively, that are comparable with previous study results. We successfully identify that some LAEs have a spatial offset between Lyα and stellar-continuum emission peaks, δ{sub Lyα}, by ∼2.5-4 kpc beyond our statistical errors. To uncover the physical origin of strong Lyα emission found in LAEs, we investigate the Lyα equivalent width (EW) dependences of three structural parameters: merger fraction, δ{sub Lyα}, and ellipticity of stellar distribution in the range of EW(Lyα) = 20-250 Å. Contrary to expectations, we find that the merger fraction does not significantly increase with Lyα EW. We reveal an anti-correlation between δ{sub Lyα} and EW(Lyα) using a Kolmogorov-Smirnov (K-S) test. There is a trend that the LAEs with a large Lyα EW have a small ellipticity. This is consistent with the recent theoretical claims that Lyα photons can more easily escape from face-on disks having a small ellipticity, due to less inter-stellar gas along the line of sight, although our K-S test indicates that this trend is not statistically significant. Our results of Lyα-EW dependence generally support the idea that an H I column density is a key quantity determining Lyα emissivity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

MO-G-BRE-04: Automatic Verification of Daily Treatment Deliveries and Generation of Daily Treatment Reports for a MR Image-Guided Treatment Machine

Purpose: Two aims of this work were to develop a method to automatically verify treatment delivery accuracy immediately after patient treatment and to develop a comprehensive daily treatment report to provide all required information for daily MR-IGRT review. Methods: After systematically analyzing the requirements for treatment delivery verification and understanding the available information from a novel MR-IGRT treatment machine, we designed a method to use 1) treatment plan files, 2) delivery log files, and 3) dosimetric calibration information to verify the accuracy and completeness of daily treatment deliveries. The method verifies the correctness of delivered treatment plans and beams, beam segments, and for each segment, the beam-on time and MLC leaf positions. Composite primary fluence maps are calculated from the MLC leaf positions and the beam-on time. Error statistics are calculated on the fluence difference maps between the plan and the delivery. We also designed the daily treatment delivery report by including all required information for MR-IGRT and physics weekly review - the plan and treatment fraction information, dose verification information, daily patient setup screen captures, and the treatment delivery verification results. Results: The parameters in the log files (e.g. MLC positions) were independently verified and deemed accurate and trustable. A computer program was developed to implement the automatic delivery verification and daily report generation. The program was tested and clinically commissioned with sufficient IMRT and 3D treatment delivery data. The final version has been integrated into a commercial MR-IGRT treatment delivery system. Conclusion: A method was developed to automatically verify MR-IGRT treatment deliveries and generate daily treatment reports. Already in clinical use since December 2013, the system is able to facilitate delivery error detection, and expedite physician daily IGRT review and physicist weekly chart review.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Photo-z outlier self-calibration in weak lensing surveys

Calibrating photometric redshift errors in weak lensing surveys with external data is extremely challenging. In this work, we show that both Gaussian and outlier photo-z parameters can be self-calibrated from the data alone. This comes at no cost for the neutrino masses, curvature and dark energy equation of state w0, but with a 65% degradation when both w0 and wa are varied. We perform a realistic forecast for the Vera Rubin Observatory (VRO) Legacy Survey of Space and Time (LSST) 3× 2 analysis, combining cosmic shear, projected galaxy clustering and galaxy - galaxy lensing. We confirm the importance of marginalizing over photo-z outliers. We examine a subset of internal cross-correlations, dubbed "null correlations", which are usually ignored in 3× 2 analyses. Despite contributing only ~ 10% of the total signal-to-noise, these null correlations improve the constraints on photo-z parameters by up to an order of magnitude. Using the same galaxy sample as sources and lenses dramatically improves the photo-z uncertainties too. Together, these methods add robustness to any claim of detected new Physics, and reduce the statistical errors on cosmology by 15% and 10% respectively. Finally, including CMB lensing from an experiment like Simons Observatory or CMB-S4 improves the cosmological and photo-z posterior constraints by about 10%, and further improves the robustness to systematics. To give intuition on the Fisher forecasts, we examine in detail several toy models that explain the origin of the photo-z self-calibration. Our Fisher code LaSSI (Large-Scale Structure Information), which includes the effect of Gaussian and outlier photo-z, shear multiplicative bias, linear galaxy bias, and extensions to LCDM, is publicly available at \href{https://github.com/EmmanuelSchaan/LaSSI}{https://github.com/EmmanuelSchaan/LaSSI}.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Galaxy blending effects in deep imaging cosmic shear probes of cosmology

ABSTRACT Upcoming deep imaging surveys such as the Vera C. Rubin Observatory Legacy Survey of Space and Time will be confronted with challenges that come with increased depth. One of the leading systematic errors in deep surveys is the blending of objects due to higher surface density in the more crowded images; a considerable fraction of the galaxies which we hope to use for cosmology analyses will overlap each other on the observed sky. In order to investigate these challenges, we emulate blending in a mock catalogue consisting of galaxies at a depth equivalent to 1.3 yr of the full 10-yr Rubin Observatory that includes effects due to weak lensing, ground-based seeing, and the uncertainties due to extraction of catalogues from imaging data. The emulated catalogue indicates that approximately 12 per cent of the observed galaxies are ‘unrecognized’ blends that contain two or more objects but are detected as one. Using the positions and shears of half a billion distant galaxies, we compute shear–shear correlation functions after selecting tomographic samples in terms of both spectroscopic and photometric redshift bins. We examine the sensitivity of the cosmological parameter estimation to unrecognized blending employing both jackknife and analytical Gaussian covariance estimators. An ∼0.025 decrease in the derived structure growth parameter S8 = σ8(Ωm/0.3)0.5 is seen due to unrecognized blending in both tomographies with a slight additional bias for the photo-z-based tomography. This bias is greater than the 2σ statistical error in measuring S8.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

State densities of heavy nuclei in the static-path plus random-phase approximation

We report that nuclear state densities are important inputs to statistical models of compound-nucleus reactions. State densities are often calculated with self-consistent mean-field approximations that do not include important correlations and must be augmented with empirical collective enhancement factors. Here, we benchmark the static-path plus random-phase approximation (SPA + RPA) to the state density in a chain of samarium isotopes 148–155 Sm against exact results (up to statistical errors) obtained with the shell-model Monte Carlo (SMMC) method. The SPA + RPA method incorporates all static fluctuations beyond the mean field together with small-amplitude quantal fluctuations around each static fluctuation. Using a pairing plus quadrupole interaction, we show that the SPA + RPA state densities agree well with the exact SMMC densities for both the even- and odd-mass isotopes. For the even-mass isotopes, we also compare our results with mean-field state densities calculated with the finite-temperature Hartree-Fock-Bogoliubov (HFB) approximation. We find that the SPA + RPA repairs the deficiencies of the mean-field approximation associated with broken rotational symmetry in deformed nuclei and with the violation of particle-number conservation in the pairing condensate. In particular, in deformed nuclei the SPA + RPA reproduces the rotational enhancement of the state density relative to the mean-field state density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

DESI DR2 Galaxy luminosity functions

We present galaxy luminosity functions (LFs) for the Dark Energy Spectroscopic Instrument (DESI) DR2 Bright Galaxy Survey (BGS) in the g, r, z, and $w1$ bands over $0.002\lt z\lt 0.6$. Our analysis uses updated k-corrections and evolutionary corrections, including new polynomial kcorrection fits derived from BGS Year 1 data that supersede earlier GAMA-based prescriptions. Exploiting the statistical power of DESI, we measure LFs to very faint magnitudes, reaching $^{0.1}M_r-5\log h\sim -10$. Independent measurements from the North and South survey regions show excellent agreement around the LF knee, but the very small statistical uncertainties reveal that simple analytic forms fail to capture the full LF shape. The bright end departs from a pure exponential decline, while the faint end exhibits complex, non-powerlaw behaviour, including a pronounced upturn at $^{0.1}M_r-5\log h\gtrsim -15$, which is stronger for red galaxies than for blue. We show that our LFs are largely complete for galaxies with surface brightness $\mu _{50}\lt 25$, and that an apparent steepening fainter than $-13$ is driven primarily by local overdensity and fragmentation of large galaxies. A systematic North–South offset at the brightest magnitudes is traced to red galaxies and may reflect shallower North photometry underestimating extended earlytype profiles, although this remains inconclusive. We therefore also provide LFs based on model Petrosian magnitudes. Redshift splitting reveals small but significant residuals, indicating limitations of a simple global evolutionary model. Using the redshift limits of J. Loveday et al. (2012), we find excellent agreement with GAMA, with substantially reduced statistical errors. These measurements provide a precise reference for studies of environmental and population-dependent LFs and for testing galaxy formation models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Stochastic density functional theory

Linear-scaling implementations of density functional theory (DFT) reach their intended efficiency regime only when applied to systems having a physical size larger than the range of their Kohn–Sham density matrix (DM). This causes a problem since many types of large systems of interest have a rather broad DM range and are therefore not amenable to analysis using DFT methods. For this reason, the recently proposed stochastic DFT (sDFT), avoiding exhaustive DM evaluations, is emerging as an attractive alternative linear-scaling approach. This review develops a general formulation of sDFT in terms of a (non)orthogonal basis representation and offers an analysis of the statistical errors (SEs) involved in the calculation. Using a new Gaussian-type basis-set implementation of sDFT, applied to water clusters and silicon nanocrystals, it demonstrates and explains how the standard deviation and the bias depend on the sampling rate and the system size in various types of calculations. We also develop a basis-set embedded-fragments theory, demonstrating its utility for reducing the SEs for energy, density of states and nuclear force calculations. Finally, we discuss the algorithmic complexity of sDFT, showing it has CPU wall-time linear-scaling. Lastly, the method parallelizes well over distributed processors with good scalability and therefore may find use in the upcoming exascale computing architectures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Theoretical Prerequisites for PIK Reactor Subcriticality Control Aided by Spectral Noise Analysis

The possibility of measuring and controlling PIK reactor subcriticality on the basis of spectral analysis of reactor noise is examined. A simple lumped-parameter working model of the reactor and an equivalent noise source in accord with Schottky’s relation are used. Attention is focused on the statistical error of measurement of the spectral density of reactor noise owing to the finite measurement time and the random nature of neutron detection in a detector. It is shown that in connection with this error the measurements will require a long time and high-sensitivity detectors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2016

Uncertainty and feasibility of dynamical downscaling for modeling tropical cyclones for storm surge simulation

This paper presents a modeling study conducted to evaluate the uncertainty of a regional model in simulating hurricane wind and pressure fields, and the feasibility of driving coastal storm surge simulation using an ensemble of region model outputs produced by 18 combinations of three convection schemes and six microphysics parameterizations, using Hurricane Katrina as a test case. Simulated wind and pressure fields were compared to observed H*Wind data for Hurricane Katrina and simulated storm surge was compared to observed high-water marks on the northern coast of the Gulf of Mexico. The ensemble modeling analysis demonstrated that the regional model was able to reproduce the characteristics of Hurricane Katrina with reasonable accuracy and can be used to drive the coastal ocean model for simulating coastal storm surge. Results indicated that the regional model is sensitive to both convection and microphysics parameterizations that simulate moist processes closely linked to the tropical cyclone dynamics that influence hurricane development and intensification. The Zhang and McFarlane (ZM) convection scheme and the Lim and Hong (WDM6) microphysics parameterization are the most skillful in simulating Hurricane Katrina maximum wind speed and central pressure, among the three convection and the six microphysics parameterizations. Error statistics of simulated maximum water levels were calculated for a baseline simulation with H*Wind forcing and the 18 ensemble simulations driven by the regional model outputs. The storm surge model produced the overall best results in simulating the maximum water levels using wind and pressure fields generated with the ZM convection scheme and the WDM6 microphysics parameterization.

Yang, Zhaoqing↗