Search NASA⌕ Search

SEARCH · Search NASA

Results for “sample average approximation”

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 361 records · Page 20

Estimating Long Term Surface Soil Moisture in the GCIP Area From Satellite Microwave Observations

Soil moisture is an important component of the water and energy balances of the Earth's surface. Furthermore, it has been identified as a parameter of significant potential for improving the accuracy of large-scale land surface-atmosphere interaction models. However, accurate estimates of surface soil moisture are often difficult to make, especially at large spatial scales. Soil moisture is a highly variable land surface parameter, and while point measurements are usually accurate, they are representative only of the immediate site which was sampled. Simple averaging of point values to obtain spatial means often leads to substantial errors. Since remotely sensed observations are already a spatially averaged or areally integrated value, they are ideally suited for measuring land surface parameters, and as such, are a logical input to regional or larger scale land process models. A nine-year database of surface soil moisture is being developed for the Central United States from satellite microwave observations. This region forms much of the GCIP study area, and contains most of the Mississippi, Rio Grande, and Red River drainages. Daytime and nighttime microwave brightness temperatures were observed at a frequency of 6.6 GHz, by the Scanning Multichannel Microwave Radiometer (SMMR), onboard the Nimbus 7 satellite. The life of the SMMR instrument spanned from Nov. 1978 to Aug. 1987. At 6.6 GHz, the instrument provided a spatial resolution of approximately 150 km, and an orbital frequency over any pixel-sized area of about 2 daytime and 2 nighttime passes per week. Ground measurements of surface soil moisture from various locations throughout the study area are used to calibrate the microwave observations. Because ground measurements are usually only single point values, and since the time of satellite coverage does not always coincide with the ground measurements, the soil moisture data were used to calibrate a regional water balance for the top 1, 5, and 10 cm surface layers in order to interpolate daily surface moisture values. Such a climate-based approach is often more appropriate for estimating large-area spatially averaged soil moisture because meteorological data are generally more spatially representative than isolated point measurements of soil moisture. Vegetation radiative transfer characteristics, such as the canopy transmissivity, were estimated from vegetation indices such as the Normalized Difference Vegetation Index (NDVI) and the 37 GHz Microwave Polarization Difference Index (MPDI). Passive microwave remote sensing presents the greatest potential for providing regular spatially representative estimates of surface soil moisture at global scales. Real time estimates should improve weather and climate modelling efforts, while the development of historical data sets will provide necessary information for simulation and validation of long-term climate and global change studies.

Owe, Manfred↗

Atmospheric parameters of field subdwarf B stars

High signal-to-noise ratio optical spectrophotometry of a sample of field subluminous B stars drawn largely from the Palomar Green ultraviolet-excess survey is analyzed with a new grid of model atmospheres and synthetic spectra. The stellar effective temperatures, surface gravities, and photospheric helium abundances are determined simultaneously from a detailed analysis of hydrogen and helium absorption line profiles. The derived temperatures and gravities place the subluminous B stars in the theoroetical Hertzsprung-Russell (H-R) diagram along and bounded below by theoretical sequences of the zero-age extended horizontal branch, lending strong support to the hypothesis that these stars are composed of helium-burning cores of approximately 0.5 solar mass overlain by very thin layers of hydrogen (approximately less than 0.02 solar mass). Various scenarios for their past evolutionary history are examined in the context of their probable future evolution into white dwarfs of lower than average mass.

Saffer, Rex A.↗

Chemical studies of H chondrites. 4: New data and comparison of Antarctic suites

