Search

Eulerian and Lagrangian numerical moisture tracking models, which are primarily used to quantify moisture contributions from global sources to specific regions, play a crucial role in hydrology and (paleo)climatology studies on the Tibetan Plateau (TP). Despite their widespread applications in the TP region, potential discrepancies in their moisture tracking results and their underlying causes remain unexplored. In this study, we compare the most widely used Eulerian and Lagrangian moisture tracking models over the TP, i.e., WAM2layers (the Water Accounting Model – 2 layers) and FLEXPART-WaterSip (the FLEXible PARTicle dispersion model coupled with the “WaterSip” moisture source diagnostic method), specifically focusing on a basin governed by the Indian summer monsoon (Yarlung Zangbo River basin, YB) and a westerly-dominated basin (upper Tarim River basin, UTB). Compared to the bias-corrected FLEXPART-WaterSip, WAM2layers generally estimates higher moisture contributions from westerly-dominated and distant sources but lower contributions from local recycling and nearby sources downwind of the westerlies. These differences become smaller with higher spatial and temporal resolutions of forcing data in WAM2layers. A notable advantage of WAM2layers over FLEXPART-WaterSip is its closer alignment of estimated moisture sources with actual evaporation, particularly in source regions with complex land–sea distributions. However, the evaporation biases in FLEXPART-WaterSip can be partly corrected through calibration with actual surface fluxes. For moisture tracking over the TP, we recommend using high-resolution forcing datasets, prioritizing temporal resolution over spatial resolution for WAM2layers, while for FLEXPART-WaterSip, we suggest applying bias corrections to optimize the filtering of precipitation particles and adjust evaporation estimates.

Using the regionally refined mesh (RRM) configuration of the US Department of Energy's Simple Cloud-Resolving Energy Exascale Earth System Model (E3SM) Atmosphere Model (SCREAM), we simulate and evaluate four meteorologically distinct atmospheric river events over California. We test five different RRM configurations, each differing in terms of the areal extent of the refined mesh and the resolution (ranging from 800 m to 3.25 km). We find that SCREAM RRM generally has a good representation of the AR-generated precipitation in CA, even for the control simulation which has a very small 3 km refined patch, and is able to capture the fine-scale regional distributions that are controlled largely by the fine-scale topography of the state. It is found that SCREAM generally has a wet bias over topography, most prominently over the Sierra Nevada mountain range, with a corresponding dry bias on the lee side. We find that refining the resolution beyond 3 km (specifically 1.6 km and 800 m) has virtually no benefit towards reducing systematic precipitation biases but that improvements can be found when increasing the areal extent of the upstream refined mesh. However, these improvements are relatively modest and only realized if the size of the refined mesh is expanded to the scale where employing RRM no longer achieves the substantial cost benefit it was intended for.

Accurate and timely global evapotranspiration (ET) data are crucial for agriculture, water resource management, and drought forecasting. Although numerous satellite-based ET products are available, few offer near-real-time data. For instance, products like NASA's ECOsystem Spaceborne Thermal Radiometer Experiment mounted on the International Space Station (ECOSTRESS) and MOD16 face challenges such as uneven coverage and delays exceeding 1 week in data availability. In this study, we refined the Variation of the Standard Evapotranspiration Algorithm (VISEA) by fully integrating satellite-based data, e.g., European Centre for Medium-Range Weather Forecasts ERA5-Land shortwave radiation (which includes satellite remote sensing data within its assimilation system) and MODIS land surface data (which include surface reflectance, temperature and/or emissivity, land cover, vegetation indices, and albedo as inputs). This enables VISEA to provide near-real-time global daily ET estimates with a maximum delay of 1 week at a resolution of 0.05°. Its accuracy was assessed globally using observation data from 149 flux towers across 12 land cover types and comparing them with five other satellite-based ET products and Global Precipitation Climatology Centre (GPCC) data. The results indicate that VISEA provides accurate ET estimates that are comparable to existing products, achieving a mean correlation coefficient (R) of about 0.6 and an RMSE of 1.4 mm d−1. Furthermore, we demonstrated VISEA's utility in drought monitoring during a drought event in the Yangtze River basin in 2022 in which ET changes correlated with precipitation. The near-real-time capability of VISEA is, thus, especially valuable in meteorological and hydrological applications for coordinating drought relief efforts. The VISEA ET dataset is available at https://doi.org/10.11888/Terre.tpdc.300782 (Huang, 2023a).

Relict subtropical coniferous forests in China face severe fragmentation, resulting in declining populations, and some are under significant threat from invasive alien species. Despite the crucial importance of understanding forest dynamics, knowledge gaps persist, particularly regarding the impact of invasive plants on vulnerable natives like Keteleeria evelyniana. In this study, we investigated the impact of invasive plants on the regeneration of forests dominated by K. evelyniana, a subtropical relict species in southwestern China. For this purpose, we characterized forest dynamics of 160 forest plots featuring K. evelyniana as the primary dominant species and determined whether the presence of invasive plants was correlated with regeneration of K. evelyniana. We identified four distinct forest types in which K. evelyniana was dominant. We found that radial growth of K. evelyniana trees is faster in younger age-classes today than it was for older trees at the same age. The population structure of K. evelyniana in each forest type exhibited a multimodal age-class distribution. However, three forest types lacked established saplings younger than 10 years old, a situation attributed to the dense coverage of the invasive alien Ageratina adenophora. This invasive species resulted in a reduction of understory species diversity. Additionally, our analysis uncovered a significant negative correlation in phylogenetic relatedness (net relatedness index) between native and invasive alien plant species in eastern Yunnan. This suggests closely related invasive species face heightened competition, hindering successful invasion. Taken together, our findings indicate that successful establishment and habitat restoration of K. evelyniana seedling/saplings require effective measures to control invasive plants.

The spatial heterogeneity related to complex topography in California demands high-resolution (< 5 km) modeling, but global convection-permitting climate models are computationally too expensive to run multi-decadal simulations. We developed a 3.25 km California climate modeling framework by leveraging regional mesh refinement (CARRM) using the U.S. Department of Energy (DOE)'s global Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) version 0. Four 5-year time periods (2015–2020, 2029–2034, 2044–2049, and 2094–2099) were simulated by nudging CARRM outside California to 1° coupled simulation of E3SMv1 under the Shared Socioeconomic Pathways (SSP)5-8.5 future scenario. The 3.25 km grid spacing adds considerable value to the prediction of the California climate changes, including more realistic high temperatures in the Central Valley and much improved spatial distributions of precipitation and snowpack in the Sierra Nevada and coastal stratocumulus. Under the SSP5-8.5 scenario, CARRM simulation predicts widespread warming of 6–10 °C over most of California, a 38 % increase in statewide average 30 d winter–spring precipitation, a near-complete loss of the alpine snowpack, and a sharp reduction in shortwave cloud radiative forcing associated with marine stratocumulus by the end of the 21st century. We note a climatological wet precipitation bias for the CARRM and discuss possible reasons. We conclude that SCREAM RRM is a technically feasible and scientifically valid tool for climate simulations in regions of interest, providing an excellent bridge to global convection-permitting simulations.

This article describes a modular ensemble-based data assimilation (DA) system which is developed for an integrated surface–subsurface hydrological model. The software environment for DA is the Parallel Data Assimilation Framework (PDAF), which provides various assimilation algorithms like the ensemble Kalman filters, non-linear filters, 3D-Var and combinations among them. The integrated surface–subsurface hydrological model is HydroGeoSphere (HGS), a physically based modelling software for the simulation of surface and variably saturated subsurface flow, as well as heat and mass transport. The coupling and capabilities of the modular DA system are described and demonstrated using an idealised model of a geologically heterogeneous alluvial river–aquifer system with drinking water production via riverbank filtration. To demonstrate its modularity and adaptability, both single and multivariate assimilations of hydraulic head and soil moisture observations are demonstrated in combination with individual and joint updating of multiple simulated states (i.e. hydraulic heads and water saturation) and model parameters (i.e. hydraulic conductivity). With the integrated model and this modular DA framework, we have essentially developed the hydrologically and DA-wise robust toolbox for developing the basic model for operational management of coupled surface water–groundwater resources.

