Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “significant quantity”

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,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Influence of Vane Trailing Edge Flow on the Formation of Cavity Cells and Rim Sealing

Abstract Hot gas ingestion into the turbine rim seal cavity is an important concern for engine designers. To prevent ingestion, rim seals use high-pressure purge flow; however, the penalty is that excessive use of the purge flow decreases engine thermal efficiency. In this paper, a one-stage turbine operating at engine-representative conditions was used to study the effect of steady and time-resolved under-platform cavity temperatures and pressures across a range of coolant flowrates in the presence of vane trailing edge (VTE) flow. This study correlates time-resolved pressure with time-resolved temperature to identify primary frequencies driving ingestion. At certain flowrates, the time-resolved pressures are out of phase with the temperatures, indicating ingestion. Measurements from high-frequency response pressure sensors in the rim seal and vane platform were also used to determine rotational speed and quantity of large-scale structures (cells). In a parallel effort, a computational model using Unsteady Reynolds-averaged Navier–Stokes (URANS) was applied to determine swirl ratio in the rim seal cavity and time-resolved rim sealing effectiveness. The experimental results confirm that at low purge flowrates, the VTE flow influences the unsteady flow field by decreasing pressure unsteadiness in the rim seal cavity. Results show an increase in purge flow increases the number of unsteady large-scale structures in the rim seal and decreases their rotational speed. However, VTE flow was shown to not significantly change the cell speed and count in the rim seal. Simulations point to the importance of the large-scale cell structures in influencing rim sealing unsteadiness, which is not captured in current rim sealing predictive models.

Engineering↗
DOE OSTI2022

Enabling High-Energy-Density Physics Using X-Ray Free Electron Lasers

New High-Energy-Density (HED) Physics research facilities, such as x-ray free electron lasers (XFEL) are becoming available to the research community. These new tools can be used to address stockpile stewardship science questions in the HED regime with greater accuracy and in parameter ranges not previously accessible. Two important challenges that must be overcome before using the XFEL are developing simulation tools needed to design the experiment and developing large quantities of inexpensive targets from which the HED state is created. We developed two-dimensional simulations in a new geometry (planar) to determine basic plasma parameters that could be expected in a notional experiment. We developed a postprocessor that used the results of the simulations to predict the x-ray Thomson scattering (XRTS) signal that yield the plasma density and pressure at a specified location and time. Finally, we developed several possible target geometries that are inexpensive and easy to use. With these tools in hand, the risk in designing an experiment for an XFEL is significantly reduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1974

An empirical analysis of the quantitative effect of data when fitting quadratic and cubic polynomials

A study is made of the extent to which the size of the sample affects the accuracy of a quadratic or a cubic polynomial approximation of an experimentally observed quantity, and the trend with regard to improvement in the accuracy of the approximation as a function of sample size is established. The task is made possible through a simulated analysis carried out by the Monte Carlo method in which data are simulated by using several transcendental or algebraic functions as models. Contaminated data of varying amounts are fitted to either quadratic or cubic polynomials, and the behavior of the mean-squared error of the residual variance is determined as a function of sample size. Results indicate that the effect of the size of the sample is significant only for relatively small sizes and diminishes drastically for moderate and large amounts of experimental data.

Canavos, G. C.↗
DOE OSTI2022

Preconditioned least‐squares Petrov–Galerkin reduced order models