We report data for the trace elements Au, Co, Sb, Ga, Rb, Ag, Se, Cs, Te, Zn, Cd, Bi, Ti, and In (ordered by putative volatility during nebular condensation and accretion) determined by neutron activation analysis in 13 H5 chondrites from Victoria Land and 20 H4-6 chondrites from Queen Maud Land, Antarctica. These and earlier results provide Antarctic sample suites of 34 chondrites from Victoria Land and 25 from Queen Maud Land. Treatment of data for the most volatile 10 elements (Rb to In) in these studies by multivariate statistical techniques more robust, as well as more conservative, than conventional linear discriminant analysis and logistic regression demonstrates that compositions differ at marginally significant levels. This difference cannot be explained by trivial (terrestrial) causes and becomes more significant, despite the smaller size of the database, when comparisons are limited to data from a single analyst and when all upper limits are eliminated from consideration. The Victoria Land and Queen Maud Land suites have different mean terrestrial ages (approximately 300 kyr and approximately 100 kyr, respectively) and age distributions, suggesting that a time-dependent variation of chondritic sources with different thermal histories is responsible. As a result, these two Antarctic suites are, on average, chemically distinguishable from each other. Since H chondrites serve as a paradigm for other meteorite classes, these results indicate that the near-Earth populations of planetary materials varied with time on the 10(exp 5)-year timescale.

Wolf, Stephen F.↗

Comparison of Morphologies of Apollo 17 Dust Particles with Lunar Simulant, JSC-1

Lunar dust (< 20 microns) makes up approx.20 wt.% of the lunar soil. Because of the abrasive and adhering nature of lunar soil, a detailed knowledge of the morphology (size, shape and abundance) of lunar dust is important for dust mitigation on the Moon. This represents a critical step towards the establishment of long-term human presence on the Moon (Taylor et al. 2005). Machinery design for in-situ resource utilization (ISRU) on the Moon also requires detailed information on dust morphology and general physical/chemical characteristics. Here, we report a morphological study of Apollo 17 dust sample 70051 and compare it to lunar soil stimulant, JSC-1. W e have obtained SEM images of dust grains from sample 70051 soil (Fig. 1). The dust grains imaged are composed of fragments of minerals, rocks, agglutinates and glass. Most particles consist largely of agglutinitic impact glass with their typical vesicular textures (fine bubbles). All grains show sub-angular to angular shapes, commonly with sharp edges, common for crushed glass fragments. There are mainly four textures: (1) ropey-textured pieces (typical for agglutinates), (2) angular shards, (3) blocky bits, and (4) Swiss-cheese grains. This last type with its high concentration of submicron bubbles, occurs on all scales. Submicron cracks are also present in most grains. Dust-sized grains of lunar soil simulant, JSC-1, were also studied. JSC-1 is a basaltic tuff with relatively high glass content (approx.50%; McKay et al. 1994). It was initially chosen in the early 90s to approximate the geotechnical properties of the average lunar soil (Klosky et al. 1996). JSC-1 dust grains also show angular blocky and shard textures (Fig. 2), similar to those of lunar dust. However, the JSC-1 grains lack the Swiss-cheese textured particles, as well as submicron cracks and bubbles in most grains.

Liu, Yang↗

Polar stratospheric optical depth observed between 1978 and 1985

Observations of the stratospheric optical depth at 1.0 micron obtained for high latitudes are presented for a 7-year period. Weekly averaged data determined from measurements made by the Stratospheric Aerosol Measurement experiment from October 1978-1985 show that the overall yearly values in both polar regions are controlled by volcanic perturbations, with most volcanic effects being experienced in Arctic latitudes. Conservatively, peak values found in the Antarctic region were approximately 0.02 and in the Arctic region about 0.55. Probable values for these regions are estimated to be 0.26 and 0.11, respectively. The weekly averaged data also show the seasonal fluctuations due to microphysical and dynamical processes. Comparison of the optical depth record with a weekly averaged 50-mbar temperature record indicates that polar stratospheric clouds are present in the southern high latitudes each year near this level from early June to early September. A depression observed in the optical depth record each austral spring season is believed to be the result of the downward displacement of particles caused by subsidence and sedimentation during the course of winter. Following the breakup of the vortex, optical depth values increase as aerosol is transported poleward. These features are noted to be present in the Arctic region as well, but on a smaller scale because of the satellite sampling methodology and the averaging scheme employed.

Mccormick, M. P.↗

FIRE extended time/limited area observations at Palisades, New York