This study documents clouds simulated by the Energy Exascale Earth System Model (E3SM) version 2 (E3SMv2) and attempts to understand what causes the model behavior change in clouds relative to E3SMv1. This is done by analyzing the last 30-year (1985–2014) data from the 165-year historical simulations using E3SMv1 and v2 and four sensitivity tests to isolate the impact of changes in model parameter choices in its turbulence, shallow convection, and cloud macrophysics parameterization (Cloud Layers Unified By Binormals, CLUBB); microphysical parameterization (MG2); and deep-convection scheme (ZM), as well as model physics changes in convective triggering. It is shown that E3SMv2 significantly improves the simulation of subtropical coastal stratocumulus clouds and clouds with optical depth larger than 3.6 over the stratocumulus-to-cumulus transition regimes, where the shortwave cloud radiative effect (SWCRE) is also improved, and the Southern Ocean (SO) while seeing an overall slight degradation in low clouds over other tropical and subtropical oceans. The better performance in E3SMv1 over those regions is partially due to error compensation between its simulated optically thin and intermediate low clouds for which E3SMv2 actually improves simulation of optically intermediate low clouds. Sensitivity tests indicate that the changes in low clouds are primarily due to the tuning done in CLUBB. The impact of the ZM tuning is mainly on optically intermediate and thick high clouds, contributing to an improved SWCRE and longwave cloud radiative effect (LWCRE). The impact of the MG2 tuning and the new convective trigger is primarily on the high latitudes and the SO. They have a relatively smaller impact on clouds than CLUBB tuning and ZM tuning do. This study offers additional insights into clouds simulated in E3SMv2 by utilizing multiple data sets and the Cloud Feedback Model Intercomparison Project (CFMIP) Observation Simulator Package (COSP) diagnostic tool as well as sensitivity tests. The improved understanding will benefit future E3SM developments.

The survival of relict Taxus wallichiana var. wallichiana (Yunnan yew) is threatened by overexploitation for its quality wood and medicinal properties, particularly for taxol extraction. Understanding the current status of its communities and populations is crucial for protecting existing natural forest resources. We established 53 vegetation plots in Yunnan, southwestern China, where T. wallichiana var. wallichiana is the primary dominant species. These plots were classified into four forest types. The forests were multi-stratified, with T. wallichiana var. wallichiana frequently dominating the subcanopy and shrub layer. Species diversity indices did not significantly differ among the four forest types. The age structure of T. wallichiana var. wallichiana exhibited a multi-modal pattern, with a maximum age of 1165 years. Growth was slow, with an average radial growth rate of 0.78 mm/year. Despite its strong sprouting ability, the species had a poor seedling/sapling bank and suffered from inadequate regeneration. Its seedlings/saplings are shade-intolerant. This study provides a scientific basis for effective conservation strategies, emphasizing the need for in situ regeneration to ensure the survival of T. wallichiana var. wallichiana and its contributions to biodiversity and ecosystem services.

Background:. The coronavirus disease 2019 (COVID-19) pandemic required an unprecedented transformation of medical education, shifting from traditional, in-person learning to distanced, online learning. This study aimed to review changes to medical education and describe the advantages and disadvantages of virtual medical education experienced by medical students during the pandemic. Methods:. An online survey study was conducted at two medical schools, University of Michigan Medical School in the United States and Koc University School of Medicine in Turkey. Medical students completed questionnaires regarding their educational experience before and during the pandemic. Survey instruments were designed to assess differences in the educational curriculum, study methods, clinical skills self-evaluations, perceptions of the quality of in-person and online learning, and overall satisfaction. Results:. A total of 184 medical students completed the survey. There was an increase in the use of online study tools since the pandemic. There was no statistically significant difference in self-reported assessments of clinical preparedness and overall clinical competence during surgical clerkship. The percentage of students interested in pursuing a career in surgery has nearly doubled from 34% to 63%. A majority of students (83%) believed that the time available for self-study and research increased during the pandemic. Fifty-two percent of students believed that online education is less efficacious than in-person education, but 86% of students still preferred a blended approach. Conclusions:. Medical schools have continued to update their curricula following the COVID-19 pandemic. This study illustrates the transformations in medical education to ensure that the most effective and suitable teaching is delivered.

Abstract Developing large‐scale monitoring strategies for marine habitats with appropriate temporal and spatial resolution is challenging logistically and economically. We compare potential DNA metabarcoding infaunal identification methods with traditional morphological assessment of marine intertidal biodiversity across two protected sites of contrasting anthropogenic impact on the south coast of England. We show that all the methodologies are effective at distinguishing the distinct communities present at each site. While morphological methods provide community abundance and biomass information, data derived from 18S metabarcoding of sediment scrapes showed the strongest discriminations between sites with contrasting anthropogenic pressures. This difference is due to the inclusion of more taxa from a wider spectrum of biodiversity in this dataset. However, 18S OTU abundance was a poor predictor of morphological taxa abundance and biomass. Examination of the presence or absence of taxa at the more or less impacted sites allows us to identify potential indicator taxa for future surveys, such as families in the phyla Ascidiacea, Cnidaria, Nematoda, and Platyhelminthes; taxa which are not traditionally incorporated into morphology‐based assessments due to the difficulty and expense of identification. DNA metabarcoding tools provide a more comprehensive snapshot of marine biodiversity compared to morphological surveys, and therefore an opportunity to assess the responses of more organisms to environmental stressors and develop novel metrics for habitat assessment. This could greatly benefit future monitoring programs and the assessment of management impacts in marine habitats.

This paper provides an overview of the United States (US) Department of Energy's (DOE's) Energy Exascale Earth System Model version 2 (E3SMv2) fully coupled regionally refined model (RRM) and documents the overall atmosphere, land, and river results from the Coupled Model Intercomparison Project 6 (CMIP6) DECK (Diagnosis, Evaluation, and Characterization of Klima) and historical simulations – a first-of-its-kind set of climate production simulations using RRM. The North American (NA) RRM (NARRM) is developed as the high-resolution configuration of E3SMv2 with the primary goal of more explicitly addressing DOE's mission needs regarding impacts to the US energy sector facing Earth system changes. The NARRM features finer horizontal resolution grids centered over NA, consisting of 25→100 km atmosphere and land, a 0.125∘ river-routing model, and 14→60 km ocean and sea ice. By design, the computational cost of NARRM is ∼3× of the uniform low-resolution (LR) model at 100 km but only ∼ 10 %–20 % of a globally uniform high-resolution model at 25 km. A novel hybrid time step strategy for the atmosphere is key for NARRM to achieve improved climate simulation fidelity within the high-resolution patch without sacrificing the overall global performance. The global climate, including climatology, time series, sensitivity, and feedback, is confirmed to be largely identical between NARRM and LR as quantified with typical climate metrics. Over the refined NA area, NARRM is generally superior to LR, including for precipitation and clouds over the contiguous US (CONUS), summertime marine stratocumulus clouds off the coast of California, liquid and ice phase clouds near the North Pole region, extratropical cyclones, and spatial variability in land hydrological processes. The improvements over land are related to the better-resolved topography in NARRM, whereas those over ocean are attributable to the improved air–sea interactions with finer grids for both atmosphere and ocean and sea ice. Some features appear insensitive to the resolution change analyzed here, for instance the diurnal propagation of organized mesoscale convective systems over CONUS and the warm-season land–atmosphere coupling at the southern Great Plains. In summary, our study presents a realistically efficient approach to leverage the fully coupled RRM framework for a standard Earth system model release and high-resolution climate production simulations.

