Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “SATELLITES”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Exomoons in Systems with a Strong Perturber: Applications to α Cen AB

The presence of a stellar companion can place constraints on occurrence and orbital evolution of satellites orbiting exoplanets, i.e., exomoons. In this work we revise earlier orbital stability limits for retrograde orbits in the case of a three-body system consisting of a star, planet, and satellite. The revised limit reads a{sub sat}{sup crit}≈0.668(1-1.236e{sub p}) for e {sub p} ≤ 0.8 in units of the Hill Radius and represents the lower critical orbit as a function of the planetary eccentricity e {sub p}. A similar formula is determined for exomoons hosted by planets in binary star systems, where e {sub p} is replaced with the components of free and forced eccentricity from secular orbit evolution theory. By exploring the dynamics of putative exomoons in α Centauri AB we find that the outer stability limit can be much less than half the Hill Radius due to oscillations in the planetary orbital eccentricity caused by the gravitational interaction with the binary star. We show, furthermore, how the resulting truncation of the outer stability limit can affect the outward tidal migration and potential observability of exomoons through transit-timing variations (TTVs). Typical TTV (rms) amplitudes induced by exomoons in binary systems are ≲10 minutes and appear more likely for planets orbiting the less massive stellar component.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

VISSTPX2DRECTG16V4MINNIS (c1)

Satellite-based retrievals of cloud and radiation properties are available in this value-added product provided by Bill Smith’s group at NASA/Langley using the VISST (Visible Infrared Solar-Infrared Split Window Technique) algorithm. This product contains VISST-derived pixel-level products from the GOES-16 satellite, rectilinear projection 2-D array.

broadband LW flux↗
DOE OSTI2024

Quantifying Error in Photovoltaic Installation Metadata: Preprint

In this research, we quantify the level of metadata error for a fleet of 2860 photovoltaic (PV) systems, using metadata values provided by fleet owners. Using satellite imagery and time series analysis techniques available in open-source Python packages Panel-Segmentation and PVAnalytics, respectively, we evaluate the accuracy of PV system metadata such as location, azimuth, tilt, and mounting configuration (fixed tilt vs. tracking). We find that approximately 75% of provided latitude-longitude coordinates are within 190 meters of the actual solar installation. We were unable to link 7.8% of latitude-longitude coordinates to any solar installation via satellite imagery analysis. We evaluate the level of error in owner-provided mounting configuration (fixed tilt vs. single-axis tracking), finding only 8 systems with an incorrect mounting configuration. When evaluating azimuth and tilt parameters, we find that approximately 64% of the data is correct, with data for 860 systems (approximately 30%) not provided by system owners. To illustrate the importance of having correct solar metadata, we evaluate how incorrect metadata affects solar performance estimates by modeling system AC energy output at ground-truth vs. incorrect latitude-longitude coordinates, mounting configurations, and azimuth-tilt configurations. Energy output estimates can vary significantly if incorrect metadata parameters are used, with incorrect mounting configuration leading to the largest discrepancy with over 20% variation in expected energy output.

azimuth↗
DOE OSTI2026

A Global Methane Observation System to Reduce Uncertainty for Anthropogenic and Natural Sources and Sinks for Detecting and Attributing Climate Feedbacks

Atmospheric methane (CH4) concentrations are accelerating global warming as net emissions increase. Observing systems that quantify sources remain too sparse and fragmented to detect trends—especially in remote regions where climate‐driven natural emissions may be rising. We provide a framework for quantifying uncertainty reductions through the implementation of a global ecosystem‐methane observing system designed to: (i) substantially lower uncertainty in sectoral and regional emissions, (ii) separate co‐occurring anthropogenic and natural fluxes, and (iii) trend detection at regional scales to verify mitigation progress and provide early warning of natural feedbacks. Using bottom‐up inventories and process‐model ensembles for 2014–2023, we show that anthropogenic emissions remain uncertain by ∼32% globally, while natural sources—tropical and boreal‐arctic wetlands, fires, and inland waters—carry far larger uncertainties (+ 70%) and trend uncertainties reaching ∼200%. Additional observations must match spatial emission structure to increase observability of emissions: high‐resolution satellite constellations for point sources combined with expanded flux networks and wetland mapping for diffuse sources, and denser ground‐based atmospheric column measurements to restore observability in under‐sampled tropics and high latitudes. Notional analyses indicate that targeted additions of flux towers and ∼20 in situ atmospheric column concentration instruments per key tropical region could reduce continental‐scale uncertainties at modest cost. Conceptual illustration of a Global Ecosystem Methane Observing System (GEM‐OS) integrating satellites, aircraft, atmospheric networks, and ecosystem measurements to quantify methane emissions from anthropogenic and natural sources. The multi‐scale observing framework improves source attribution, reduces uncertainty in regional methane budgets, and enables early detection of climate‐driven feedbacks from wetlands, fires, permafrost, agriculture, and fossil‐fuel emissions.

