Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear Regression”

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,117 records · Page 62

Friction Velocity Estimation Using Dual-Frequency Altimeter Data

Global estimation of wind speeds near the sea-surface can be obtained using normalized radar cross sections sigma(o) from satellite-borne single-frequency altimeters. Indeed recent validation studies of Geosat altimeter-derived wind speeds show that several often cited semiempirical algorithms yield linear trends between wind-speeds from altimeter data and wind speeds from (a) ships of opportunity (Etcheto 1992), (b) buoys (Carter 1992), and (c) atmospheric models (Guillaume 1992). Although the linear trends are striking, there is considerable scatter since the range of standard deviations is from 1.5 to 3 m/s. In fact, the scatter of some particular measurements about the regressions is striking because they depart from the 'standard' by as much as +/- 5 m/s, which is large - even for full gale winds of 20 m/s. These anomalies can cause serious problems for gas exchange, heat, and aerosol estimates. They can also induce considerable errors in wind-driven oceanic circulation models and sea-surface wave models. This paper explores the possibility of improving the accuracy of wind speed estimation using dual-frequency altimeters.

Bliven, Larry F.↗

Full-Field Strain Measurement Integrated with Two Dimension Regression Analysis to Evaluate the Bi-Modulus Elastic Properties of Isotropic and Transversely Isotropic Materials

Background: Measuring the physical properties of shale is critical for optimizing engineering activities such as geothermal energy generation and hydraulic fracturing. Shale is a transversely isotropic material. Furthermore, this material can also include micro and macro cracks at different locations and orientations that cause it to behave differently under tensile or compressive loading. Objective: In this work, a combined experimental–numerical approach is proposed to evaluate the bi-modulus elastic properties of isotropic and transversely isotropic materials. Methods: Full-field strain measurements for a circular disk under diametral compression are integrated with a regression analysis technique to evaluate the elastic properties of bi-modulus materials subjected to tensile and compressive loads using two loading configurations on the same specimen. Digital Image Correlation (DIC) is used to measure the full-field strains. Subsequently, in the case of an isotropic material, a linear least-squares approach is utilized to process the experimentally determined strains in conjunction with analytical expressions of the stress fields (in terms of far-field loading) to determine the elastic modulus E, the shear modulus G, and the Poisson’s ratio $v$. In the case of a transversely isotopic material, such as shale, a finite element model is implemented to determine the stress fields (again in terms of far-field loading), which is followed by repeating the previous regression analysis in an iterative process to estimate the elastic parameters. Results: The results show that the proposed technique successfully provides a complete set of elastic properties as a function of both the loading condition and the principal material directions. The technique is validated by measurements on a known isotropic material and then applied to determine the properties of shale. Conclusion: In this work, the proposed approach is successfully used to calculate the bi-modulus elastic response of poly(methyl meth- acrylate) (PMMA) and shale. As expected, PMMA exhibits an isotropic response with no bi-modulus effect, however, shale exhibits both transverse isotropy and a bi-modulus effect. Therefore, this approach holds promise for investigating the elastic properties of materials like rocks and fiber-reinforced composite laminates as functions of the principal material directions and the loading conditions.

42 ENGINEERING↗

Quantile-specific heritability of sibling leptin concentrations and its implications for gene-environment interactions

“Quantile-dependent expressivity” occurs when the effect size of a genetic variant depends upon whether the phenotype (e.g., leptin) is high or low relative to its distribution. Leptin concentrations are strongly related to adiposity, whose heritability is quantile dependent. Whether inheritance of leptin concentrations is quantile dependent, and whether this explains the greater heritability in women than men in accordance with their greater adiposity, and explains other gene-environment interactions, remains to be determined. Therefore, leptin and leptin receptor concentrations from 3068 siblings in 1133 sibships from the Framingham Heart Study Third Generation Cohort were analyzed. Free leptin index (FLI) was calculated as the ratio of leptin to soluble leptin receptor concentrations. Full-sib (β FS ) regression slopes were robustly estimated by quantile regression with nonparametric significance assigned from 1000 bootstrap samples. The analyses showed β FS increased significantly with increasing percentiles of the offspring’s age- and sex-adjusted leptin distribution (P linear = 0.0001), which was accelerated at the higher concentrations (P quadratic = 0.0003). β FS at the 90th percentile (0.418 ± 0.066) was 4.7-fold greater than at the 10th percentile (0.089 ± 0.032, P difference = 3.6 × 10 –6 ). Consistent with quantile-dependent expressivity, the β FS was greater in female sibs, which was attributable to their higher leptin concentrations. Reported gene-environment interactions involving adiposity and LEP , LEPR , MnSOD , PPARγ , PPARγ2 , and IRS-1 polymorphisms were consistent with quantile-dependent expressivity of leptin concentrations. β FS for leptin receptor concentrations and free leptin index also increased significantly with increasing percentiles of their distributions (P linear = 0.04 and P linear = 8.5 × 10 –6 , respectively). In conclusion, inherited genetic and shared environmental effects on leptin concentrations were quantile dependent, which likely explains male–female differences in heritability and some gene-environment interactions.