The S1 stream surface inverse method for a single-stage axial transonic compressor was developed and studied under the guidance of the quasi-three-dimensional viscous design theory, based on computational fluid dynamics (CFD), in order to establish and solve governing equations. The re--adjusted load distribution of the S1 stream surface was imposed to obtain a special profile for the rotor, which is known as S-shaping. This changed the structure and weaken the strength of the shock due to the pre-compression effect. In order to match the inlet flow angle of the downstream stator, inverse modification for multiple S1 stream surfaces of the stator blade was conducted in the present study at the same time. As it is known, the S-shaped profile is commonly applied to supersonic flows. Therefore, NASA stage 37 was selected as the design case to verify whether the present inverse method is effective and reliable. Stage 37 was re-designed by stacking five S1 stream surfaces. The profiles of these surfaces were renewed through the inverse design. The results revealed that, compared to the prototype, the aerodynamic performance of the re-designed one was apparently promoted within the stable working range. The adiabatic efficiency at the design point increased by approximately 0.4%, and the pressure ratio improved by approximately 4.5%. In addition, analysis result for the characteristic line revealed that the performance of the redesigned blades at off-design points significantly improved.

As the adverse impacts of hydrological extremes increase in many regions of the world, a better understanding of the drivers of changes in risk and impacts is essential for effective flood and drought risk management and climate adaptation. However, there is currently a lack of comprehensive, empirical data about the processes, interactions, and feedbacks in complex human–water systems leading to flood and drought impacts. Here we present a benchmark dataset containing socio-hydrological data of paired events, i.e. two floods or two droughts that occurred in the same area. The 45 paired events occurred in 42 different study areas and cover a wide range of socio-economic and hydro-climatic conditions. The dataset is unique in covering both floods and droughts, in the number of cases assessed and in the quantity of socio-hydrological data. The benchmark dataset comprises (1) detailed review-style reports about the events and key processes between the two events of a pair; (2) the key data table containing variables that assess the indicators which characterize management shortcomings, hazard, exposure, vulnerability, and impacts of all events; and (3) a table of the indicators of change that indicate the differences between the first and second event of a pair. The advantages of the dataset are that it enables comparative analyses across all the paired events based on the indicators of change and allows for detailed context- and location-specific assessments based on the extensive data and reports of the individual study areas. The dataset can be used by the scientific community for exploratory data analyses, e.g. focused on causal links between risk management; changes in hazard, exposure and vulnerability; and flood or drought impacts. The data can also be used for the development, calibration, and validation of socio-hydrological models. The dataset is available to the public through the GFZ Data Services (Kreibich et al., 2023, https://doi.org/10.5880/GFZ.4.4.2023.001).

Surface structure is used to interfere with the turbulent boundary layer in the groove drag reduction, which is important to the endurance and stability of high-speed and ultrahigh-speed aircraft. The size of the groove structure directly affects the flow in the turbulent boundary layer and changes the drag reduction effect. The drag reduction characteristics of bionic triangular (V-groove) riblets were studied through Particle Image Velocimetry (PIV) experiment and Finite Volume Method (FVM) simulation. Triangular riblets with adjacent height ratios (AHR) of 1.00, 1.74, and 3.02 were considered in this research, and the influence of these groove structures on the flow characteristics of turbulence near the wall is compared with those of the smooth plates. The distribution of time-averaged velocity, turbulence intensity, and coherent structures of turbulent boundary layer on the riblet surface is analyzed to document the effects of the geometric parameters of various groove structures on drag reduction rates. Results showed that the best drag reduction is obtained using the V-groove riblets with adjacent height ratio of 1:1 under the low free-stream velocity. The results can be used as a reference for further optimization of drag reduction structures with surface grooves.

In this paper we develop a novel framework aimed to significantly reduce biases related to marine stratocumulus clouds in general circulation models (GCMs) while circumventing excessive computational cost requirements. Our strategy is to increase the horizontal resolution using a regionally refined mesh (RRM) over our region of interest in addition to using the Framework for Improvement by Vertical Enhancement (FIVE) to increase the vertical resolution only for specific physical processes that are important for stratocumulus. We apply the RRM off the coast of Peru in the southeastern Pacific, a region that climatologically contains the most marine stratocumulus in the subtropics. We find that our new modeling framework is able to replicate the results of our high-resolution benchmark simulation with much fidelity, while reducing the computational cost by several orders of magnitude. In addition, this framework is able to greatly reduce the long-standing biases associated with marine stratocumulus in GCMs when compared to the standard-resolution control simulation.

The E3SM Diagnostics Package (E3SM Diags) is a modern, Python-based Earth system model (ESM) evaluation tool (with Python module name e3sm_diags), developed to support the Department of Energy (DOE) Energy Exascale Earth System Model (E3SM). E3SM Diags provides a wide suite of tools for evaluating native E3SM output, as well as ESM data on regular latitude–longitude grids, including output from Coupled Model Intercomparison Project (CMIP) class models. E3SM Diags is modeled after the National Center for Atmospheric Research (NCAR) Atmosphere Model Working Group (AMWG, 2022) diagnostics package. In its version 1 release, E3SM Diags included a set of core essential diagnostics to evaluate the mean physical climate from model simulations. As of version 2.7, more process-oriented and phenomenon-based evaluation diagnostics have been implemented, such as analysis of the quasi-biennial oscillation (QBO), the El Niño–Southern Oscillation (ENSO), streamflow, the diurnal cycle of precipitation, tropical cyclones, ozone and aerosol properties. An in situ dataset from DOE's Atmospheric Radiation Measurement (ARM) program has been integrated into the package for evaluating the representation of simulated cloud and precipitation processes. This tool is designed with enough flexibility to allow for the addition of new observational datasets and new diagnostic algorithms. Additional features include customizable figures; streamlined installation, configuration and execution; and multiprocessing for fast computation. The package uses an up-to-date observational data repository maintained by its developers, where recent datasets are added to the repository as they become available. Finally, several applications for the E3SM Diags module were introduced to fit a diverse set of use cases from the scientific community.

Engineered chromodomains are promising probes to analyze histone methylation. Here the authors designed high-affinity chromodomains for enhanced genome-wide binding analysis, live-cell imaging and ultra-potent CRISPR-based gene repression.

The naturalized streamflow, i.e., streamflow without water management effects, in the Yellow River basin (YRB) has been significantly decreased at a rate of -3.71×108 m3 yr−1 during 1982–2018, although annual precipitation experienced an insignificantly positive trend. Explicit detection and attribution of naturalized streamflow are critical to manage limited water resources for the sustainable development of ecosystems and socio-economic systems. The effects from temporally explicit changes of climate variables and underlying surfaces on the streamflow trend were assessed using the variable infiltration capacity (VIC) model prescribed with continuously dynamic leaf area index (LAI) and land cover. The results show a sharp increase of the LAI trend and land use change as a conversion of cropland into forest grass in the basin. The decrease in naturalized streamflow can primarily be attributed to the vegetation changes including an interannual LAI increase and intra-annual LAI temporal pattern change, which account for the streamflow reduction of 1.99×108 and 0.45×108 m3 yr−1, respectively. The impacts of the LAI change are largest at the subregion of Longmen–Huayuankou where the LAI increasing trend is high and land use change is substantial. Attribution based on simulations with multiyear average LAI changes obviously underestimates the impacts of the interannual LAI change and intra-annual LAI temporal change on the natural streamflow trend. Overall, the effect of climate variation on streamflow is slight because the positive effect from precipitation and wind speed changes was offset by the negative effect from increasing temperature. Although climate variation is decisive for streamflow change, this study suggests that change in underlying surfaces has imposed a substantial trend on naturalized streamflow. This study improves the understanding of the spatiotemporal patterns and the underlying mechanisms of natural streamflow reduction across the YRB between 1982 and 2018.

The effective radiative forcing of anthropogenic aerosols (ERFaer) is an important measure of the anthropogenic aerosol effects simulated by a global climate model. Here we analyze ERFaer simulated by the E3SM version 1 (E3SMv1) atmospheric model using both century-long free-running atmosphere–land simulations and short nudged simulations. We relate the simulated ERFaer to characteristics of the aerosol composition and optical properties, and we evaluate the relationships between key aerosol and cloud properties. In terms of historical changes from the year 1870 to 2014, our results show that the global mean anthropogenic aerosol burden and optical depth increase during the simulation period as expected, but the regional averages show large differences in the temporal evolution. The largest regional differences are found in the emission-induced evolution of the burden and optical depth of the sulfate aerosol: a strong decreasing trend is seen in the Northern Hemisphere high-latitude region after around 1970, while a continued increase is simulated in the tropics. The relationships between key aerosol and cloud properties (relative changes between pre-industrial and present-day conditions) also show evident changes after 1970, diverging from the linear relationships exhibited for the period of 1870–1969. In addition to the regional differences in the simulated relationships, a reduced sensitivity in cloud droplet number and other cloud properties to aerosol perturbations is seen when the aerosol perturbation is large. Consequently, the global annual mean ERFaer magnitude does not increase after around 1970. The ERFaer in E3SMv1 is relatively large compared to the recently published multi-model estimates; the primary reason is the large indirect aerosol effect (i.e., through aerosol–cloud interactions). Compared to other models, E3SMv1 features large relative changes in the cloud droplet effective radius in response to aerosol perturbations. Large sensitivity is also seen in the liquid cloud optical depth, which is determined by changes in both the effective radius and liquid water path. Aerosol-induced changes in liquid and ice cloud properties in E3SMv1 are found to have a strong correlation, as the evolution of anthropogenic sulfate aerosols affects both the liquid cloud formation and the homogeneous ice nucleation in cirrus clouds (that causes a large effect on longwave ERFaer). As suggested by a previous study, the large ERFaer appears to be one of the reasons why the model cannot reproduce the observed global mean temperature evolution in the second half of the 20th century. Sensitivity simulations are performed to understand which parameterization and/or parameter changes have a large impact on the simulated ERFaer. The ERFaer estimates in E3SMv1 for the shortwave and longwave components are sensitive to the parameterization changes in both liquid and ice cloud processes. When the parameterization of ice cloud processes is modified, the top-of-model forcing changes in the shortwave and longwave components largely offset each other, so the net effect is negligible. This suggests that, to reduce the magnitude of the net ERFaer, it would be more effective to reduce the anthropogenic aerosol effect through liquid or mixed-phase clouds.

In order to improve the performance of single-stage turbines, blade profiling optimization was conducted for the guide vane and rotor under design condition. Support vector regression (SVR) and non-dominated sorting genetic algorithm-II (NSGA-II) were used to execute the optimization, with the objective of maximizing the efficiency and total pressure ratio of single-stage turbines. The gas turbine chosen for the initial study was the KJ66, which is one of the most robust and primitive small gas turbine designs available. The influence mechanism of the stator and rotor profiling on flow field and performance was discussed. The results revealed that compared with the prototype, the adiabatic efficiency increased by 5.95% and the total pressure ratio increased by 0.9%. Furthermore, the matching of flows between the stator blade and rotor blade obviously improved. The optimized guide vane suppressed the flow separation by increasing the leading edge and improving the distribution of the inlet angle of attack. The load distribution of rotors with a 50% spanwise position changed from the original "C" loaded to post-loaded. The leading load obviously decreased, and the angle of attack was smaller than that of the prototype, which effectively weakened the flow separation at the leading edge of the rotor. Compared with the original rotor, the higher lean angle and pressure ratio of the turbine stage also improved. However, the leakage loss near the shroud of the rotor increased, which led to decreased efficiency.

This paper documents the experimental setup and general features of the coupled historical and future climate simulations with the first version of the US Department of Energy (DOE) Energy Exascale Earth System Model (E3SMv1.0). The future projected climate characteristics of E3SMv1.0 at the highest emission scenario (SSP5-8.5) designed in the Scenario Model Intercomparison Project (ScenarioMIP) and the SSP5-8.5 greenhouse gas (GHG) only forcing experiment are analyzed with a focus on regional responses of atmosphere, ocean, sea ice, and land. Due to its high equilibrium climate sensitivity (ECS of 5.3 K), E3SMv1.0 is one of the Coupled Model Intercomparison Project phase 6 (CMIP6) models with the largest surface warming by the end of the 21st century under the high-emission SSP5-8.5 scenario. The global mean precipitation change is highly correlated with the global temperature change, while the spatial pattern of the change in runoff is consistent with the precipitation changes. The oceanic mixed layer generally shoals throughout the global ocean. The annual mean Atlantic meridional overturning circulation (AMOC) is overly weak with a slower change from ∼ 11 to ∼ 6 Sv (Sverdrup) relative to other CMIP6 models. The sea ice, especially in the Northern Hemisphere, decreases rapidly with large seasonal variability. We detect a significant polar amplification in E3SMv1.0 from the atmosphere, ocean, and sea ice. Comparing the SSP5-8.5 all-forcing experiment with the GHG-only experiment, we find that the unmasking of the aerosol effects due to the decline of the aerosol loading in the future projection period causes transient accelerated warming in the all-forcing experiment in the first half of the 21st century. While the oceanic climate response is mainly controlled by the GHG forcing, the land runoff response is impacted primarily by forcings other than GHG over certain regions, e.g., southern North America, southern Africa, central Africa, and eastern Asia. However, the importance of the GHG forcing on the land runoff changes grows in the future climate projection period compared to the historical period.

The inherently dry summer climate of the Iberian Peninsula (IP) is undergoing drought exacerbated by more intense warming and reduced precipitation. Although many studies have studied changes in summer climate factors, it is still unclear how the changes in moisture contribution from the sources lead to the decrease in summer precipitation. This study investigates the differences in the IP precipitationshed from 1980 to 1997 and 1998 to 2019 using the Water Accounting Model-2layers with ERA5 data, and assesses the role of local moisture recycling and external moisture in reducing summer precipitation. Our findings indicate that the moisture contributions from the local IP, and from the west and the east of the precipitationshed, contributed 1.7, 3.6 and 1.1 mm per month less precipitation after 1997 than before 1997, accounting for 26 %, 57 % and 17 % of the main source supply reduction, respectively. The significant downward trend of the IP local moisture recycling closely links to the disappearance of the wet years after 1997 as well as the decrease in local contribution in the dry years. Moreover, the feedback between the weakened local moisture recycling and the drier land surface can exacerbate the local moisture scarcity and summer drought.

Realistic simulation of the Earth's mean-state climate remains a major challenge, and yet it is crucial for predicting the climate system in transition. Deficiencies in models' process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections and incorrect attribution of the physical mechanisms governing past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and long-standing biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when it is run as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback, as well as the fidelity of the simulated base climate state, are important for constraining the climate in the past and future.

Tibet, which is characterized by collisional orogens, has undergone the process of delamination or convective removal. The lower crust and mantle lithosphere appear to have been removed through delamination during orogenic development. Numerical and analog experiments demonstrate that the metamorphic eclogitized oceanic subduction slab or lower crust may promote gravitational instability due to increased density. The eclogitized oceanic subduction slab or crustal root is believed to be denser than the underlying mantle and tends to sink. However, the density of eclogite under high-pressure and high-temperature conditions as well as density differences from the surrounding mantle are not preciously constrained. Here, we offer new insights into the derivation of eclogite density with a single experiment to constrain delamination in Tibet. Using in situ synchrotron X-ray diffraction combined with a diamond anvil cell, experiments focused on minerals (garnet, omphacite, and epidote) of eclogite are conducted under simultaneous high-pressure and high-temperature conditions, which avoids systematic errors. Fitting the pressure–temperature–volume data with the third-order Birch–Murnaghan equation of state, the thermal equation of state (EoS) parameters, including the bulk modulus (KT0), its pressure derivative (KT0′), and the thermal expansion coefficient (α0), are derived. The densities of rock-forming minerals and eclogite are modeled along with the geotherms of two types of delamination. The delamination processes of subduction slab break-off and the removal of the eclogitized lower crust in Tibet are discussed. The Tibetan eclogite, which contains 40 vol %–60 vol % garnet and 44 %–70 % eclogitization, can promote the delamination of slab break-off in Tibet. Our results indicate that eclogite is a major controlling factor in the initiation of delamination. A high abundance of garnet, a high Fe content, and a high degree of eclogitization are more conducive to instigating the delamination.

Trehalose is a nonreducing disaccharide that is involved in the regulation of plant responses to a variety of environmental stresses. Trehalose 6-phosphate synthase (TPS) and trehalose 6-phosphate phosphatase (TPP) are two key enzymes in trehalose synthesis and they are widely distributed in higher plants. At present, TPS family genes have been systematically identified and analyzed in many plant species, but the TPP family genes have been rarely studied. In this study, ten TPS and six TPP genes in cannabis (Cannabis sativa L.) were identified at the genomic level. The phylogenetic tree of TPS and TPP family members in cannabis, Arabidopsis, and rice was constructed, and all the genes were divided into three subgroups: Class I, Class II, and Class III. The number of exons and motif types among Class I members was exactly the same, as were Class II members, but the gene structure and motif types of Class III members were slightly different. There were four pairs of CsTPSs and CsTPPs that had gene duplication, indicating that gene duplication events played an important role in the amplification of TPS and TPP families in cannabis. The results of expression analysis under abiotic stresses showed that 68.75 % of CsTPS and CsTPP genes were significantly induced by at least one abiotic stress. Among these genes, the expression of CsTPS1, CsTPS9, and CsTPPA was highest under at least one abiotic stress. These three genes may play a key role in abiotic stress responses. Most of the CsTPS and CsTPP genes that are closely located in the evolutionary tree have the same or similar functions. To our knowledge, this is the first paper that systematically reports the TPS and TPP gene families in cannabis.

Measurement(s) climate stability index Technology Type(s) general circulation models of climate change Factor Type(s) temporal interval Sample Characteristic - Environment climate system • climate change Sample Characteristic - Location global Machine-accessible metadata file describing the reported data: https://doi.org/10.6084/m9.figshare.19050569

Protein turnover, a highly regulated process governed by the ubiquitin–proteasome system (UPS), is essential for maintaining cellular homeostasis. Dysregulation of the UPS has been implicated in various diseases, including viral infections and cancer, making the proteins in the UPS attractive targets for therapeutic intervention. However, the functional and structural redundancies of UPS enzymes present challenges in identifying precise drug targets and achieving target selectivity. Consequently, only 26S proteasome inhibitors have successfully advanced to clinical use thus far. To overcome these obstacles, engineered peptides and proteins, particularly engineered ubiquitin, have emerged as promising alternatives. In this review, we examine the impact of engineered ubiquitin on UPS and non-UPS proteins, as well as on viral enzymes. Furthermore, we explore their potential to guide the development of small molecules targeting novel surfaces, thereby expanding the range of druggable targets.

In order to improve the efficiency of identifying parameters using the maximum likelihood method and to avoid the sensitivity of initial values, a proposed method that combines the micro-genetic algorithm with the advance and retreat method is presented in order to identify the eccentricity of the spindle-tool system with random input and output parameters, which obey a certain probability distribution. Eccentricity without prior information is determined through an iterative procedure. The initial value starts from zero, and the interval is determined by the advance and retreat method. Then, the optimal value is searched in the corresponding interval, utilizing the micro-genetic algorithm. The initial value and interval at each of iterations are changed to ensure a fast and stable convergence. Eventually, a numerical example with three kinds of random deviations verifies the feasibility and validity of the proposed method.

Subseasonal-to-seasonal (S2S) prediction, especially the prediction of extreme hydroclimate events such as droughts and floods, is not only scientifically challenging, but also has substantial societal impacts. Motivated by preliminary studies, the Global Energy and Water Exchanges (GEWEX)/Global Atmospheric System Study (GASS) has launched a new initiative called “Impact of Initialized Land Surface Temperature and Snowpack on Subseasonal to Seasonal Prediction” (LS4P) as the first international grass-roots effort to introduce spring land surface temperature (LST)/subsurface temperature (SUBT) anomalies over high mountain areas as a crucial factor that can lead to significant improvement in precipitation prediction through the remote effects of land–atmosphere interactions. LS4P focuses on process understanding and predictability, and hence it is different from, and complements, other international projects that focus on the operational S2S prediction. More than 40 groups worldwide have participated in this effort, including 21 Earth system models, 9 regional climate models, and 7 data groups. This paper provides an overview of the history and objectives of LS4P, provides the first-phase experimental protocol (LS4P-I) which focuses on the remote effect of the Tibetan Plateau, discusses the LST/SUBT initialization, and presents the preliminary results. Multi-model ensemble experiments and analyses of observational data have revealed that the hydroclimatic effect of the spring LST on the Tibetan Plateau is not limited to the Yangtze River basin but may have a significant large-scale impact on summer precipitation beyond East Asia and its S2S prediction. Preliminary studies and analysis have also shown that LS4P models are unable to preserve the initialized LST anomalies in producing the observed anomalies largely for two main reasons: (i) inadequacies in the land models arising from total soil depths which are too shallow and the use of simplified parameterizations, which both tend to limit the soil memory; (ii) reanalysis data, which are used for initial conditions, have large discrepancies from the observed mean state and anomalies of LST over the Tibetan Plateau. Innovative approaches have been developed to largely overcome these problems.

The analysis of environmental DNA (eDNA) by high throughput sequencing (HTS) is proving to be a promising tool for freshwater fish biodiversity assessment in Europe within the Water Framework Directive (WFD, 2000/60/EC), especially for large rivers and lakes where current fish monitoring techniques have known shortcomings. These new biomonitoring methods based on eDNA show several advantages compared to classical morphological methods. The sampling procedures are easier and cheaper and eDNA metabarcoding is non-invasive and very sensitive, allowing for the detection of traces of DNA. However, eDNA metabarcoding methods need careful standardization to make the results of different surveys comparable. The aim of the EU project Eco-AlpsWater is to test and validate molecular biodiversity monitoring tools for aquatic ecosystems (i.e., eDNA metabarcoding) to improve the traditional WFD monitoring approaches in Alpine waterbodies. To this end, an inter-calibration test was performed using fish mock community samples containing either tissue-extracted DNA, eDNA collected from aquaculture tanks and eDNA samples collected from Lake Bourget (France). Samples were analysed using a DNA metabarcoding approach, relying on the amplification and HTS of a 12S rDNA marker, in two separate laboratories, to evaluate if different laboratory and bioinformatic protocols can provide a reliable and comparable description of the fish communities in both mock and natural samples. Our results highlight good replicability of the molecular laboratory protocols for HTS and good amplification success of selected primers, providing essential information concerning the taxonomic resolution of the 12S mitochondrial marker in describing the Alpine fish communities. Interestingly, different concentrations of species DNA in the mock samples were well represented by the relative DNA reads abundance. These tests confirm the reproducibility of eDNA metabarcoding analyses for the biomonitoring of freshwater fish inhabiting Alpine and peri-Alpine lakes and rivers.

Abstract Submarine mega‐slides involving hundreds of cubic kilometers of slope material pose a major threat due to their potential to destroy offshore infrastructure and trigger devastating tsunamis. The Sahara Slide Complex affected about 50,000 km2 of the northwestern (NW) African continental margin. Previous studies focused either on its distal depositional zone or the uppermost headwall area, but failed in reconstructing the succession of individual slide events within the entire headwall area. New hydroacoustic data reveal a complex slide morphology including three main acoustic facies, large scale slide blocks, linear troughs, multiple glide planes and three major headwall scarps (the upper, southern and lower headwall). The evacuated slide scar hosts chaotic slide deposits that cover stratified sediments in the upper and southern headwall area, but are vertically stacked onto older slide deposits in the lower headwall area. Based on these observations, and dating of recently collected sediment samples, we reconstructed the history of slope failures that led to the formation of the structurally and morphologically complex headwall area of the Sahara Slide. Slope instability initiated when the lower headwall failed at ∼60 kyr, followed by the failure of the northeastern upper headwall at ∼14 kyr. Around 6 kyr, a major slide within the upper headwall area took place, followed by a series of smaller events within the southern and most‐proximal upper headwall area. The youngest of these slides occurred around 2 kyr. This scenario suggests a long‐lasting history of successive slope failures for the Sahara Slide Complex along the NW African continental slope.

Abstract Quantification of dust aerosols in Earth System Models (ESMs) has important implications for water cycle and biogeochemistry studies. This study examines the global life cycle and direct radiative effects (DREs) of dust in the U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1), and the impact of increasing model resolution both horizontally and vertically. The default 1° E3SMv1 captures the spatial and temporal variability in the observed dust aerosol optical depth (DAOD) reasonably well, but overpredicts dust absorption in the shortwave (SW). Simulations underestimate the dust vertical and long‐range transport, compared with the satellite dust extinction profiles. After updating dust refractive indices and correcting for a bias in partitioning size‐segregated emissions, both SW cooling and longwave (LW) warming of dust simulated by E3SMv1 are increased and agree better with other recent studies. The estimated net dust DRE of −0.42 Wm−2 represents a stronger cooling effect than the observationally based estimate −0.2 Wm−2 (−0.48 to +0.2), due to a smaller LW warming. Constrained by a global mean DAOD, model sensitivity studies of increasing horizontal and vertical resolution show strong influences on the simulated global dust burden and lifetime primarily through the change of dust dry deposition rate; there are also remarkable differences in simulated spatial distributions of DAOD, DRE, and deposition fluxes. Thus, constraining the global DAOD is insufficient for accurate representation of dust climate effects, especially in transitioning to higher‐ or variable‐resolution ESMs. Better observational constraints of dust vertical profiles, dry deposition, size, and LW properties are needed.