Ciais, P↗
DOE OSTI2024

Using phenology to unravel differential soil water use and productivity in a semiarid savanna

Savannas are water-limited ecosystems characterized by two dominant plant types: trees and an understory primarily made up grass. Different phenology and root structures of these plant types complicate how savanna primary productivity responds to changes in water availability. We tested the hypothesis that productivity in savannas is controlled by the temporal and vertical distribution of soil water content (SWC) and differences in growing season length of understory and tree plant functional types. To quantify the relationship between tree, understory, and savanna-wide phenology and productivity, we used PhenoCam and satellite observations surrounding an eddy covariance tower at a semiarid savanna site in Arizona, USA. We distinguished between SWC across two different depth intervals (shallow, <0–30 cm and deep, >30–100 cm). We found that tree greenness increased with SWC at both depths, while understory greenness was only sensitive to the shallower SWC measurements. Onset of ecosystem dormancy, estimated from satellite observations close to the eddy covariance tower, explained more variability in annual gross primary productivity (GPP) than in other phenometrics. Higher SWC led to an extended growing season, caused by delayed dormancy in trees, but the understory showed no evidence of delayed dormancy in wetter periods. We infer that the timing of ecosystem scale dormancy, driven by trees, is important in understanding changes in a savanna's GPP. These findings highlight the important effects of rainfall during the winter. These findings suggest that savanna GPP is conditional on different responses to moisture availability in each of the dominant vegetation components.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A Nonstationary and Non-Gaussian Moving Average Model for Solar Irradiance

Historically, power has flowed from large power plants to customers. Increasing penetration of distributed energy resources such as solar power from rooftop photovoltaic has made the distribution network a two-way-street with power being generated at the customer level. The incorporation of renewables introduces additional uncertainty and variability into the power grid. Distribution network operation studies are being adapted to include renewables; however, such studies require high quality solar irradiance data that adequately reflect realistic meteorological variability. Data from satellite-based products are spatially complete, but temporally coarse, whereas solar irradiances exhibit high frequency variation at very fine timescales. We propose a new stochastic method for temporally downscaling global horizontal irradiance (GHI) to 1 min resolution, but we do not consider the spatial aspect due to limited availability of the in situ irradiance measurements. Solar irradiance's first and second-order structures vary diurnally and seasonally, and our model adapts to such nonstationarity. Empirical irradiance data exhibits highly non-Gaussian behavior; we develop a nonstationary and non-Gaussian moving average model that is shown to capture realistic solar variability at multiple timescales. We also propose a new estimation scheme based on Cholesky factors of empirical autocovariance matrices, bypassing difficult and inaccessible likelihood-based approaches. The model is demonstrated for a case study of three locations that are located in diverse climates through the United States. The model is compared against competitors from the literature and is shown to provide better uncertainty and variability quantification on testing data.

Cholesky factor↗
DOE OSTI2023

Insights into the relation between vertical cloud structure and dynamics of three polar lows: Observations from COMBLE

The vertical structure of three polar lows is described using profiling radar, lidar, and passive remote sensors deployed at a coastal site in northern Norway. The data were collected as part of the Cold-air Outbreaks in the Marine Boundary Layer Experiment (COMBLE). These data are examined in the context of satellite imagery, operational weather radar reflectivity imagery, and output from the AROME-Arctic model. A comparison to satellite images shows that AROME-Arctic captured the polar lows well. All three polar lows form in marine cold-air outbreaks and are surrounded by open-cellular convection typical of such outbreaks. Two of the lows form in a forward shear environment, and one under reverse shear. Initially, the three polar lows are comma-shaped; two of them transition to be more spirali-form at their mature stage and have mostly cloud-free warm cores. The warm cores are the result of a warm-seclusion process. All three lows have stratiform precipitation bands, marked by little cloud liquid water, and rather high surface precipitation rate. The vertical drafts and turbulence in these stratiform clouds are generally weak. In conclusion, all three lows also have convective clouds, which have stronger vertical drafts, stronger turbulence, and pockets of high liquid water content.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Opportunities for Earth Observation to Inform Risk Management for Ocean Tipping Points