Downwelling shortwave and longwave irradiation are being continuously monitored at Palisades, New York as part of the First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE) Extended Time/Limited Area Initiative. In addition, fisheye (180 degree) sky photographs are taken at the times of NOAA 9 and LANDSAT satellite overpasses on select days, particularly when cirrus clouds are present. Measurements of incoming shortwave (0.28 to 2.80 microns) hemispheric and diffuse, hemispheric near infrared (0.7 to 2.80 microns), and downwelling hemispheric infrared (4.0 to 50.0 microns) irradiation have been made from a rooftop location on the grounds of the Lamont-Doherty Geological Observatory since December 1986. The three Eppley Precision Spectral Pyranometers and the Eppley Pyrgeometer used to measure these variables were calibrated with Colorado State University instruments at Madison, Wisconsin as part of the FIRE Intensive Laboratory. Pyrgeometer output contains an adjustment for body temperature but not for dome temperature. Data are transmitted to a Campbell CR-21 Digital Recorder, where one minute averages of ten second samples are stored and subsequently dumped to a cassette recorder. Using a Campbell C-20 Cassette Interface, these data are transferred to an Apple Macintosh computer for analysis and for archiving on floppy disks. In addition to the raw irradiances collected, variables derived from these data are generated and stored. These include: the ratio of near infrared irradiation to visible irradiation and the fraction of the full shortwave irradiation which is diffuse; and will soon include: shortwave transmissivity and optical depth in the shortwave. Sky photographs are taken with an Olympus OM2-N 35 mm camera and are timed to be coincident with overpassing NOAA and LANDSAT satellites. Palisades is within the field of view of the NOAA 9 daily in the middle to late afternoon. The satellite viewing angle is within 45 degrees of nadir over Palisades on approximately half of the passes.

Robinson, David A.↗

Comparison of Marine Boundary Layer Cloud Properties from CERES-MODIS Edition 4 and DOE ARM AMF Measurements at the Azores

Marine boundary layer (MBL) cloud properties derived from the NASA Clouds and the Earth's Radiant Energy System (CERES) project using Terra and Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) data are compared with observations taken at the Department of Energy Atmospheric Radiation Measurement (ARM) Mobile Facility at the Azores (AMF-Azores) site from June 2009 through December 2010. Cloud properties derived from ARM ground-based observations were averaged over a 1 h interval centered at the satellite overpass time, while the CERES-MODIS (CM) results were averaged within a 30 km×30 km grid box centered over the Azores site. A total of 63 daytime and 92 nighttime single-layered overcast MBL cloud cases were selected from 19 months of ARM radar-lidar and satellite observations. The CM cloud top/base heights (Htop/Hbase) were determined from cloud top/base temperatures (Ttop/Tbase) using a regional boundary layer lapse rate method. For daytime comparisons, the CM-derived Htop (Hbase), on average, is 0.063 km (0.068 km) higher (lower) than its ARM radar-lidar-observed counterpart, and the CM-derived Ttop and Tbase are 0.9 K less and 2.5 K greater than the surface values with high correlations (R(sup 2) = 0.82 and 0.84, respectively). In general, the cloud top comparisons agree better than the cloud base comparisons, because the CM cloud base temperatures and heights are secondary products determined from cloud top temperatures and heights. No significant day-night difference was found in the analyses. The comparisons of MBL cloud microphysical properties reveal that when averaged over a 30 km× 30 km area, the CM-retrieved cloud droplet effective radius (re) at 3.7 micrometers is 1.3 micrometers larger than that from the ARM retrievals (12.8 micrometers), while the CM-retrieved cloud liquid water path (LWP) is 13.5 gm( exp -2) less than its ARM counterpart (114.2 gm( exp-2) due to its small optical depth (9.6 versus 13.7). The differences are reduced by 50% when the CM averages are computed only using the MODIS pixel nearest the AMF site. Using the effective radius retrieved using 2.1 micrometers channel to calculate LWP can reduce the difference between the CM and ARM microwave radiometer retrievals from 13.7 to 2.1 gm2. The 10% differences between the ARM and CERES-MODIS LWP and r(sub e) retrievals are within the uncertainties of the ARM LWP (approximately 20gm( exp -2)) and r(sub e) (approximately 10%) retrievals; however, the 30% difference in optical depth is significant. Possible reasons contributing to this discrepancy are increased sensitivities in optical depth from both surface retrievals when τ is approximately 10 and topography. The τ differences vary with wind direction and are consistent with the island orography.Much better agreement in τ is obtained when using only those data taken when the wind is from the northeast, where topographical effects on the sampled clouds are minimal.

