Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time resolution data”

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 1,279 records · Page 71

Project Description and Publications List for UAH CMMR

This research combines a state of the art X-ray Transmission Microscope, XTM, with a specially designed x-ray transparent horizontal Bridgman furnace to image (with resolutions up to 3 micrometers) the solidification of metal alloys in real-time. The objective is to obtain real-time dynamic data to provide direct measure of the solute profile in the liquid, phase coalescence and growth in the liquid, and the detailed interface morphology (e,g., dendrites and cells) during solidification. We are also enhancing the XTM data with precise solid-liquid interfacial temperature and the thermal gradient measurement techniques, and working on the application of this technology to the study of the fundamentals of solidification in microgravity. Over the last several years we have successfully imaged in real-time: interfacial-morphologies, phase growth, coalescence, incorporation of phases into the growing interface, and the solute boundary layer in the liquid at the solid-liquid interface. We have also measured true local growth rates and can evaluate segregation structures in the solid. Interfacial undercoolings are being measured either with a special Seebeck furnace or with micro-thermocouple arrays we are developing. These later techniques are presently being incorporated with the XTM furnace. This last year emphasized the investigation of the solute layer in the melt during solidification. Methods were developed to quantify the solute concentrations using x-ray absorption and to compare to predictions from simulations. In addition, work is being completed on a brass-board portable XTM that incorporates a vertical Bridgman furnace.

Kaukler, William F.↗

Validation of MODIS Aerosol Optical Depth Retrieval Over Land

Aerosol optical depths are derived operationally for the first time over land in the visible wavelengths by MODIS (Moderate Resolution Imaging Spectroradiometer) onboard the EOSTerra spacecraft. More than 300 Sun photometer data points from more than 30 AERONET (Aerosol Robotic Network) sites globally were used in validating the aerosol optical depths obtained during July - September 2000. Excellent agreement is found with retrieval errors within (Delta)tau=+/- 0.05 +/- 0.20 tau, as predicted, over (partially) vegetated surfaces, consistent with pre-launch theoretical analysis and aircraft field experiments. In coastal and semi-arid regions larger errors are caused predominantly by the uncertainty in evaluating the surface reflectance. The excellent fit was achieved despite the ongoing improvements in instrument characterization and calibration. This results show that MODIS-derived aerosol optical depths can be used quantitatively in many applications with cautions for residual clouds, snow/ice, and water contamination.

Chu, D. A.↗

Further Studies on the Physical and Biogeochemical Causes for Large Interannual Changes in the Patagonian Shelf Spring-Summer Phytoplankton Bloom Biomass

A very strong and persistent phytoplankton bloom was observed by ocean color satellites during September - December 2003 along the northern Patagonian shelf. The 2003 bloom had the highest extent and chlorophyll a (Chl-a) concentrations of the entire Sea-viewing Wide Field-of-view Sensor (SeaWiFS) period (1997 to present). SeaWiFS-derived Chl-a exceeded 20 mg/cu m in November at the bloom center. The bloom was most extensive in December when it spanned more than 300 km across the shelf and nearly 900 km north-south (35degS to 43degS). The northward reach and the deep penetration on the shelf of the 2003 bloom were quite anomalous when compared with other years, which showed the bloom more confined to the Patagonian shelf break (PSB). The PSB bloom is a conspicuous austral spring-summer feature detected by ocean color satellites and its timing can be explained using the Sverdrup critical depth theory. Based on high-resolution numerical simulations, in situ and remote sensing data, we provide some suggestions for the probable mechanisms responsible for that large interannual change of biomass as seen by ocean color satellites. Potential sources of macro and micro (e.g., Fe) nutrients that sustain the high phytoplankton productivity of the Patagonian shelf waters are identified, and the most likely physical processes that maintain the nutrient balance in the region are discussed.

Signorini, Sergio R.↗

The Geospace Dynamics Observatory; A Paradigm Changing Geospace Mission

