Characterization of individual spores of two biological insecticides, Bacillus thuringiensis and Lysinibacillus sphaericus, in response to glutaraldehyde using single-cell optical approaches.

Bacillus thuringiensis (Bt) and Lysinibacillus sphaericus (Ls) are the most widely used microbial insecticides. Both encounter unfavorable environmental factors and pesticides in the field. Here, the responses of Bt and Ls spores to glutaraldehyde were characterized using Raman spectroscopy and differential interference contrast imaging at the single-cell level. Bt spores were more sensitive to glutaraldehyde than Ls spores under prolonged exposure: <1.0% of Bt spores were viable after 10 min of 0.5% (v/v) glutaraldehyde treatment, compared to ~ 20% of Ls spores. The Raman spectra of glutaraldehyde-treated Bt and Ls spores were almost identical to those of untreated spores; however, the germination process of individual spores was significantly altered. The time to onset of germination, the period of rapid Ca2+-2,6-pyridinedicarboxylic acid (CaDPA) release, and the period of cortex hydrolysis of treated Bt spores were significantly longer than those of untreated spores, with dodecylamine germination being particularly affected. Similarly, the germination of treated Ls spores was significantly prolonged, although the prolongation was less than that of Bt spores. Although the interiors of Bt and Ls spores were undamaged and CaDPA did not leak, proteins and structures involved in spore germination could be severely damaged, resulting in slower and significantly prolonged germination. This study provides insights into the impact of glutaraldehyde on bacterial spores at the single cell level and the variability in spore response to glutaraldehyde across species and populations.

Novosphingobium album sp. nov., Novosphingobium organovorum sp. nov. and Novosphingobium mangrovi sp. nov. with the organophosphorus pesticides degrading ability isolated from mangrove sediments.

Members of the genus Novosphingobium were frequently isolated from polluted environments and possess great bioremediation potential. Here, three species, designated B2637T, B2580T and B1949T, were isolated from mangrove sediments and might represent novel species in the genus Novosphingobium based on a polyphasic taxonomy study. Phylogenomic analysis revealed that strains B2580T, B1949T and B2637T clustered with Novosphingobium naphthalenivorans NBRC 102051T, 'N. profundi' F72 and N. decolorationis 502str22T, respectively. The average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values between isolates and their closely related species were less than 94 and 54 %, respectively, all below the threshold of species discrimination. The sizes of the genomes of isolates B2580T, B2637T and B1949T ranged from 4.4 to 4.6 Mb, containing 63.3-66.4 % G+C content. Analysis of their genomic sequences identified genes related to pesticide degradation, heavy-metal resistance, nitrogen fixation, antibiotic resistance and sulphur metabolism, revealing the biotechnology potential of these isolates. Except for B2637T, B1949T and B2580T were able to grow in the presence of quinalphos. Results from these polyphasic taxonomic analyses support the affiliation of these strains to three novel species within the genus Novosphingobium, for which we propose the name Novosphingobium album sp. nov. B2580T (=KCTC 72967T=MCCC 1K04555T), Novosphingobium organovorum sp. nov. B1949T (=KCTC 92158T=MCCC 1K03763T) and Novosphingobium mangrovi sp. nov. B2637T (KCTC 72969T=MCCC 1K04460T).

Taxonomy and anticancer potential of Streptomyces niphimycinicus sp. nov. against nasopharyngeal carcinoma cells.

Streptomyces species are ubiquitous, Gram-positive, spore-forming bacteria with the ability to produce various clinically relevant compounds. The strain 4503 T was isolated from mangrove sediments, showing morphological and chemical properties which were consistent with those of members of the genus Streptomyces. Phylogenetic analysis based on 16S rRNA gene sequences revealed that the isolate was primarily identified as members of the genus Streptomyces, sharing more than 99% sequence identity to Streptomyces yatensis DSM 41771 T, S. antimycoticus NBRC 12839 T, and S. melanosporofaciens NBRC 13061 T. Average nucleotide identities (ANI) and digital DNA-DNA hybridization (dDDH) values between strain 4503 T and its close relatives were all below 95-96% and 75% of the novel species threshold, respectively. Results from phylogenetic, genomic, phenotypic, and chemotaxonomic characteristics analyses confirmed that the isolate represented a novel species of the genus Streptomyces, for which the name Streptomyces niphimycinicus sp. nov. 4503 T (= MCCC 1K04557T = JCM 34996 T) is proposed. The bioassay-guided fractionation of the extract of strain 4503 T resulted in the isolation of a known compound niphimycin C, which showed cytotoxic activity against nasopharyngeal carcinoma (NPC) cell lines TW03 and 5-8F with half maximal inhibitory concentration (IC50) values of 12.24 µg/mL and 9.44 µg/mL, respectively. Further experiments revealed that niphimycin C not only exhibited the capacity of anti-proliferation, anti-metastasis, induction of cell cycle arrest, and apoptosis, but was also able to increase the reactive oxygen species (ROS) production and regulate several signaling pathways in NPC cells. KEY POINTS: • Strain 4503 T was classified as a novel species of Streptomyces. • Niphimycin C correlates with the cytotoxic effect of strain 4503 T against NPC cells. • Niphimycin C induces apoptosis, autophagic flux disruption and cell cycle arrest.