Xi, Baike↗

A Possible Explanation for the Z-R Parameter Inconsistency when Comparing Stratiform and Convective Rainfall

The well-known Z-R power law Z A R(sup b) uses two parameters, A and b, in order to relate rainfall rate R to measured weather radar reflectivity Z. A common method used by researchers is to compute Z and R from disdrometer data and then extract the A-b parameter pair from a log-linear line fit to a scatter plot of Z-R pairs. Even though it may seem far more truthful to extract the parameter pair from a fit of radar Z(sub R) versus gauge rainfall rate R(sub G), the extreme difference in spatial and temporal sampling volumes between radar and rain gauge creates a slew of problems that can generally only be solved by using rain gauge arrays and long sampling averages. Disdrometer derived A-b parameters are easily obtained and can provide information for the study of stratiform versus convective rainfall. However, an inconsistency appears when comparing averaged A-b pairs from various researchers. Values of b range from 1.26 to 1.51 for both stratiform and convective events. Paradoxically the values of A fall into three groups: 150 to 200 for convective; 200 to 400 for stratiform; and 400 to 500 again for convective. This apparent inconsistency can be explained by computing the A-b pair using the gamma DSD coupled with a modified drop terminal velocity model, v(D) alpha D(sup beta) - w, where w is a somewhat artificial constant vertical velocity of the air above the disdrometer. This model predicts three regions of A, corresponding to w less than 0, w = 0, and w greater than 0, which approximately matches observed data.

gamma DSD↗

Non-intrusive reduced-order modeling for dynamical systems with spatially localized features

This work presents a non-intrusive reduced-order modeling framework for dynamical systems with spatially localized features characterized by slow singular value decay. The proposed approach builds upon two existing methodologies for reduced and full-order non-intrusive modeling, namely Operator Inference (OpInf) and sparse Full-Order Model (sFOM) inference. We decompose the domain into two complementary subdomains that exhibit fast and slow singular value decay. The dynamics of the subdomain exhibiting slow singular value decay are learned with sFOM while the dynamics with intrinsically low dimensionality on the complementary subdomain are learned with OpInf. The resulting, coupled OpInf-sFOM formulation leverages the computational efficiency of OpInf and the high resolution of sFOM, and thus enables fast non-intrusive predictions for conditions beyond those sampled in the training data set. A novel regularization technique with a closed-form solution based on the Gershgorin disk theorem is introduced to promote stable sFOM and OpInf models. We also provide a data-driven indicator for subdomain selection and ensure solution smoothness over the interface via a post-processing interpolation step. We evaluate the efficiency of the approach in terms of offline and online speedup through a quantitative, parametric computational cost analysis. We demonstrate the coupled OpInf-sFOM formulation for two test cases: a one-dimensional Burgers’ model for which accurate predictions beyond the span of the training snapshots are presented, and a two-dimensional parametric model for the Pine Island Glacier ice thickness dynamics, for which the OpInf-sFOM model achieves an average prediction error on the order of 1% with an online speedup factor of approximately 8$\times$ compared to the numerical simulation.

42 ENGINEERING↗

Vitrification of Hanford Tank Wastes for Condensate Recycle and Feed Composition Changeover Testing