59 BASIC BIOLOGICAL SCIENCES↗

In-Situ Measurement of Hall Thruster Erosion Using a Fiber Optic Regression Probe

One potential life-limiting mechanism in a Hall thruster is the erosion of the ceramic material comprising the discharge channel. This is especially true for missions that require long thrusting periods and can be problematic for lifetime qualification, especially when attempting to qualify a thruster by analysis rather than a test lasting the full duration of the mission. In addition to lifetime, several analytical and numerical models include electrode erosion as a mechanism contributing to enhanced transport properties. However, there is still a great deal of dispute over the importance of erosion to transport in Hall thrusters. The capability to perform an in-situ measurement of discharge channel erosion is useful in addressing both the lifetime and transport concerns. An in-situ measurement would allow for real-time data regarding the erosion rates at different operating points, providing a quick method for empirically anchoring any analysis geared towards lifetime qualification. Erosion rate data over a thruster's operating envelope would also be useful in the modeling of the detailed physics inside the discharge chamber. A recent fundamental sensor development effort has led to a novel regression, erosion, and ablation sensor technology (REAST). The REAST sensor allows for measurement of real-time surface erosion rates at a discrete surface location. The sensor was tested using a linear Hall thruster geometry, which served as a means of producing plasma erosion of a ceramic discharge chamber. The mass flow rate, discharge voltage, and applied magnetic field strength could be varied, allowing for erosion measurements over a broad thruster operating envelope. Results are presented demonstrating the ability of the REAST sensor to capture not only the insulator erosion rates but also changes in these rates as a function of the discharge parameters.

Polzink, Kurt A.↗

In-Situ Measurement of Hall Thruster Erosion Using a Fiber Optic Regression Probe

One potential life-limiting mechanism in a Hall thruster is the erosion of the ceramic material comprising the discharge channel. This is especially true for missions that require long thrusting periods and can be problematic for lifetime qualification, especially when attempting to qualify a thruster by analysis rather than a test lasting the full duration of the mission. In addition to lifetime, several analytical and numerical models include electrode erosion as a mechanism contributing to enhanced transport properties. However, there is still a great deal of dispute over the importance of erosion to transport in Hall thrusters. The capability to perform an in-situ measurement of discharge channel erosion is useful in addressing both the lifetime and transport concerns. An in-situ measurement would allow for real-time data regarding the erosion rates at different operating points, providing a quick method for empirically anchoring any analysis geared towards lifetime qualification. Erosion rate data over a thruster s operating envelope would also be useful in the modeling of the detailed physics inside the discharge chamber. There are many different sensors and techniques that have been employed to quantify discharge channel erosion in Hall thrusters. Snapshots of the wear pattern can be obtained at regular shutdown intervals using laser profilometry. Many non-intrusive techniques of varying complexity and sensitivity have been employed to detect the time-varying presence of erosion products in the thruster plume. These include the use quartz crystal microbalances, emission spectroscopy, laser induced flourescence, and cavity ring-down spectroscopy. While these techniques can provide a very accurate picture of the level of eroded material in the thruster plume, it is more difficult to use them to determine the location from which the material was eroded. Furthermore, none of the methods cited provide a true in-situ measure of erosion at the channel surface while the thruster is in operation (i.e. none yield a continuous channel erosion measurement). A recent fundamental sensor development effort has led to a novel regression, erosion, and ablation sensor technology (REAST). The REAST sensor allows for measurement of real-time surface erosion rates at a discrete surface location. The sensor was tested using a linear Hall thruster geometry (see Fig. 1), which served as a means of producing plasma erosion of a ceramic discharge chamber. The mass flow rate, discharge voltage, and applied magnetic field strength could be varied, allowing for erosion measurements over a broad thruster operating envelope. Results are presented demonstrating the ability of the REAST sensor to capture not only the insulator erosion rates but also changes in these rates as a function of the discharge parameters.