Single-Acquisition 2-D Multifocal Raman Spectroscopy Using Compressive Sensing.

Confocal Raman microscopy is a powerful method for nondestructive and noninvasive detection of chemicals with high spatial resolution, but its long acquisition time hinders its applications in large-scale monitoring of fast dynamics. Here, we report the development of a compressive sensing technique for single-acquisition multifocal Raman spectroscopy, which is capable of improving the speed of conventional confocal Raman spectroscopy by 2-3 orders of magnitude. A sample is excited with a 2-D multifocus pattern, and the Raman scatterings from the multiple foci were projected onto the spectrometer's entrance in a 2-D array. The superimposed spectra within each row of the array were processed with an algorithm such that the spectra from the individual foci were retrieved in a single acquisition and with reduced noise. The performances of the developed technique were demonstrated by parallel Raman spectroscopy of multiple individual particles as well as by single-acquisition confocal Raman imaging of a large scale with high spatial resolution when combined with spatially sparse sampling. The technique is expected to find wide applications in investigating fast dynamics in large-scale biological systems.

Parallel micro-Raman spectroscopy of multiple cells in a single acquisition using hierarchical sparsity.

Parallel micro-Raman spectroscopy can significantly expand the analytical capacity of single biological cells. By positioning the Raman spectra of multiple trapped cells on a detector array along the grating dispersion direction, the throughput of single-cell analysis can be improved by orders of magnitude. However, accurate retrieval of the individual spectra from the superimposed spectrum in a single acquisition presents great challenges. In this work, we developed a hierarchical sparsity method under a compressive sensing framework. The method combined a group-selection strategy with in-group sparsity for spectral reconstruction. The performances of the developed method were demonstrated with both simulated and experimental data, and the Raman spectra of the individual trapped cells were retrieved with both high accuracy and low noises; especially, with a group-selection mechanism, the developed method successfully avoided wrong selection of the eigenspectra for spectral reconstruction. The technique is expected to find wide applications in simultaneous monitoring of long biological processes of multiple cells by Raman spectroscopy.

Single-cell analysis reveals the effects of glutaraldehyde and formaldehyde on individual Nosema bombycis spores.

Nosema bombycis (Nb) is the pathogen that causes pebrine in silkworms. Aldehydes are effective disinfectants commonly used in sericulture. However, the precise mechanism of their action on Nb spores remains unclear. Here, we used laser tweezers Raman spectroscopy to investigate the effects of glutaraldehyde and formaldehyde on individual Nb spores, as well as phase contrast microscopy imaging to monitor the germination dynamics of individual treated spores, to acquire a deeper understanding of the mechanism of action of aldehydes and to provide a theoretical reference for establishing an effective strategy for disease control in sericulture. The positions of the Raman peaks remained constant during treatment. The Raman intensity was enhanced and the germination rate of the spores significantly decreased with treatment time. Tlag, the time when individual spores begin to germinate, and Tgerm, the time for complete germination, increased with enhanced treatment. The germination time (ΔTgerm) showed no significant difference from that for untreated spores. Heterogeneity was shown, which is relevant to the resistance of Nb spores to aldehydes. The results indicate that glutaraldehyde and formaldehyde do not destroy the spore wall and plasma membrane, do not cause the leakage of intracellular components, and might not damage the extrusion apparatus. The effects of aldehydes on Nb spores are mainly on the spore coat. They may block the external factors that stimulate spore germination. Single-cell analysis based on novel optical techniques reveals the action of chemical sporicides on microsporidia spores in real time and explains the heterogeneity of cell stress resistance. These applications of new techniques offer new insight into traditional disinfectants.