The Geospace Dynamics Observatory (GDO) mission observes the near-Earth region in space called Geospace with unprecedented resolution, scale and sensitivity. At a distance of 60 Earth Radii (Re) in a near-polar circular orbit and a approx. 27-day period, GDO images the earth's full disk with (1) a three-channel far ultraviolet imager, (2) an extreme ultraviolet imager of the plasmasphere, and (3) a spectrometer in the near to far ultraviolet range that probes any portion of the disk and simultaneously observes the limb. The exceptional capabilities of the GDO mission include (1) unprecedented improvement in signal to noise for globalscale imaging of Earth's space environment that enable changes in the Earth's space environment to be resolved with orders of magnitude higher in temporal and spatial resolution compared to existing data and other approaches, and (2) unrivaled capability for resolving the temporal evolution, over many days, in local time or latitude with a continuous view of Earth's global-scale evolution while simultaneously capturing the changes at scales smaller than are possible with other methods. This combination of new capabilities is a proven path to major scientific advances and discoveries. The GDO mission (1) has the first full disk imagery of the density and composition variability that exist during disturbed "storm" periods and the circulation systems of the upper atmosphere, (2) is able to image the ionosphere on a global and long time scale basis, (3) is able to probe the mechanisms that control the evolution of planetary atmospheres, and (4) is able to test our understanding of how the Earth is connected to the Sun. This paper explores the optical and technical aspects of the GDO mission and the implementation strategy. Additionally, the case will be made that GDO addresses a significant portion of the priority mission science articulated in the recent Solar and Space Physics Decadal Survey.

Spann, James↗

Measurement Report :Emission Factors of NH3 and NHx for Wildfires and Agricultural Fires in the United States

During the 2019 Fire Influence on Regional to Global Environments and Air Quality (FIREX-AQ) study, the NASA DC-8 carried out in situ chemical measurements in smoke plumes emitted from wildfires and agricultural fires in the contiguous United States. The DC-8 payload included a modified proton-transfer-reaction time-of-flight mass spectrometer (PTR-ToF-MS) for the fast measurement of gaseous ammonia (NH 3 ) and a high-resolution time-of-flight aerosol mass spectrometer (AMS) for the fast measurement of submicron particulate ammonium (NH 4 + ). We herein report data collected in smoke plumes emitted from 6 wildfires in the Western United States, 2 prescribed grassland fires in the Central United States, 1 prescribed forest fire in the Southern United States, and 66 small agricultural fires in the Southeastern United States. Smoke plumes contained double to triple digit ppb levels of NH 3 . In the wildfire plumes, a significant fraction of NH 3 had already been converted to NH 4 + at the time of sampling (≥2 h after emission). Substantial amounts of NH 4 + were also detected in freshly emitted smoke from corn and rice field fires. We herein present a comprehensive set of emission factors of NH 3 and NH x , with NH x = NH 3 + NH 4 + . Average NH 3 and NH x emission factors for wildfires in the Western United States were 1.86±0.75 g kg −1 and 2.47±0.80 g kg −1 of fuel burned, respectively. Average NH 3 and NH x emission factors for agricultural fires in the Southeastern United States were 0.89±0.58 and 1.74±0.92 g kg −1 , respectively. Our data show no clear inverse correlation between modified combustion efficiency (MCE) and NH 3 emissions. The observed NH 3 emissions were significantly higher than measured in previous laboratory experiments in the FIREX FireLab 2016 study.

Laura Tomsche↗

Measurement of the muon spin precession frequency using the straw tracking detectors at the Fermilab Muon g-2 experiment