During the vitrification of Hanford Site nuclear waste at the Waste Treatment and Immobilization Plant (WTP), which is a part of the safe and efficient retrieval, treatment, and disposal mission of the U.S. Department of Energy Office of River Protection, the offgas condensate generated from the waste-to-glass conversion is currently planned to be concentrated by evaporation in the Effluent Management Facility (EMF). This concentrated condensate can then be recycled back to the incoming waste and vitrified. To test the recycle process, an apparatus was designed and built to mimic the EMF evaporator and was then used to concentrate a volume of condensate produced during the vitrification of a sample of Hanford tank 241-AP-107 (referred to herein as AP-107) waste in a continuous laboratory-scale melter (CLSM). The concentrated condensate was added to an additional sample of AP-107 waste, to mimic one round of the recycle process, and the combined solution was vitrified, producing a second round of recycle condensate. In the current study, the EMF test apparatus was used to concentrate the second-round recycle condensate under evaporation conditions (at 45 °C and 1.4 psia) designed to emulate EMF operation. The condensate was successfully concentrated by a factor of ~10 while retaining over 95 % of the technetium-99 ( 99 Tc), Cs, and I inventories in the concentrate. Another portion of AP-107 waste was retrieved by Washington River Protection Solutions, LLC (WRPS) and transferred to Pacific Northwest National Laboratory (PNNL), where it was pretreated and then combined with the second-round recycle AP-107 condensate concentrate and glass-forming chemicals (GFCs) to form the two-time recycle AP-107 melter feed, approximating a second round to the recycling action to be performed at the WTP. A portion of AP-105 waste was also retrieved by WRPS and provided to PNNL for pretreatment and combining with GFCs to form AP-105 melter feed. The two-time recycle AP-107 and AP-105 melter feeds were processed consecutively in the CLSM. The CLSM run proceeded for 13.63 hours, producing 9.70 kg of glass for an average glass production rate of 1464 kg m 2 d -1 during the two-time recycle AP-107 feed charging and 1568 kg m 2 d -1 during the AP-105 feed charging. The rate during AP-107 charging was essentially equivalent to the rate when processing no-recycle AP-107 feed and lower than that achieved when processing one-time recycle AP-107 feed. However, all rates were within the potential range of variability when processing melter feeds with similar composition in the CLSM. Likewise, the rate during AP-105 charging was higher than the previous rate processing AP-105, but within the potential CLSM range. The cold-cap characteristics changed from the typically thin AP-107 cold cap to a foamy-edged cold cap as previously seen with AP-105 shortly after transitioning to the AP-105 melter feed. The glass produced during the CLSM run was within 10 % of its target composition for the primary glass components. The CaO and Li 2 O targets varied by more than 1 wt% between the two-time recycle AP-107 and AP-105 glass targets and it took about 2 turnovers of the CLSM glass inventory to reach a relative chemical steady state in the glass for CaO and Li 2 O after the melter feed inputs were switched. The 99 Tc and total cesium content in the melter feeds were maintained at concentrations expected to be experienced at the WTP. During the CLSM run, while processing the two-time recycle AP-107 melter feed at a relative chemical steady state, the 99 Tc/Cs ratio was 10, and 34% of 99 Tc and 74% of Cs were retained in the glass. These values were higher than those measured in the CLSM run with one-time recycle AP-107 melter feed. After the transition to processing the AP-105 melter feed, when the production reached a relative chemical steady state, the 99 Tc/Cs ratio was 77 while 44% of 99 Tc was retained in the glass. The C's retention during this time frame reached 200% due to the excess C's in the glass after the target content decreased to 15% of its initial level in the two-time recycle AP-107 melter feed to the lower target in the Ap-105 melter feed. While iodine was below inductively coupled plasma mass spectrometry analytical reporting limits in the melter feed and glass samples, it was detected in quantities above the analytical reporting limits in the liquid and filter samples collected from the CLSM offgas treatment system. The behavior of iodine in the CLSM offgas treatment system followed a similar pattern to those of 99 Tc and Cs.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Fiscal Year 2022. Filtration of Hanford Tank 241-AP-101 Supernatant at 16 °C

