<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Archiving and Interchange DTD v2.3 20070202//EN" "archivearticle.dtd">
<article article-type="methods-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Remote Sens.</journal-id>
<journal-title>Frontiers in Remote Sensing</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Remote Sens.</abbrev-journal-title>
<issn pub-type="epub">2673-6187</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1449854</article-id>
<article-id pub-id-type="doi">10.3389/frsen.2024.1449854</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Remote Sensing</subject>
<subj-group>
<subject>Methods</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Retrieval of aerosol particle size distribution using an improved l&#xe9;vy flight and circle chaos sparrow search algorithm</article-title>
<alt-title alt-title-type="left-running-head">Xun et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/frsen.2024.1449854">10.3389/frsen.2024.1449854</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Xun</surname>
<given-names>Lina</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2768548/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chu</surname>
<given-names>Yuting</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zeng</surname>
<given-names>Hao</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Siyu</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Yan</surname>
<given-names>Qing</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Jingjing</given-names>
</name>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff>
<institution>The Key Laboratory of Intelligent Computing and Signal Processing of Ministry of Education</institution>, <institution>School of Electrical Engineering and Automation</institution>, <institution>Anhui University</institution>, <addr-line>Hefei</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1011993/overview">Quan Wang</ext-link>, Shizuoka University, Japan</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1055158/overview">Zhenzhu Wang</ext-link>, CAS, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1802346/overview">Lei Liu</ext-link>, National University of Defense Technology, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Lina Xun, <email>ahuxln@ahu.edu.cn</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>30</day>
<month>08</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>5</volume>
<elocation-id>1449854</elocation-id>
<history>
<date date-type="received">
<day>16</day>
<month>06</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>21</day>
<month>08</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Xun, Chu, Zeng, Wang, Yan and Zhang.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Xun, Chu, Zeng, Wang, Yan and Zhang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>The aerosol size distribution is a crucial metric for evaluating aerosol optical characteristics and determines the direct and indirect radiative forcings of aerosols. Traditional methods for inversion of aerosol size distribution often suffer from ambiguities and limitations. To address these shortcomings, this paper introduces a method for inferring aerosol volume size distribution utilizing an improved L&#xe9;vy flight and chaotic sparrow search algorithm (ILCSSA). The algorithm incorporates Circle chaotic mapping onto the basic Sparrow Search Algorithm (SSA) to obtain a high-quality initial population and employs a L&#xe9;vy flight strategy to enhance population diversity. Additionally, the process of updating the population&#x2019;s positions is optimized to improve algorithm accuracy. To validate the feasibility of the proposed method, the measured aerosol optical depth (AOD) data obtained from a Precision Solar Radiometer (PSR) sun photometer in Shouxian are utilized and a series of comparisons were conducted among the Sparrow Search Algorithm, Standard Particle Swarm Algorithm, Improved Particle Swarm Algorithm, and Spider Wasp Optimization Algorithm. The results demonstrate a significant performance advantage for the ILCSSA, evidenced by an average reduction of 50<inline-formula id="inf1">
<mml:math id="m1">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> in Sum of Squared Errors (SSE) and 36<inline-formula id="inf2">
<mml:math id="m2">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> in Root Mean Squared Error (RMSE) when compared to the other four algorithms. Additionally, the AOD obtained by ILCSSA had a correlation coefficient of 0.9748 with the original AOD data. Furthermore, we analyzed the aerosol volume size distribution in Shouxian under conditions of good air quality, moderate pollution, and mild pollution. The proposed method holds significant reference value in the field of aerosol volume spectrum inversion.</p>
</abstract>
<kwd-group>
<kwd>sparrow search algorithm</kwd>
<kwd>aerosol size distribution</kwd>
<kwd>aerosol optical depth</kwd>
<kwd>sun photometer</kwd>
<kwd>L&#xe9;vy flight</kwd>
<kwd>circle chaos</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Multi- and Hyper-Spectral Imaging</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Aerosols are defined as minute particles suspended within a specific volume of gas, existing in the atmosphere in solid, liquid, or a combination of both states, thus forming a fundamental component of the atmospheric system. These microscopic aerosol particles can impact the balance of the atmosphere-Earth system through direct and indirect radiation <xref ref-type="bibr" rid="B9">Goldsmith et al. (1998)</xref>, consequently affecting weather conditions and climate change, as well as posing potential threats to human health. Atmospheric aerosols play a crucial role in climate, ecology, and human health <xref ref-type="bibr" rid="B4">Carslaw et al. (2010)</xref>.</p>
<p>Aerosols are characterized by a range of optical parameters, such as aerosol optical depth (AOD), scattering coefficient, phase function, and particle size distribution. AOD, recognized as a fundamental parameter, provides information on the combined impacts of aerosol light absorption and scattering across the atmosphere, thus serving as a crucial tool for assessing the degree of air pollution <xref ref-type="bibr" rid="B18">Xun et al. (2021)</xref>. Moreover, particle size distribution characterizes the distribution of aerosol particles with varying radii in the atmosphere, affecting other optical parameters and thereby influencing climate radiation effects. Accurate calculation of these optical parameters holds crucial significance for comprehensive research into the impacts of aerosols on human health and climate change <xref ref-type="bibr" rid="B13">Mao and Li (2014)</xref>.</p>
<p>Since the 20th century, researchers in the field of engineering science have encountered the challenge of inferring original images from measured values, a problem commonly known as the mathematical-physical inverse problem <xref ref-type="bibr" rid="B7">Engl et al. (1996)</xref>. AOD can be directly obtained through direct solar irradiance, and its relationship with particle size distribution follows the first type of Fredholm integral equation. The process of solving for the particle size distribution entails resolving the unknown variables in this integral equation. Nevertheless, because of the oscillations of the Mie kernel function, the aforementioned equation is ill-posed. The first type of Fredholm integral equation was originally proposed by Philips and Tikhonov and further developed by <xref ref-type="bibr" rid="B15">Twomey (1963)</xref>. This theoretical framework forms the basis for most inversions of particle size distribution using sun photometers. For example, King et al. utilized this theory to conduct linearly constrained inversion for aerosol particle size distribution. However, the inversion results exhibited significant errors because the actual aerosol particle size distribution may not necessarily conform to a smooth function. Furthermore, due to the nonlinearity of the problem and limitations in solutions, the uniqueness of the solution remains uncertain <xref ref-type="bibr" rid="B12">Lienert et al. (2001)</xref>. Hence, when the assumption of a smooth function is unreasonable, the Philips-Twomeny method becomes ineffective. Currently, there are numerous methods for aerosol particle size distribution inversion. These encompass regularization algorithms for multi-wavelength lidar data <xref ref-type="bibr" rid="B7">Engl et al. (1996)</xref>, improved particle swarm optimization algorithms <xref ref-type="bibr" rid="B14">Mao and Li (2015)</xref>, spectral extinction fitting techniques for Microtops II sun photometer data <xref ref-type="bibr" rid="B10">Lee et al. (2015)</xref>, genetic algorithms <xref ref-type="bibr" rid="B22">Ye et al. (1999)</xref>, and ant colony optimization algorithms <xref ref-type="bibr" rid="B9">Goldsmith et al. (1998)</xref>. The utilization of these methods provides researchers with additional options and opportunities to conduct comprehensive investigations into aerosol particle size distribution.</p>
<p>The Sparrow Search Algorithm is a novel swarm intelligence optimization algorithm introduced by Xue and Shen in 2020, inspired by the foraging and feeding behaviors of sparrows <xref ref-type="bibr" rid="B17">Xue and Shen (2020)</xref>. SSA is widely adopted in various domains owing to its straightforward structure, limited parameters, and rapid convergence. For instance, Xingjia Li employed SSA to improve the control precision of robot manipulators <xref ref-type="bibr" rid="B11">Li et al. (2022)</xref>. Shenghao Bi utilized SSA to design efficient workshop layouts <xref ref-type="bibr" rid="B2">Bi et al. (2024)</xref>; Ahmed Fathy applied SSA to optimize the operation of microgrids <xref ref-type="bibr" rid="B8">Fathy et al. (2022)</xref>. However, like other swarm intelligence algorithms, SSA faces challenges such as blind initialization of the population, limited global search capabilities, and susceptibility to local optima. To address these issues, we propose an improved algorithm called the Improved Levy Flight and Circle Chaotic Sparrow Search Algorithm (ILCSSA), which outperforms the standard SSA and other classical swarm intelligence algorithms in terms of computational accuracy and global optimization capability.</p>
<p>This study utilizes AOD data collected by a PSR sun photometer and applies ILCSSA for the inversion of aerosol volume size distribution. The paper provides a detailed description of the principles underlying aerosol particle size distribution based on ILCSSA. To validate the efficacy of ILCSSA and determine its feasibility in inverting volume distribution under varying air quality conditions, a series of experiments were conducted, and their outcomes were thoroughly compared with several other methods.</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>2 Materials and methods</title>
<sec id="s2-1">
<title>2.1 Aerosol volume size distribution</title>
<p>According to Mie scattering theory, assuming aerosol particles are spherical, the relationship between aerosol optical depth and aerosol size distribution is as follows:<disp-formula id="e1">
<mml:math id="m3">
<mml:mtable class="aligned">
<mml:mtr>
<mml:mtd columnalign="right"/>
<mml:mtd columnalign="left">
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mo>&#x222b;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>z</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:msubsup>
<mml:mrow>
<mml:mo>&#x222b;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3b7;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mi>&#x3c0;</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:mi>n</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">z</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">ext</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="normal">m</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>d</mml:mi>
<mml:mi>z</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="right"/>
<mml:mtd columnalign="left">
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mo>&#x222b;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3b7;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mi>&#x3c0;</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
<mml:msub>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">ext</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mi>n</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
<label>(1)</label>
</disp-formula>Where Qext(<inline-formula id="inf3">
<mml:math id="m4">
<mml:mi>&#x3bb;</mml:mi>
</mml:math>
</inline-formula>,r,m) is the extinction efficiency factor, <inline-formula id="inf4">
<mml:math id="m5">
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> and <inline-formula id="inf5">
<mml:math id="m6">
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> are the upper and lower limits of aerosol particle radius, respectively, and n(r, z) is the aerosol particle size distribution at a certain altitude. In addition, n(r, z)dz is the number of aerosol particles within the radius range from r to r &#x2b; dr at altitude z and <inline-formula id="inf6">
<mml:math id="m7">
<mml:msub>
<mml:mrow>
<mml:mi>z</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> is the height of the atmospheric column top.</p>
<p>A conversion relationship exists between aerosol volume size distribution and particle size distribution <xref ref-type="bibr" rid="B6">Dubovik and King (2000)</xref>. <xref ref-type="disp-formula" rid="e1">Equation 1</xref> can be rewritten as:<disp-formula id="e2">
<mml:math id="m8">
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mo>&#x222b;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mi>K</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mi>v</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:math>
<label>(2)</label>
</disp-formula>Additionally,<disp-formula id="e3">
<mml:math id="m9">
<mml:mi>K</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>4</mml:mn>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:msub>
<mml:mrow>
<mml:mi>Q</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">ext</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(3)</label>
</disp-formula>Where K(<inline-formula id="inf7">
<mml:math id="m10">
<mml:mi>&#x3bb;</mml:mi>
</mml:math>
</inline-formula>,r,m) denotes the efficiency kernel function, which signifies the sensitivity of AOD to variations in aerosol particle radius, as described in <xref ref-type="disp-formula" rid="e3">Equation 3</xref>. The relationship between its distribution across different wavelengths and radii is depicted in <xref ref-type="fig" rid="F1">Figure 1</xref>. According to Mie scattering theory, the extinction efficiency factor is intricately linked to the complex refractive index of aerosols, where the real component relates to the speed of light propagation and the imaginary component to light absorption <xref ref-type="bibr" rid="B3">Bohren and Huffman (2008)</xref>. Its values vary at different wavelengths, exhibiting oscillatory decay as the aerosol particle radius increases and gradually approaches a constant of 2.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>The curve of extinction efficiency kernel function corresponding to different wavelengths.</p>
</caption>
<graphic xlink:href="frsen-05-1449854-g001.tif"/>
</fig>
</sec>
<sec id="s2-2">
<title>2.2 Aerosol size distribution</title>
<p>In aerosol particle size distribution studies, distribution functions such as the single log-normal distribution and multiple log-normal distributions are commonly used to describe the characteristics of aerosols across different scales. These functions accurately reflect the variation in aerosol particle sizes in the atmosphere and are suitable for analyzing and modeling various types of aerosols.</p>
<sec id="s2-2-1">
<title>2.2.1 Single-lognormal spectrum</title>
<p>The single-lognormal distribution is a statistically significant aerosol particle size distribution function, expressed as follows:<disp-formula id="e4">
<mml:math id="m11">
<mml:mi>n</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>n</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
<mml:mi>exp</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="normal">l</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">l</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msup>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(4)</label>
</disp-formula>where N is the normalized particle number constant, <inline-formula id="inf8">
<mml:math id="m12">
<mml:mi>&#x3c3;</mml:mi>
</mml:math>
</inline-formula> is the standard deviation, and <inline-formula id="inf9">
<mml:math id="m13">
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> is the mean particle radius, as described in <xref ref-type="disp-formula" rid="e4">Equation 4</xref>.</p>
</sec>
<sec id="s2-2-2">
<title>2.2.2 Multi-lognormal distribution</title>
<p>In reality, aerosol patterns in the atmosphere are not singular but rather multifaceted. The multiple-peak lognormal distribution provides a better description of multi-source, multi-scale aerosol spectrum distribution. The expression for the number concentration spectrum is <xref ref-type="bibr" rid="B16">Whitby (1978)</xref>:<disp-formula id="e5">
<mml:math id="m14">
<mml:mi>n</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="bold">r</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>n</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="bold">r</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="bold">r</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mrow>
<mml:mo>&#x2211;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:munderover>
</mml:mstyle>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
<mml:mi>exp</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>ln</mml:mi>
<mml:mi mathvariant="bold">r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="bold">r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msubsup>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(5)</label>
</disp-formula>Where <inline-formula id="inf10">
<mml:math id="m15">
<mml:msub>
<mml:mrow>
<mml:mi>N</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>, <inline-formula id="inf11">
<mml:math id="m16">
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> and <inline-formula id="inf12">
<mml:math id="m17">
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> represent the aerosol particle number concentration, geometric standard deviation, and median radius, respectively.</p>
<p>While analyzing aerosol particle size distribution, number density spectra might not always capture the true particle distribution accurately. In contrast, volume size distribution offers substantial advantages. The expression for the volume spectrum is written as <xref ref-type="bibr" rid="B5">Dubovik et al. (2011)</xref>:<disp-formula id="e6">
<mml:math id="m18">
<mml:mi>&#x3bd;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="bold">r</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mi>&#x3bd;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="bold">r</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi>d</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi mathvariant="bold">r</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mrow>
<mml:mo>&#x2211;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:munderover>
</mml:mstyle>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>V</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
<mml:mi>exp</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">l</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="normal">l</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msubsup>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(6)</label>
</disp-formula>Where <inline-formula id="inf13">
<mml:math id="m19">
<mml:msub>
<mml:mrow>
<mml:mi>V</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>, <inline-formula id="inf14">
<mml:math id="m20">
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> and <inline-formula id="inf15">
<mml:math id="m21">
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> represent the volume concentration, standard deviation, and geometric mean radius, respectively. The symbol &#x2018;i&#x2019; in <xref ref-type="disp-formula" rid="e5">Equation 5</xref> and <xref ref-type="disp-formula" rid="e6">Equation 6</xref> indicates the aerosol mode, distinguishing between the coarse mode and fine mode. The result of substituting <xref ref-type="disp-formula" rid="e6">Equation 6</xref> into <xref ref-type="disp-formula" rid="e2">Equation 2</xref> is given in <xref ref-type="disp-formula" rid="e7">Equation 7</xref>.<disp-formula id="e7">
<mml:math id="m22">
<mml:mi>&#x3c4;</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mo>&#x222b;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:msubsup>
<mml:mi>K</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mrow>
<mml:mo>&#x2211;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:munderover>
</mml:mstyle>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>V</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msqrt>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:mfrac>
<mml:mi>exp</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>ln</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msubsup>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
<mml:mi>d</mml:mi>
<mml:mi>r</mml:mi>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
</sec>
</sec>
<sec id="s2-3">
<title>2.3 Principle of SSA</title>
<p>The Sparrow Search Algorithm <xref ref-type="bibr" rid="B17">Xue and Shen (2020)</xref> primarily achieves position optimization by mimicking the foraging behavior and reciprocity behavior of sparrows. Research has indicated that in captive sparrow populations, two distinct types exist: discoverers and followers. Discoverers possess the richest resources, responsible for locating food and providing direction and area foraging for the followers. The followers exploit the discoverer to obtain food. Additionally, when faced with external attacks, they use a danger alert mechanism. The mathematical model of SSA is presented as follows.</p>
<p>Assuming the search space is D-dimensional and there are n sparrows, the spatial positions of the population X can be represented as follows, as described in <xref ref-type="disp-formula" rid="e8">Equation 8</xref>:<disp-formula id="e8">
<mml:math id="m23">
<mml:mi>X</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfenced open="[" close="]">
<mml:mrow>
<mml:mtable class="matrix">
<mml:mtr>
<mml:mtd columnalign="center">
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1,1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd columnalign="center">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="bold">x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1,2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd columnalign="center">
<mml:mo>&#x22ef;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="bold">x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="center">
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2,1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd columnalign="center">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="bold">x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2,2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd columnalign="center">
<mml:mo>&#x22ef;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="bold">x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="center">
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
<mml:mtd columnalign="center">
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
<mml:mtd columnalign="center">
<mml:mo>&#x22ee;</mml:mo>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="center">
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd columnalign="center">
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="bold">x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mtd>
<mml:mtd columnalign="center">
<mml:mo>&#x22ef;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="bold">x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>d</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>Based on the randomly initialized positions, the fitness value of each sparrow can be calculated. Those with superior fitness are chosen as discoverers, and the remaining are considered followers. The updated formula for the positions of the discoverer is as follows:<disp-formula id="e9">
<mml:math id="m24">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>i,j</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mfenced open="{" close="">
<mml:mrow>
<mml:mtable class="cases">
<mml:mtr>
<mml:mtd columnalign="left">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>i,j</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>t</mml:mtext>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>max</mml:mtext>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
<mml:mspace width="1em"/>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>T</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="left">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>i,j</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mtext>t</mml:mtext>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="normal">Q</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi mathvariant="normal">L</mml:mi>
<mml:mspace width="1em"/>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>T</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(9)</label>
</disp-formula>
</p>
<p>Where t is the current iteration number, j &#x3d; 1,2,3, &#x2026; ,d. <inline-formula id="inf16">
<mml:math id="m25">
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> denotes the maximum iteration times, and <inline-formula id="inf17">
<mml:math id="m26">
<mml:mi>&#x3b1;</mml:mi>
</mml:math>
</inline-formula> is a random number between 0 and 1. R2 (R2<inline-formula id="inf18">
<mml:math id="m27">
<mml:mo>&#x2208;</mml:mo>
</mml:math>
</inline-formula>[0,1]) and ST (ST<inline-formula id="inf19">
<mml:math id="m28">
<mml:mo>&#x2208;</mml:mo>
</mml:math>
</inline-formula>[0.5,1.0]) respectively represent the warning and safety values. Q is a random number following a normal distribution. L represents a matrix of size 1&#x2a;D, where each element in the matrix is 1. If <inline-formula id="inf20">
<mml:math id="m29">
<mml:msub>
<mml:mrow>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mspace width="0.3333em"/>
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula> ST, it indicates that there are no predators nearby, and the discoverer can conduct an extensive search. Conversely, if <inline-formula id="inf21">
<mml:math id="m30">
<mml:msub>
<mml:mrow>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mspace width="0.3333em"/>
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula> ST, it means that some sparrows have detected predators, and the group needs to immediately evacuate to a safe area.</p>
<p>The update formula for the follower positions is given by <xref ref-type="disp-formula" rid="e10">Equation 10</xref>:<disp-formula id="e10">
<mml:math id="m31">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mfenced open="{" close="">
<mml:mrow>
<mml:mtable class="cases">
<mml:mtr>
<mml:mtd columnalign="left">
<mml:mi>Q</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>exp</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">worst</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
<mml:mspace width="1em"/>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mo>&#x3e;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="left">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mfenced open="|" close="|">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x22c5;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:mrow>
</mml:msup>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>L</mml:mi>
<mml:mspace width="1em"/>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mi>o</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(10)</label>
</disp-formula>
</p>
<p>Where <inline-formula id="inf22">
<mml:math id="m32">
<mml:msub>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> is the current best position of the discoverer, and <inline-formula id="inf23">
<mml:math id="m33">
<mml:msub>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">worst</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> denotes the current globally worst position. A represents a matrix of size 1&#x2a;D, where each element is randomly assigned either 1 or -1, <inline-formula id="inf24">
<mml:math id="m34">
<mml:msup>
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:mrow>
</mml:msup>
</mml:math>
</inline-formula> &#x3d; <inline-formula id="inf25">
<mml:math id="m35">
<mml:msup>
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mspace width="0.3333em"/>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
</mml:math>
</inline-formula>. When (i <inline-formula id="inf26">
<mml:math id="m36">
<mml:mo>&#x3e;</mml:mo>
</mml:math>
</inline-formula> n/2), it indicates that the <inline-formula id="inf27">
<mml:math id="m37">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula>th follower has not obtained food, with very low energy, and needs to fly elsewhere to forage for energy.</p>
<p>Opting for a minority of individuals from the population to function as sentinels, the formula for updating vigilantes&#x2019; positions is as follows:<disp-formula id="e11">
<mml:math id="m38">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mfenced open="{" close="">
<mml:mrow>
<mml:mtable class="matrix">
<mml:mtr>
<mml:mtd columnalign="left">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">best</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:mfenced open="|" close="|">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">best</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mtd>
<mml:mtd columnalign="center">
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3e;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="center">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>K</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mfenced open="|" close="|">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">worst</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b5;</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mtd>
<mml:mtd columnalign="center">
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(11)</label>
</disp-formula>Where <inline-formula id="inf28">
<mml:math id="m39">
<mml:msub>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">best</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> denotes the current global best solution. <inline-formula id="inf29">
<mml:math id="m40">
<mml:mi>&#x3b2;</mml:mi>
</mml:math>
</inline-formula> is a random number following the standard normal distribution. K is a random number within the interval [-1, 1], and <inline-formula id="inf30">
<mml:math id="m41">
<mml:mi>&#x3b5;</mml:mi>
</mml:math>
</inline-formula> serves to prevent division by zero. <inline-formula id="inf31">
<mml:math id="m42">
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> and <inline-formula id="inf32">
<mml:math id="m43">
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>w</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> respectively denote the current global best and worst fitness values.</p>
<p>In brief, if <inline-formula id="inf33">
<mml:math id="m44">
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mspace width="0.3333em"/>
<mml:mo>&#x3e;</mml:mo>
<mml:mspace width="0.3333em"/>
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>, it indicates that the sparrow at this position is situated at the edge of the group, making it susceptible to predator attacks. Conversely, if <inline-formula id="inf34">
<mml:math id="m45">
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> &#x3d; <inline-formula id="inf35">
<mml:math id="m46">
<mml:mi>f</mml:mi>
<mml:mi>i</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>, it signifies that the sparrow positioned in the middle of the group has recognized the danger and needs to move closer to neighboring sparrows to decrease the risk of being captured. Here, K represents the direction and stride of the sparrow&#x2019;s movement.</p>
</sec>
<sec id="s2-4">
<title>2.4 Principle of ILCSSA</title>
<p>During the initialization phase of the population, random positions are generated for various dimensions of sparrow individuals. This often leads to the local clustering of sparrow individuals. Hence, the introduction of Circle chaotic mapping is incorporated to enhance the distribution of the initial population. This allows the sparrow population to explore the entire space more extensively, thereby improving the optimization capability of the algorithm. The mathematical model is written as <xref ref-type="disp-formula" rid="e12">Equation 12</xref>:<disp-formula id="e12">
<mml:math id="m47">
<mml:msub>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="normal">m</mml:mi>
<mml:mi mathvariant="normal">o</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mfenced open="{" close="}">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mfenced open="[" close="]">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mi>&#x3c0;</mml:mi>
</mml:mrow>
</mml:mfrac>
<mml:mi>sin</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mi>&#x3c0;</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(12)</label>
</disp-formula>Where n is the dimensionality of the solution. <inline-formula id="inf36">
<mml:math id="m48">
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>, <inline-formula id="inf37">
<mml:math id="m49">
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>, and <inline-formula id="inf38">
<mml:math id="m50">
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> are constants whose values can be adjusted to alter the resulting outcomes. For the experiment, <inline-formula id="inf39">
<mml:math id="m51">
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> is chosen as 2, <inline-formula id="inf40">
<mml:math id="m52">
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> as 0.2, and <inline-formula id="inf41">
<mml:math id="m53">
<mml:msub>
<mml:mrow>
<mml:mi>d</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>3</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> as 0.5.</p>
<p>L&#xe9;vy flight <xref ref-type="bibr" rid="B20">Yang and Deb (2009)</xref> is a type of random walk model that simulates long-distance movement with random step lengths and directions. A key characteristic of this model is that its step lengths follow a L&#xe9;vy distribution, which has a long tail. This means there is a probability, albeit low, of encountering extreme values with large step lengths. This allows for both small and large step lengths, enhancing the likelihood of long-distance jumps <xref ref-type="bibr" rid="B21">Yang and Deb (2010)</xref>. Consequently, L&#xe9;vy flight is widely used in swarm intelligence optimization algorithms, effectively broadening the exploration range of the search space, enhancing the randomness of the swarm&#x2019;s search, and aiding in avoiding local optima <xref ref-type="bibr" rid="B19">Yang (2010)</xref>.</p>
<p>The formula for updating positions using L&#xe9;vy flight is given in <xref ref-type="disp-formula" rid="e13">Equation 13</xref>:<disp-formula id="e13">
<mml:math id="m54">
<mml:msubsup>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b4;</mml:mi>
<mml:mo>&#x2295;</mml:mo>
<mml:mi>L</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mi>y</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(13)</label>
</disp-formula>Where <inline-formula id="inf42">
<mml:math id="m55">
<mml:msubsup>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:math>
</inline-formula> indicates the position of <inline-formula id="inf43">
<mml:math id="m56">
<mml:msub>
<mml:mrow>
<mml:mi>x</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> in the <inline-formula id="inf44">
<mml:math id="m57">
<mml:mi>t</mml:mi>
</mml:math>
</inline-formula>th generation; <inline-formula id="inf45">
<mml:math id="m58">
<mml:mi>&#x3b4;</mml:mi>
</mml:math>
</inline-formula> represents the step size, with <inline-formula id="inf46">
<mml:math id="m59">
<mml:mi>&#x3b4;</mml:mi>
</mml:math>
</inline-formula> &#x3e; 0; <inline-formula id="inf47">
<mml:math id="m60">
<mml:mo>&#x2295;</mml:mo>
</mml:math>
</inline-formula> denotes point-by-point multiplication; L&#xe9;vy<inline-formula id="inf48">
<mml:math id="m61">
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> is the random search path, where the step lengths follow a L&#xe9;vy distribution. The mathematical model is given in <xref ref-type="disp-formula" rid="e14">Equation 14</xref>:<disp-formula id="e14">
<mml:math id="m62">
<mml:mi>L</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>v</mml:mi>
<mml:mi>y</mml:mi>
<mml:mo>&#x223c;</mml:mo>
<mml:mi>u</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x3c;</mml:mo>
<mml:mi>&#x3bb;</mml:mi>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>3</mml:mn>
</mml:math>
<label>(14)</label>
</disp-formula>
</p>
<p>L&#xe9;vy flight involves random step lengths that conform to a L&#xe9;vy distribution, which is very complex and has not been fully implemented. Therefore, the Mantegna algorithm is commonly used for simulation. The mathematical representation of this algorithm is given in <xref ref-type="disp-formula" rid="e15">Equation 15</xref>:<disp-formula id="e15">
<mml:math id="m63">
<mml:msub>
<mml:mrow>
<mml:mi>s</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="|" close="|">
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:math>
<label>(15)</label>
</disp-formula>
<disp-formula id="e16">
<mml:math id="m64">
<mml:mi>u</mml:mi>
<mml:mo>&#x223c;</mml:mo>
<mml:mi>N</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>u</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
<mml:mo>,</mml:mo>
<mml:mi>v</mml:mi>
<mml:mo>&#x223c;</mml:mo>
<mml:mi>N</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>,</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>v</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(16)</label>
</disp-formula>
<disp-formula id="e17">
<mml:math id="m65">
<mml:msub>
<mml:mrow>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>u</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mfenced open="{" close="}">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mi>sin</mml:mi>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>&#x3c0;</mml:mi>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">&#x393;</mml:mi>
<mml:mfenced open="[" close="]">
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfenced>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo>&#x2217;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfenced>
<mml:mo>/</mml:mo>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>/</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:msup>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3b4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>&#x3bd;</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:math>
<label>(17)</label>
</disp-formula>In <xref ref-type="disp-formula" rid="e16">Equation 16</xref>, u and v follow a normal distribution, and <inline-formula id="inf49">
<mml:math id="m66">
<mml:mi>&#x3b2;</mml:mi>
</mml:math>
</inline-formula> is typically set to 1.5, as shown in <xref ref-type="disp-formula" rid="e17">Equation 17</xref>.</p>
<p>According to <xref ref-type="disp-formula" rid="e9">Equation 9</xref>, when <inline-formula id="inf50">
<mml:math id="m67">
<mml:msub>
<mml:mrow>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mspace width="0.3333em"/>
<mml:mo>&#x3c;</mml:mo>
</mml:math>
</inline-formula> ST, the value of each dimension of the sparrow individual decreases with the number of iterations. Obviously, this is not an optimal strategy for global optimization. Therefore, the position update for discoverers in the SSA is:<disp-formula id="e18">
<mml:math id="m68">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mtext>i,j</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mfenced open="{" close="">
<mml:mrow>
<mml:mtable class="cases">
<mml:mtr>
<mml:mtd columnalign="left">
<mml:msubsup>
<mml:mrow>
<mml:mi mathvariant="bold">X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold">i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold">j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold">t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x22c5;</mml:mo>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold">Q</mml:mi>
</mml:mrow>
</mml:mfenced>
<mml:mspace width="1em"/>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn mathvariant="normal">2</mml:mn>
</mml:mrow>
</mml:msub>
<mml:mo>&#x3c;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>T</mml:mi>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="left">
<mml:msubsup>
<mml:mrow>
<mml:mi mathvariant="bold">X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold">i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold">j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold">t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="bold">Q</mml:mi>
<mml:mspace width="1em"/>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi mathvariant="bold">R</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold">2</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>T</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(18)</label>
</disp-formula>
</p>
<p>In each iteration, the sparrow population is uniformly distributed. When followers distance themselves from poor positions in all dimensions and move closer to the optimal position, they tend to converge quickly but are susceptible to local optima. Therefore, the position update formula for followers in the SSA is:<disp-formula id="e19">
<mml:math id="m69">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mfenced open="{" close="">
<mml:mrow>
<mml:mtable class="cases">
<mml:mtr>
<mml:mtd columnalign="left">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>r</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>d</mml:mi>
<mml:mo>&#x22c5;</mml:mo>
<mml:mfenced open="|" close="|">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
<mml:mspace width="1em"/>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mi>i</mml:mi>
<mml:mi>f</mml:mi>
<mml:mtext>i</mml:mtext>
<mml:mo>&#x3e;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd columnalign="left">
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:mfenced open="|" close="|">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:mfenced>
<mml:mo>&#x22c5;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mo>&#x2b;</mml:mo>
</mml:mrow>
</mml:msup>
<mml:mo>&#x22c5;</mml:mo>
<mml:mi>L</mml:mi>
<mml:mspace width="1em"/>
</mml:mtd>
<mml:mtd columnalign="left">
<mml:mi>o</mml:mi>
<mml:mi>t</mml:mi>
<mml:mi>h</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>w</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>e</mml:mi>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:math>
<label>(19)</label>
</disp-formula>
</p>
</sec>
<sec id="s2-5">
<title>2.5 Experimental inversion of aerosol particle size distribution using ILCSSA</title>
<p>Initially, the population is initialized using the Circle chaotic mapping. In <xref ref-type="disp-formula" rid="e6">Equation 6</xref>, <inline-formula id="inf51">
<mml:math id="m70">
<mml:msub>
<mml:mrow>
<mml:mi>V</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>, <inline-formula id="inf52">
<mml:math id="m71">
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> and <inline-formula id="inf53">
<mml:math id="m72">
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> represent unknown parameters, necessitating the determination of six position parameters. This implies that the dimensionality of the sparrow population is 6. <xref ref-type="table" rid="T1">Table 1</xref> outlines the range of each dimension for the particles, enabling the adjustment of the positions of the sparrow individuals.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Range of parameters for different Aerosol mode particles.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Aerosol modes</th>
<th align="left">
<inline-formula id="inf54">
<mml:math id="m73">
<mml:msub>
<mml:mrow>
<mml:mi>V</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
</th>
<th align="left">
<inline-formula id="inf55">
<mml:math id="m74">
<mml:msub>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
</th>
<th align="left">
<inline-formula id="inf56">
<mml:math id="m75">
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c3;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Fine Mode</td>
<td align="left">1&#x2013;0.001</td>
<td align="left">0.01&#x2013;1</td>
<td align="left">0.2&#x2013;0.8</td>
</tr>
<tr>
<td align="left">Coarse Mode</td>
<td align="left">1&#x2013;0.001</td>
<td align="left">1&#x2013;5</td>
<td align="left">0.2&#x2013;0.8</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In this study, the AOD measurements at six different wavelengths obtained from the PSR sun photometer are used as the true values, denoted by <inline-formula id="inf57">
<mml:math id="m76">
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>. The AOD values derived from <xref ref-type="disp-formula" rid="e7">Equation 7</xref> are referred to as the calculated values, <inline-formula id="inf58">
<mml:math id="m77">
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>. Since the goal of ILCSSA is to minimize the error between these values, the root mean square error (RMSE) is chosen as the fitness function value, <inline-formula id="inf59">
<mml:math id="m78">
<mml:mi>F</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
</mml:math>
</inline-formula>, which is calculated as shown in <xref ref-type="disp-formula" rid="e20">Equation 20</xref>:<disp-formula id="e20">
<mml:math id="m79">
<mml:mi>F</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>t</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>6</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mrow>
<mml:mo>&#x2211;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>6</mml:mn>
</mml:mrow>
</mml:munderover>
</mml:mstyle>
<mml:msubsup>
<mml:mrow>
<mml:mi>e</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msubsup>
</mml:mrow>
</mml:msqrt>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>6</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mrow>
<mml:mo>&#x2211;</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn>6</mml:mn>
</mml:mrow>
</mml:munderover>
</mml:mstyle>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>r</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c4;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
</mml:msqrt>
<mml:mtext>i</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1,2,3</mml:mn>
<mml:mo>,</mml:mo>
<mml:mo>&#x2026;</mml:mo>
<mml:mn>6</mml:mn>
<mml:mo>.</mml:mo>
</mml:math>
<label>(20)</label>
</disp-formula>Where <inline-formula id="inf60">
<mml:math id="m80">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> represents the measurement wavelength of the PSR sun photometer. <xref ref-type="fig" rid="F2">Figure 2</xref> shows the comprehensive flowchart of the aerosol volume size distribution inversion process using ILCSSA.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Flowchart illustrating the inversion of aerosol particle size distribution using the enhanced Circle chaotic mapping sparrow search algorithm.</p>
</caption>
<graphic xlink:href="frsen-05-1449854-g002.tif"/>
</fig>
<p>Step 1: Utilize the Circle chaotic mapping to generate unknown parameters, which will be employed in the subsequent inversion process of the aerosol volume size distribution.</p>
<p>Step 2: Employ <xref ref-type="disp-formula" rid="e6">Equation 6</xref> to invert the aerosol volume size distribution based on the generated parameters, and compute the corresponding AOD values using <xref ref-type="disp-formula" rid="e7">Equation 7</xref>.</p>
<p>Step 3: Compute the fitness function value (Fit) for each sparrow individual and rank them accordingly. Based on the sorted results, designate 20<inline-formula id="inf61">
<mml:math id="m81">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> of the individuals as discoverers, 70<inline-formula id="inf62">
<mml:math id="m82">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> as followers, and 10<inline-formula id="inf63">
<mml:math id="m83">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> as vigilantes, assigning roles for the subsequent optimization process.</p>
<p>Step 4: Retrieve the current individual positions <inline-formula id="inf64">
<mml:math id="m84">
<mml:msub>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> and the global best position <inline-formula id="inf65">
<mml:math id="m85">
<mml:msub>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> for each sparrow, and record their corresponding fitness function values <inline-formula id="inf66">
<mml:math id="m86">
<mml:msub>
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula> and the global best fitness value <inline-formula id="inf67">
<mml:math id="m87">
<mml:msub>
<mml:mrow>
<mml:mi>G</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">best</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>.</p>
<p>Step 5: Update the positions of discoverers, followers, and vigilantes using <xref ref-type="disp-formula" rid="e18">Equation 18</xref>, <xref ref-type="disp-formula" rid="e19">Equation 19</xref>, and <xref ref-type="disp-formula" rid="e11">Equation 11</xref>, respectively, and recalculate the fitness function values. Repeat the process from Step 2 until the set number of iterations is achieved.</p>
<p>Step 6: After completing the iterations, output the global best position <inline-formula id="inf68">
<mml:math id="m88">
<mml:msub>
<mml:mrow>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi>g</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>, and substitute it into <xref ref-type="disp-formula" rid="e6">Equation 6</xref> to obtain the final aerosol volume size distribution.</p>
</sec>
</sec>
<sec sec-type="results|discussion" id="s3">
<title>3 Results and discussion</title>
<sec id="s3-1">
<title>3.1 The comparison and analysis of multiple algorithms</title>
<p>The research employed measurement data from the PSR sun photometer for analysis. As shown in <xref ref-type="fig" rid="F3">Figure 3</xref>, the PSR sun photometer is located at the Shouxian National Meteorological Observatory (32.56&#xb0;N, 116.78&#xb0;E). This device features broad spectral coverage, with each scan covering 9 bands, comprising 340&#xa0;nm and 380&#xa0;nm in the ultraviolet spectrum, 440&#xa0;nm, 500&#xa0;nm, and 670&#xa0;nm in the visible spectrum, and 870&#xa0;nm, 940&#xa0;nm, 1020&#xa0;nm, and 1640&#xa0;nm in the infrared spectrum.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>PSR sun photometer device utilized in the experiment.</p>
</caption>
<graphic xlink:href="frsen-05-1449854-g003.tif"/>
</fig>
<p>The sun photometer is engineered with complete automation, enabling direct irradiance observation through sun tracking and sky radiance measurement. Developed by Anhui Zhongke Puruida Optoelectronics Co., Ltd., detailed instrument parameters and other relevant information can be referenced from <xref ref-type="bibr" rid="B18">Xun et al. (2021)</xref>.</p>
<p>The aerosol size distribution inversion method based on the Improved Particle Swarm Optimization Algorithm (IPSO) proposed by Mao et al <xref ref-type="bibr" rid="B14">Mao and Li (2015)</xref>. has received widespread recognition and application in the field of atmospheric science. This method enhances global optimization capabilities by adjusting the speed and position of particles and has effectively solved the accurate inversion of aerosol size distribution under various weather conditions. Mao et al.&#x2018;s approach has demonstrated exceptional performance, particularly in dealing with nonlinear and ill-posed problems, which is crucial for aerosol size distribution inversion.</p>
<p>Mao et al.&#x2019;s research not only provides a theoretical foundation for subsequent studies but has also validated its effectiveness in several practical applications. For instance, comparisons with traditional methods show that their approach offers higher accuracy and robustness under dust and haze conditions. Consequently, Mao et al.&#x2019;s method is regarded as an important reference standard for aerosol size distribution inversion.</p>
<p>In this context, to evaluate the feasibility of the ILCSSA algorithm, we have compared it with SSA, PSO, IPSO, and SMO <xref ref-type="bibr" rid="B1">Abdel-Basset et al. (2023)</xref> algorithms. <xref ref-type="table" rid="T2">Table 2</xref> provides a detailed list of the experimental parameters of these algorithms for further analysis.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Range of parameters for different Aerosol mode particles.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Inversion method</th>
<th align="left">ILCSSA</th>
<th align="left">PSO</th>
<th align="left">SSA</th>
<th align="left">IPSO</th>
<th align="left">SMO</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Population Size</td>
<td align="left">30</td>
<td align="left">30</td>
<td align="left">30</td>
<td align="left">30</td>
<td align="left">30</td>
</tr>
<tr>
<td align="left">Convergence Threshold</td>
<td align="left">0.8</td>
<td align="left"/>
<td align="left">0.8</td>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="left">Maximum Iterations</td>
<td align="left">25</td>
<td align="left">50</td>
<td align="left">50</td>
<td align="left">50</td>
<td align="left">50</td>
</tr>
<tr>
<td align="left">Initial Inertia Weight <inline-formula id="inf69">
<mml:math id="m89">
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">ini</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
</td>
<td align="left"/>
<td align="left">0.9</td>
<td align="left"/>
<td align="left">0.9</td>
<td align="left"/>
</tr>
<tr>
<td align="left">Final Inertia Weight <inline-formula id="inf70">
<mml:math id="m90">
<mml:msub>
<mml:mrow>
<mml:mi>&#x3c9;</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="italic">end</mml:mi>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
</td>
<td align="left"/>
<td align="left">0.4</td>
<td align="left"/>
<td align="left">0.4</td>
<td align="left"/>
</tr>
<tr>
<td align="left">Cognitive Parameters <inline-formula id="inf71">
<mml:math id="m91">
<mml:msub>
<mml:mrow>
<mml:mi>c</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
</td>
<td align="left"/>
<td align="left">1.8</td>
<td align="left"/>
<td align="left">1.8</td>
<td align="left"/>
</tr>
<tr>
<td align="left">Social Parameters <inline-formula id="inf72">
<mml:math id="m92">
<mml:msub>
<mml:mrow>
<mml:mi>c</mml:mi>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:msub>
</mml:math>
</inline-formula>
</td>
<td align="left"/>
<td align="left">1.8</td>
<td align="left"/>
<td align="left">1.8</td>
<td align="left"/>
</tr>
<tr>
<td align="left">Discoverers Ratio</td>
<td align="left">
<inline-formula id="inf73">
<mml:math id="m93">
<mml:mn>20</mml:mn>
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula>
</td>
<td align="left"/>
<td align="left">
<inline-formula id="inf74">
<mml:math id="m94">
<mml:mn>20</mml:mn>
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula>
</td>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="left">Follower Ratio</td>
<td align="left">
<inline-formula id="inf75">
<mml:math id="m95">
<mml:mn>70</mml:mn>
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula>
</td>
<td align="left"/>
<td align="left">
<inline-formula id="inf76">
<mml:math id="m96">
<mml:mn>70</mml:mn>
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula>
</td>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="left">Sentinels Ratio</td>
<td align="left">
<inline-formula id="inf77">
<mml:math id="m97">
<mml:mn>10</mml:mn>
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula>
</td>
<td align="left"/>
<td align="left">
<inline-formula id="inf78">
<mml:math id="m98">
<mml:mn>10</mml:mn>
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula>
</td>
<td align="left"/>
<td align="left"/>
</tr>
<tr>
<td align="left">Dimension d</td>
<td align="left">6</td>
<td align="left">6</td>
<td align="left">6</td>
<td align="left">6</td>
<td align="left">6</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>
<xref ref-type="fig" rid="F4">Figure 4</xref> shows the variations in fitness values over different numbers of iterations for the ILCSSA, SSA, PSO, IPSO, and SMO methods. Despite initially demonstrating faster convergence rates, both the ILCSSA and PSO algorithms show that after 10 iterations, the ILCSSA&#x2019;s fitness value decreases more rapidly than that of PSO. Furthermore, compared to the SSA, IPSO, and SMO methods, the ILCSSA significantly addresses premature convergence issues while achieving superior results in terms of both convergence speed and accuracy.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>The variation of fitness function values with evolutionary generations for the ILCSSA, SSA, PSO, IPSO, and SMO methods.</p>
</caption>
<graphic xlink:href="frsen-05-1449854-g004.tif"/>
</fig>
<p>The aerosol volume size distributions retrieved by the ILCSSA, SSA, PSO, IPSO, and SMO algorithms are presented in <xref ref-type="fig" rid="F5">Figure 5</xref>. To further evaluate the inversion performance of these five algorithmic sets, we also randomly selected the sum of squared errors and the mean square error between the calculated and true values of eight sets of AOD at six wavelengths, as depicted in <xref ref-type="fig" rid="F6">Figure 6</xref>. Notably, the superiority of the ILCSSA algorithm is evident.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Aerosol volume size distributions retrieved by various algorithms.</p>
</caption>
<graphic xlink:href="frsen-05-1449854-g005.tif"/>
</fig>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Using the Sparrow Search Algorithm and the Circle Chaotic Sparrow Search Algorithm based on improved Levy flight for solving the sum of squared errors <bold>(A)</bold> and mean square errors <bold>(B)</bold>.</p>
</caption>
<graphic xlink:href="frsen-05-1449854-g006.tif"/>
</fig>
<p>The correlation between the calculated AOD values obtained through the ILCSSA algorithm and the measured values is presented in <xref ref-type="fig" rid="F7">Figure 7</xref>. We randomly selected 24 data sets to assess their correlation, achieving a maximum correlation coefficient of 0.9748. This result further confirms the superiority of the ILCSSA algorithm in aerosol inversion.</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>The correlation between the original and calculated values of aerosol optical depth using the Circle Chaotic Sparrow Search Algorithm based on improved Levy flight.</p>
</caption>
<graphic xlink:href="frsen-05-1449854-g007.tif"/>
</fig>
</sec>
<sec id="s3-2">
<title>3.2 Inversion and analysis of volume size distribution</title>
<p>To assess the practical application of ILCSSA, we describe the characteristics of atmospheric aerosol volume size distribution in the Shouxian. Experimental data were selected from periods of excellent, good, and mild pollution air quality in Shouxian.</p>
<p>
<xref ref-type="fig" rid="F8">Figure 8</xref> shows the diurnal variations in AOD and volume size distribution based on the ILCSSA algorithm, during a period of excellent air quality in Shouxian County on 1 August 2020. It is evident from <xref ref-type="fig" rid="F8">Figure 8A</xref> that there is a pronounced similarity between the wavelength-dependent trend of AOD and its diurnal variation. At the same moment, the relationship between AOD values and wavelength exhibits an inverse correlation: longer wavelengths correspond to smaller AOD values and trends, whereas shorter wavelengths correspond to larger AOD values and trends. The temporal variations in AOD values across different wavelengths also demonstrate similarities. Following sunrise in the morning, rapid heating of the Earth&#x2019;s surface induces a gradual increase in atmospheric temperature, leading to intense convection. This convective process transports surface pollutants upwards to higher altitudes, leading to an increase in AOD during the morning hours. As human activities intensify around noon, emissions from vehicular exhausts and industrial pollutants escalate, coupled with the introduction of ground-level dust into the atmosphere. These factors collectively contribute to a rise in AOD, peaking at 12:54 PM. Subsequently, AOD values start declining in the afternoon due to the decrease in surface temperature and reduced human activities, facilitating the settling of particulate matter back to the ground. As depicted in <xref ref-type="fig" rid="F8">Figure 8B</xref>, the volume size distributions at four different time points exhibit bimodal spectra, mirroring the trends observed in AOD values, with higher aerosol volume concentrations in the morning and lower concentrations in the afternoon. Influenced by various factors, the amplitude of the aerosol volume size distribution curves for the day is substantial, indicating a predominance of anthropogenic fine particulate matter in the atmospheric pollution of Shouxian County.</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Aerosol optical depth <bold>(A)</bold> and volume size distribution <bold>(B)</bold> based on improved L&#xe9;vy Flight and Circle Chaos Sparrow Search Algorithm for excellent air day in Shouxian on 1 August 2020.</p>
</caption>
<graphic xlink:href="frsen-05-1449854-g008.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F9">Figure 9</xref> presents AOD and aerosol volume size distribution in Shouxian on 14 March 2021, during a period of favorable air quality, as determined through the ILCSSA algorithm. Despite the overcast weather conditions, AOD values remained relatively steady until 11:42 a.m. Subsequently, over, for the next hour, there was a slight decrease followed by an increase back to the initial values, peaking abruptly at 2:51 p.m., after which AOD values began declining across all spectral bands. Analysis of the volume spectra at different time intervals reveals marginal discrepancies in volume concentration throughout the day, yet the primary and secondary peaks exhibit consistent values, suggesting a stable overall structure of atmospheric aerosols in Shouxian County for the day. The retrieved aerosol volume size distribution demonstrates minimal variability, consistent with the observed good air quality conditions.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Aerosol optical depth <bold>(A)</bold> and volume size distribution <bold>(B)</bold> based on improved L&#xe9;vy Flight and Circle Chaos Sparrow Search Algorithm for good air day in Shouxian on 14 March 2021.</p>
</caption>
<graphic xlink:href="frsen-05-1449854-g009.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F10">Figure 10</xref> illustrates the AOD and aerosol volume size distribution in Shouxian County on 4 December 2020, during a period of mild pollution, as derived using the ILCSSA algorithm. In <xref ref-type="fig" rid="F10">Figure 10A</xref>, before 12:40 p.m., AOD values exhibit significant fluctuations at wavelengths of 340&#xa0;nm, 440&#xa0;nm, 500&#xa0;nm, and 675&#xa0;nm, while remaining relatively stable at 870&#xa0;nm and 1020&#xa0;nm. Subsequently, AOD values commence an ascent, reaching their peak at 2:51 p.m., indicating a pronounced influence of fine particles on afternoon AOD levels. <xref ref-type="fig" rid="F10">Figure 10B</xref> portrays a comparable trend in volume size distribution at 9:18 a.m. and 11:33 a.m., characterized by relatively low volume concentrations, followed by a marked increase in particles smaller than 1 micron in volume concentration at 1:30 p.m. and 2:33 p.m. as aerosols accumulate. Relative to the volume size distribution of aerosols at various time points on 1 August 2020, and 14 March 2021, the aerosol volume concentration on this day is elevated, indicative of poor air quality conditions consistent with the prevailing weather. Furthermore, the elevated presence of fine particles suggests that anthropogenic emissions of fine particulate matter primarily influenced the air quality in Shouxian County on that day.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Aerosol optical depth <bold>(A)</bold> and volume size distribution <bold>(B)</bold> based on improved L&#xe9;vy Flight and Circle Chaos Sparrow Search Algorithm for mild pollution day in Shouxian on 4 December 2020.</p>
</caption>
<graphic xlink:href="frsen-05-1449854-g010.tif"/>
</fig>
<p>The AOD values for 14 March 2021, were slightly lower or comparable to those for 4 December 2020, but exhibited instability. On 4 December 2020, the air quality index (AQI) was 104, classified as lightly polluted, while on 14 March 2021, the AQI was 94, classified as good. Notably, the AOD values for the mild pollution air day were higher than those for the good air quality day and this phenomenon is further analysed.</p>
<p>This discrepancy may be attributed to differences in meteorological conditions and ozone concentrations. 4 December 2020, was a clear day, whereas 14 March 2021, was cloudy. Cloud cover likely reduced solar radiation and increased atmospheric humidity, which could have led to aerosol particle condensation or growth, affecting the accuracy of AOD measurements and resulting in lower AOD values. Additionally, the higher ozone concentration on 14 March 2021 (128 <inline-formula id="inf79">
<mml:math id="m99">
<mml:mi mathvariant="normal">&#x3bc;</mml:mi>
</mml:math>
</inline-formula>g/m3 compared to 64 <inline-formula id="inf80">
<mml:math id="m100">
<mml:mi mathvariant="normal">&#x3bc;</mml:mi>
</mml:math>
</inline-formula>g/m3), may have intensified the absorption in the ozone layer, thereby reducing the scattering and absorption of sunlight by aerosols, further influencing the stability and accuracy of the AOD retrieval.</p>
</sec>
</sec>
<sec sec-type="conclusions" id="s4">
<title>4 Conclusions</title>
<p>Aerosols constitute a vital component of the atmosphere, exerting significant influence on atmospheric radiation balance, air quality, and human health. The aerosol size distribution is a crucial parameter for assessing aerosol optical properties. In this study, the ILCSSA algorithm was successfully employed to overcome the limitations of traditional methods in solving particle spectra. This approach initializes the population using a typical Circle chaotic mapping, optimizes the positions of explorers and followers through optimization, and introduces Levy flights to prevent convergence to local optima.</p>
<p>Compared to SSA, PSO, IPSO, and SMO algorithms, the ILCSSA algorithm demonstrates superior performance in terms of computational time and accuracy. Specifically, ILCSSA achieves an average reduction of 60<inline-formula id="inf81">
<mml:math id="m101">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> in SSE and 43<inline-formula id="inf82">
<mml:math id="m102">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> in RMSE compared to SSA; a 30<inline-formula id="inf83">
<mml:math id="m103">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> reduction in SSE and a 17<inline-formula id="inf84">
<mml:math id="m104">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> reduction in RMSE compared to PSO; a 52<inline-formula id="inf85">
<mml:math id="m105">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> reduction in SSE and a 34<inline-formula id="inf86">
<mml:math id="m106">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> reduction in RMSE compared to IPSO; and a 57<inline-formula id="inf87">
<mml:math id="m107">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> reduction in SSE and a 51<inline-formula id="inf88">
<mml:math id="m108">
<mml:mi>%</mml:mi>
</mml:math>
</inline-formula> reduction in RMSE compared to SMO. Additionally, the correlation coefficient between the AOD obtained using ILCSSA and the original AOD data is 0.974, reflecting high accuracy and reliability. By utilizing AOD data obtained from a PSR sun photometer, we successfully inverted aerosol volume spectra under different air quality conditions in Shouxian County, validating the feasibility of the ILCSSA method. The proposed method holds significant promise in the field of aerosol volume spectrum inversion.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s5">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6">
<title>Author contributions</title>
<p>LX: Conceptualization, Investigation, Methodology, Writing&#x2013;original draft. YC: Methodology, Software, Writing&#x2013;original draft. HZ: Data curation, Writing&#x2013;original draft. SW: Writing&#x2013;review and editing. QY: Validation, Writing&#x2013;review and editing. JZ: Supervision, Writing&#x2013;review and editing.</p>
</sec>
<sec sec-type="funding-information" id="s7">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article. This work was supported in part by the National Key R &#x26; D Plan of China under Grant 2018YFB0504604.</p>
</sec>
<ack>
<p>Authors are grateful to the referees and handling editor for their constructive comments.</p>
</ack>
<sec sec-type="COI-statement" id="s8">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s9">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abdel-Basset</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Mohamed</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Jameel</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Abouhawwash</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Spider wasp optimizer: a novel meta-heuristic optimization algorithm</article-title>. <source>Artif. Intell. Rev.</source> <volume>56</volume>, <fpage>11675</fpage>&#x2013;<lpage>11738</lpage>. <pub-id pub-id-type="doi">10.1007/s10462-023-10446-y</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bi</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Shao</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>R.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Workshop layout optimization method based on sparrow search algorithm: a new approach</article-title>. <source>J. Industrial Prod. Eng.</source> <volume>41</volume>, <fpage>324</fpage>&#x2013;<lpage>343</lpage>. <pub-id pub-id-type="doi">10.1080/21681015.2024.2302630</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Bohren</surname>
<given-names>C. F.</given-names>
</name>
<name>
<surname>Huffman</surname>
<given-names>D. R.</given-names>
</name>
</person-group> (<year>2008</year>). <source>Absorption and scattering of light by small particles</source>. <publisher-name>John Wiley and Sons</publisher-name>.</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Carslaw</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Boucher</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Spracklen</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Mann</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Rae</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Woodward</surname>
<given-names>S.</given-names>
</name>
<etal/>
</person-group> (<year>2010</year>). <article-title>A review of natural aerosol interactions and feedbacks within the earth system</article-title>. <source>Atmos. Chem. Phys.</source> <volume>10</volume>, <fpage>1701</fpage>&#x2013;<lpage>1737</lpage>. <pub-id pub-id-type="doi">10.5194/acp-10-1701-2010</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dubovik</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Herman</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Holdak</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Lapyonok</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Tanr&#xe9;</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Deuz&#xe9;</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2011</year>). <article-title>Statistically optimized inversion algorithm for enhanced retrieval of aerosol properties from spectral multi-angle polarimetric satellite observations</article-title>. <source>Atmos. Meas. Tech.</source> <volume>4</volume>, <fpage>975</fpage>&#x2013;<lpage>1018</lpage>. <pub-id pub-id-type="doi">10.5194/amt-4-975-2011</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dubovik</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>King</surname>
<given-names>M. D.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>A flexible inversion algorithm for retrieval of aerosol optical properties from sun and sky radiance measurements</article-title>. <source>J. Geophys. Res. Atmos.</source> <volume>105</volume>, <fpage>20673</fpage>&#x2013;<lpage>20696</lpage>. <pub-id pub-id-type="doi">10.1029/2000jd900282</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Engl</surname>
<given-names>H. W.</given-names>
</name>
<name>
<surname>Hanke</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Neubauer</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>1996</year>). <source>Regularization of inverse problems</source>, <volume>375</volume>. <publisher-name>Springer Science and Business Media</publisher-name>.</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fathy</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Alanazi</surname>
<given-names>T. M.</given-names>
</name>
<name>
<surname>Rezk</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Yousri</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Optimal energy management of micro-grid using sparrow search algorithm</article-title>. <source>Energy Rep.</source> <volume>8</volume>, <fpage>758</fpage>&#x2013;<lpage>773</lpage>. <pub-id pub-id-type="doi">10.1016/j.egyr.2021.12.022</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Goldsmith</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Blair</surname>
<given-names>F. H.</given-names>
</name>
<name>
<surname>Bisson</surname>
<given-names>S. E.</given-names>
</name>
<name>
<surname>Turner</surname>
<given-names>D. D.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>Turn-key Raman lidar for profiling atmospheric water vapor, clouds, and aerosols</article-title>. <source>Appl. Opt.</source> <volume>37</volume>, <fpage>4979</fpage>&#x2013;<lpage>4990</lpage>. <pub-id pub-id-type="doi">10.1364/ao.37.004979</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lee</surname>
<given-names>K. H.</given-names>
</name>
<name>
<surname>Wong</surname>
<given-names>M. S.</given-names>
</name>
<name>
<surname>Nichol</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Chan</surname>
<given-names>P. W.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Retrieval of aerosol size distribution from microtops ii sunphotometer in Hong Kong</article-title>. <source>Aerosol Air Qual. Res.</source> <volume>15</volume>, <fpage>1712</fpage>&#x2013;<lpage>1719</lpage>. <pub-id pub-id-type="doi">10.4209/aaqr.2015.01.0048</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Gu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Parameter identification of robot manipulators with unknown payloads using an improved chaotic sparrow search algorithm</article-title>. <source>Appl. Intell.</source> <volume>52</volume>, <fpage>10341</fpage>&#x2013;<lpage>10351</lpage>. <pub-id pub-id-type="doi">10.1007/s10489-021-02972-5</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lienert</surname>
<given-names>B. R.</given-names>
</name>
<name>
<surname>Porter</surname>
<given-names>J. N.</given-names>
</name>
<name>
<surname>Sharma</surname>
<given-names>S. K.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Repetitive genetic inversion of optical extinction data</article-title>. <source>Appl. Opt.</source> <volume>40</volume>, <fpage>3476</fpage>&#x2013;<lpage>3482</lpage>. <pub-id pub-id-type="doi">10.1364/ao.40.003476</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Dust particle size distribution inversion based on the multi population genetic algorithm</article-title>. <source>TAO Terr. Atmos. Ocean. Sci.</source> <volume>25</volume>, <fpage>791</fpage>. <pub-id pub-id-type="doi">10.3319/tao.2014.06.12.01(a)</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Retrieval of particle size distribution from aerosol optical thickness using an improved particle swarm optimization algorithm</article-title>. <source>Opt. Rev.</source> <volume>22</volume>, <fpage>809</fpage>&#x2013;<lpage>818</lpage>. <pub-id pub-id-type="doi">10.1007/s10043-015-0130-5</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Twomey</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>1963</year>). <article-title>On the numerical solution of Fredholm integral equations of the first kind by the inversion of the linear system produced by quadrature</article-title>. <source>J. ACM (JACM)</source> <volume>10</volume>, <fpage>97</fpage>&#x2013;<lpage>101</lpage>. <pub-id pub-id-type="doi">10.1145/321150.321157</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Whitby</surname>
<given-names>K. T.</given-names>
</name>
</person-group> (<year>1978</year>). &#x201c;<article-title>The physical characteristics of sulfur aerosols</article-title>,&#x201d; in <source>Sulfur in the atmosphere</source> (<publisher-name>Elsevier</publisher-name>), <fpage>135</fpage>&#x2013;<lpage>159</lpage>.</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xue</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>A novel swarm intelligence optimization approach: sparrow search algorithm</article-title>. <source>Syst. Sci. and control Eng.</source> <volume>8</volume>, <fpage>22</fpage>&#x2013;<lpage>34</lpage>. <pub-id pub-id-type="doi">10.1080/21642583.2019.1708830</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xun</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Qian</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Lyu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Analysis of aerosol optical depth from sun photometer at shouxian, China</article-title>. <source>Atmosphere</source> <volume>12</volume>, <fpage>1226</fpage>. <pub-id pub-id-type="doi">10.3390/atmos12091226</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>X.-S.</given-names>
</name>
</person-group> (<year>2010</year>). <source>Nature-inspired metaheuristic algorithms</source>. <publisher-loc>Frome, United Kingdom</publisher-loc>: <publisher-name>Luniver press</publisher-name>.</citation>
</ref>
<ref id="B20">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>X.-S.</given-names>
</name>
<name>
<surname>Deb</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2009</year>). &#x201c;<article-title>Cuckoo search via l&#xe9;vy flights</article-title>,&#x201d; in <source>2009 World congress on nature and biologically inspired computing (NaBIC)</source> (<publisher-loc>Ieee</publisher-loc>), <fpage>210</fpage>&#x2013;<lpage>214</lpage>.</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>X.-S.</given-names>
</name>
<name>
<surname>Deb</surname>
<given-names>S.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Engineering optimisation by cuckoo search</article-title>. <source>Int. J. Math. Model. Numer. Optimisation</source> <volume>1</volume>, <fpage>330</fpage>&#x2013;<lpage>343</lpage>. <pub-id pub-id-type="doi">10.1504/ijmmno.2010.035430</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ye</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Inversion of particle-size distribution from angular light-scattering data with genetic algorithms</article-title>. <source>Appl. Opt.</source> <volume>38</volume>, <fpage>2677</fpage>&#x2013;<lpage>2685</lpage>. <pub-id pub-id-type="doi">10.1364/ao.38.002677</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>