Abstract In this article, we introduce a methodology for improving the accuracy and efficiency of reduced order models (ROMs) constructed using the least‐squares Petrov–Galerkin (LSPG) projection method through the introduction of preconditioning. Unlike prior related work, which focuses on preconditioning the linear systems arising within the ROM numerical solution procedure to improve linear solver performance, our approach leverages a preconditioning matrix directly within the minimization problem underlying the LSPG formulation. Applying preconditioning in this way has the potential to improve ROM accuracy for several reasons. First, preconditioning the LSPG formulation changes the norm defining the residual minimization, which can improve the residual‐based stability constant bounding the ROM solution's error. The incorporation of a preconditioner into the LSPG formulation can have the additional effect of scaling the components of the residual being minimized to make them roughly of the same magnitude, which can be beneficial when applying the LSPG method to problems with disparate scales (e.g., dimensional equations, multi‐physics problems). Importantly, we demonstrate that an “ideal preconditioned” LSPG ROM (a ROM in which the preconditioner is the inverse of the Jacobian of its corresponding full order model) emulates projection of the full order model solution increment onto the reduced basis. This quantity defines a lower bound on the error of a ROM solution for a given reduced basis. By designing preconditioners that approximate the Jacobian inverse—as is common in designing preconditioners for solving linear systems—it is possible to obtain a ROM whose error approaches this lower bound. The proposed approach is evaluated on several mechanical and thermo‐mechanical problems implemented within the Albany HPC code and run in the predictive regime, with prediction across material parameter space. We demonstrate numerically that the introduction of simple Jacobi, Gauss‐Seidel, and ILU preconditioners into the proper orthogonal decomposition/LSPG formulation reduces significantly the ROM solution error, the reduced Jacobian condition number, the number of nonlinear iterations required to reach convergence, and the wall time (thereby improving efficiency). Moreover, our numerical results reveal that the introduction of preconditioning can deliver a robust and accurate solution for test cases in which the unpreconditioned LSPG method fails to converge.

Lindsay, Payton↗
DOE OSTI2025

Exploring diversion-pathway analysis of a generic molten-salt fast reactor using multiphysics informed signatures

Molten salt reactors are being explored by multiple commercial ventures due to their inherent safety features, flexibility in fuel sources, and high fuel utilization and thermal efficiency. The continual flow of fuel salt, large fissile quantities present, and ability to add or divert material due to the liquid nature introduces new challenges for international safeguards. To understand how international safeguards should be applied, it is important to capture the inherent multi-physics nature of a molten salt reactor. This work examines a generic molten salt fast reactor to understand how potential diversion scenarios would affect the concentration of radionuclides in the primary and auxiliary systems. Three types of diversion were examined: a slow drip of fuel salt, gaseous plutonium extraction, and uranium metal plating. The analysis determined that several key isotopes become statistically significant once diversion begins, indicating that detection of such diversion cases would be possible through measuring specific signatures such as gamma spectra.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A simple soil mass correction for a more accurate determination of soil carbon stock changes

Abstract Agricultural soils can act as a sink for large quantities of soil organic carbon (SOC) but can also be sources of carbon to the atmosphere. The international standard for assessing SOC stock and measuring stock change stipulates fixed depth sampling to at least 30 cm. The tendency of bulk density (BD) to decrease with decreasing disturbance and increasing SOC concentration and the assumption of constant SOC and BD within this depth profile promotes error in the estimates of SOC stock. A hypothetical but realistic change in BD from 1.5 to 1.1 g cm −3 from successive fixed depth sampling to 30 cm underestimates SOC stock change by 17%. Significant effort has been made to evaluate and reduce this fixed depth error by using the equivalent soil mass (ESM) approach, but with limited adoption. We evaluate the error in SOC stock assessment and change generated from fixed depth measurements over time relative to the ESM approach and propose a correction that can be readily adopted under current sampling and analytical methods. Our approach provides a more accurate estimate of SOC stock accumulation or loss that will help incentivize management practice changes that reduce the environmental impacts of agriculture and further legitimize the accounting practices used by the emerging carbon market and organizations that have pledged to reduce their supply chain greenhouse gas (GHG) footprints.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2018

Elliptic-Cylindrical Analytical Flux Rope Model for Magnetic Clouds

In this paper, we present the elliptic-cylindrical analytical flux rope model, which constitutes the first level of complexity above that of a circular-cylindrical geometry. The framework of this series of models was established by Nieves-Chinchilla et al. with the circular-cylindrical analytical flux rope model. The model describes the magnetic flux rope topology with distorted cross section as a possible consequence of the flux rope interaction with the solar wind. In this model, for the first time, a flux rope is completely described by a nonorthogonal geometry. The Maxwell equations can be consistently solved using tensorial analysis, and relevant physical quantities can be derived, such as magnetic fluxes, number of turns, or Lorentz force distribution. The model is generalized in termsof the radial dependence of the poloidal and axial current density components. The circular-cylindrical reconstruction technique has been adapted to the new geometry for a specific case of the model and tested againstan interplanetary coronal mass ejection observed by the Wind spacecraft on 2005 June 12. In this specific case,from the comparative analysis between the circular-cylindrical and elliptic-cylindrical models, the inclusion ofthe cross-section distortion in the 3D reconstruction results in significant changes in the derived axis orientation,size, central magnetic field, magnetic fluxes, and force-freeness. The case studied in this paper exemplifies the useof the model and reconstruction technique developed. Furthermore, the novel mathematical formulation to modelflux ropes in heliophysics paves the way to the inclusion of more complex magnetic field configurations.