Single-cell Raman spectroscopy reveals microsporidia spore heterogeneity in various insect hosts.

Single-cell Raman spectroscopy was used to analyze the spore heterogeneity of 16 microsporidia strains from various insect hosts in order to better understand the basic biology of microsporidia. The Raman spectrum of a single spore revealed basic spore composition, and microsporidia spores in various hosts were found to be rich in trehalose. Principal component analysis and Raman intensity showed obvious heterogeneity in the trehalose, nucleic acid, and protein content of various spores; however, there was no correlation between various spore groups and host type. Trehalose content correlated with spore infectivity on Bombyx mori. Raman spectroscopy is an excellent tool for label-free investigation of intercellular molecular constituents, providing insight into the heterogeneity of microsporidia spores.

Probing the germination kinetics of ethanol-treated Bacillus thuringiensis spores.

Bacillus thuringiensis (Bt) is the most widely used microbial insecticide. To clarify the mechanism of bacterial resistance to ethanol toxicity, the present study investigated the effects of 70% (v/v) ethanol at a moderate temperature (65°C) on Bt spore germination by single-cell Raman spectroscopy and differential interference contrast microscopy. We found that over 80% of Bt spores were inviable after 30 min of treatment. Moreover, ethanol treatment affected spore germination; the time for initiation of rapid calcium dipicolinate (CaDPA) release (i.e., lag time, Tlag), time taken for rapid CaDPA release (i.e., ΔTrelease), and time required for complete hydrolysis of the peptidoglycan cortex of spores (i.e., ΔTlys) were increased with longer treatment times. Alanine-initiated germination upon ethanol treatment for 30-90 min showed a 2- to 4-fold longer Tlag, 2- to 3.5-fold longer ΔTrelease, and ∼2-fold longer ΔTlys relative to the control. Dodecylamine-initiated germination treated for 15-30 min had 3- to 5-fold longer Tlag and 1.4- to 1.7-fold longer ΔTrelease than the control. Germination induced by exogenous CaDPA was observed only in a small fraction of spores treated with ethanol for 5 min. Single-cell Raman spectroscopy revealed that more than 52% of spores lost CaDPA after 30 min of ethanol treatment; these showed reductions in the intensity of 1280 and 1652 cm-1 bands (corresponding to protein α-helical structure) and increases in that of 1245 and 1665 cm-1 bands (attributed to irregularities in protein structure). These results indicate that CaDPA in the core of Bt spores confers resistance to ethanol, and that damage to the spore inner membrane by ethanol treatment results in CaDPA leakage. Additionally, moderate-temperature ethanol treatment and consequent denaturation of germination-related proteins affected spore germination, specifically by inactivating the cortex-lytic enzyme CwlJ. Our findings provide a theoretical basis for the development of more effective methods for killing spore-forming bacteria; microscopy imaging and Raman spectroscopy can provide novel insight into the effects of chemical agents on microbial cells.

Probing the Kinetic Anabolism of Poly-Beta-Hydroxybutyrate in Cupriavidus necator H16 Using Single-Cell Raman Spectroscopy.

Poly-beta-hydroxybutyrate (PHB) can be formed in large amounts in Cupriavidus necator and is important for the industrial production of biodegradable plastics. In this investigation, laser tweezers Raman spectroscopy (LTRS) was used to characterize dynamic changes in PHB content-as well as in the contents of other common biomolecule-in C. necator during batch growth at both the population and single-cell levels. PHB accumulation began in the early stages of bacterial growth, and the maximum PHB production rate occurred in the early and middle exponential phases. The active biosynthesis of DNA, RNA, and proteins occurred in the lag and early exponential phases, whereas the levels of these molecules decreased continuously during the remaining fermentation process until the minimum values were reached. The PHB content inside single cells was relatively homogenous in the middle stage of fermentation; during the late growth stage, the variation in PHB levels between cells increased. In addition, bacterial cells in various growth phases could be clearly discriminated when principle component analysis was performed on the spectral data. These results suggest that LTRS is a valuable single-cell analysis tool that can provide more comprehensive information about the physiological state of a growing microbial population.

Compressive Raman imaging by combining scattering-projection interleaving with context-aware excitation.