Bench-scale filtration testing of ~9 liters of supernatant from Hanford waste tank 241-AP-101, chilled to 16 °C, was conducted using a backpulse dead-end filter (BDEF) filtration system equipped with a feed vessel and a Mott inline filter Model 6610 (Media Grade 5) in the hot cells of the Radiochemical Processing Laboratory at Pacific Northwest National Laboratory. This was done to assess the performance of the anticipated third feed to the Tank Side Cesium Removal (TSCR) system. The as-received samples were diluted to the target sodium concentration and transferred to 1.5-liter polyethylene bottles and held at 16 °C for approximately 1 week prior to filtration. The feed was filtered through the BDEF system at a targeted flux of 0.065 gpm/ft 2 to match the prototypic operation of the TSCR system. During filtration, the differential pressure required to effect filtration at 0.065 gpm/ft 2 increased little over the filtration campaign and never reached 2 psid (the TSCR action limit). This indicates that the TSCR filter should perform well when processing AP-101 supernatant. After completing filtration of the AP-101 feed, the filter was cleaned. Solids concentrated from the backpulse solutions displayed sodium nitrate-type phases, aluminum and silicon phases reported as cancrinite or nitrate-cancrinite, a mixed chromium-aluminum oxide, iron oxides, and Ca-bearing phases (calcite). Scanning electron microscopy analyses showed that the average particle size was 0.5 micron.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Numerical simulations of cometary dust

Most observations of comets are done photometrically or spectrophotometrically. The interpretation of the aperture-averaged flux is relatively simple for an isotropic, radially expanding coma of infinite extent - the canonical model. However, the interpretation of the observations is not so clear when the motion of the dust is affected by radiation pressure, or when the emission is time-varying and anisotropic. For example, in a sample of CCD images of 10 comets, Jewitt and Meech (1987, Ap.J. 317, 992) found that the photometric profiles of only three comets were consistent, within the observational errors, with the profiles predicted form the canonical model. Photometric observations with large apertures, however, seem to suggest that the canonical model may be quite adequate (c.f. Osip, Schleicher, and Millis, 1992, Icarus 98, 115). The dust itself is characterized by a size distribution, with size dependencies on the expansion velocity, the scattered and thermal radiation, the response to radiation pressure, and probably the density. How good then are the approximations normally used in determining the production rates of the dust when these effects are present? As part of a program to better understand the dynamics of cometary dust and gas, a computer program has been developed which numerically simulates the emission of both dust and gas from a tilted rotating nucleus.

Lien, David J.↗

The role of thermal shock in cyclic oxidation

The effect of thermal shock on the spalling of oxides from the surfaces of several commercial alloys was determined. The average cooling rate was varied from approximately 240 C per second to less than 1.0 C per second during cyclic oxidation tests in air. The tests consisted of one hundred cycles of one hour at the maximum temperature (1100 or 1200 C). The alloys were HOS-875, TD-Ni, TD-NiCrAl, IN-601, IN-702, and B-1900 plus Hf. All of these alloys exhibited partial spalling within the oxide rather than total oxide loss down to bare metal. Thermal shock resulted in deformation of the metal which in turn resulted, in most cases, in changing the oxide failure mode from compressive to tensile. Tensile failures were characterized by cracking of the oxide and little loss, while compressive failures were characterized by explosive loss of platelets of oxide. This behavior was confirmed by examination of mechanically stressed oxide scales. The thermally shocked oxides spalled less than the slow cooled samples with the exception of TD-NiCrAl. This material failed in a brittle manner rather than by plastic deformation.

Lowell, C. E.↗

Burn propagation rates of metals and alloys in gaseous oxygen

The average burn rates of several metals and alloys were determined at oxygen pressures between 3.45 and 68.91 MPa (500 and 10,000 psig) and ambient temperature. Several materials were tested at elevated sample temperatures. The test materials were fabricated into solid cylindrical rods and mounted vertically in the test chamber. A magnesium igniter was positioned at the bottom end of each test specimen to promote upward burn propagations. Nickel 200 and copper 102 could not be ignited at all oxygen pressures tested whereas Monel 400 appeared to ignite but quickly self-extinguished. The other materials tested burned the entire length of the test sample. Aluminum 6061 exhibited the fastest burn propagation rate. Inconel 718 burned slower than aluminum but faster than the stainless steels (types 304 and 316). Increasing oxygen pressure generally increased the burn propagation rate of the materials. Increasing the ambient temperature of the test specimens for several materials to approximately 850 K (1070 F) had little effect upon the ignition or burn properties of nickel 200 or Monel 400. Type 316 stainless steel exhibited an increase in its burn propagation rate at this higher temperature.