Abstract As climate change continues, the likelihood of passing critical thresholds or tipping points increases. Hence, there is a need to advance the science for detecting such thresholds. In this paper, we assess the needs and opportunities for Earth Observation (EO, here understood to refer to satellite observations) to inform society in responding to the risks associated with ten potential large-scale ocean tipping elements: Atlantic Meridional Overturning Circulation; Atlantic Subpolar Gyre; Beaufort Gyre; Arctic halocline; Kuroshio Large Meander; deoxygenation; phytoplankton; zooplankton; higher level ecosystems (including fisheries); and marine biodiversity. We review current scientific understanding and identify specific EO and related modelling needs for each of these tipping elements. We draw out some generic points that apply across several of the elements. These common points include the importance of maintaining long-term, consistent time series; the need to combine EO data consistently with in situ data types (including subsurface), for example through data assimilation; and the need to reduce or work with current mismatches in resolution (in both directions) between climate models and EO datasets. Our analysis shows that developing EO, modelling and prediction systems together, with understanding of the strengths and limitations of each, provides many promising paths towards monitoring and early warning systems for tipping, and towards the development of the next generation of climate models.

Wood, Richard A. (ORCID:0000000239609513)↗
DOE OSTI2025

SUDA: A SUrface Dust Analyser for Compositional Mapping of the Galilean Moon Europa

Abstract The Surface Dust Analyser (SUDA) is a mass spectrometer onboard the Europa Clipper mission for investigating the surface composition of the Galilean moon Europa. Atmosphereless planetary moons such as the Galilean satellites are wrapped into a ballistic dust exosphere populated by tiny samples from the moon’s surface produced by impacts of fast micrometeoroids. SUDA will measure the composition of such surface ejecta during close flybys of Europa to obtain key chemical signatures for revealing the satellite’s composition such as organic molecules and salts, history, and geological evolution. Because of their ballistic orbits, detected ejecta can be traced back to the surface with a spatial resolution roughly equal to the instantaneous altitude of the spacecraft. SUDA is a Time-Of-Flight (TOF), reflectron-type impact mass spectrometer, optimized for a high mass resolution which only weakly depends on the impact location. The instrument will measure the mass, speed, charge, elemental, molecular, and isotopic composition of impacting grains. The instrument’s small size of$268 ~\mathrm {mm} \times 250 ~\mathrm {mm} \times 171$ 268 mm × 250 mm × 171 $~\mathrm {mm}$ mm , radiation-hard design, and rather large sensitive area of 220 cm 2 matches well the challenging demands of the Clipper mission.

Astronomy & Astrophysics↗
DOE OSTI2023

DownScaleBench for developing and applying a deep learning based urban climate downscaling- first results for high-resolution urban precipitation climatology over Austin, Texas

Abstract Cities need climate information to develop resilient infrastructure and for adaptation decisions. The information desired is at the order of magnitudes finer scales relative to what is typically available from climate analysis and future projections. Urban downscaling refers to developing such climate information at the city (order of 1 – 10 km) and neighborhood (order of 0.1 – 1 km) resolutions from coarser climate products. Developing these higher resolution (finer grid spacing) data needed for assessments typically covering multiyear climatology of past data and future projections is complex and computationally expensive for traditional physics-based dynamical models. In this study, we develop and adopt a novel approach for urban downscaling by generating a general-purpose operator using deep learning. This ‘DownScaleBench’ tool can aid the process of downscaling to any location. The DownScaleBench has been generalized for both in situ (ground- based) and satellite or reanalysis gridded data. The algorithm employs an iterative super-resolution convolutional neural network (Iterative SRCNN) over the city. We apply this for the development of a high-resolution gridded precipitation product (300 m) from a relatively coarse (10 km) satellite-based product (JAXA GsMAP). The high-resolution gridded precipitation datasets is compared against insitu observations for past heavy rain events over Austin, Texas, and shows marked improvement relative to the coarser datasets relative to cubic interpolation as a baseline. The creation of this Downscaling Bench has implications for generating high-resolution gridded urban meteorological datasets and aiding the planning process for climate-ready cities.