Exciting an object with a laser-focus array and randomly interleaving its scattering projection has been proved to be an effective strategy for speeding up Raman imaging. The so-called scattering interleaved Raman imaging (SIRI) method allows Raman hyperspectral imaging with a single snapshot and exhibits excellent reconstruction fidelity and signal-to-noise ratios (SNRs). Here, we show that the performance of SIRI is significantly improved when combined with context-aware excitation. The experiments on micro-plastics demonstrate that the restriction of Raman excitation within a smaller region of interest as guided by bright-field microscopy improves the signal intensity and the SNR, and it is surprising that the spectral resolution is also significantly improved. The context-aware SIRI method is successfully used for imaging of lipid-producing yeast cells, suggesting that it is a promising analytical tool for studying live cells or tissues.

Gaining insights into the responses of individual yeast cells to ethanol fermentation using Raman tweezers and chemometrics.

Saccharomyces cerevisiae is the most common microbe used for the industrial production of bioethanol, and it encounters various stresses that inhibit cell growth and metabolism during fermentation. However, little is currently known about the physiological changes that occur in individual yeast cells during ethanol fermentation. Therefore, in this work, Raman spectroscopy and chemometric techniques were employed to monitor the metabolic changes of individual yeast cells at distinct stages during high gravity ethanol fermentation. Raman tweezers was used to acquire the Raman spectra of individual yeast cells. Multivariate curve resolution-alternating least squares (MCR-ALS) and principal component analysis were employed to analyze the Raman spectra dataset. MCR-ALS extracted the spectra of proteins, phospholipids, and triacylglycerols and their relative contents in individual cells. Changes in intracellular biomolecules showed that yeast cells undergo three distinct physiological stages during fermentation. In addition, heterogeneity among yeast cells significantly increased in the late fermentation period, and different yeast cells may respond to ethanol stress via different mechanisms. Our findings suggest that the combination of Raman tweezers and chemometrics approaches allows for characterizing the dynamics of molecular components within individual cells. This approach can serve as a valuable tool in investigating the resistance mechanism and metabolic heterogeneity of yeast cells during ethanol fermentation.

Rapid and accurate identification of marine bacteria spores at a single-cell resolution by laser tweezers Raman spectroscopy and deep learning.

Marine bacteria have been considered as important participants in revealing various carbon/sulfur/nitrogen cycles of marine ecosystem. Thus, how to accurately identify rare marine bacteria without a culture process is significant and valuable. In this work, we constructed a single-cell Raman spectra dataset from five living bacteria spores and utilized convolutional neural network to rapidly, accurately, nondestructively identify bacteria spores. The optimal CNN architecture can provide a prediction accuracy of five bacteria spore as high as 94.93% ± 1.78%. To evaluate the classification weight of extracted spectra features, we proposed a novel algorithm by occluding fingerprint Raman bands. Based on the relative classification weight arranged from large to small, four Raman bands located at 1518, 1397, 1666, and 1017 cm-1 mostly contribute to producing such high prediction accuracy. It can be foreseen that, LTRS combined with CNN approach have great potential for identifying marine bacteria, which cannot be cultured under normal condition.

[Baseline correction of Raman spectrum based on piecewise linear fitting].

The baseline drifts of Raman spectra occur in many types of instrumental measurements. It is an important part and a routine step to correct the baseline drift for the data preprocessing. In the present work, the limitations of the baseline correction method based on polynomial fitting were highlighted and a modified polynomial fitting method, i. e. piecewise linear fitting method, was proposed. Combined with the computer, this method could eliminate the baseline automatically. A series of Raman spectra of single polystyrene bead, red blood cell or yeast cell acquired by laser tweezers Raman spectroscopy were preprocessed by this method and its efficiency was verified. The results demonstrated that piecewise linear fitting can correct the baseline shifts effectively and provides more accurate information for further data analysis. It is a feasible method for correction of Raman spectrum.

Copula-Based Assessment and Regionalization of Drought Risk in China.