Benz, F. J.↗

Observational and Numerical Diagnostics of Galaxy Cluster Outer Regions

Aims. We present the analysis of a local (z = 0.04 - 0.2) sample of 31 galaxy clusters with the aim of measuring the density of the X-ray emitting gas in cluster outskirts. We compare our results with numerical simulations to set constraints on the azimuthal symmetry and gas clumping in the outer regions of galaxy clusters. Methods. We exploit the large field-of-view and low instrumental background of ROSAT/PSPC to trace the density of the intracluster gas out to the virial radius. We perform a stacking of the density profiles to detect a signal beyond r(sub 200) and measure the typical density and scatter in cluster outskirts. We also compute the azimuthal scatter of the profiles with respect to the mean value to look for deviations from spherical symmetry. Finally, we compare our average density and scatter profiles with the results of numerical simulations. Results. As opposed to several recent results, we observe a steepening of the density profiles beyond approximately 0.3r(sub 500). Comparing our density profiles with simulations, we find that non-radiative runs predict too steep density profiles, whereas runs including additional physics and/or gas clumping are in better agreement with the observed gas distribution. We note a systematic difference between cool-core and non-cool core clusters beyond approximately 0.3r(sub 200), which we explain by a different distribution of the gas in the two classes. Beyond approximately r(sub 500), galaxy clusters deviate significantly from spherical symmetry, with only little differences between relaxed and disturbed systems. We find good agreement between the observed and predicted scatter profiles, but only when the 1% densest clumps are filtered out in the simulations. Conclusions. The general trend of steepening density around the virial radius indicates that the shallow density profiles found in several recent works were probably obtained along particular directions (e.g., filaments) and are not representative of the typical behavior of clusters. Comparing our results with numerical simulations, we find that non-radiative simulations fail to reproduce the gas distribution, even well outside cluster cores. Therefore, a detailed treatment of gas cooling, star formation, clumping, and AGN feedback is required to construct realistic models of cluster outer regions.

Eckert, D.↗

Canaveral National Seashore Water Quality and Aquatic Resource Inventory

Mosquito Lagoon is a shallow, bar-built estuary located on the east central Florida Coast, primarily within the KSC boundary. The lagoon and watershed cover approximately 327 sq km (79422 acres) .The Lagoon occupies 159 sq km (37853 acres). Water depths average approximately 1m. The lagoon volume is approximately 1.6 x 10(exp 8)cu m. Water quality in Mosquito Lagoon is good. Salinity data typically range between 20 ppt and 35 ppt. The lowest value recorded was 4.5 ppt and the highest value was 37 ppt. Water temperatures fluctuate 2 - 3 C over a 24 h period. Cold front passage can rapidly alter water temperatures by 5 - 10 C or more in a short period of time. The highest temperature was 33.4 C and the lowest temperature was 8.8 C after a winter storm. Dissolved oxygen concentrations ranged from a low of 0.4 mg/l to a high of 15.3 mg/l. Extended periods of measurements below the Florida Department of Environmental Protection criteria of 4.0 mg/l were observed in fall and spring months suggesting high system respiration and oxygen demand. Metals such as antimony, arsenic, molybdenum and mercury were report as below detection limits for all samples. Cadmium, copper, chromium, silver, and zinc were found to be periodically above the Florida Department of Environmental Protection criteria for Class II and Class III surface waters.

Hall, C. R.↗

A hybrid 3D/2D field response calculation for liquid argon detectors with PCB based anode plane

