Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Ecosystem functioning”

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 937 records · Page 52

BSEC VPRM 10m Hourly Biogenic Fluxes in Baltimore (2021)

Model outputs from the Vegetation Photosynthesis and Respiration Model (VPRM: version from Horne et al. in prep). Model remote sensing inputs come from Sential 2-derived EVI and LSWI. Model meteorological inputs for two-meter air temperature and shortwave incoming come from the BSEC WRF 2021 Control Run (Foust, W. 2023). Plant functional Types (PFTs) are spatially classified using the Chesapeake Bay Program 2018 land use land cover product. The final biogenic flux (µmol CO2 m^-2 s^-1) outputs of NEE, RESP, and GEE are a weighted average based on the portion of PFTs within the cell. Individual PFT outputs are saved inside PFT directories (e.g., Crops, Grass, etc.) inside the specific month directory. Model outputs are denoted as a negative flux into the land system (i.e., photosynthesis) and a positive flux as a net release into the overlying atmosphere. Respiration (RESP) fluxes are positive and combine heterotrophic (only soil) and autotrophic sources. Gross ecosystem exchange (GEE) is a negative flux driven by only photosynthetic activity from vegetation, and the Net ecosystem exchange (NEE) is the sum of the two (i.e., NEE=RESP+GEE). Data Characteristics Spatial Resolution: 10m Temporal Resolution: Hourly File Format: VPRM_ _BSEC. .tif (Hour is in UTC) For more information on the model results, please email Jason Horne (jph6488@psu.edu). References: Foust, W. (2023). BSEC WRF 2021 Control Run Output (v0.1.0) [Data set]. MSD-LIVE Data Repository. https://data.msdlive.org/records/m0e6m-vvq17

Baltimore↗

Impacts of Forest Fire Ash on Aquatic Mercury Cycling

Mercury (Hg) is a ubiquitous contaminant in the environment and its methylated form, methylmercury (MeHg), poses a worldwide health concern for humans and wildlife, primarily through fish consumption. Global production of forest fire ash, derived from wildfires and prescribed burns, is rapidly increasing due to a warming climate, but their interactions with aqueous and sedimentary Hg are poorly understood. In this study, we compared the differences of wildfire ash with activated carbon and biochar on the sorption of aqueous inorganic Hg and sedimentary Hg methylation. Sorption of aqueous inorganic Hg was greatest for wildfire ash materials (up to 0.21 μg g ₋1 or 2.2 μg g ₋1 C) among all of the solid sorbents evaluated. A similar Hg adsorption mechanism for activated carbon, biochar made of walnut, and wildfire ash was found that involves the formation of complexes between Hg and oxygen-containing functional groups, especially the ₋COO group. Notably, increasing dissolved organic matter from 2.4 to 70 mg C L ₋1 remarkably reduced Hg sorption (up to 40% reduction) and increased the time required to reach Hg–sorbent pseudo-equilibrium. Surprisingly, biochar and wildfire ash, but not activated carbon, stimulated MeHg production during anoxic sediment incubation, possibly due to the release of labile organic matter. Overall, our study indicates that while wildfire ash can sequester aqueous Hg, the leaching of its labile organic matter may promote production of toxic MeHg in anoxic sediments, which has an important implication for potential MeHg contamination in downstream aquatic ecosystems after wildfires.

54 ENVIRONMENTAL SCIENCES↗

Non-growing Season Plant Nutrient Uptake Controls Arctic Tundra Vegetation Composition: Modeling Archive

This Modeling Archive is in support of a NGEE-Arctic publication: Riley et al. (2021) “Non-growing season plant nutrient uptake controls Arctic tundra vegetation composition under future climate” (Riley et al., 2021). The dataset contains “ecosys” model outputs reported in Riley et al. (2021) that examines the role of non-growing season (NGS) nutrient dynamics on 21st century vegetation composition. The study highlighted the importance of these nutrient dynamics on the emergent plant functional type distributions, focusing on their role in shrub expansion. The study showed that ignoring NGS nutrient dynamics led to a dramatic under-prediction of shrub expansion by year 2100, motivating the need to include these processes in Earth System Models. Included are modeled net primary production (NPP) and nitrogen (N) uptake for two scenarios: (1) baseline and (2) no NGS N uptake. The output is provided (1) at 25 km resolution across the North American tundra; (2) for years 2012 through 2100; and (3) PFT-specific. Dataset contains 20 NetCDF files in a zip file and one pdf.The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy’s Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy’s Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