Droughts are widespread in China and have brought considerable losses to the economy and society. Droughts are intricate, stochastic processes with multi-attributes (e.g., duration, severity, intensity, and return period). However, most drought assessments tend to focus on univariate drought characteristics, which are inadequate to describe the intrinsic characteristics of droughts due to the existence of correlations between drought attributes. In this study, we employed the standardized precipitation index to identify drought events using China's monthly gridded precipitation dataset from 1961 to 2020. Univariate and copula-based bivariate methods were then used to examine drought duration and severity on 3-, 6-, and 12-month time scales. Finally, we used the hierarchical cluster method to identify drought-prone regions in mainland China at various return periods. Results revealed that time scale played an essential role in the spatial heterogeneity of drought behaviors, such as average characteristics, joint probability, and risk regionalization. The main findings were as follows: (1) 3- and 6-month time scales yielded comparable regional drought features, but not 12-month time scales; (2) higher drought severity was associated with longer drought duration; (3) drought risk was higher in the northern Xinjiang, western Qinghai, southern Tibet, southwest China, and the middle and lower reaches of the Yangtze River, and lower in the southeastern coastal areas of China, the Changbai Mountains, and the Greater Khingan Mountains; (4) mainland China was divided into six subregions according to joint probabilities of drought duration and severity. Our study is expected to contribute to better drought risk assessment in mainland China.

New insights into the effects of UV light on individual Nosema bombycis spores, as determined using single-cell optical approaches.

Nosema bombycis (Nb) is a pathogen causing pebrine in sericulture. Ultraviolet (UV) light exposure is a common physical disinfection method, but the mechanisms underlying UV-based disinfection have only been studied at the population level. In this study, changes in and germination of UV-irradiated spores were observed using Raman tweezers and phase-contrast imaging to evaluate the effects of UV radiation on Nb spores at the single-cell level. We found that irradiation caused the complete leakage of trehalose from individual spores. We also found that more spores leaked as the UV dose increased. There was no significant loss of intracellular biomacromolecules and no marked changes in the peaks associated with protein secondary structures. Low-dose radiation promoted spore germination and high-dose radiation decreased the germination rate, while the germination time did not undergo significant alterations. These results suggest that UV radiation disrupts the permeability of the inner membrane and alters the spore wall, thereby affecting the ability of the spore to sense and respond to extracellular stimuli, which further triggers germination and reduces or stops spore germination. This study provides new insights into the molecular mechanisms underlying conventional disinfection measures on microsporidian spores.

Determining optimal probability distributions for gridded precipitation data based on L-moments.

Understanding the probability distributions of precipitation is crucial for predicting climatic events and constructing hydraulic facilities. To overcome the inadequacy of precipitation data, regional frequency analysis was commonly used by "trading space for time". However, with the increasing availability of gridded precipitation datasets with high spatial and temporal resolutions, the probability distributions of precipitation for these datasets have been less explored. We used L-moments and goodness-of-fit criteria to identify the probability distributions of annual, seasonal, and monthly precipitation for a 0.5° × 0.5° dataset across the Loess Plateau (LP). We examined five 3-parameter distributions, namely General Extreme Value (GEV), Generalized Logistic (GLO), Generalized Pareto (GPA), Generalized Normal (GNO), and Pearson type III (PE3), and evaluated the accuracy of estimated rainfall using the leave-one-out method. We also presented pixel-wise fit-parameters and quantiles of precipitation as supplements. Our findings indicated that precipitation probability distributions vary by location and time scale, and the fitted probability distribution functions are reliable for estimating precipitation under various return periods. Specifically, for annual precipitation, GLO was prevalent in humid and semi-humid areas, GEV in semi-arid and arid areas, and PE3 in cold-arid areas. For seasonal precipitation, spring precipitation mainly conforms to GLO distribution, summer precipitation around the 400 mm isohyet prevalently follows GEV distribution, autumn precipitation primarily meets GPA and PE3 distributions, and winter precipitation in the northwest, south, and east of the LP mainly conforms to GPA, PE3 and GEV distributions, respectively. Regarding monthly precipitation, the common distribution functions are PE3 and GPA for the less-precipitation months, whereas the distribution functions of precipitation for more-precipitation months vary substantially across different regions of the LP. Our study contributes to a better understanding of precipitation probability distributions in the LP and provides insights for future studies on gridded precipitation datasets using robust statistical methods.

CCUS development in China and forecast its contribution to emission reduction.