Liquid Argon Time Projection Chamber (LArTPC) technology is commonly utilized in neutrino detector designs. It enables detailed reconstruction of neutrino events with high spatial precision and low energy threshold. Its field response (FR) model describes the time-dependent electric currents induced in the anode-plane electrodes when ionization electrons drift nearby. An accurate and precise FR is a crucial input to LArTPC detector simulations and charge reconstruction. Established LArTPC designs have been based on parallel wire planes. It allows accurate and computationally economic two-dimensional (2D) FR models utilizing the translational symmetry along the direction of the wires. Recently, novel LArTPC designs utilize electrodes formed on printed circuit board (PCB) in the shape of strips with through holes. The translational symmetry is no longer a good approximation near the electrodes and a new FR calculation that employs regions with three dimensions (3D) has been developed. Extending the 2D models to 3D would be computationally expensive. Fortuitously, the nature of strips with through holes allows for a computationally economic approach based on the finite-difference method (FDM). In this paper, we present a new software package pochoir that calculates LArTPC field response for these new strip-based anode designs. This package combines 3D calculations in the volume near the electrodes with 2D far-field solutions to achieve fast and precise field response computation. We apply the resulting FR to simulate and reconstruct samples of cosmic-ray muons and 39 Ar decays from a Vertical Drift (VD) detector prototype operated at CERN. We find the difference between real and simulated data within 5%. Current state-of-the-art LArTPC software requires a 2D FR which we provide by averaging over one dimension and estimate that variations lost in this average are smaller than 7%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Correcting for filter-based aerosol light absorption biases at the Atmospheric Radiation Measurement program's Southern Great Plains site using photoacoustic measurements and machine learning

Abstract. Measurement of light absorption of solar radiation by aerosols is vital for assessing direct aerosol radiative forcing, which affects local and global climate. Low-cost and easy-to-operate filter-based instruments, such as the Particle Soot Absorption Photometer (PSAP), that collect aerosols on a filter and measure light attenuation through the filter are widely used to infer aerosol light absorption. However, filter-based absorption measurements are subject to artifacts that are difficult to quantify. These artifacts are associated with the presence of the filter medium and the complex interactions between the filter fibers and accumulated aerosols. Various correction algorithms have been introduced to correct for the filter-based absorption coefficient measurements toward predicting the particle-phase absorption coefficient (Babs). However, the inability of these algorithms to incorporate into their formulations the complex matrix of influencing parameters such as particle asymmetry parameter, particle size, and particle penetration depth results in prediction of particle-phase absorption coefficients with relatively low accuracy. The analytical forms of corrections also suffer from a lack of universal applicability: different corrections are required for rural and urban sites across the world. In this study, we analyzed and compared 3 months of high-time-resolution ambient aerosol absorption data collected synchronously using a three-wavelength photoacoustic absorption spectrometer (PASS) and PSAP. Both instruments were operated on the same sampling inlet at the Department of Energy's Atmospheric Radiation Measurement program's Southern Great Plains (SGP) user facility in Oklahoma. We implemented the two most commonly used analytical correction algorithms, namely, Virkkula (2010) and the average of Virkkula (2010) and Ogren (2010)–Bond et al. (1999) as well as a random forest regression (RFR) machine learning algorithm to predict Babs values from the PSAP's filter-based measurements. The predicted Babs was compared against the reference Babs measured by the PASS. The RFR algorithm performed the best by yielding the lowest root mean square error of prediction. The algorithm was trained using input datasets from the PSAP (transmission and uncorrected absorption coefficient), a co-located nephelometer (scattering coefficients), and the Aerosol Chemical Speciation Monitor (mass concentration of non-refractory aerosol particles). A revised form of the Virkkula (2010) algorithm suitable for the SGP site has been proposed; however, its performance yields approximately 2-fold errors when compared to the RFR algorithm. To generalize the accuracy and applicability of our proposed RFR algorithm, we trained and tested it on a dataset of laboratory measurements of combustion aerosols. Input variables to the algorithm included the aerosol number size distribution from the Scanning Mobility Particle Sizer, absorption coefficients from the filter-based Tricolor Absorption Photometer, and scattering coefficients from a multiwavelength nephelometer. The RFR algorithm predicted Babs values within 5 % of the reference Babs measured by the multiwavelength PASS during the laboratory experiments. Thus, we show that machine learning approaches offer a promising path to correct for biases in long-term filter-based absorption datasets and accurately quantify their variability and trends needed for robust radiative forcing determination.

54 ENVIRONMENTAL SCIENCES↗