CMEs↗
DOE OSTI2023

Order Parameter Engineering for Random Systems

The chemical short-range order (CSRO) in crystalline materials influences the properties, and its effect is significant in the context of multicomponent materials. Here, we propose a scheme for the CSRO parameter or Δ-parameter in terms of the number of like and unlike bonds in the multicomponent systems. The OPERA or Order Parameter Engineering for RAndom Systems scheme for semi-canonical and canonical ensembles is proposed. The proposed framework of Δ-parameter with OPERA framework can generate the single-phase supercell with desired CSRO without explicit energy calculations and provides a computationally efficient scheme for exploration of the CSRO. We demonstrate the applicability of the Δ-parameter as a scalar quantity for describing the CSRO in multicomponent alloys and oxides (FCC-CoCrNi, BCC-MoNbTaW, and (CoCuMgNiZn)O).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Developments on two lithium vapor-box linear test-stand experiments

The lithium vapor-box divertor is a possible fusion power exhaust solution.It uses condensation pumping to create a gradient of vapor density in a divertor slot; this should allow a stable detachment front without active feedback.As initial explorations of the concept, two test stands which take the form of three connected cylindrical stainless steel boxes are being developed: one without plasma at PPPL, to test models of lithium evaporation and flow; and one for the linear plasma device Magnum-PSI (at DIFFER in Eindhoven, The Netherlands) to test the ability of a lithium vapor cloud to induce volumetric detachment and redistribute the plasma power.The first experiment uses boxes with diameters of 6 cm, joined by apertures with diameters of 2.2 cm. Up to 1 g of Li is placed in one box, which is heated to up to 600 degrees C. The Li evaporates, then flows to and condenses in the two other, cooler boxes over several minutes. The quantity of Li transported is assessed by weighing the boxes before and after the heating cycle, and is compared to the quantity predicted to flow for the box at its measured temperature using a Direct Simulation Monte Carlo code, SPARTA. With good experimental conditions, the two values agree to within 15%.The experiment on Magnum-PSI is in the conceptual design stage.The design is assessed by simulations using the code B2.5-Eunomia.They show that when the hydrogen-ion plasma beam, with n_e = 4e20 per cubic meter, T_e = 1.5 eV, and r = 1 cm, is passed through a 16 cm long, 12 Pa, 625 degree C Li vapor cloud, the plasma heat flux and pressure on the target are significantly reduced compared to the case without Li.With the Li present, the plasma is cooled by excitation of Li neutrals followed by radiation until it volumetrically recombines, lowering the heat flux from 3.7 MW/m^2 to 0.13 MW/m^2, and the pressure is reduced by 93%, largely by collisions of hydrogen ions with neutral Li.

divertor↗
NASA NTRS1982

A study of dirigibles for use in the Peruvian Selva Central Region

The potential for dirigibles as transports in the Selva Central region of Peru was evaluated by means of a mission and economic analysis. A total requirement to transport over 19 million tons (t) of agricultural produce, lumber, and meat was projected by the year 2004. A primary route involving zones for loading and delivering this cargo was identified. Although dirigibles are capable of short field operation, all existing airfields must be enlarged in width to allow for all conditions of wind and weather and to provide space for overnight mooring. A maintenance base and operations headquarters, complete with hangar and other service facilities would be required. The quantities of cargo to be carried establish requirements for fleets of dirigibles of various sizes and capacities. Cargo capacities of 5 to 100 tons were identified. Fleet sizes up to 106 dirigibles (in 20 t capacities) would be required. Dirigibles were assumed to be of the nonrigid type except in the 100 t category for which rigid characteristics were assumed. A method of determining dirigible costs was developed. The values derived were then applied to an economic analysis to determine initial investment and operating costs. It was found that larger dirigibles of approximately 20 t capacities or higher could offer significant cost benefits over airplanes, provided cruise speeds were higher than 100 km/hr.