The measurement of the anomalous magnetic dipole moment of the muon ($a_{\mu}$) has long stood as an excellent precision test of the Standard Model (SM). The Fermilab Muon g-2 experiment has recently finished data-taking and in July 2023 published its latest determination of $a_\mu$ with a world-leading precision of 0.2\,ppm. In this publication, it surpassed the systematic uncertainty goal defined in the TDR. The analyses of a dataset approximately four times larger than this recent publication is now underway. The principle measurement of the Muon g-2 experiment measures $a_{\mu}$ by taking the ratio of two frequencies; the anomalous precession frequency ($\omega_a$) and the muon-weighted magnetic field of the experiment's storage ring measured from the precession frequency of protons in water using nuclear magnetic resonance (NMR) probes. In all publications to date, $\omega_a$ has been determined using energy deposits in the 24 calorimeters. However, the Fermilab experiment has t wo straw tracker detectors measuring the time and momentum of charged particles which can in principle also be used to to measure $\omega_a$ and such a measurement can provide an invaluable cross-check of the calorimeter result with different, and reduced, systematic uncertainties. This thesis presents the first (blinded) determination of $\omega_a$ using just charged tracks from the straw tracking detectors as opposed to calorimeter energy deposits. This analysis was undertaken using the Run-2/3 dataset which represents approximately 25\% of the final dataset. A total uncertainty of 2.19\,ppm on $\omega_a$ was obtained which is dominated by the statistical uncertainty of 2.16\,ppm. Additionally two new methodologies important to the analysis of the straw tracking data have been developed: one to better determine the track arrival time ($t_0$) and one to determine the level of pileup in the tracking detectors. The new $t_0$ algorithm which incorporates angular information improves t he resolution on the determination of the $t_0$ by a factor of two and results in 19\% more tracks being successfully reconstructed. The data from the trackers is also used to determine the beam profile that weights the magnetic field in the determination of $a_\mu$ and in determining several of the systematic uncertainties in the calorimeter-based $\omega_a$ analysis. A detailed study of the impact of the internal alignment of the tracker, the $t_0$ and pileup on the determination of the beam position was undertaken and propagated through to an uncertainty in the $\omega_a$ determination. These uncertainties were used in the Fermilab Muon g-2 experiment's recent publication in Phys. Rev. Lett.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Recovering the phase and amplitude of X-ray FEL pulses using neural networks and differentiable models

Dynamics experiments are an important use-case for X-ray free-electron lasers (XFELs), but time-domain measurements of the X-ray pulses themselves remain a challenge. Shot-by-shot X-ray diagnostics could enable a new class of simpler and potentially higher-resolution pump-probe experiments. Here, we report training neural networks to combine low-resolution measurements in both the time and frequency domains to recover X-ray pulses at high-resolution. Critically, we also recover the phase, opening the door to coherent-control experiments with XFELs. The model-based generative neural-network architecture can be trained directly on unlabeled experimental data and is fast enough for real-time analysis on the new generation of MHz XFELs.

42 ENGINEERING↗

Deforestation and Secondary Growth in Rondonia, Brazil from SIR-C SAR and Landsat.SPOT data

Covers problems with existing data collected with high-resolution optical sensors. They say active microwave sensors could complement other sensors in getting through things like cloud cover. They analyzed SIR-C data in combination with Landsat TM data, a 9-year time series of SPOT XS data, and a preliminary field survey. They report findings and draw conclusions, including that SARs operating at long radar wavelengths, with both like and cross-polarizations, are needed for tropical deforestation studies.

high-resolution optical sensors SIR-C SAR SPOT↗

A study of marine stratocumulus using lidar and other FIRE aircraft observations