The Scenario Model Intercomparison Project (ScenarioMIP) defines and coordinates the main set of future climate projections, based on concentration-driven simulations, within the Coupled Model Intercomparison Project phase 6 (CMIP6). This paper presents a range of its outcomes by synthesizing results from the participating global coupled Earth system models. We limit our scope to the analysis of strictly geophysical outcomes: mainly global averages and spatial patterns of change for surface air temperature and precipitation. We also compare CMIP6 projections to CMIP5 results, especially for those scenarios that were designed to provide continuity across the CMIP phases, at the same time highlighting important differences in forcing composition, as well as in results. The range of future temperature and precipitation changes by the end of the century (2081–2100) encompassing the Tier 1 experiments based on the Shared Socioeconomic Pathway (SSP) scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0 and SSP5-8.5) and SSP1-1.9 spans a larger range of outcomes compared to CMIP5, due to higher warming (by close to 1.5 ∘C) reached at the upper end of the 5 %–95 % envelope of the highest scenario (SSP5-8.5). This is due to both the wider range of radiative forcing that the new scenarios cover and the higher climate sensitivities in some of the new models compared to their CMIP5 predecessors. Spatial patterns of change for temperature and precipitation averaged over models and scenarios have familiar features, and an analysis of their variations confirms model structural differences to be the dominant source of uncertainty. Models also differ with respect to the size and evolution of internal variability as measured by individual models' initial condition ensemble spreads, according to a set of initial condition ensemble simulations available under SSP3-7.0. These experiments suggest a tendency for internal variability to decrease along the course of the century in this scenario, a result that will benefit from further analysis over a larger set of models. Benefits of mitigation, all else being equal in terms of societal drivers, appear clearly when comparing scenarios developed under the same SSP but to which different degrees of mitigation have been applied. It is also found that a mild overshoot in temperature of a few decades around mid-century, as represented in SSP5-3.4OS, does not affect the end outcome of temperature and precipitation changes by 2100, which return to the same levels as those reached by the gradually increasing SSP4-3.4 (not erasing the possibility, however, that other aspects of the system may not be as easily reversible). Central estimates of the time at which the ensemble means of the different scenarios reach a given warming level might be biased by the inclusion of models that have shown faster warming in the historical period than the observed. Those estimates show all scenarios reaching 1.5 ∘C of warming compared to the 1850–1900 baseline in the second half of the current decade, with the time span between slow and fast warming covering between 20 and 27 years from present. The warming level of 2 ∘C of warming is reached as early as 2039 by the ensemble mean under SSP5-8.5 but as late as the mid-2060s under SSP1-2.6. The highest warming level considered (5 ∘C) is reached by the ensemble mean only under SSP5-8.5 and not until the mid-2090s.