Mayer, N. J.↗
DOE OSTI2020

Iron Heterogeneity in Early Active Multiple Sclerosis Lesions

Objective: Multiple sclerosis (MS) is a heterogeneous inflammatory demyelinating disease. Iron distribution is altered in MS patients' brains, suggesting iron liberation within active lesions amplifies demyelination and neurodegeneration. Whether the amount and distribution of iron are similar or different among different MS immunopatterns is currently unknown. Methods: We used synchrotron X-ray fluorescence imaging, histology, and immunohistochemistry to compare the iron quantity and distribution between immunopattern II and III early active MS lesions. We analyzed archival autopsy and biopsy tissue from 21 MS patients. Results: Immunopattern II early active lesions contain 64% more iron (95% confidence interval [CI] = 17-127, p = 0.004) than immunopattern III lesions, and 30% more iron than the surrounding periplaque white matter (95% CI = 3–64%, p = 0.03). Iron in immunopattern III lesions is 28% lower than in the periplaque white matter (95% CI = -40 to -14%, p < 0.001). When normalizing the iron content of early active lesions to that of surrounding periplaque white matter, the ratio is significantly higher in immunopattern II (p < 0.001). Microfocused X-ray fluorescence imaging shows that iron in immunopattern II lesions localizes to macrophages, whereas macrophages in immunopattern III lesions contain little iron. Interpretation: Iron distribution and content are heterogeneous in early active MS lesions. Iron accumulates in macrophages in immunopattern II, but not immunopattern III lesions. This heterogeneity in the two most common MS immunopatterns may be explained by different macrophage polarization, origin, or different demyelination mechanisms, and paves the way for developing new or using existing iron-sensitive magnetic resonance imaging techniques to differentiate among immunopatterns in the general nonbiopsied MS patient population.

60 APPLIED LIFE SCIENCES↗
NASA NTRS2010

Three Methods of Detection of Hydrazines

Three proposed methods for measuring trace quantities of hydrazines involve ionization and detection of hydrazine derivatives. These methods are intended to overcome the limitations of prior hydrazine- detection methods. Hydrazine (Hz), monomethylhydrazine (MMH), and unsymmetrical dimethylhydrazine (UDMH) are hypergolic fuels and are highly reactive, toxic, and corrosive. A capability to measure concentrations of hydrazines is desirable for detecting leaks and ensuring safety in aerospace settings and in some industrial settings in which these compounds are used. One of the properties (high reactivity) that make it desirable to detect trace amounts of hydrazines also makes it difficult to detect hydrazines and measure their concentrations accurately using prior methods: significant amounts are lost to thermal and catalytic decomposition prior to detection. Further complications arise from the sticky nature of hydrazines: Sample hydrazine molecules tend to become irreversibly adsorbed onto solid surfaces with which they come into contact during transport to detectors, giving rise to drift in detector responses. In each proposed method, the reactive, sticky nature of hydrazines would be turned to advantage by providing a suitably doped substrate surface with which the hydrazines would react. The resulting hydrazine derivatives would be sufficiently less sticky and sufficiently more stable so that fewer molecules would be lost to decomposition or adsorption during transport. Consequently, it would be possible to measure concentration with more sensitivity and less error than in prior techniques. The first proposed method calls for the use of a recently developed technique known as desorption electrospray ionization (DESI), in which a pneumatically assisted micro -electrospray at ambient pressure is directed at a surface of interest. In this case, the surface of interest would be that of a substrate described above.

Griffin, Timothy↗
NASA NTRS2010

Psychrophilic Biomass Producers in the Trophic Chain of the Microbial Community of Lake Untersee, Antarctica