Singh, Manmeet (ORCID:0000000233747149)↗
DOE OSTI2021

Early season prediction of within-field crop yield variability by assimilating CubeSat data into a crop model

Accurate early season predictions of crop yield at the within-field scale can be used to address a range of crop production, management, and precision agricultural challenges. While the remote sensing of within-field insights has been a research goal for many years, it is only recently that observations with the required spatio-temporal resolutions, together with efficient assimilation methods to integrate these into modeling frameworks, have become available to advance yield prediction efforts. Here we explore a yield prediction approach that combines daily high-resolution CubeSat imagery with the APSIM crop model. The approach employs APSIM to train a linear regression that relates simulated yield to simulated leaf area index (LAI). That relationship is then used to identify the optimal regression date at which the LAI provides the best prediction of yield: in this case, approximately 14 weeks prior to harvest. Instead of applying the regression on satellite imagery that is coincident, or closest to, the regression date, our method implements a particle filter that integrates CubeSat-based LAI into APSIM to provide end-of-season high-resolution (3 m) yield maps weeks before the optimal regression date. The approach is demonstrated on a rainfed maize field located in Nebraska, USA, where suitable collections of both imagery and in-situ data were available for assessment. The procedure does not require in-field data to calibrate the regression model, with results showing that even with a single assimilation step, it is possible to provide yield estimates with good accuracy up to 21 days before the optimal regression date. Yield spatial variability was reproduced reasonably well, with a strong correlation to independently collected measurements (R 2 = 0.73 and rRMSE = 12%). When the field averaged yield was compared, our approach reduced yield prediction error from 1 Mg/ha (control case based on a calibrated APSIM model), to 0.5 Mg/ha (using satellite imagery alone), and then to 0.2 Mg/ha (results with assimilation up to three weeks prior to the optimal regression date). Such a capacity to provide spatially explicit yield predictions early in the season has considerable potential to enhance digital agricultural goals and improve end-of-season yield predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Can upscaling ground nadir SIF to eddy covariance footprint improve the relationship between SIF and GPP in croplands?

Ground solar-induced chlorophyll fluorescence (SIF) is important for the mechanistic understanding of the dynamics of vegetation gross primary production (GPP) at fine spatiotemporal scales. However, eddy covariance (EC) observations generally cover larger footprint areas than ground SIF observations (a bare fiber with nadir), and this footprint mismatch between nadir SIF and GPP could complicate the canopy SIF-GPP relationships. Here, we upscaled nadir SIF observations to EC footprint and investigated the change in SIF-GPP relationships after the upscaling in cropland. We included 13 site-years data in our study, with seven site-years corn, four site-years soybeans, and two site-years miscanthus, all located in the US Corn Belt. All sites’ crop nadir SIF observations collected from the automated FluoSpec2 system (a hemispheric-nadir system) were upscaled to the GPP footprint-based SIF using vegetation indices (VIs) calculated from high spatiotemporal satellite reflectance data. We found that SIF-GPP relationships were not substantially changed after upscaling nadir SIF to GPP footprint at our crop sites planted with corn, soybean, and miscanthus, with R 2 change after the upscaling ranging from -0.007 to 0.051 and root mean square error (RMSE) difference from -0.658 to 0.095 umol m –2 s –1 relative to original nadir SIF-GPP relationships across all the site-years. The variation of the SIF-GPP relationship within each species across different site-years was similar between the original nadir SIF and upscaled SIF. Different VIs, EC footprint models, and satellite data led to marginal differences in the SIF-GPP relationships when upscaling nadir SIF to EC footprint. Furthermore, our study provided a methodological framework to correct this spatial mismatch between ground nadir SIF and GPP observations for croplands and potentially for other ecosystems. Our results also demonstrated that the spatial mismatch between ground nadir SIF and GPP might not significantly affect the SIF-GPP relationship in cropland that are largely homogeneous.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Interannual variability of spring and summer monsoon growing season carbon exchange at a semiarid savanna over nearly two decades