Stratospheric ozone affects climate directly as the predominant heat source in the stratosphere and indirectly through chemical reactions controlling other greenhouse gases. The U.S. Department of Energy's Energy Exascale Earth System Model version 1 (E3SMv1) implemented a new ozone chemistry module that improves the simulation of the sharp tropopause gradients, replacing a version based partly on long-term average climatologies that poorly represented heating rates in the lowermost stratosphere. The new O3v2 module extends seamlessly into the troposphere and preserves the naturally sharp cross-tropopause gradient, with 20 %–40 % less ozone in this region. Additionally, O3v2 enables the diagnosis of stratosphere–troposphere exchange flux of ozone, a key budget term lacking in E3SMv1. Here, we evaluate key features in ozone abundance and other closely related quantities in atmosphere-only E3SMv1 simulations driven by observed sea surface temperatures (SSTs, years 1990–2014), comparing them with satellite observations of ozone and also with the University of California, Irvine chemistry transport model (UCI CTM) using the same stratospheric chemistry scheme but driven by European Centre forecast fields for the same period. In terms of stratospheric column ozone, O3v2 shows reduced mean bias and improved northern midlatitude variability, but it is not quite as good as the UCI CTM. As expected, SST-forced E3SMv1 simulations cannot synchronize with observed quasi-biennial oscillations (QBOs), but they do show the typical QBO pattern seen in column ozone. This new O3v2 E3SMv1 model mostly retains the same climate state and climate sensitivity as the previous version, and we recommend its use for other climate models that still use ozone climatologies.