The National Center for Atmospheric Research (NCAR) airborne infrared lidar system (NAILS) used in the 1987 First ISCCP Regional Experiment (FIRE) off the coast of California is a 10.6 microns wavelength carbon dioxide lidar system constructed by Ron Schwiesow and co-workers at NCAR. The lidar is particularly well suited for detailed observations of cloud shapes; i.e., height of cloud top (when flying above cloud and looking down) and cloud base (when flying below cloud and looking up) along the flight path. A brief summary of the lidar design characteristics is given. The lidar height resolution of plus or minus 3 m allows for the distance between the aircraft and cloud edge to be determined with this accuracy; however, the duration of the emitted pulse is approximately 3 microseconds, which corresponds to a 500 m pulse length. Therefore, variations in backscatter intensities within the clouds can normally not be resolved. Hence the main parameter obtainable from the lidar is distance to cloud; in some cases the cloud depth can also be determined. During FIRE the lidar was operational on 7 of the 10 Electra flights, and data were taken when the distance between cloud and aircraft (minimum range) was at least 500 m. The lidar was usually operated at 8 Hz, which at a flight speed of 100 m s(-1) translates into a horizontal resolution of about 12 m. The backscatter as function of time (equivalent to distance) for each laser pulse is stored in digital form on magnetic tape. Currently, three independent variables are available to the investigators on the FIRE Electra data tapes: lidar range to cloud, strength of return (relative power), and pulse width of return, which is related to penetration depth.

Jensen, Jorgen B.↗

Measuring the neutrino-oxygen neutral current quasielastic cross section using the accelerator neutrino neutron interaction experiment

The Accelerator Neutrino Neutron Interaction Experiment (ANNIE) is a 26-ton gadolinium-doped water Cherenkov detector located on-axis to Fermilab’s Booster Neutrino Beam (BNB). ANNIE is uniquely positioned to perform high-statistics measurements of neutrino-nucleus interactions in water, benefiting from a large neutrino flux due to a short (100-meter) baseline. A central focus of ANNIE’s physics program is the measurement of both charged current (CC) and neutral current (NC) cross sections on water, including neutral current quasielastic (NCQE) and CC-inclusive channels. The NCQE measurement is particularly critical for constraining uncertainties in rare-event searches such as the Diffuse Supernova Neutrino Background (DSNB), where atmospheric $\nu$NCQE interactions constitute a significant and poorly constrained background. This dissertation presents a measurement of the flux-averaged neutrino-oxygen neutral current quasielastic ($\nu$NCQE) cross section using $2.573 \times 10^{20}$~POT of BNB exposure from the 2022 and 2023 beam years. The $\nu$NCQE interaction is identified through the primary $\gamma$-rays produced by nuclear de-excitation of the residual $^{15}$N$^*$ or $^{15}$O$^*$ nucleus following nucleon knockout from $^{16}$O. A dedicated Monte Carlo (MC) re-tuning campaign was conducted using an americium-beryllium (AmBe) calibration source, Michel electrons from stopped muons, and throughgoing dirt muons originating upstream of the detector. This multi-sample approach provided a wide-ranging $\mathcal{O}(\text{MeV})$--$\mathcal{O}(\text{GeV})$ dataset for tuning the simulated detector response, which was subsequently validated against AmBe neutron and Michel electron data for use in the $\nu$NCQE analysis. A dedicated laser calibration campaign was carried out to reduce timing uncertainties across the PMT system, enabling reconstruction of the BNB bunch substructure with sufficient resolution to serve as a background rejection tool. By selecting events in-time with individual neutrino bunches, beam-correlated $\nu$NCQE events are separated from diffuse and accelerator-induced backgrounds, notably skyshine neutrons and externally-originating events, that would otherwise dominate traditional charge-based selections within a small-scale, surface-level, short-baseline detector. A data-driven estimation of the skyshine neutron and external background rates was performed and incorporated into the systematic uncertainty budget. The flux-averaged $\nu$NCQE cross section on oxygen is measured to be $1.57 \pm 0.06\,(\text{stat.})$ $^{+0.91}_{-0.67}\,(\text{syst.})$ $\times 10^{-38}\ \text{cm}^{2}$. A full systematic budget is constructed by propagating uncertainties in the secondary hadronic interaction modeling, background cross section normalizations, detector response, neutrino flux, and the primary $\gamma$-ray emission probabilities from oxygen nuclear de-excitation. An idealized de-excitation model, constructed from existing measurements in the literature is developed to benchmark the predictions of the \textsc{GENIE} event generator. A comparison reveals that \textsc{GENIE} systematically overpredicts the primary $\gamma$-ray emission probability from oxygen de-excitation by a factor of $1.49\times$ for $E_\gamma > 6$~MeV and $3.07\times$ in the $3$--$6$~MeV band. This comparison motivates the dominant systematic uncertainty in this analysis, where a conservative uncertainty of $^{+39.9\%}_{-0\%}$ on the primary $\gamma$-ray signal prediction is assigned. The ANNIE result is consistent with and complementary to existing flux-averaged $\nu$NCQE cross section measurements from T2K and Super-Kamiokande, providing an independent measurement with a different detector, neutrino beam, and analysis methodology. Looking ahead, an upgrade to the ANNIE DAQ infrastructure enabling continuous extended readout will allow a complementary $\nu$NCQE neutron multiplicity measurement, directly relevant to constraining the NCQE background in DSNB searches, competitive with the recent T2K measurement at SK-Gd. The planned Super-SANDI upgrade, deploying a large Water-based Liquid Scintillator (WbLS) volume, will further extend ANNIE's reach to hadronic final states and exclusive NC channels, and enable joint measurements with liquid argon detectors sharing the BNB beamline ahead of DUNE and Hyper-Kamiokande.