UAS remote sensing (DJI Phantom 4 RTK platform): RGB orthomosaic, digital surface and canopy height models, plant functional type map, Seward Peninsula, Alaska, 2019

Airborne remote sensing data collected using a DJI Phantom 4 RTK unoccupied aerial system (UAS) ? operated by the University of Maine Wheatland Geospatial Lab (https://wheatlandlab.org/). This package includes data from 17 flights flown over the NGEE-Arctic Council, Kougarok, Kougarok Mile80, Teller, and Teller Mile32 sites in July 2019. The Phantom 4 RTK is a drone platform that collects very high spatial resolution optical red/green/blue (RGB) imagery. Derived image products include point cloud, ortho-mosaiced RGB, a digital surface model (DSM) using the structure from motion (SfM) technique, and a canopy height model (CHM). Unprocessed and processed data products (1,000+ files) are included in this package (processing levels 0-3). Data and metadata are provided as text (*.txt, *.json, *hdr,), ENVI image file (*.dat), point cloud (*.laz) and image (*.jpg, *.tif, *png) formats. The Next-Generation Ecosystem Experiments: Arctic (NGEE Arctic), was a research effort to reduce uncertainty in Earth System Models by developing a predictive understanding of carbon-rich Arctic ecosystems and feedbacks to climate. NGEE Arctic was supported by the Department of Energy's Office of Biological and Environmental Research. The NGEE Arctic project had two field research sites: 1) located within the Arctic polygonal tundra coastal region on the Barrow Environmental Observatory (BEO) and the North Slope near Utqiagvik (Barrow), Alaska and 2) multiple areas on the discontinuous permafrost region of the Seward Peninsula north of Nome, Alaska. Through observations, experiments, and synthesis with existing datasets, NGEE Arctic provided an enhanced knowledge base for multi-scale modeling and contributed to improved process representation at global pan-Arctic scales within the Department of Energy's Earth system Model (the Energy Exascale Earth System Model, or E3SM), and specifically within the E3SM Land Model component (ELM).

54 ENVIRONMENTAL SCIENCES↗

Synoptic NO3 in Slate River Watershed, Colorado (2022)

The synoptic nitrate (NO3) dataset in the Slate River Watershed, Colorado consists of NO3 data collected at 19 locations three times during the summer of 2022. Stream samples were collected in early summer (early July), mid summer (late August), and late summer (late September). The samples include mainstem, tributary, and point source input water samples. These data were collected to evaluate spatiotemporal variability in stream NO3 during the summer, and evaluate anthropogenic controls on stream NO3 dynamics. This data package contains: (1) a csv of all NO3 samples and (2) a csv of locations for each sampling site. The dataset additionally includes a file-level metadata (flmd.csv) file that lists each file contained in the dataset with associated metadata; and a data dictionary (dd.csv) file that contains column/row headers used throughout the files along with a definition, units, and data type. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > BIOSPHERE > ECOSYSTEMS > FRESHWATE↗

Understanding SYCL Portability for Pseudorandom Number Generation: a Case Study with Gene-Expression Connectivity Mapping

Towards the goal of improving functional and performance portability of SYCL, we study a bioinformatics application that has been accelerated with CUDA and fast pseudorandom number generation on a GPU. We describe the experience of migrating pseudorandom number generation from CUDA to SYCL, evaluate the performance of pseudorandom number generators using the CUDA random number generation library, suggest the support of the XORWOW pseudorandom number generator in the oneAPI math kernel library (oneMKL) interface for performance portability, and identify the performance gap using the MKL interface in SYCL that supports pseudorandom number generation with third-party libraries. We hope that the results are valuable for the development of the SYCL ecosystem.

Jin, Zheming↗

Melanization slows the rapid movement of fungal necromass carbon and nitrogen into both bacterial and fungal decomposer communities and soils

Microbial necromass contributes significantly to both soil carbon (C) persistence and ecosystem nitrogen (N) availability, but quantitative estimates of C and N movement from necromass into soils and decomposer communities are lacking. Additionally, while melanin is known to slow fungal necromass decomposition, how it influences microbial C and N acquisition as well as elemental release into soils remains unclear. Here, we tracked decomposition of isotopically labeled low and high melanin fungal necromass and measured 13 C and 15 N accumulation in surrounding soils and microbial communities over 77 d in a temperate forest in Minnesota, USA. Mass loss was significantly higher from low melanin necromass, corresponding with greater 13 C and 15 N soil inputs. A taxonomically and functionally diverse array of bacteria and fungi was enriched in 13 C and/or 15 N at all sampling points, with enrichment being consistently higher on low melanin necromass and earlier in decomposition. Similar patterns of preferential C and N enrichment of many bacterial and fungal genera early in decomposition suggest that both microbial groups co-contribute to the rapid assimilation of resource-rich soil organic matter inputs. While overall richness of taxa enriched in C was higher than in N for both bacteria and fungi, there was a significant positive relationship between C and N in co-enriched taxa. Collectively, our results demonstrate that melanization acts as a key ecological trait mediating not only fungal necromass decomposition rate but also necromass C and N release and that both elements are rapidly co-utilized by diverse bacterial and fungal decomposers in natural settings.

54 ENVIRONMENTAL SCIENCES↗

A Modified Vegetation Photosynthesis and Respiration Model (VPRM) for the Eastern USA and Canada, Evaluated With Comparison to Atmospheric Observations and Other Biospheric Models

Atmospheric CO 2 measurements from a dense surface network can help to evaluate terrestrial biosphere model (TBM) simulations of Net Ecosystem Exchange (NEE) with two key benefits. First, gridded CO 2 flux estimates can be evaluated over regional scales, not possible using flux tower observations at discrete locations for model evaluation. Second, TBM ability to explain atmospheric CO 2 fluctuations due to the biosphere can be directly tested, an important objective for anthropogenic emissions monitoring using atmospheric observations. Here, we customize the Vegetation Photosynthesis and Respiration Model (VPRM) for an eastern North American domain with strong biological activity upwind of urban areas. Parameters are optimized using flux tower observations from a historical database with sites in (and near) the domain. In addition, the respiration model (originally a linear function of temperature) is modified to account for impacts of changing foliage, non-linear temperature, and water stress. Flux estimates from VPRM, the Carnegie-Ames-Stanford Approach (CASA) model and the Simple Biosphere Model v4 (SiB4), are convolved with footprints from atmospheric transport models for evaluation with CO 2 observations at 21 towers in the domain, with roughly half of the towers used here for the first time. Results show that the new respiration model in VPRM helps to correct a growing season sink bias in the atmosphere associated with underestimated summertime respiration using the original model with annual parameters. The new VPRM also better explains fine-scale atmospheric CO 2 variability compared to other TBMs, due to higher resolution diagnostic phenology, the new respiration model, domain-specific parameters, and high-quality input data sets.

54 ENVIRONMENTAL SCIENCES↗

Disentangling top-down drivers of mortality underlying diel population dynamics of Prochlorococcus in the North Pacific Subtropical Gyre

Photosynthesis fuels primary production at the base of marine food webs. Yet, in many surface ocean ecosystems, diel-driven primary production is tightly coupled to daily loss. This tight coupling raises the question: which top-down drivers predominate in maintaining persistently stable picocyanobacterial populations over longer time scales? Motivated by high-frequency surface water measurements taken in the North Pacific Subtropical Gyre (NPSG), we developed multitrophic models to investigate bottom-up and top-down mechanisms underlying the balanced control of Prochlorococcus populations. Here, we find that incorporating photosynthetic growth with viral- and predator-induced mortality is sufficient to recapitulate daily oscillations of Prochlorococcus abundances with baseline community abundances. In doing so, we infer that grazers in this environment function as the predominant top-down factor despite high standing viral particle densities. The model-data fits also reveal the ecological relevance of light-dependent viral traits and non-canonical factors to cellular loss. Finally, we leverage sensitivity analyses to demonstrate how variation in life history traits across distinct oceanic contexts, including variation in viral adsorption and grazer clearance rates, can transform the quantitative and even qualitative importance of top-down controls in shaping Prochlorococcus population dynamics.

54 ENVIRONMENTAL SCIENCES↗

Environmental impacts of dam reservoir filling in the East Amazon

Mitigating the environmental impacts caused by hydroelectric dams is a worldwide challenge. The aquatic ecosystem is the most impacted during the reservoir filling phase, yet biogeochemical dynamics at this stage are not well-studied. Here we evaluate water quality and hydraulic parameters in the Araguari River (Amapá/Brazil) during the filling of the Ferreira Gomes Hydroelectric Power Plant reservoir (UHEFG). Five field campaigns were performed from July 2014 to August 2015 across nine sample sites within the reservoir (P1, P2, …P9) and only one downstream (P10). The following key variables were monitored: Trophic State Index (TSI), Total Coliforms (TC), E. coli (EC) and Chlorophyll-a (Chl-a), along with physical and chemical variables (Temperature, Suspended Solids, Total Dissolved Solids, Electrical Conductivity, Turbidity, Color, pH, DO, Al3+, NH4+, Cl-, Mg1+, Ca2+, total Phosphorous, NH3, NO3-, and SO42). Besides that, hydraulic-operational variables were also investigated: inflows (Qa), outflows (Qd), and variation in reservoir volume (Vol%). Multiple Regression Analyses showed that the key parameters were both significantly influenced by physio-chemical and hydraulic variables (0.46=R2adj=0.99, p<0.05). The DO showed significant spatial variation, being influenced by the turbulence from the Coaracy Nunes dam (UHECN) upstream and the UHEFG dam downstream. The Vol% influenced the TSI, which ranged from oligotrophic to hypertrophic and eventually stabilized at mesotrophic. The levels of TSI, TC, and Chl-a decreased and the level of E. coli increased (p<0.05) as a function of Vol%. A Cluster Analysis showed the formation of three spatial groups - two inside the reservoir and one downstream (P10). This suggests that in the rainy season or transition season, the hydraulic residence time in the reservoir is very low (16 = thr = 36 h) when hydrodynamic processes are dominant. In the dry season (thr ˜ 1 month), biogeochemical and hydrodynamic processes occur at similar time scales to the nutrient and microbial abundance. These results confirm the hypothesis that the filling phase has significant impacts on the key parameters of the water quality (p<0.05). We conclude that the filling phase of the UHEFG reservoir generated significant environmental impacts, which have repercussions even retrospectively.

Amazon, biogeochemistry, ecosystem, extreme event,↗

A Mixed Length Scale Model for Migrating Fluvial Bedforms

With the expansion of hydropower, in-stream converters, flood-protection infrastructures, and growing concerns on deltas fragile ecosystems, there is a pressing need to evaluate and monitor bedform sediment mass flux. It is critical to estimate real-time bedform size and migration velocity and provide a theoretical framework to convert easily accessible time histories of bed elevations into spatially evolving patterns. In this study, we collected spatiotemporally resolved bathymetries from laboratory flumes and the Colorado River in statistically steady, homogeneous, subcritical flow conditions. Wave number and frequency spectra of bed elevations show compelling evidence of scale-dependent velocity for the hierarchy of migrating bedforms observed in the laboratory and field. New scaling laws were applied to describe the full range of migration velocities as function of two dimensionless groups based on the bed shear velocity, sediment diameter, and water depth. Further simplification resulted in a mixed length scale model estimating scale-dependent migration velocities, without requiring bedform classification or identification.

58 GEOSCIENCES↗

Organizing principles for vegetation dynamics

Understanding vegetation dynamics is very challenging because of the multitude of contributing processes at 64widely different spatial and temporal scales. In this Perspective we propose that understanding of vegetation dynamics can be improved, permitting better predictions, based on organizing principles that constrain plant and ecosystem be haviour: natural selection, self-organization, and entropy maximization. Although these ideas are increasingly used,a limited common understanding of their theoretical basis has prevented their full potential to be realized. We explain the power of natural selection-based optimality to predict photosynthesis and carbon allocation responses to multiple environmental drivers, and how individual plasticity leads to the predictable self-organization of forest canopies. We show how models of natural selection acting on a few key traits can generate realistic plant communities, and how entropy maximization can distinguish between stochastic and deterministic drivers of vegetation patterns. In combination with empirical exploration of patterns in plant functional variation, these principles can accelerate the development of dynamic vegetation models as well as trait-based ecology resting on strengthened theoretical and empirical foundations.

Geosciences↗

A Near-Real-Time Model for Predicting Electricity Disruptions in Texas During Winter Storms

There has been an increase in extreme weather events, posing a threat to power grid systems, potentially influenced by factors such as population growth, changes in ecosystems, land cover, and land use in the service area, as well as the growth of certain vegetation types. This research seeks to develop a predictive model to mitigate potential damages caused by future winter storms. This research utilizes the Light Gradient Boosting Machine (LightGBM), incorporating the number of power outages experienced at the county level, geographic details, weather information, and lagged outage and lagged weather data. The developed models were broadly divided into two groups, with six models in each group - one group without optimization and another with optimization, totaling 12 trained models. For model optimization, Bayesian optimization was employed using Root Mean Squared Error (RMSE) as the objective function. In results, when comparing Group 2 (the optimized group) with Group 1 (the non-optimized group), it was found that optimization did not always lead to a reduction in RMSE and Mean Absolute Error (MAE). However, in terms of Mean Directional Accuracy (MDA), while all results in Group 1 were below the baseline accuracy of 0.33, all results in Group 2 exceeded 0.33, with some cases showing an increase of more than three times the baseline. The results indicated that, in the optimized model group, Population and Pressure were the most influential factors when using current weather data and geographical information. When using lagged data, lagged recorded outages and lagged Pressure emerged as the most significant factors. Among the 12 developed models, the L-1-2-O model showed the lowest RMSE and MAE, as well as the highest accuracy, with values of 390.62 households and 168.13 households, respectively. To normalize the RMSE and MAE values, each metric was divided by the average number of households among the counties in Texas. For the L-1-2-O model, the scaled RMSE was 0.88% and the scaled MAE was 0.38%. In terms of MDA, which indicates the accuracy of the prediction direction, the L-1-O model achieved the highest score of 0.41. Although this study focused on Texas, which suffered the greatest impact from the winter storms in 2021, with additional validation, the methodology used in this research could be applied to other regions.

Lee, Jangjae [Texas A & M Univ., College Station, ↗

Bioenergy cropping systems shape ant community composition and functional roles

R code and partial data for Haan, Helms, & Landis 2023, Bioenergy cropping systems shape ant community composition and functional roles, to be published in Frontiers in Conservation Science. The ant species list we use to generate the species matrix was uploaded previously (10.5281/zenodo.8215581) in association with another manuscript (Haan et al. 2023, Science Advances, Contrasting effects of bioenergy crops on biodiversity). Therefore here we include a small table with trait information for each species along with the R code used specifically for this manuscript.

ants↗

Large-scale genome sequencing of mycorrhizal fungi provides insights into the early evolution of symbiotic traits

Mycorrhizal fungi are mutualists that play crucial roles in nutrient acquisition in terrestrial ecosystems. Mycorrhizal symbioses arose repeatedly across multiple lineages of Mucoromycotina, Ascomycota, and Basidiomycota. Considerable variation exists in the capacity of mycorrhizal fungi to acquire carbon from soil organic matter. Here, we present a combined analysis of 135 fungal genomes from 73 saprotrophic, endophytic and pathogenic species, and 62 mycorrhizal species, including 29 new mycorrhizal genomes. This study samples ecologically dominant fungal guilds for which there were previously no symbiotic genomes available, including ectomycorrhizal Russulales, Thelephorales and Cantharellales. Our analyses show that transitions from saprotrophy to symbiosis involve (1) widespread losses of degrading enzymes acting on lignin and cellulose, (2) co-option of genes present in saprotrophic ancestors to fulfill new symbiotic functions, (3) diversification of novel, lineage-specific symbiosis-induced genes, (4) proliferation of transposable elements and (5) divergent genetic innovations underlying the convergent origins of the ectomycorrhizal guild.

59 BASIC BIOLOGICAL SCIENCES↗

Assessing CLUBB PDF Closure Assumptions for a Continental Shallow-to-Deep Convective Transition Case Over Multiple Spatial Scales

Assumed-PDF (probability density function) higher-order turbulence closures (APHOCs) are now widely used for parameterizing boundary layer turbulence and shallow convection in Earth system models (ESMs). A better understanding of the resolution-dependent behavior of APHOCs is essential for improving the performance of next-generation ESMs with intended horizontal resolutions finer than 10 km. In this study, we evaluate the PDF family of Analytic double-Gaussian 1 implemented in Cloud Layers Unified By Binormals (CLUBB) over a range of spatial scales (Dx) from 2 to 100 km. A 120-km-wide large eddy simulation (LES) for a continental convection case during 2016 Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign serves as benchmark to evaluate the PDF closure using an off-line approach. We find during the shallow convection period, the CLUBB PDF closure tends to produce positive biases of cloud properties and liquid water flux near cloud base for all scales of analysis. It produces negative biases for these variables near cloud top that are more severe for Dx larger than 25 km. Results show that replacing the CLUBB-parameterized moisture and temperature skewnesses with LES-derived ones can fix most of the biases if clipping of input moments is allowed to prevent the occurrence of unrealizable solutions. Overall, the performance of the PDF closure is better for smaller Dx = 2–5 km than for larger Dx = 50–100 km; for a given grid spacing, it is better when the convective clouds become deeper in the late afternoon. Likely causes for the resolution dependence and implications for improving the PDF closure are discussed.

54 ENVIRONMENTAL SCIENCES↗

Decrypting bacterial polyphenol metabolism in an anoxic wetland soil

Abstract Microorganisms play vital roles in modulating organic matter decomposition and nutrient cycling in soil ecosystems. The enzyme latch paradigm posits microbial degradation of polyphenols is hindered in anoxic peat leading to polyphenol accumulation, and consequently diminished microbial activity. This model assumes that polyphenols are microbially unavailable under anoxia, a supposition that has not been thoroughly investigated in any soil type. Here, we use anoxic soil reactors amended with and without a chemically defined polyphenol to test this hypothesis, employing metabolomics and genome-resolved metaproteomics to interrogate soil microbial polyphenol metabolism. Challenging the idea that polyphenols are not bioavailable under anoxia, we provide metabolite evidence that polyphenols are depolymerized, resulting in monomer accumulation, followed by the generation of small phenolic degradation products. Further, we show that soil microbiome function is maintained, and possibly enhanced, with polyphenol addition. In summary, this study provides chemical and enzymatic evidence that some soil microbiota can degrade polyphenols under anoxia and subvert the assumed polyphenol lock on soil microbial metabolism.

59 BASIC BIOLOGICAL SCIENCES↗

A chromosome-level genome assembly of the Chinese cork oak (Quercus variabilis)

Quercus variabilis (Fagaceae) is an ecologically and economically important deciduous broadleaved tree species native to and widespread in East Asia. It is a valuable woody species and an indicator of local forest health, and occupies a dominant position in forest ecosystems in East Asia. However, genomic resources from Q. variabilis are still lacking. Here, we present a high-quality Q. variabilis genome generated by PacBio HiFi and Hi-C sequencing. The assembled genome size is 787 Mb, with a contig N50 of 26.04 Mb and scaffold N50 of 64.86 Mb, comprising 12 pseudo-chromosomes. The repetitive sequences constitute 67.6% of the genome, of which the majority are long terminal repeats, accounting for 46.62% of the genome. We used ab initio , RNA sequence-based and homology-based predictions to identify protein-coding genes. A total of 32,466 protein-coding genes were identified, of which 95.11% could be functionally annotated. Evolutionary analysis showed that Q. variabilis was more closely related to Q. suber than to Q. lobata or Q. robur. We found no evidence for species-specific whole genome duplications in Quercus after the species had diverged. This study provides the first genome assembly and the first gene annotation data for Q. variabilis. These resources will inform the design of further breeding strategies, and will be valuable in the study of genome editing and comparative genomics in oak species.

Han, Biao↗