Search NASA⌕ Search

SEARCH · Search NASA

Results for “higher accuracy”

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 991 records · Page 55

Pre- and postprocessing techniques for determining goodness of computational meshes

Research in error estimation, mesh conditioning, and solution enhancement for finite element, finite difference, and finite volume methods has been incorporated into AUDITOR, a modern, user-friendly code, which operates on 2D and 3D unstructured neutral files to improve the accuracy and reliability of computational results. Residual error estimation capabilities provide local and global estimates of solution error in the energy norm. Higher order results for derived quantities may be extracted from initial solutions. Within the X-MOTIF graphical user interface, extensive visualization capabilities support critical evaluation of results in linear elasticity, steady state heat transfer, and both compressible and incompressible fluid dynamics.

Oden, J. Tinsley↗

A demonstration of TDRS orbit determination using differential tracking observables from GPS ground receivers

A Global Positioning System (GPS) flight receiver provides a means of precisely determining orbits for satellites in low to moderate altitude orbits. Above 5000 km altitude, however, relatively few GPS satellites are visible. We present a new approach to orbit determination for satellites at higher altitudes. Modification of GPS ground receivers enables a beacon from the orbiter to be tracked simultaneously with GPS data. The orbit accuracy expected from this GPS-like tracking (GLT) technique in principle could be comparable to accuracies of GPS orbits. Present-day GPS orbit quality from a daily semi-automated analysis system at the Jet Propulsion Laboratory is at the 30 - 50 cm level: expected accuracies for orbiters with GLT are calculated to be in the few-meter range for altitudes up to 100000 km. For geosynchronous satellites, however, there are unique challenges due to geometrical limitations and to the lack of strong dynamical signature in tracking data. We examine two approaches for tracking Tracking and Data Relay Satellites (TDRS) geostationary orbiters. One uses GLT with a global network; the other relies on a small `connected element' ground network with a distributed clock for short baseline differential carrier phase (SB Delta Phi). In priciple, both could meet 50-m TDRS operational requirements. However there are practical difficulties with either pure GLT or SB Delta Phi tracking schemes for present-day TDRS satellites. We describe an experiment planned for late 1993 which will combine aspects of both GLT and SB Delta Phi to demonstrate a new approach for tracking TDRS which offers a number of operationally convenient and attractive features. The TDRS demo will in effect ba a proof of concept experiment for both the GLT and SB Delta Phi approaches to tracking spacecraft.

Lichten, Stephen M.↗

Orion Entry Performance-Based Center-of-Gravity Box

The Orion capsule is designed both for Low Earth Orbit missions to the ISS and for missions to the moon. For ISS class missions, the capsule will use an Apollo-style direct entry. For lunar return missions, depending on the timing of the mission, the capsule could perform a direct entry or a skip entry of up to 4800 n.mi. in order to land in the coastal waters of California. The physics of atmospheric re-entry determine the capability of the Orion vehicle. For a given vehicle mass and shape, physics tells us that the driving parameters for an entry vehicle are the hypersonic lift-to-drag ratio (L/D) and the flight path angle at entry interface (gamma(sub EI)). The design of the Orion atmospheric re-entry must meet constraints during both nominal and dispersed flight conditions on landing accuracy, heating rate, total heat load, sensed acceleration, and proper disposal of the Service Module. These constraints define an entry corridor in the space of L/D-gamma(sub EI); if the vehicle falls within this corridor, then all constraints are met. The gamma(sub EI) dimension of the corridor can be further constrained by the gloads experienced during emergency entries. Thus, the entry performance for the Orion vehicle can be described completely by the L/D. Bounds on the hypersonic L/D necessary to achieve all the mission requirements can be defined for the given entry corridor. Landing accuracy performance drives the lower limit on L/D. In order to achieve the desired landing accuracy, a minimum L/D must be ensured. The design of the Thermal Protection System (TPS) drives the upper limit on L/D. A higher L/D can drive mass into the design of the TPS. Conversely, once the TPS is designed, the L/D must be ensured to stay below a certain limit in order for the TPS to stay within its design envelop. The L/D must stay within its upper and lower bounds during dispersed flight conditions. L/D is a function of both the aerodynamics and the center-of-gravity (CG) of the vehicle. The aerodynamics of the vehicle are determined by Computational Fluid Mechanics (CFD) and wind tunnel tests. However, the aerodynamics are not known precisely. Instead, an aerodynamic database has been developed where the aerodynamic coefficients are known to fall within a probabilistic band defined by upper and lower bounds. It is expected that the probabilistic band will shrink after the first missions are flown and real-world data is collected. Until that time, the Orion must be designed to the current aerodynamic database. Thus, for a given aerodynamic database with given uncertainties, the allowable range in L/D can be mapped to an allowable box for the CG location. The CG box is used to set requirements on the dispersions allowed for vehicle packaging and cargo storage. As the aerodynamic uncertainties decrease, the size of the CG box can increase. This paper discusses the technique used to map the minimum and maximum L/D bounds set by the entry performance requirements to the allowable dispersions in CG while accounting for aerodynamic uncertainties. The L/D is defined as the ratio of the lift force to the drag force. It is equivalent to the ratio of lift coefficient (C(sub L)) over drag coefficient (C(sub D)). C(sub L) and C(sub D) are functions of Mach number (M) and angle of attack (alpha). A Mach number of 25 is used as a measuring point of the hypersonic L/D. Variations in C(sub L), C(sub D) and alpha cause variations in L/D. Equation (1) shows the three contributions to the variation in L/D.

Rea, Jeremy R.↗

The Need for Technology Maturity of Any Advanced Capability to Achieve Better Life Cycle Cost (LCC)

Programs such as space transportation systems are developed and deployed only rarely, and they have long development schedules and large development and life cycle costs (LCC). They have not historically had their LCC predicted well and have only had an effort to control the DDT&E phase of the programs. One of the factors driving the predictability, and thus control, of the LCC of a program is the maturity of the technologies incorporated in the program. If the technologies incorporated are less mature (as measured by their Technology Readiness Level - TRL), then the LCC not only increases but the degree of increase is difficult to predict. Consequently, new programs avoid incorporating technologies unless they are quite mature, generally TRL greater than or equal to 7 (system prototype demonstrated in a space environment) to allow better predictability of the DDT&E phase costs unless there is no alternative. On the other hand, technology development programs rarely develop technologies beyond TRL 6 (system/subsystem model or prototype demonstrated in a relevant environment). Currently the lack of development funds beyond TRL 6 and the major funding required for full scale development leave little or no funding available to prototype TRL 6 concepts so that hardware would be in the ready mode for safe, reliable and cost effective incorporation. The net effect is that each new program either incorporates little new technology or has longer development schedules and costs, and higher LCC, than planned. This paper presents methods to ensure that advanced technologies are incorporated into future programs while providing a greater accuracy of predicting their LCC. One method is having a dedicated organization to develop X-series vehicles or separate prototypes carried on other vehicles. The question of whether such an organization should be independent of NASA and/or have an independent funding source is discussed. Other methods are also discussed. How to make the choice of which technologies to pursue to the prototype level is also discussed since, to achieve better LCC, first the selection of the appropriate technologies.

Robinson, John W.↗

Goddard Laser for Absolute Measurement of Radiance

High accuracy achieved with temperature stabilized unfiltered trap detectors illuminated by monochromatic light. Tunable lasers and optical parametric oscillators provide orders of magnitude higher spectral radiance than blackbody or other broadband sources; calibrate at high signal levels.

McAndrew, Brendan↗

On-Orbit Tracking of Sub-Sample Gain Differences in SNPP and NOAA-20 VIIRS Imagery Bands

The VIIRS instruments on board the SNPP and NOAA-20 (N20) satellites have 14 reflective solar bands covering a spectral range from 412 nm to 2250 nm. Three of these are imaging bands (I bands) with a nadir spatial resolution of 375 m and 11 are moderate resolution bands (M bands) with a resolution of 750 m. The higher resolution in the I bands is achieved by a combination of more detectors, with the I bands having twice as many detectors of half the size for every M band detector, and a higher data rate, with the I bands having two sub-samples for every sample of M band data. To ensure calibration accuracy, any systematic difference in the response of the two sub-samples needs to be monitored and corrected in the calibrated products. In this paper, we use the solar diffuser calibrations to monitor the gain differences between the two sub-samples of the I bands both as a function of time and signal level. We find gain differences of about 0.1% for I1, 0.3% for I2, and <0.1% for I3 that are mostly constant over the range of signal values available in the SD calibration. These values are mostly consistent throughout the mission for both instruments. The results are remarkably similar for the two VIIRS instruments, including a slightly out-of-family behavior seen in a few detectors. We discuss possible causes for the difference and the impact on the aggregated Earth view images.

On-orbit↗

Evaluation of Voronoi Meshes for Large Eddy Simulations of High Lift Aerodynamics

Numerical sensitivity to 3 different Voronoi seeding methods is investigated for Large Eddy Simulations (LES). A second order accurate, non-dissipative finite volume discretization is used to systematically investigate the effects of different polyhedral Voronoi mesh types using a sequence of three canonical problems with increasing complexity. First, inviscid isentropic vortex propagation is studied to demonstrate the substantial reduction of errors for rhombic dodecahedron and truncated octahedral cell types over Cartesian hexagonal cells of identical spacing. Furthermore, the reduction of errors induced at cell-size transitions (grid-coarsening interfaces) due to Lloyd smoothing iterations is quantified. It is shown that by utilizing an appropriate viscous flux discretization, a constant coefficient subgrid scale model is sufficient for non-linear stability at the 2:1 cell-size transitions on polyhedral grids, although some further error reduction does occur when smoothing is utilized. Next, forced homogeneous isotropic turbulence at an asymptotically large Reynolds number is studied to demonstrate the non-dissipative character of the inviscid flux discretization, and the non-linear robustness and accuracy offered by the viscous flux discretization using a subgrid scale model for all three Voronoi grid types. Finally,Wall-Modeled Large Eddy Simulations (WMLES) are performed to study the high-lift aerodynamics on the McDonnell Douglas 30P30N multi-element airfoil at two distinct grid levels and for two distinct Voronoi cell types. The formulation is shown to predict the aerodynamic loading with high accuracy at all angles of attack when sufficient resolution is reached, and the hexagonal prism grid topology, while computationally more expensive, has higher effective resolution compared to the Cartesian grid topology with the same spacing.

TTT↗

Evaluation of Voronoi Meshes for Large Eddy Simulations of High Lift Aerodynamics

Numerical sensitivity to 3 different Voronoi seeding methods is investigated for Large Eddy Simulations (LES). A second order accurate, non-dissipative finite volume discretization is used to systematically investigate the effects of different polyhedral Voronoi mesh types using a sequence of three canonical problems with increasing complexity. First, inviscid isentropic vortex propagation is studied to demonstrate the substantial reduction of errors for rhombic dodecahedron and truncated octahedral cell types over Cartesian hexagonal cells of identical spacing. Furthermore, the reduction of errors induced at cell-size transitions (grid-coarsening interfaces) due to Lloyd smoothing iterations is quantified. It is shown that by utilizing an appropriate viscous flux discretization, a constant coefficient subgrid scale model is sufficient for non-linear stability at the 2:1 cell-size transitions on polyhedral grids, although some further error reduction does occur when smoothing is utilized. Next, forced homogeneous isotropic turbulence at an asymptotically large Reynolds number is studied to demonstrate the non-dissipative character of the inviscid flux discretization, and the non-linear robustness and accuracy offered by the viscous flux discretization using a subgrid scale model for all three Voronoi grid types. Finally, Wall-Modeled Large Eddy Simulations (WMLES) are performed to study the high-lift aerodynamics on the McDonnell Douglas 30P30N multi-element airfoil at two distinct grid levels and for two distinct Voronoi cell types. The formulation is shown to predict the aerodynamic loading with high accuracy at all angles of attack when sufficient resolution is reached, and the hexagonal prism grid topology, while computationally more expensive, has higher effective resolution compared to the Cartesian grid topology with the same spacing.

TTT↗

Forecasted Dynamic Line Ratings for Regional Transmission Lines

This report was prepared for the Wind Energy Technology Office for the FY 2022, quarter 3 deliverable. This details the use of forecasted dynamic line rating as an improvement over static rating for three transmission lines in three regions of interest. A selection of wind plants across Idaho, the Columbia River Gorge and offshore wind in the NYSERDA territory were modeled for their concurrent cooling effects for dynamic line rating. Regional transmission lines for the wind plant interconnections were determined and a CFD domain was created to determine local wind flows over the lines. A dynamic line rating for each of the lines was calculated for both weather observation data and HRRR model forecast data. The comparison of peaks in the wind power production were compared to the DLR peaks. It was determined that for shorter regional transmission lines, the effect of concurrent cooling that commonly occur with gen-tie lines were still valid. Concurrent cooling effects drop as the length of the regional transmission lines increase in distance from the wind plants. In addition, the accuracy of the forecasted ampacity to the weather observation data for ampacity was assessed. In general, the accuracy was best for the INL region with the NOAA-maintained sites. The Columbia River Gorge and Long Island sites showed similar levels of accuracy for the forecasted DLR. The accuracy levels were around 11% for the INL region, 15% for Columbia River Gorge and 7% for New York Long Island. The higher error for the Columbia Gorge region is likely due to the sparse availability of weather stations available, whereas both the INL site and Long Island are well instrumented.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Comparison of surface wind stress measurements - Airborne radar scatterometer versus sonic anemometer

Sea surface wind stress measurements recorded by a sonic anemometer are correlated with airborne scatterometer measurements of ocean roughness (cross section of radar backscatter) to establish the accuracy of remotely sensed data and assist in the definition of geophysical algorithms for the scatterometer sensor aboard Seasat A. Results of this investigation are as follows: Comparison of scatterometer and sonic anemometer wind stress measurements are good for the majority of cases; however, a tendency exists for scatterometer wind stress to be somewhat high for higher wind conditions experienced in this experiment (6-9 m/s). The scatterometer wind speed algorithm tends to overcompute the higher wind speeds by approximately 0.5 m/s. This is a direct result of the scatterometer overestimate of wind stress from which wind speeds are derived. Algorithmic derivations of wind speed and direction are, in most comparisons, within accuracies defined by Seasat A scatterometer sensor specifications.

Brucks, J. T.↗

Higher order accurate difference solutions of fluid mechanics problems by a compact differencing technique

A general fourth order differencing scheme is developed and applied to three viscous test problems to verify the accuracy and applicability of the technique. The procedure is atypical since only three nodes are necessary to obtain the desired fourth order accuracy. This is accomplished by a differencing technique which considers the function and all necessary derivatives as unknowns. The relations for these derivatives yield simple tridiagonal equations which can be easily solved. Comparisons of the fourth order results with those computed using second order methods are presented for the test problems and clearly indicate that the accuracy achieved by these fourth order computations is always significantly better than current second order procedures.

Hirsh, R. S.↗

Sensor Calibration Impacts on Dust Detection Based On MODIS And VIIRS Thermal Emissive Bands

Dust detection using remotely sensed measurements has been one of the challenging problems encountered by atmospheric scientists. MODerate Resolution Imaging Spectroradiometer (MODIS) on the Terra (T)and Aqua (A) platforms have been a versatile sensor for well over 21 and 18 years respectively and have been extremely useful in the retrieval of aerosol information over the entire globe. The MODIS radiances from the Level1B in general are expected to be within 5% accuracy in the reflective wavelengths and within 1% in the thermal emissive wavelengths. In this paper, we evaluate the sensitivity of previously developed dust detection technique based on thermal emissive wavelengths, which correspond to MODIS bands 20, 29, 31, and 32 respectively. The Thermal Emissive Dust Index (TEDI) performed very comparably to the traditional Aerosol Optical Thickness (AOT) retrievals by MODIS reflective channels. Since the MODIS Thermal Emissive Bands (TEB) are well calibrated on-orbit using a BlackBody (BB) source, the calibration of these long wave infrared bands is quite robust. As A-MODIS continues to perform well beyond its designed lifetime of 6 years, the instrument has undergone various levels of degradation during its mission time. As a consequence, it is imperative to check the impacts of calibration on the higher-level retrievals. In this paper, we rigorously analyze the sensitivity of TEDI due to the impact of calibration by the afore-mentioned TEB. The perturbation of the dominant (linear) calibration term demonstrated the following: first, there was a correlation in the sensitivity of the TEDI due to the uncertainty in the linear calibration term. Based on a perturbation in the linear calibration term for all aforementioned bands over a range of ± 5 % yielded the TEDI sensitivity to vary from approximately -3.2 % to about -3.6 %. When considering the uncertainty in each individual band significant changes were observed. The least change was observed for the perturbation in the calibration of band 20 with the TEDI sensitivity and the largest sensitivity in TEDI was observed in the perturbation of band 31 calibration. Thus, in the case of TEDI, noticeable sensitivity due to calibration uncertainty was observed in bands 29, 31, and 32, reiterating the importance of the TEB calibration in these bands. Also, the dust detection scheme based on A-MODIS was successfully transferred to the follow-on sensors such as Suomi (SNPP) and NOAA 20 (N20) VIIRS. The results presented in this paper would be extremely helpful in understanding impacts of calibration on the higher-level products for both current and future missions based on the MODIS heritage. Finally, the work also identifies the importance of radiometric fidelity in maintaining the accuracy of the dust detection. Results presented will show drastic improvement of the Saharan dust detection after the reduction of the electronic crosstalk in the 8.5 μm channel of T-MODIS.

Sriharsha Madhavan↗

Developing a Consistent Travel-Time Framework for Comparing Three-Dimensional Velocity Models for Seismic Location Accuracy

Abstract Location algorithms have historically relied on simple, one-dimensional (1D) velocity models for fast seismic event locations. 1D models are generally used as travel-time lookup tables, one for each seismic phase, with travel-times pre-calculated for event distance and depth. These travel-time lookup tables are extremely fast to use and this fast computational speed makes them the preferred type of velocity model for operational needs. Higher-dimensional (i.e., three-dimensional—3D) seismic velocity models are becoming readily available and provide more accurate event locations over 1D models. The computational requirements of these 3D models tend to make their operational use prohibitive. Additionally, comparing location accuracy for 3D seismic velocity models tends to be problematic, as each model is determined using different ray-tracing algorithms. Attempting to use a different algorithm than the one used to develop a model usually results in poor travel-time prediction. We demonstrate and test a framework to create first-P and first-S 3D travel-time correction surfaces using an open-source framework ( PCalc + GeoTess , https://www.sandia.gov/salsa3d/software/geotess ) that easily stores 3D travel-time and uncertainty data. This framework produces fast travel-time and uncertainty predictions and overcomes the ray-tracing algorithm hurdle because the lookup tables can be generated using the exact ray-tracing algorithm that is preferred for a model.

3D velocity models↗

Evaluation of GEOS Precipitation Flagging for SMAP Soil Moisture Retrieval Accuracy

The precipitation flag in the Soil Moisture Active Passive (SMAP) Level 2 passive soil moisture (L2SMP) retrieval product indicates the presence or absence of heavy precipitation at the time of the SMAP overpass. The flag is based on precipitation estimates from the Goddard Earth Observing System (GEOS) Forward Processing numerical weather prediction system. An error in flagging during an active or recent precipitation event can produce either 1) an overestimation of soil moisture due to short-term surface wetting of vegetation and/or surface ponding (if soil moisture retrieval was attempted in the presence of rain) or 2) an unnecessary nonretrieval of soil moisture and loss of data (if retrieval is flagged due to an erroneous indication of rain). Satellite precipitation estimates from the Integrated Multisatellite Retrievals for GPM (IMERG), version 06, Early Run (latency of ~4 h) precipitationCal product are used here to evaluate the GEOS-based precipitation flag in the L2SMP product for both the 1800 local time (LT) ascending and 0600 LT descending SMAP overpasses over the first five years of the mission (2015–20). Consisting of blended precipitation measurements from the Global Precipitation Mission (GPM) satellite constellation, IMERG is treated as the “truth” when comparing to the GEOS model forecasts of precipitation used by SMAP. Key results include (i) IMERG measurements generally show higher spatial variability than the GEOS forecast precipitation, (ii) the IMERG product has a higher frequency of light precipitation amounts, and (iii) the effect of incorporating IMERG rainfall measurements in lieu of GEOS precipitation forecasts are minimal on the L2SMP retrieval accuracy (determined vs in situ soil moisture measurements at core validation sites). Our results indicate that L2SMP retrievals continue to meet the mission’s accuracy requirement [standard deviation of the unbiased RMSE (ubRMSE) less than 0.04 cu. m/cu. m].

IMERG-precipitationCal↗

Coordinate parameterisation and spectral method optimisation for Beltrami field solver in stellarator geometry

The numerical solution of the stepped pressure equilibrium (Hudson et al 2012 Phys. Plasmas 19 112502) requires a fast and robust solver to obtain the Beltrami field in three-dimensional geometry such as stellarators. The spectral method implemented in the stepped pressure equilibrium code (SPEC) is efficient when the domain is a hollow torus, but ill-conditioning of the discretised linear equations occurs in the (solid) toroid due to the artificially singular coordinate parameterisation near the axis. Here, we propose an improved choice for the reference axis to prevent coordinates surfaces from overlapping. Then, we examine the parity and asymptotics of the magnetic vector potential near the axis and suggest the use of recombined and rescaled Zernike radial basis functions. The maximum relative error in the magnetic field of the Wendelstein 7-X geometry is shown to reach 10 –9 at high resolution in a series of convergence tests and benchmarks against the boundary integral equation solver for Taylor states. The new method is also reported to significantly improve the accuracy of multi-volume SPEC calculations. A comparison between free-boundary SPEC and the analytical Dommaschk potential is presented with higher-than-usual Fourier resolution. It is illustrated that we are able to resolve low amplitude current sheets when an interface is placed where there is no flux surface in the analytic solution. This was previously concealed because of insufficient numerical resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Improvements to the Blade Element Momentum Formulation of OpenFAST for Skewed Inflows

In this work, we modify the blade element momentum algorithm of OpenFAST to improve its predictions under large skewed inflow conditions. We use the well-known Glauert's skew correction and introduce continuous extension of the model for high-thrust conditions. We present the rationale behind Glauert's empirical model and discuss the different conventions possible for the axial induction factor. We verify the model against the higher-fidelity lifting-line vortex method and blade-resolved computational fluid dynamics, and we observe that the new implementation enhances the accuracy and reliability of OpenFAST's aerodynamic modeling capabilities in conditions involving large skew angles. For the parametric studies run using the different codes, we find that the power changes with the skew angle as cos 1.7 (θ skew ) and the thrust as cos 0.65 (θ skew ). An analysis of the azimuthal variation of the induced velocities in the rotor plane reveals that current redistribution models used in blade element momentum codes may need to be refined.

17 WIND ENERGY↗

Machine-learning prediction for quasiparton distribution function matrix elements

There have been rapid developments in the direct calculation in lattice QCD (LQCD) of the Bjorken-x dependence of hadron structure through large-momentum effective theory (LaMET). LaMET overcomes the previous limitation of LQCD to moments (that is, integrals over Bjorken x) of hadron structure, allowing LQCD to directly provide the kinematic regions where the experimental values are least known. LaMET requires large-momentum hadron states to minimize its systematics and allow us to reach small- x reliably. This means that very fine lattice spacing to minimize lattice artifacts at order $(P_za)^n$ will become crucial for next-generation LaMET-like structure calculations. Furthermore, such calculations require operators with long Wilson-link displacements, especially in finer lattice units, increasing the communication costs relative to that of the propagator inversion. In this work, we explore whether machine-learning algorithms can make predictions of correlators to reduce the computational cost of these LQCD calculations. We consider two algorithms, gradient-boosting decision tree and linear models, applied to LaMET data, the matrix elements needed to determine the kaon and $η_s$ unpolarized parton distribution functions (PDFs), meson distribution amplitude (DA), and the nucleon gluon PDF. We find that both algorithms can reliably predict the target observables with different prediction accuracy and systematic errors. The predictions from smaller displacement z to larger ones work better than those for momentum p due to the higher correlation among the data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Upscaling Methane Flux From Plot Level to Eddy Covariance Tower Domains in Five Alaskan Tundra Ecosystems

Spatial heterogeneity in methane (CH 4 ) flux requires a reliable upscaling approach to reach accurate regional CH 4 budgets in the Arctic tundra. In this study, we combined the CLM-Microbe model with three footprint algorithms to scale up CH 4 flux from a plot level to eddy covariance (EC) tower domains (200 m × 200 m) in the Alaska North Slope, for three sites in Utqiaġvik (US-Beo, US-Bes, and US-Brw), one in Atqasuk (US-Atq) and one in Ivotuk (US-Ivo), for a period of 2013–2015. Three footprint algorithms were the homogenous footprint (HF) that assumes even contribution of all grid cells, the gradient footprint (GF) that assumes gradually declining contribution from center grid cells to edges, and the dynamic footprint (DF) that considers the impacts of wind and heterogeneity of land surface. Simulated annual CH 4 flux was highly consistent with the EC measurements at US-Beo and US-Bes. In contrast, flux was overestimated at US-Brw, US-Atq, and US-Ivo due to the higher simulated CH 4 flux in early growing seasons. The simulated monthly CH 4 flux was consistent with EC measurements but with different accuracies among footprint algorithms. At US-Bes in September 2013, RMSE and NNSE were 0.002 μmol m −2 s −1 and 0.782 using the DF algorithm, but 0.007 μmol m −2 s −1 and 0.758 using HF and 0.007 μmol m −2 s −1 and 0.765 using GF, respectively. DF algorithm performed better than the HF and GF algorithms in capturing the temporal variation in daily CH 4 flux each month, while the model accuracy was similar among the three algorithms due to flat landscapes. Temporal variations in CH 4 flux during 2013–2015 were predominately explained by air temperature (67–74%), followed by precipitation (22–36%). Spatial heterogeneities in vegetation fraction and elevation dominated the spatial variations in CH 4 flux for all five tower domains despite relatively weak differences in simulated CH 4 flux among three footprint algorithms. The CLM-Microbe model can simulate CH 4 flux at both plot and landscape scales at a high temporal resolution, which should be applied to other landscapes. Integrating land surface models with an appropriate algorithm provides a powerful tool for upscaling CH 4 flux in terrestrial ecosystems.

54 ENVIRONMENTAL SCIENCES↗