Polzin, Kurt↗

Materials representation and transfer learning for multi-property prediction

The adoption of machine learning in materials science has rapidly transformed materials property prediction. Hurdles limiting full capitalization of recent advancements in machine learning include the limited development of methods to learn the underlying interactions of multiple elements as well as the relationships among multiple properties to facilitate property prediction in new composition spaces. To address these issues, we introduce the Hierarchical Correlation Learning for Multi-property Prediction (H-CLMP) framework that seamlessly integrates: (i) prediction using only a material's composition, (ii) learning and exploitation of correlations among target properties in multi-target regression, and (iii) leveraging training data from tangential domains via generative transfer learning. The model is demonstrated for prediction of spectral optical absorption of complex metal oxides spanning 69 three-cation metal oxide composition spaces. H-CLMP accurately predicts non-linear composition-property relationships in composition spaces for which no training data are available, which broadens the purview of machine learning to the discovery of materials with exceptional properties. This achievement results from the principled integration of latent embedding learning, property correlation learning, generative transfer learning, and attention models. The best performance is obtained using H-CLMP with transfer learning [H-CLMP(T)] wherein a generative adversarial network is trained on computational density of states data and deployed in the target domain to augment prediction of optical absorption from composition. H-CLMP(T) aggregates multiple knowledge sources with a framework that is well suited for multi-target regression across the physical sciences.

36 MATERIALS SCIENCE↗

Quantitative Non-Destructive Evaluation of Fatigue Damage Based on Multi-Sensor Fusion

Based on sensor fusion and machine learning, this project developed a novel non-destructive evaluation (NDE) methodology, which consists of a remaining useful life (RUL) prediction framework and regression models for predicting residual stress and full width at half maximum (FWHM). A series of fatigue testing experiments were conducted using 5052-H32 aluminum alloy specimens. All specimens were measured using linear ultrasonic (LU) and nonlinear ultrasonic (NLU) testing methods non-destructively. Machine learning models were developed to use LU and NLU measurements to predict loading condition, fatigue level, residual stress, and FWHM. It was demonstrated that the developed methodology could distinguish new and fatigue specimens with an accuracy of 97.53%. Also, the prediction errors for residual stress and FWHM were as low as 4.73% and 1.62%, respectively. An interactive database was created to publicly share the data generated from the project. It is envisioned that the developed NDE technology will equip manufacturers with a responsive screening system for incoming used metallic components, and potentially lead to a significant increase in using used metallic components for remanufacturing.

42 ENGINEERING↗

Geographical Insights into Suicide Mortality Through Spatial Machine Learning

Suicide mortality is a leading cause of death in the United States, with an upward trend that emphasizes its significance as a public health issue. Previous research has employed global models like ordinary least squares (OLS) regression and local models such as geographically weighted regression (GWR). While local models are useful for analyzing spatial variations in suicide mortality, they share limitations with traditional global models, particularly about their inability to handle multi-collinearity and non-linear relationships. Machine learning approaches, like random forests (RF), can address some of these limitations but often fail to account for spatial variability. This gap highlights the need for spatial ML models specifically designed to tackle suicide mortality. This research seeks to fill this void by using a geographically weighted random forest model (GWRF) to examine the associations between county-level suicide mortality in the U.S. from 2010 to 2020 and various social and environmental determinants of health. A key aspect of our methodology is disciplined feature selection, which reduces the pool of explanatory variables by about 90%. This refinement enhances the explanatory power of both global (R2 improved from 0.59 to 0.67) and local (R2 improved from 0.64 to 0.67) RF models while reducing their run times. An analysis of the importance scores for these selected features reveals that the drivers of suicide mortality vary by context. Thus, to effectively address regional disparities and inform targeted public health interventions, a holistic approach that incorporates multiple county-level characteristics is essential.