Doran, Steven [Iowa State U.]↗

Humming Trains in Seismology: An Opportune Source for Probing the Shallow Crust

Abstract Seismologists are eagerly seeking new and preferably low-cost ways to map and track changes in the complex structure of the top few kilometers of the crust. By understanding it better, they can build on what is known regarding important, practical issues. These include telling us whether imminent earthquakes and volcanic eruptions are generating telltale underground signs of hazard, about mitigation of induced seismicity such as from deep injection of wastewater, how the Earth and its atmosphere couple, and where accessible natural resources are. Passive seismic imaging usually relies on blind correlations within extended recordings of Earth’s ceaseless “hum” or coda of well-mixed, small vibrations. In this article, we propose a complementary approach. It is seismic interferometry using opportune sources—specifically ones not stationary in time and moving in a well-understood configuration. Its interpretation relies on an accurate understanding of how these sources radiate seismic waves, precise timing, careful placement of pairs of listening stations, and seismic phase differentiation (surface and body waves). Massive freight trains were only recently recognized as such a persistent, powerful cultural (human activity-caused) seismic source. One train passage may generate a tremor with an energy output of a magnitude 1 earthquake and be detectable for up to 100 km from the track. We discuss the source mechanisms of train tremors and review the basic theory on sources. Finally, we present case studies of body- and surface-wave retrieval as an aid to mineral exploration in Canada and to monitoring of a southern California fault zone. We believe noise recovery from this new signal source, together with dense data acquisition technologies such as nodes or distributed acoustic sensing, will deeply transform our ability to monitor activity in the shallow crust at sharpened resolution in time and space.

Geochemistry & Geophysics↗

Augmented signal processing in Liquid Argon Time Projection Chambers with a deep neural network

The Liquid Argon Time Projection Chamber (LArTPC) is an advanced neutrino detector technology widely used in recent and upcoming accelerator neutrino experiments. It features a low energy threshold and high spatial resolution that allow for comprehensive reconstruction of event topologies. In current-generation LArTPCs, the recorded data consist of digitized waveforms on wires produced by induced signal on wires of drifting ionization electrons, which can also be viewed as two-dimensional (2D) (time versus wire) projection images of charged-particle trajectories. For such an imaging detector, one critical step is the signal processing that reconstructs the original charge projections from the recorded 2D images. For the first time, we introduce a deep neural network in LArTPC signal processing to improve the signal region of interest detection. By combining domain knowledge (e.g., matching information from multiple wire planes) and deep learning, this method shows significant improvements over traditional methods. This work details the method, software tools, and performance evaluated with realistic detector simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Fusion Approach for Remotely-Sensed Mapping of Agriculture (FARMA): A Scalable Open Source Method for Land Cover Monitoring Using Data Fusion