Nowadays environmental issues have been of great concern to the world, among which the problem of global warming caused by greenhouse gas emissions is particularly prominent. All countries in the Kyoto Protocol and the Paris Agreement have committed to control greenhouse gas emissions, and China, as the largest carbon emitter, has assumed a heavier burden. China has been striving to develop low-carbon technologies such as hydrogen, nuclear, wind, and solar energy, but the most attention should be paid to CCUS, which many scholars have high expectations that CCUS can help China reduce emissions to some extent. Therefore, this paper presents a prediction that CCUS can reduce 3.8% of carbon emissions for China in 2040 when CCUS emission reductions increase at a rate of 30%. The power and chemical industries could reduce carbon emissions by 2.3% and 17.3%, respectively.

Large Improvement of Thermoelectric Performance by Magnetism in Co-Based Full-Heusler Alloys.

Full-Heusler alloys (fHAs) exhibit high mechanical strength with earth-abundant elements, but their metallic properties tend to display small electron diffusion thermopower, limiting potential applications as excellent thermoelectric (TE) materials. Here, it is demonstrated that the Co-based fHAs Co2 XAl (X = Ti, V, Nb) exhibit relatively high thermoelectric performance due to spin and charge coupling. Thermopower contributions from different magnetic mechanisms, including spin fluctuation and magnon drag are extracted. A significant contribution to thermopower from magnetism compared to that from electron diffusion is demonstrated. In Co2 TiAl, the contribution to thermopower from spin fluctuation is higher than that from electron diffusion, resulting in an increment of 280 µW m-1  K-2 in the power factor value. Interestingly, the thermopower contribution from magnon drag can reach up to -47 µV K-1 , which is over 2400% larger than the electron diffusion thermopower. The power factor of Co2 TiAl can reach 4000 µW m-1  K-2 which is comparable to that of conventional semiconducting TE materials. Moreover, the corresponding figure of merit zT can reach ≈0.1 at room temperature, which is significantly larger than that of traditional metallic materials. The work shows a promising unconventional way to create and optimize TE materials by introducing magnetism.

Analysis of the effects of a gerP mutation on the germination of spores of Bacillus subtilis.

As previously reported, gerP Bacillus subtilis spores were defective in nutrient germination triggered via various germinant receptors (GRs), and the defect was eliminated by severe spore coat defects. The gerP spores' GR-dependent germination had a longer lag time between addition of germinants and initiation of rapid release of spores' dipicolinic acid (DPA), but times for release of >90% of DPA from individual spores were identical for wild-type and gerP spores. The gerP spores were also defective in GR-independent germination by DPA with its associated Ca(2+) divalent cation (CaDPA) but germinated better than wild-type spores with the GR-independent germinant dodecylamine. The gerP spores exhibited no increased sensitivity to hypochlorite, suggesting that these spores have no significant coat defect. Overexpression of GRs in gerP spores did lead to faster germination via the overexpressed GR, but this was still slower than germination of comparable gerP(+) spores. Unlike wild-type spores, for which maximal nutrient germinant concentrations were between 500 µM and 2 mM for l-alanine and ≤10 mM for l-valine, rates of gerP spore germination increased up to between 200 mM and 1 M l-alanine and 100 mM l-valine, and at 1 M l-alanine, the rates of germination of wild-type and gerP spores with or without all alanine racemases were almost identical. A high pressure of 150 MPa that triggers spore germination by activating GRs also triggered germination of wild-type and gerP spores identically. All these results support the suggestion that GerP proteins facilitate access of nutrient germinants to their cognate GRs in spores' inner membrane.

Application of L-moment method for regional frequency analysis of meteorological drought across the Loess Plateau, China.

Water shortages have always been the primary bottleneck for the healthy and sustainable development of the ecological environment on the Loess Plateau (LP). Proper water resource management requires knowledge of the spatiotemporal characteristics of precipitation frequency. This paper employed the gridded precipitation dataset obtained from the China Meteorological Data Service Centre to present a spatially explicit characterization of precipitation frequencies in tandem with their return periods on the LP based on the L-moment method. The 60% and 80% of the mean annual precipitation from 1981 to 2010 were synonymous with severe and moderate droughts, respectively. Droughts occurred more frequently in the northwest than in the southeast of the LP. Moreover, the frequencies of moderate drought showed a slight difference throughout the area, while those of severe droughts demonstrated considerable differences between the northwestern arid zone and the southeastern semi-humid zone. The maps associated with various return periods of precipitation deficits can be used to produce drought risk maps together with drought vulnerability maps. These findings could also provide useful information for drought management, water resource management and the development of food security policies.