The study of photosynthetic microorganisms from the Lake Untersee samples showed dispersed distribution of phototrophs within ~80 m water column. Lake Untersee represents a unique ecosystem that experienced complete isolation: sealed by the Anuchin Glacier for many millennia. Consequently, its biocenosis has evolved over a significant period of time without exchange or external interaction with species from other environments. The major producers of organic matter in Lake Untersee are represented by phototrophic and chemolithotrophic microorganisms. This is the traditional trophic scheme for lacustrine ecosystems on Earth. Among the phototrophs, diatoms were not found, which differentiates this lake from other known ecosystems. The dominant species among phototrophs was Chlamydomonas sp. with typical morphostructure: green chloroplasts, bright red round spot, and two polar flagella near the opening. As expected, the physiology of studied phototrophs was limited by low temperature, which defined them as obligate psychrophilic microorganisms. By the quantity estimation of methanogenesis in this lake, the litho-autotrophic production of organic matter is competitive with phototrophic production. However, pure cultures of methanogens have not yet been obtained. We discuss the primary producers of organic matter and the participation of our novel psychrophilic homoacetogen into the litho-autotrophic link of biomass production in Lake Untersee.

Pikuta, Elena V.↗
DOE OSTI2022

A hybrid nodal-staggered pseudo-spectral electromagnetic particle-in-cell method with finite-order centering

Electromagnetic particle-in-cell (PIC) codes are widely used to perform computer simulations of a variety of physical systems, including fusion plasmas, astrophysical plasmas, plasma wakefield particle accelerators, and secondary photon sources driven by ultra-intense lasers. In a PIC code, Maxwell's equations are solved on a grid with a numerical method of choice. This article focuses on pseudo-spectral analytical time-domain (PSATD) algorithms and presents a novel hybrid PSATD PIC scheme that combines the respective advantages of standard nodal and staggered methods. The novelty of the hybrid scheme consists in using finite-order centering of grid quantities between nodal and staggered grids, in order to combine the solution of Maxwell's equations on a staggered grid with the deposition of charges and currents and the gathering of electromagnetic forces on a nodal grid. The correctness and performance of the novel hybrid scheme are assessed by means of numerical tests that employ different classes of PSATD equations in a variety of physical scenarios, ranging from the modeling of electron-positron pair creation in vacuum to the simulation of laser-driven and particle beam-driven plasma wakefield acceleration. It is shown that the novel hybrid scheme offers significant numerical and computational advantages, compared to purely nodal or staggered methods, for all the test cases presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Impact of operating reserve rules on electricity prices with high penetrations of renewable energy

In competitive wholesale electricity markets, significant effort is devoted to designing markets that set efficient prices for maintaining supply-demand balance. One factor that can impact prices is administratively-set scarcity pricing, which sets prices to a preset level when the market is not able to meet operating reserve or energy requirements. When energy and operating reserves are co-optimized, assumptions surrounding operating reserve requirements and scarcity pricing can impact system-wide price outcomes for both operating reserves and energy. This study uses production cost modeling of an ERCOT-like system to evaluate the impact of operating reserve eligibility, scarcity pricing, and quantity rules on electricity prices, and therefore also on generator revenues. Results reveal economic and operational benefits with allowing open participation in reserve markets, as this enables greater access to the full set of capable resources at lowest cost. Furthermore, both energy and reserve prices are strongly impacted by reserve scarcity pricing events, which reveals that reserve scarcity pricing assumptions can impact price outcomes even for units not providing reserves. This study highlights the importance of operating reserve scarcity pricing rules because of the strong coupling between energy and reserve prices and because these rules serve as proxies for true price responsive demand.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Modeling Global Vegetation Gross Primary Productivity, Transpiration and Hyperspectral Canopy Radiative Transfer Simultaneously Using a Next Generation Land Surface Model—CliMA Land

Recent progress in satellite observations has provided unprecedented opportunities to monitor vegetation activity at global scale. However, a major challenge in fully utilizing remotely sensed data to constrain land surface models (LSMs) lies in inconsistencies between simulated and observed quantities. For example, gross primary productivity (GPP) and transpiration (T) that traditional LSMs simulate are not directly measurable from space, although they can be inferred from spaceborne observations using assumptions that are inconsistent with those LSMs. In comparison, canopy reflectance and fluorescence spectra that satellites can detect are not modeled by traditional LSMs. To bridge these quantities, we presented an overview of the next generation land model developed within the Climate Modeling Alliance (CliMA), and simulated global GPP, T, and hyperspectral canopy radiative transfer (RT; 400–2,500 nm for reflectance, 640–850 nm for fluorescence) at hourly time step and 1° spatial resolution using CliMA Land. CliMA Land predicts vegetation indices and outgoing radiances, including solar-induced chlorophyll fluorescence (SIF), normalized difference vegetation index (NDVI), enhanced vegetation index (EVI), and near infrared reflectance of vegetation (NIRv) for any given sun-sensor geometry. The spatial patterns of modeled GPP, T, SIF, NDVI, EVI, and NIRv correlate significantly with existing data-driven products (mean R 2 = 0.777 for 9 products). CliMA Land would be also useful in high temporal resolution simulations, for example, providing insights into when GPP, SIF, and NIRv diverge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A 13‐Year Record Indicates Differences in the Duration and Depth of Soil Carbon Accrual Among Potential Bioenergy Crops