Lebakula, Viswadeep [ORNL] (ORCID:0000000152935914↗

Horizontal vestibuloocular reflex evoked by high-acceleration rotations in the squirrel monkey. I. Normal responses

The horizontal angular vestibuloocular reflex (VOR) evoked by high-frequency, high-acceleration rotations was studied in five squirrel monkeys with intact vestibular function. The VOR evoked by steps of acceleration in darkness (3,000 degrees /s(2) reaching a velocity of 150 degrees /s) began after a latency of 7.3 +/- 1.5 ms (mean +/- SD). Gain of the reflex during the acceleration was 14.2 +/- 5.2% greater than that measured once the plateau head velocity had been reached. A polynomial regression was used to analyze the trajectory of the responses to steps of acceleration. A better representation of the data was obtained from a polynomial that included a cubic term in contrast to an exclusively linear fit. For sinusoidal rotations of 0.5-15 Hz with a peak velocity of 20 degrees /s, the VOR gain measured 0.83 +/- 0.06 and did not vary across frequencies or animals. The phase of these responses was close to compensatory except at 15 Hz where a lag of 5.0 +/- 0.9 degrees was noted. The VOR gain did not vary with head velocity at 0.5 Hz but increased with velocity for rotations at frequencies of >/=4 Hz (0. 85 +/- 0.04 at 4 Hz, 20 degrees /s; 1.01 +/- 0.05 at 100 degrees /s, P < 0.0001). No responses to these rotations were noted in two animals that had undergone bilateral labyrinthectomy indicating that inertia of the eye had a negligible effect for these stimuli. We developed a mathematical model of VOR dynamics to account for these findings. The inputs to the reflex come from linear and nonlinear pathways. The linear pathway is responsible for the constant gain across frequencies at peak head velocity of 20 degrees /s and also for the phase lag at higher frequencies being less than that expected based on the reflex delay. The frequency- and velocity-dependent nonlinearity in VOR gain is accounted for by the dynamics of the nonlinear pathway. A transfer function that increases the gain of this pathway with frequency and a term related to the third power of head velocity are used to represent the dynamics of this pathway. This model accounts for the experimental findings and provides a method for interpreting responses to these stimuli after vestibular lesions.

NASA Program Biomedical Research and Countermeasur↗

Winter Wheat Yield Assessment from Landsat 8 and Sentinel-2 Data: Incorporating Surface Reflectance, Through Phenological Fitting, into Regression Yield Models

A combination of Landsat 8 and Sentinel-2 offers a high frequency of observations (3–5 days) at moderate spatial resolution (10–30 m), which is essential for crop yield studies. Existing methods traditionally apply vegetation indices (VIs) that incorporate surface reflectances (SRs) in two or more spectral bands into a single variable, and rarely address the incorporation of SRs into empirical regression models of crop yield. In this work, we address these issues by normalizing satellite data (both VIs and SRs) derived from NASA’s Harmonized Landsat Sentinel-2 (HLS) product, through a phenological fitting. We apply a quadratic function to fit VIs or SRs against accumulated growing degree days (AGDDs), which affects the rate of crop development. The derived phenological metrics for VIs and SRs, namely peak, area under curve (AUC), and fitting coefficients from a quadratic function, were used to build empirical regression winter wheat models at a regional scale in Ukraine for three years, 2016–2018. The best results were achieved for the model with near infrared (NIR) and red spectral bands and derived AUC, constant, linear, and quadratic coefficients of the quadratic model. The best model yielded a root mean square error (RMSE) of 0.201 t/ha (5.4%) and coefficient of determination R2 = 0.73 on cross-validation.

phenological fitting↗

Integrating Reservoirs into the Dissolved Organic Matter Versus Primary Production Paradigm: How Does Chlorophyll-$a$ Change Across Dissolved Organic Carbon Concentrations in Reservoirs?

Primary production in freshwater ecosystems is largely a function of light and nutrient availability, both of which have been changing in many lakes and reservoirs in response to anthropogenic pressures. Recent studies focusing on natural lakes have found a hump-shaped response of primary production (sometimes measured as chlorophyll-$a$) to dissolved organic matter (DOM, measured as dissolved organic carbon, DOC), which has both light-absorbing chromophoric properties and DOM-bound nutrients. We used the United States National Lakes Assessment dataset to integrate reservoirs into this paradigm in comparison with natural lakes and assessed the relative differences in the predicted response’s model structure, regression parameter values, and drivers of the chlorophyll-$a$ residuals. We found that chlorophyll-$a$ in reservoirs exhibited a hump-shaped response to DOC, while natural lakes from this dataset were better fit with a linear response, differing from previous studies focused on boreal lakes. Despite this, reservoirs had a greater maximum chlorophyll-a response compared to natural lakes in this study (45.5 versus 33.8 μg L -1 ), which occurred at a lower DOC concentration threshold (18.3 versus 26.4 mg L -1 ) when compared using quadratic models. Reservoirs had lower median light:nutrient values compared to natural lakes, and greater median surface area and total phosphorus (TP), that can all influence the light environment and the peak chlorophyll-a responses. In both reservoirs and natural lakes, chlorophyll-$a$ residuals were most strongly influenced by TP, where TP < 25-30 µg L -1 suppressed chlorophyll-a residuals and higher TP amplified them. Light:nutrient values were somewhat important predictors, and patterns with chlorophyll-$a$ residuals supported previous work showing low light:nutrient values amplified chlorophyll-$a$ responses and higher values suppressed them. In conclusion, quantifying the shape of the response of primary production to DOM quantity and quality as well as the drivers of the residuals, namely TP for lakes and reservoirs in this dataset, will be important for understanding the effects that changes in water quality may have on primary production and freshwater ecosystem processes.

59 BASIC BIOLOGICAL SCIENCES↗

MODIS Tree Cover Validation for the Circumpolar Taiga-Tundra Transition Zone

A validation of the 2005 500m MODIS vegetation continuous fields (VCF) tree cover product in the circumpolar taiga-tundra ecotone was performed using high resolution Quickbird imagery. Assessing the VCF's performance near the northern limits of the boreal forest can help quantify the accuracy of the product within this vegetation transition area. The circumpolar region was divided into longitudinal zones and validation sites were selected in areas of varying tree cover where Quickbird imagery is available in Google Earth. Each site was linked to the corresponding VCF pixel and overlaid with a regular dot grid within the VCF pixel's boundary to estimate percent tree crown cover in the area. Percent tree crown cover was estimated using Quickbird imagery for 396 sites throughout the circumpolar region and related to the VCF's estimates of canopy cover for 2000-2005. Regression results of VCF inter-annual comparisons (2000-2005) and VCF-Quickbird image-interpreted estimates indicate that: (1) Pixel-level, inter-annual comparisons of VCF estimates of percent canopy cover were linearly related (mean R(sup 2) = 0.77) and exhibited an average root mean square error (RMSE) of 10.1 % and an average root mean square difference (RMSD) of 7.3%. (2) A comparison of image-interpreted percent tree crown cover estimates based on dot counts on Quickbird color images by two different interpreters were more variable (R(sup 2) = 0.73, RMSE = 14.8%, RMSD = 18.7%) than VCF inter-annual comparisons. (3) Across the circumpolar boreal region, 2005 VCF-Quickbird comparisons were linearly related, with an R(sup 2) = 0.57, a RMSE = 13.4% and a RMSD = 21.3%, with a tendency to over-estimate areas of low percent tree cover and anomalous VCF results in Scandinavia. The relationship of the VCF estimates and ground reference indicate to potential users that the VCF's tree cover values for individual pixels, particularly those below 20% tree cover, may not be precise enough to monitor 500m pixel-level tree cover in the taiga-tundra transition zone.

Montesano, P. M.↗

An ecologic study of dietary links to prostate cancer

BACKGROUND: The etiology of prostate cancer has not been fully resolved in the scientific and medical literature, although the non-fat portion of milk and calcium are emerging as leading dietary risk factors, with lycopene (found in tomatoes) and vitamin D apparently being risk reduction factors. METHODS: The ecologic (multi-country statistical) approach is used to study dietary links to prostate cancer. Mortality data from 1986 for various age groups in 41 countries are compared with national consumer macronutrient supply values for 1983 and tomato supply values for 1985. RESULTS: For 28 countries with more than five Kcal/day of tomatoes in the consumer supply, a linear combination of non-fat milk (risk factor) and tomatoes (risk reduction factor) was found to have the highest statistical association with prostate cancer mortality rates for men over the age of 35, with the Pearson regression coefficient (R2) for those aged 65-74 years = 0.67 and p < 0.001. For the 13 countries with fewer than six Kcal/day of tomatoes, non-fat milk had the highest association (R2 = 0.92, p < 0.001 for men aged 65-74 years). For 41 countries combined, the non-fat portion of milk had the highest association with prostate cancer mortality rates (R2 = 0.73, p < 0.001 for men aged 65-74 years). CONCLUSIONS: These results support the results of several cohort studies which found the non-fat portion of milk to have the highest association with prostate cancer, likely due to the calcium, and tomatoes to reduce the risk of prostate cancer, most likely due to lycopene.

Diet/adverse effects↗

MIONet: Learning Multiple-Input Operators via Tensor Product

As an emerging paradigm in scientific machine learning, neural operators aim to learn operators, via neural networks, that map between infinite-dimensional function spaces. Several neural operators have been recently developed. However, all the existing neural operators are only designed to learn operators defined on a single Banach space; i.e., the input of the operator is a single function. Here, for the first time, we study the operator regression via neural networks for multiple-input operators defined on the product of Banach spaces. We first prove a universal approximation theorem of continuous multiple-input operators. We also provide a detailed theoretical analysis including the approximation error, which provides guidance for the design of the network architecture. Based on our theory and a low-rank approximation, we propose a novel neural operator, MIONet, to learn multiple-input operators. MIONet consists of several branch nets for encoding the input functions and a trunk net for encoding the domain of the output function. Here, we demonstrate that MIONet can learn solution operators involving systems governed by ordinary and partial differential equations. In our computational examples, we also show that we can endow MIONet with prior knowledge of the underlying system, such as linearity and periodicity, to further improve accuracy.

97 MATHEMATICS AND COMPUTING↗

Data Driven Optimization Framework for Volt-age Regulation in Distribution Systems

This letter proposes a data-driven optimization framework for voltage regulation problems to address the challenge of model inaccuracy and parameter varying. To achieve online voltage optimization, the recursive kernel regression and interior point methods are integrated. The IEEE 123-Bus system and EPRI Ckt5 feeder are selected to validate the effectiveness of the proposed data-driven optimization framework. The proposed method is also compared with a linear function based method.

Hong, Tianqi↗

Heuristic Computational Model for Predicting Lignin Solubility in Tailored Organic Solvents

Lignin is a random heteropolymer that has been extensively studied as a renewable source of aromatic precursors for high-value chemicals, biofuels, and bioplastics. A key challenge in lignin valorization is the structural and compositional heterogeneity of lignin feedstocks. Solvent-based approaches are commonly used to fractionate lignin to reduce this heterogeneity, but solvent selection can be challenging due to variability in lignin composition. In this work, we developed computational methods to predict good and poor organic solvents as a function of lignin composition. We analyzed 28 different linear pentamer structures, 18 from known libraries and 10 hypothetical polymers, and calculated their activity coefficients in 50 different organic solvents by using the conductor-like screening model for realistic solvents. We used these data to train a regression model that enabled the extensive investigation of the impact of solvent and monolignol compositions on predicted lignin solubility. The exhaustive exploration of solubility trends using model predictions revealed sets of solvents, identified using Kamlet–Taft parameters, that are predicted to promote lignin dissolution regardless of lignin composition. We further identified solvents expected to selectively isolate lignin fractions enriched in certain subunits. Furthermore, these results establish heuristic guidelines for solvent selection that can be used to tailor fractionation processes for lignin feedstocks of distinct composition or to design new processes that isolate fractions with higher proportions of selected subunits.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Survey of Constrained Gaussian Process: Approaches and Implementation Challenges

Gaussian process regression is a popular Bayesian framework for surrogate modeling of expensive data sources. As part of a larger effort in scientific machine learning, many recent works have incorporated physical constraints or other a priori information within Gaussian process regression to supplement limited data and regularize the behavior of the model. We provide an overview and survey of several classes of Gaussian process constraints, including positivity or bound constraints, monotonicity and convexity constraints, differential equation constraints provided by linear PDEs, and boundary condition constraints. We compare the strategies behind each approach as well as the differences in implementation, concluding with a discussion of the computational challenges introduced by constraints.

97 MATHEMATICS AND COMPUTING↗

[Re] Drivers of evapotranspiration from boreal wildfires

Computational reproducibility is a difficult challenge across science. I attempted to use R 3.6.1 to reproduce linear model fits, done originally using v2.6.0 for a 2009 paper on the drivers of large-scale forest evapotranspiration after wildfire. Model outputs were largely identical, aside from minor formatting changes, except for one–out of 12 total– regression in which the median residual value changed very slightly (in the sixth decimal place). I suggest that this essentially successful reproducibility is due to the relative simplicity of the script, its use of only base R functions, and R’s historically conservative approach to breaking changes.

54 ENVIRONMENTAL SCIENCES↗