The increasing availability of very-high resolution (VHR; <2 m) imagery has the potential to enable agricultural monitoring at increased resolution and cadence, particularly when used in combination with widely available moderate-resolution imagery. However, scaling limitations exist at the regional level due to big data volumes and processing constraints. Here, we demonstrate the Fusion Approach for Remotely-Sensed Mapping of Agriculture (FARMA), using a suite of open source software capable of efficiently characterizing time-series field-scale statistics across large geographical areas at VHR resolution. We provide distinct implementation examples in Vietnam and Senegal to demonstrate the approach using WorldView VHR optical, Sentinel-1 Synthetic Aperture Radar, and Sentinel-2 and Sentinel-3 optical imagery. This distributed software is open source and entirely scalable, enabling large area mapping even with modest computing power. FARMA provides the ability to extract and monitor sub-hectare fields with multisensor raster signals, which previously could only be achieved at scale with large computational resources. Implementing FARMA could enhance predictive yield models by delineating boundaries and tracking productivity of smallholder fields, enabling more precise food security observations in low and lower-middle income countries.

fusion↗

Improving the SMAP Level-4 Soil Moisture Product

The NASA Soil Moisture Active Passive (SMAP) mission generates, among other data sets, the Level 4 Soil Moisture (L4 SM) product. The L4 SM algorithm is based on the assimilation of SMAP radiometer brightness temperature observations into the NASA Catchment land surface model using a spatially distributed ensemble Kalman filter. The L4 SM data are published with a mean latency of approx. 2.5 days from the time of observation and provide global, three-hourly, 9 km resolution estimates of surface and root-zone soil moisture and related land surface states and fluxes. In 2018, the product was upgraded from Version 3 to Version 4. Underlying the new version is a revised modeling system that includes improved input parameter datasets for land cover, topography, and vegetation height that are based on recent, high quality space-borne remote sensing observations. Land cover inputs were updated to the GlobCover2009 product, which is based on satellite observations from the Medium Resolution Imaging Spectrometer. Topographic statistics now rely on observations from the Shuttle Radar Topography Mission. Finally, vegetation height inputs are derived from space-borne Lidar measurements. Additionally, SMAP Level-2 soil moisture retrievals and in situ soil moisture measurements were used to calibrate a particular Catchment model parameter that governs the recharge of soil moisture from the models root-zone excess reservoir into the surface excess reservoir. Specifically, the replenishment of soil moisture near the surface from below under non-equilibrium conditions was substantially reduced, which brings the models surface soil moisture more in line with the SMAP Level 2 and in situ soil moisture. Finally, the calibration of the assimilated SMAP brightness temperatures changed substantially from Version 3 to Version 4. Considerable effort went into the version upgrade, creating an expectation that the new version is improved over the old version. Indeed, some aspects of the new version are clearly better. However, other aspects are not, and on balance, the overall improvement is modest at best. In this presentation we summarize the skill of the new and old versions vs. independent in situ measurements and in terms of data assimilation diagnostics, including, for example, the statistics of the (soil moisture) analysis increments and the observation minus forecast (brightness temperatures) residuals. We share our experience with trying to improve to the L4 SM product and the lessons learned from the effort.

Reichle, Rolf↗

An Examination of Parameters Affecting Large Eddy Simulations of Flow Past a Square Cylinder

Separated flow over a bluff body is analyzed via large eddy simulations. The turbulent flow around a square cylinder features a variety of complex flow phenomena such as highly unsteady vortical structures, reverse flow in the near wall region, and wake turbulence. The formation of spanwise vortices is often times artificially suppressed in computations by either insufficient depth or a coarse spanwise resolution. As the resolution is refined and the domain extended, the artificial turbulent energy exchange between spanwise and streamwise turbulence is eliminated within the wake region. A parametric study is performed highlighting the effects of spanwise vortices where the spanwise computational domain's resolution and depth are varied. For Re=22,000, the mean and turbulent statistics computed from the numerical large eddy simulations (NLES) are in good agreement with experimental data. Von-Karman shedding is observed in the wake of the cylinder. Mesh independence is illustrated by comparing a mesh resolution of 2 million to 16 million. Sensitivities to time stepping were minimized and sampling frequency sensitivities were nonpresent. While increasing the spanwise depth and resolution can be costly, this practice was found to be necessary to eliminating the artificial turbulent energy exchange.