A stochastic deep convection parameterization is implemented into the US Department of Energy (DOE) Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1.0 (EAMv1). This study evaluates its performance in simulating precipitation. Compared to the default model, the probability distribution function (PDF) of rainfall intensity in the new simulation is greatly improved. The well-known problem of “too much light rain and too little heavy rain” is alleviated, especially over the tropics. As a result, the contribution from different rain rates to the total precipitation amount is shifted toward heavier rain. The less frequent occurrence of convection contributes to suppressed light rain, while more intense large-scale and convective precipitation contributes to enhanced heavy total rain. The synoptic and intraseasonal variabilities of precipitation are enhanced as well to be closer to observations. The sensitivity of the rainfall intensity PDF to the model vertical resolution is examined. The relationship between precipitation and dilute convective available potential energy in the stochastic simulation agrees better with that in the Atmospheric Radiation Measurement (ARM) observations compared with the standard model simulation. The annual mean precipitation is largely unchanged with the use of the stochastic scheme except over the tropical western Pacific, where a moderate increase in precipitation represents a slight improvement. The responses of precipitation and its extremes to climate warming are similar with or without the stochastic deep convection scheme.

Cancer immunotherapy, which enhances the body's immune response and inhibits immune escape to eliminate or combat tumors, holds great promise for cancer treatment. However, it still faces many challenges, such as low response rate and immune-related adverse events, which limit its clinical applications. Ultrasound, as a safe, effective and non-invasive technique, has been widely used in diagnostics and therapeutics of tumors. In recent years, many studies have confirmed that ultrasound can also mediate tumor immunotherapy. High intensity focused ultrasound can destroy tumor tissue and lead to the in situ production of tumor fragments at the same time, which can be used as antigens of the immune system to form in situ vaccines to induce a systemic immune response. The ultrasound-targeted microbubble destruction technique with microbubbles as the carrier mediates the targeted delivery of immunomodulatory adjuvants to tumor cells or immune cells. It can not only enhance the host immune response and improve immunotherapy efficacy, but also reduce the dosage of immunoregulatory adjuvants and side effects. Sonodynamic therapy can induce immunogenic cell death to elicit an immune response and the combination with immunotherapy can achieve higher therapeutic efficacy. In this review, research advances in ultrasound-mediated tumor immunotherapy are summarized, and the potential and challenges of its application are discussed.