Six years after replacing a maize/soybean cropping system, perennial grasses miscanthus (Miscanthus × giganteus) and switchgrass (Panicum virgatum), and a 28-species restored prairie increased particulate organic carbon in surface soils without increasing soil organic carbon (SOC). To resolve potential changes in the quantity and distribution of SOC, soils were resampled after seven to thirteen years to measure bulk density, carbon (C) content, and stable C isotopes to a depth of 1 m. SOC stocks increased between 1.75 and 2.5 Mg ha −1 year −1 in all perennial crops between 2008 and 2016 (nine growing seasons). Despite relatively low litter inputs and belowground biomass, the highest rate of SOC accrual was in restored prairie (2.5 Mg ha −1 year −1 ), followed by miscanthus (2.0 Mg ha −1 year −1 ) and switchgrass (1.75 Mg ha −1 year −1 ). The change in SOC in maize/soybean was not significant. After 2016, total SOC decreased in maize/soybean and miscanthus, resulting in slower overall rates of SOC accumulation over the full sampling period for miscanthus (0.8 Mg ha −1 year −1 ). The rate of SOC accumulation was greatest below 50 cm depth for restored prairie and switchgrass but in the top 10 cm for miscanthus. Stable isotope analysis showed 13 C enrichment in all depths of switchgrass soils, an indication of new organic C accumulation, but mixed results in all other crops. Planting perennial crops on land formerly in an annual maize/soybean cropping system can slow or reverse soil carbon losses, with the greatest increases in SOC from species-rich prairie.

agriculture↗
NASA NTRS2016

The Darkening of the Greenland Ice Sheet: Trends, Drivers and Projections (1981-2100)

The surface energy balance and meltwater production of the Greenland ice sheet (GrIS) are modulated by snow and ice albedo through the amount of absorbed solar radiation. Here we show, using space-borne multispectral data collected during the 3 decades from 1981 to 2012, that summertime surface albedo over the GrIS decreased at a statistically significant (99 %) rate of 0.02 decade(sup -1) between 1996 and 2012. Over the same period, albedo modelled by the Modele Atmospherique Regionale (MAR) also shows a decrease, though at a lower rate (approximately -0.01 decade(sup -1)) than that obtained from space-borne data. We suggest that the discrepancy between modelled and measured albedo trends can be explained by the absence in the model of processes associated with the presence of light-absorbing impurities. The negative trend in observed albedo is confined to the regions of the GrIS that undergo melting in summer, with the dry snow zone showing no trend. The period 1981-1996 also showed no statistically significant trend over the whole GrIS. Analysis of MAR outputs indicates that the observed albedo decrease is attributable to the combined effects of increased near-surface air temperatures, which enhanced melt and promoted growth in snow grain size and the expansion of bare ice areas, and to trends in light-absorbing impurities (LAI) on the snow and ice surfaces. Neither aerosol models nor in situ and remote sensing observations indicate increasing trends in LAI in the atmosphere over Greenland. Similarly, an analysis of the number of fires and BC emissions from fires points to the absence of trends for such quantities. This suggests that the apparent increase of LAI in snow and ice might be related to the exposure of a "dark band" of dirty ice and to increased consolidation of LAI at the surface with melt, not to increased aerosol deposition. Albedo projections through to the end of the century under different warming scenarios consistently point to continued darkening, with albedo anomalies averaged over the whole ice sheet lower by 0.08 in 2100 than in 2000, driven solely by a warming climate. Future darkening is likely underestimated because of known underestimates in modelled melting (as seen in hindcasts) and because the model albedo scheme does not currently include the effects of LAI, which have a positive feedback on albedo decline through increased melting, grain growth, and darkening.

surface energy↗