Eddy covariance measurements of land-atmosphere energy, carbon, and water exchange now span multiple decades at some sites, supporting an improved understanding of flux interannual variability (IAV) and its ecophysiological and physical controls. Most eddy covariance IAV studies have focused on temperate forest ecosystems, where carbon fluxes are large and flux records are longest – but also where IAV is much lower than in dryland regions, which have been identified as an essential driver of the trend and variability in the global terrestrial carbon sink. In this study, we leveraged 19 years of continuous micrometeorological measurements at the AmeriFlux US-SRM mesquite savanna site in southern Arizona, USA to quantify the IAV, trends, and drivers of carbon fluxes during the distinct spring and summer growing seasons. We also assessed the ability of modern satellite and land surface models to capture the IAV of seasonal water and carbon fluxes. Annual net ecosystem production (NEP) was small and highly variable (23 +/- 64 gC m –2 yr –1 ). Precipitation and associated measures of water availability determined most of the variability in NEP, largely through their influence on annual and seasonal gross ecosystem productivity (GEP) as opposed to ecosystem respiration (ER). Root-zone soil moisture captured between 73% (spring) and 85% (summer) of GEP variability and between 73% (spring) and 58% (summer) of ER variability. Throughout the study period, soil moisture and greenness increased with associated increases in GEP, ER and NEP. These trends were strongly influenced by very productive and wet summer growing seasons during the last two years, which were characterized by abundant understory grass cover. Typically, less than half of the variability in growing season GEP and evapotranspiration was captured by satellite-based estimates and land surface model simulations with local site forcing and calibration, highlighting the ongoing utility of long-term datasets to support careful model testing and improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mapping tree height in complex terrain of northern China using ultra-high-resolution images

Tree height is a key parameter for estimating forest biomass and carbon sequestration. In recent years, notable progress has been made in mapping tree height using satellite imagery. However, existing tree height products show low accuracy in mountainous and complex terrains, and few studies typically addressed tree height estimations in mountain areas. This study examines the Mentougou district of Beijing, China, characterized by complex terrain and mountainous landscapes. We analyzed two methods for estimating tree height: one using only spectral features and another combining spectral features with topographic factors (elevation, slope, aspect). We used 3-m resolution PlanetScope 8-band multispectral imagery, with 710 field-measured individual tree heights averaged to obtain 471 pixel-level tree height values as ground-truth, to develop tree height prediction models using eXtreme Gradient Boosting (XGBoost), Random Forest (RF), and Gradient Boosting Machine (GBM) models. The results show that the XGBoost model consistently presented the highest accuracy for both methods evaluated. Specifically, the XGBoost model that combined spectral data with elevation and slope variables with an R² of 0.75 and an RMSE of 2.69 m. Using the XGBoost model, we generated the tree height map for the Mentougou area at 3 m resolution, showing tree heights ranging from 0.5 to 30.4 m, and the model’s prediction error standard deviations ranged from 2.50 to 4.71 m, indicating reliable performance across varied terrain. Additionally, we compared and evaluated the global tree height products, identifying limitations in the accuracy within complex terrains. This study demonstrates the potential for accurately predicting tree heights by combining high-resolution multispectral satellites with a terrain factor modeling approach.

Complex terrain↗
DOE OSTI2022

Assessing mass balance-based inverse modeling methods via a pseudo-observation test to constrain NO x emissions over South Korea

This study constrains NO x emissions by utilizing mass balance–based inverse modeling methods and examines the optimization of inverse modeling conditions through a pseudo-observation test for South Korea, which has complex topography. Here, we applied mass balance–based inverse modeling methods such as Basic Mass Balance (BMB), Finite Difference Mass Balance (FDMB), and Iterative Finite Difference Mass Balance (IFDMB). We performed a numerical simulation of air quality using the Community Multi-scale Air Quality (CMAQ) model to calculate the NO 2 column density required for the inverse modeling. The pseudo-observation test was performed according to season, modeling resolution, and regridding methodology of satellite observation data to identify various conditions while applying the inverse modeling for South Korea. Comparing the inverse modeling results from the BMB, FDMB, and IFDMB methods, IFDMB was the most effective method in constraining NO x emissions in the South Korean region since it minimized smearing effects (i.e., transport-induced errors) through iterative calculations. The accuracy of the constrained NO x emissions using mass balance–based inversions in South Korea was the highest in the summer due to the minimized smearing effects, and the 9 km resolution modeling was the most efficient for inverse modeling of the region. In addition, the results of inverse modeling varied depending on the regridding methods, implying the importance of using a suitable regridding method and modeling resolution. This study used pseudo-observations, but the inversions are expected to be applied based on actual satellite data in the future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Ammonia in northeast Colorado is increasing, rising most quickly in regions close to confined animal feeding operations