Large volcanic eruptions reaching the stratosphere have caused marked perturbations to the global climate including cooling at the Earth's surface, changes in large-scale circulation and precipitation patterns and marked temporary reductions in global ocean heat content. Many studies have investigated these effects using climate models; however, uncertainties remain in the modelled response to these eruptions. This is due in part to the diversity of forcing datasets that are used to prescribe the distribution of stratospheric aerosols resulting from these volcanic eruptions, as well as uncertainties in optical property derivations from these datasets. To improve this situation for the sixth phase of the Coupled Model Intercomparison Project (CMIP6), a two-step process was undertaken. First, a combined stratospheric aerosol dataset, the Global Space-based Stratospheric Aerosol Climatology (GloSSAC; 1979–2016), was constructed. Next, GloSSAC, along with information from ice cores and Sun photometers, was used to generate aerosol distributions, characteristics and optical properties to construct a more consistent stratospheric aerosol forcing dataset for models participating in CMIP6. This “version 3” of the stratospheric aerosol forcing has been endorsed for use in all contributing CMIP6 simulations. Recent updates to the underlying GloSSAC from version 1 to version 1.1 affected the 1991–1994 period and necessitated an update to the stratospheric aerosol forcing from version 3 to version 4. As version 3 remains the official CMIP6 input, quantification of the impact on radiative forcing and climate is both relevant and timely for interpreting results from experiments such as the CMIP6 historical simulations. This study uses two models, the Canadian Earth System Model version 5 (CanESM5) and the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1), to estimate the difference in instantaneous radiative forcing in simulated post-Pinatubo climate response when using version 4 instead of version 3. Differences in temperature, precipitation and radiative forcings are generally found to be small compared to internal variability. An exception to this is differences in monthly temperature anomalies near 24 km altitude in the tropics, which can be as large as 3 ∘C following the eruption of Mt. Pinatubo.

Although many studies have explored the effect of teleconnection patterns on floods, few investigations have focused on the assessment of expected flood damage under such large-scale atmospheric signals. This study aims to determine the effect of the most emblematic teleconnection, El Niño, on the expected damage due to floods with short return periods in the Kan River basin, Iran. To determine the flood damage costs, the median of annual precipitation changes (ΔP) during El Niño conditions was used, although ΔP cannot necessarily be transferred to extreme values. Then the flooded area was determined under the increased rainfall due to El Niño for 5-, 10-, and 50-year return periods. The results showed that El Niño has increased the annual precipitation by 12.2 %. Flood damage assessment using damage–depth curves showed that the relative increase in expected damage during El Niño conditions is much higher for short return period floods than that for long return period floods. In general, a 12.2 % increase in the annual precipitation would increase the damage by 1671 % and 176 %, respectively, for the return periods of 5 and 10 years. However, in the case of a 50-year flood, this increased percentile decreased to 52 %. These results indicate the importance of small flood events in flood management planning during El Niño.

Data assimilation integrates information from observational measurements with numerical models. When used with coupled models of Earth system compartments, e.g., the atmosphere and the ocean, consistent joint states can be estimated. A common approach for data assimilation is ensemble-based methods which utilize an ensemble of state realizations to estimate the state and its uncertainty. These methods are far more costly to compute than a single coupled model because of the required integration of the ensemble. However, with uncoupled models, the ensemble methods also have been shown to exhibit a particularly good scaling behavior. This study discusses an approach to augment a coupled model with data assimilation functionality provided by the Parallel Data Assimilation Framework (PDAF). Using only minimal changes in the codes of the different compartment models, a particularly efficient data assimilation system is generated that utilizes parallelization and in-memory data transfers between the models and the data assimilation functions and hence avoids most of the file reading and writing, as well as model restarts during the data assimilation process. This study explains the required modifications to the programs with the example of the coupled atmosphere–sea-ice–ocean model AWI-CM (AWI Climate Model). Using the case of the assimilation of oceanic observations shows that the data assimilation leads only to small overheads in computing time of about 15 % compared to the model without data assimilation and a very good parallel scalability. The model-agnostic structure of the assimilation software ensures a separation of concerns in which the development of data assimilation methods can be separated from the model application.

Aims: Pinus yunnanesis is commercially, culturally and economically important, but there is a lack of ecological data on its role in stand dynamics. Our aims are to clarify the structure, composition, regeneration and growth trends of primary mature P. yunnanensis forests. Study area: The Tianchi National Nature Reserve in the Xuepan Mountains, Yunlong County, northwestern Yunnan, China. Methods: We investigated forests containing P. yunnanensis, measured tree ages and analyzed the data. Results: Six forest types were identified: (1) coniferous forest: Pinus yunnanensis (Type 1); (2) mixed coniferous and evergreen broad-leaved forest: P. yunnanensis-Lithocarpus variolosus (Type 2); (3) mixed coniferous and deciduous broad-leaved forest: P. yunnanensis-Quercus griffithii (Type 3); (4) mixed evergreen broad-leaved and coniferous forest: Castanopsis orthacantha-P. yunnanensis-Schima argentea (Type 4); (5) mixed coniferous, evergreen and deciduous broad-leaved forest: Pinus yunnanensis-Schima argentea-Quercus griffithii (Type 5); (6) mixed coniferous and evergreen broad-leaved forest: Pinus armandii-Quercus rehderiana-Pinus yunnanensis (Type 6). The size- and age-structure and regeneration patterns of P. yunnanensis were highly variable within these six forest types. P. yunnanensis regeneration was well balanced in forest Type 1 as compared to the other five types. All six forest types were identified as rare and old-growth with P. yunnaensis trees reaching ages of more than 105 years (a maximum age of 165 years with a diameter 116 cm at breast height) except for the Type 4 forest (a 90-year-old stand). Growth rates of P. yunnanensis, based upon ring width measurements, were high for the first 10 years, then declined after the 10th year. In contrast, basal area increment (BAI) increased for the first 25 years, plateaued, and only declined as trees became older. Trees in the older age classes grew more quickly than younger trees at the same age, a consequence of either site quality or competitive differences. The BAI of P. yunnanensis in all age classes in the Tianchi National Nature Reserve was much higher than those of the secondary and degraded natural P. yunnanensis forests of other areas. Conclusions: The P. yunnanensis forests of the Tianchi area appear to be some of the last remnants of primeval and old-growth forests of this species. These forests are structurally diverse and contain a rich diversity of overstory, mid-story, and understory species. Taxonomic reference: Editorial Committee of Flora Republicae Popularis Sinicae (1959–2004) for vascular plants. Abbreviations: BA = basal area; BAI = basal area increment; DBH = diameter at breast height; H = height; RBA = relative basal area.

Abstract Background Schima genus of Theaceae is confined to subtropics and tropics of South, East and Southeast Asia. Thirteen species of Schima are distributed in subtropical China. Many of them appear as dominant canopy species in the subtropical forests. To date, Schima species richness distribution patterns of China have remained unknown. Meanwhile, there has been a longtime debate as to whether forests dominated by Schima species are early or late successional forests. We aim to clarify Schima species richness patterns and these species’ roles in the forest succession and regeneration dynamics of the subtropical ecosystem in Yunnan Province, China. Method We mapped Schima species richness distribution patterns in China. Based on 71 vegetation plots, we analyzed forest characteristics, population structure, and regeneration dynamics of Schima species in Yunnan. Results Yunnan was found to harbor the greatest richness and the highest rarity-weighted richness of Schima species in the subtropical regions of China. We classified five primary and six secondary forest types containing Schima species as one of dominants. Yunnan had the high floristic diversity and varying stand structure of forests containing Schima species. The Schima species studied generally had a sporadic regeneration type and a long life-span. Four species (Schima argentea, Schima villosa, Schima sinensis, Schima sericans) were shade-intolerant. But three species (Schima noronhae, Schima khasiana and Schima wallichii) were considered as bi-modal type species having shade-intolerant and shade-tolerant traits. Schima noronhae was seen to be a top dominant in late successional forests, while S. wallichii was found as a top-dominant in early or middle or late successional forests. S. khasiana, Schima villosa, Schima sinensis usually appeared as a top dominant in early or middle successional secondary forests, though they also presented as a second dominant in late-successional forests. Schima argentea and Schima sericans dominated only in the early or middle/seral successional forests. Schima species’ regeneration establishment depended mainly on forest canopy gap formation through moderate human and natural disturbances. Conclusions Yunnan has high species richness and rarity-weighted richness of Schima. Both moderate human and natural disturbances have provided regeneration niches for Schima species. Some of the Schima species studied as a second dominant (rare as the top-dominant) present in the late-successional forests. Some of them are more often as the top-dominant in early or middle successional forests, where as time goes by the dominance of Schima species would be replaced by their associated dominant taxa such as Castanopsis species.