Mankbadi, M. R.↗

Temporal resolution for crop discrimination estimated using J-M distance

A study is presented which examines techniques utilizing spectral/phenological data to classify and mensurate target crops at predetermined accuracies in order to estimate the temporal resolution requirements for a remote sensing system. An analysis is made of the physical target (crop) characteristics as they vary in time. A model is developed relating the frequency required to attain a particular probability of correct discrimination with N accumulated observations to that required with single observations. Finally, results for a multicrop spectral data set are shown to substantiate the theoretically-modeled results.

Baran, J. W.↗

The AgMIP GRIDded Crop Modeling Initiative (AgGRID) and the Global Gridded Crop Model Intercomparison (GGCMI)

Climate change is a significant risk for agricultural production. Even under optimistic scenarios for climate mitigation action, present-day agricultural areas are likely to face significant increases in temperatures in the coming decades, in addition to changes in precipitation, cloud cover, and the frequency and duration of extreme heat, drought, and flood events (IPCC, 2013). These factors will affect the agricultural system at the global scale by impacting cultivation regimes, prices, trade, and food security (Nelson et al., 2014a). Global-scale evaluation of crop productivity is a major challenge for climate impact and adaptation assessment. Rigorous global assessments that are able to inform planning and policy will benefit from consistent use of models, input data, and assumptions across regions and time that use mutually agreed protocols designed by the modeling community. To ensure this consistency, large-scale assessments are typically performed on uniform spatial grids, with spatial resolution of typically 10 to 50 km, over specified time-periods. Many distinct crop models and model types have been applied on the global scale to assess productivity and climate impacts, often with very different results (Rosenzweig et al., 2014). These models are based to a large extent on field-scale crop process or ecosystems models and they typically require resolved data on weather, environmental, and farm management conditions that are lacking in many regions (Bondeau et al., 2007; Drewniak et al., 2013; Elliott et al., 2014b; Gueneau et al., 2012; Jones et al., 2003; Liu et al., 2007; M¨uller and Robertson, 2014; Van den Hoof et al., 2011;Waha et al., 2012; Xiong et al., 2014). Due to data limitations, the requirements of consistency, and the computational and practical limitations of running models on a large scale, a variety of simplifying assumptions must generally be made regarding prevailing management strategies on the grid scale in both the baseline and future periods. Implementation differences in these and other modeling choices contribute to significant variation among global-scale crop model assessments in addition to differences in crop model implementations that also cause large differences in site-specific crop modeling (Asseng et al., 2013; Bassu et al., 2014).

climate↗

Differentiable physics-enabled closure modeling for Burgers’ turbulence

Abstract Data-driven turbulence modeling is experiencing a surge in interest following algorithmic and hardware developments in the data sciences. We discuss an approach using the differentiable physics paradigm that combines known physics with machine learning to develop closure models for Burgers’ turbulence. We consider the one-dimensional Burgers system as a prototypical test problem for modeling the unresolved terms in advection-dominated turbulence problems. We train a series of models that incorporate varying degrees of physical assumptions on an a posteriori loss function to test the efficacy of models across a range of system parameters, including viscosity, time, and grid resolution. We find that constraining models with inductive biases in the form of partial differential equations that contain known physics or existing closure approaches produces highly data-efficient, accurate, and generalizable models, outperforming state-of-the-art baselines. Addition of structure in the form of physics information also brings a level of interpretability to the models, potentially offering a stepping stone to the future of closure modeling.

97 MATHEMATICS AND COMPUTING↗