The Colorado Front Range urban corridor and nearby agricultural operations are important source regions of atmospheric ammonia (NH 3 ). Upslope flows periodically transport these emissions into Rocky Mountain National Park (RMNP), located 50 km west of the urban corridor, where wet and dry deposition of excess reactive nitrogen (N) impacts ecosystems. Here, we use a combination of in situ passive NH 3 measurements and NH 3 vertical column density retrievals from the Infrared Atmospheric Sounding Interferometer (IASI) to assess variability and changes in NH 3 across three land use categories in the northeast Colorado source region (agricultural, urban, and remote) during the period 2013-2023. A strong seasonal cycle is present across the region with increased NH 3 during summer months. Elevated NH 3 is spatially correlated with the number of permitted animal units in confined animal feeding operations (CAFOs) within 12 km. Ground-level NH 3 concentrations are strongly positively correlated with monthly gridded IASI satellite column densities. Satellite retrievals reveal an increasing trend in NH 3 column amounts of ∼3% per year in agricultural and ∼2% per year in urban sub-regions. The magnitude of the trend observed in NH 3 columns averaged over the agricultural sub-region is > 3 times larger than observed near and over Denver. The largest increases in NH 3 are closely aligned with the distribution of CAFOs. Reductions in particle sulfate associated with declining sulfur dioxide (SO 2 ) emissions could account for only ∼0.1% per year increase in gaseous NH 3 . Wildfire smoke across the region has increased but appears unlikely to explain the majority of the observed NH 3 increase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Mining metatranscriptomes reveals a vast world of viroid-like circular RNAs

Viroids and viroid-like covalently closed circular (ccc) RNAs are minimal replicators that typically encode no proteins and hijack cellular enzymes for replication. The extent and diversity of viroid-like agents are poorly understood. We developed a computational pipeline to identify viroid-like cccRNAs and applied it to 5,131 metatranscriptomes and 1,344 plant transcriptomes. The search yielded 11,378 viroid-like cccRNAs spanning 4,409 species-level clusters, a 5-fold increase compared to the previously identified viroid-like elements. Within this diverse collection, we discovered numerous putative viroids, satellite RNAs, retrozymes, and ribozy-like viruses. Diverse ribozyme combinations and unusual ribozymes within the cccRNAs were identified. Self-cleaving ribozymes were identified in ambiviruses, some mito-like viruses and capsid-encoding satellite virus-like cccRNAs. The broad presence of viroid-like cccRNAs in diverse transcriptomes and ecosystems implies that their host range is far broader than currently known, and matches to CRISPR spacers suggest that some cccRNAs replicate in prokaryotes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

An efficient computational framework for charge density estimation in twisted bilayer graphene

Electronic properties such as band structure and Fermi velocity in low-angle twisted bilayer graphene (TBG) are intrinsically dependent on the atomic structure. Rigid rotation between individual graphene layers provides an approximate description of the bilayer symmetry. Upon relaxation, in-plane displacement of the atoms in low angle TBG causes a change in the symmetry through the enlargement of the AB stacking regions and the reduction in size of AA and SP stacking regions. However, the effect of this in-plane relaxation on the charge density remains unexplored, because the necessary electronic structure calculations of such large supercells of low twist angle TBG are computationally infeasible. Therefore, we develop a computationally efficient framework that enables the exploration of the charge density symmetry of the low twist angle TBG. This framework is based on the Fourier representation of the charge density which presents high intensity Bragg peaks. Here we find that with the decrease of twist angle, low intensity satellite peaks also become apparent. Our framework incorporates these satellite peaks which reveals transformation of symmetry in the charge density distribution from high to low twist angle TBG. One striking outcome is the demonstration of the electron localization in the AA region of low twist angle TBG. Our framework helps to explain the effect of the atomistic relaxation on the charge density distribution and thus, it provides information about exotic electronic properties of low twist angle TBG at a low computational expense.

36 MATERIALS SCIENCE↗