Female weevils of the family Attelabidae (Coleoptera: Curculionoidea) possess a unique behaviour of partially cutting the branches connecting egg-bearing organs of their host plants during oviposition. However, the consequence of such behaviour remains unclear. Using Rhynchites foveipennis and its host pear (Pyrus pyrifolia), the present study tested the hypothesis that the oviposition behaviour could disarm the host plants’ defence. We compared the survival rates, growth rates, and performance of eggs and larvae under two conditions: (1) the fruit stems were naturally damaged by the females before and after oviposition, and (2) the fruit stems were artificially protected from the females. When fruit stems were protected from female damage, the survival rates of eggs and larvae were only 21.3–32.6%, respectively; and the larval weight was 3.2–4.1 mg 30 days after laying eggs. When the fruit stems were damaged, the survival rates of eggs and larvae reached 86.1–94.0%, respectively; and the larval weight reached 73.0–74.9 mg 30 days after laying eggs. The contents of tannin and flavonoids in the pears did not change significantly along with the oviposition and larval feeding, but weevil eggs were crushed and killed by the callus in the pears. Once the stunted larvae in branch-growing pears were moved into the picked-off ones, the growth and development recovered. The findings indicate that the oviposition behaviour can significantly increase the survival of the offspring. Our study suggested that the oviposition behaviour of attelabid weevils is a strategy to overcome plant defence.

Plasmacytoid dendritic cells (pDCs) constitute a rare type of immune cell with multifaceted functions, but their potential use as a cell-based immunotherapy is challenged by the scarce cell numbers that can be extracted from blood. Here, we systematically investigate culture parameters for generating pDCs from hematopoietic stem and progenitor cells (HSPCs). Using optimized conditions combined with implementation of HSPC pre-expansion, we generate an average of 465 million HSPC-derived pDCs (HSPC-pDCs) starting from 100,000 cord blood-derived HSPCs. Furthermore, we demonstrate that such protocol allows HSPC-pDC generation from whole-blood HSPCs, and these cells display a pDC phenotype and function. Using GMP-compliant medium, we observe a remarkable loss of TLR7/9 responses, which is rescued by ascorbic acid supplementation. Ascorbic acid induces transcriptional signatures associated with pDC-specific innate immune pathways, suggesting an undescribed role of ascorbic acid for pDC functionality. This constitutes the first protocol for generating pDCs from whole blood and lays the foundation for investigating HSPC-pDCs for cell-based immunotherapy.

Rapidly evolving computational techniques are making a large gap between scientific aspiration and code implementation in climate modeling. In this work, we design a simple computing library to bridge the gap and decouple the work of ocean modeling from parallel computing. This library provides 12 basic operators that feature user-friendly interfaces, effective programming, and implicit parallelism. Several state-of-the-art computing techniques, including computing graph and just-in-time compiling, are employed to parallelize the seemingly serial code and speed up the ocean models. These operator interfaces are designed using native Fortran programming language to smooth the learning curve. We further implement a highly readable and efficient ocean model that contains only 1860 lines of code but achieves a 91 % parallel efficiency in strong scaling and 99 % parallel efficiency in weak scaling with 4096 Intel CPU cores. This ocean model also exhibits excellent scalability on the heterogeneous Sunway TaihuLight supercomputer. This work presents a promising alternative tool for the development of ocean models.

Across all forest systems, the primary driver of throughfall (TF) amount is rainfall amount (Pg) though no work has addressed the sensitivity of the TF response to projected shifts in Pg due to climate change. We assessed how climate change may impact TF sensitivity to variability in Pg for eleven typical forest sites across the main climate types of Iran using a nondimensional relative sensitivity coefficient. The Coupled Model Intercomparison Project phase 5 (CMIP5) HadGEM2-ES product was used under two emission scenarios (Representative Concentration Pathway (RCP) 2.6 and 8.5) to project yearly precipitation and Pg for the measurement sites during 2020–50. There was a strong linear relationship between TF and Pg at all sites [TF = 0.66 (Pg) – 0.30; R2 = 0.91; n = 639]. The sensitivity coefficient ranged from 0.96–5.3 across the eleven forest sites. Large sensitivity coefficient differences were found between small (< mean annual Pg) and large (> mean annual Pg) storms for arid plantations. To buffer expected shifts in storm size due to climate change, it may be suitable to incorporate TF sensitivity when choosing landscaping and urban greening. Shifts in Pg and increased small storm frequency are predicted for 2020–50 per CMIP5 HadGEM2-ES low and high emission scenarios.

LIDAR and photogrammetry techniques were used to quantify the land surface morphology dynamics of ephemeral gully system based on the field investigation. LIDAR monitoring results showed that ephemeral gully occurred in the same location on the slope surface after every rainy season in the ephemeral gully system, then lateral topsoil was tilled and brought into the ephemeral gully channel, inducing 2 cm decrease around channel before rainy season of next year, which make it a cycle of erosion-tillage-erosion. During this process, imbricated landform was formed in ephemeral gully system. The photogrammetry monitoring results showed that most drop-sills distances in ephemeral gully channel were in 10 to 25 cm, while slope gradients were among 15° to 40°. The drop-sills distance and slope gradient showed negative exponent relationship. The results of this study showed that LIDAR technology can quickly acquire the topographic characteristics of the whole ephemeral gully system, while photogrammetry method could quickly acquire detailed morphology in ephemeral gully channel.

Evapotranspiration (ET) has changed greatly across the world in the past decades. There are various drivers that may affect the changes in ET. Understanding the changes in ET and the potential drivers is of great importance to water resources management in the arid regions. In this study, we examine the trends in ET using satellite-based data in the Asian arid regions from 2003 to 2017. Trends in different ET components, i.e., the evaporation from bare soil and transpiration from vegetation, are estimated to examine their respective contributions to the changes in ET. The results show that the evaporation from bare soil is the main component of ET and has contributed larger to the changes in ET from 2003 to 2017. Statistically significant increasing trend is mostly found in Southwest Pakistan, Northwest India, Southeast Iran, many areas in Northwest China, and some parts in South Kazakhstan. The relations between ET and the global sea surface temperature (SST) are also examined. It suggests that the interannual variation of ET is closely related to annual SST anomalies in most areas, especially in Kazakhstan, Turkmenistan, Northwest China, Southwest Pakistan, and some parts of eastern Iran. This analysis provides useful information for water resources management and planning with respect to forest and agricultural development in the Asian arid regions.

The Chinese water pine Glyptostrobus pensilis is the sole surviving species of the genus Glyptostrobus. It is endemic to southern China, central Vietnam, and eastern Laos, and today it is nearly extinct in the wild. Forest community characteristics and population structure of G. pensilis in China have remained unknown up to now. We investigated six swamp forest stands and analyzed their forest community characteristics (i.e. vertical stratification, species composition, and diversity) and population structure, including the frequency distribution of DBH (diameter at breast height) and age-classes as found in Fujian Province, southeastern China.The vertical stratifications of all the forest stands were rather simple. The remaining wild specimens ranged from roughly 15 to some 357 years for an average of ca. 85 years, with only a few individuals less than 20 years old. Compared with the stands and populations of G. pensilis in Vietnam, the taxonomic compositions of the stands in the two regions were different, except for the dominant species-G. pensilis. The Shannon–Wiener index showed the overstory of each stand had much lower diversity (0.26 on average) in Fujian Province than that (1.97 on average) in Vietnam, whereas the diversity indices were about the same (around 2.41) for the understories in the two regions. Furthermore, we discovered 18 G. pensilis seedlings at the study sites in Fujian Province. This discovery demonstrates that G. pensilis regeneration is extremely poor and its populations are declining, although these populations are relatively healthier than those in Vietnam. Keywords: Age-class, Fujian Province, Glyptostrobus, Regeneration, Relict plant, Vietnam