<?xml version='1.0' encoding='UTF-8'?>
<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Publishing DTD v1.1d1 20130915//EN" "JATS-journalpublishing1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta id="journal-meta-aaee4b1957f449ef9e28a1dadde08f2b">
      <journal-id journal-id-type="nlm-ta">Sciresol</journal-id>
      <journal-id journal-id-type="publisher-id">Sciresol</journal-id>
      <journal-id journal-id-type="journal_submission_guidelines">http://ugit.net/publication_fsjoaj3qdho/geographical-analysis_su-zbsigk49/</journal-id>
      <journal-title-group>
        <journal-title>Geographical Analysis</journal-title>
      </journal-title-group>
      <issn publication-format="electronic">XXXX-XXXX</issn>
      <issn publication-format="print"/>
    </journal-meta>
    <article-meta id="article-meta-9389c5c00e944c14a18817fb62174ef6">
      <article-id pub-id-type="doi">10.53989/bu.ga.v13i2.18</article-id>
      <article-categories>
        <subj-group>
          <subject>Research Article</subject>
        </subj-group>
      </article-categories>
      <title-group>
        <article-title id="article-title-b9bf4ef4d47746b2b8afaa643a7b891a">
          <bold id="strong-427f99fee8f74457984bff827c3774f1">Comprehensive Analysis of Land Surface Temperature Changes in </bold>
          <bold id="strong-7bdaa8803de34becb773bb33ee83f1cf">Chikmagalur Taluk Using Landsat 8 Level 1 Data</bold>
        </article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <name id="name-4e6e5d55b31542bc82cbb74f99f8aa3b">
            <surname>Soundarya</surname>
            <given-names>R</given-names>
          </name>
          <xref id="xref-f07374a3794249cd9c3cee3a256f7f32" rid="aff-98e6f5274fe6474da4d3687362c9f3cf" ref-type="aff">1</xref>
        </contrib>
        <contrib contrib-type="author" corresp="yes">
          <name id="name-f3eeacd7c39f42ad804223a5fd32ba77">
            <surname>Anil</surname>
            <given-names>Sawant Sushant</given-names>
          </name>
          <email>geo_sushant@jssuni.edu.in</email>
          <xref id="xref-ec594861cbe445299c8865464b190e95" rid="aff-98e6f5274fe6474da4d3687362c9f3cf" ref-type="aff">1</xref>
        </contrib>
        <aff id="aff-98e6f5274fe6474da4d3687362c9f3cf">
          <institution>School of Life Sciences, JSS Academy of Higher Education and Research</institution>
          <addr-line>Mysuru, Karnataka</addr-line>
          <country country="IN">India</country>
        </aff>
      </contrib-group>
      <volume>13</volume>
      <issue>2</issue>
      <fpage>10</fpage>
      <permissions>
        <copyright-year>2024</copyright-year>
      </permissions>
      <abstract id="abstract-abstract-title-7477c19e6faa4bb0b49c85a27c00338f">
        <title id="abstract-title-7477c19e6faa4bb0b49c85a27c00338f">Abstract</title>
        <p id="paragraph-3be20bcbd06c496b9babb6fd377b30ae">This comprehensive report presents a meticulous geospatial analysis of land surface temperature (LST) changes in Chikmagalur Taluk, Karnataka, India, leveraging Landsat 8 Level 1 satellite imagery. The study aims to understand the temporal dynamics of LST and its environmental implications. We detail a comprehensive methodology encompassing radiometric calibration, spectral indices, and advanced spatial modeling. Notably, the results reveal a diverse pattern of LST changes, with the eastern part of the Taluk showing localized decreases, while the hill regions in the west and central areas exhibit notable temperature increases. These findings suggest a potential link between recent environmental incidents and LST fluctuations in the region.</p>
      </abstract>
      <kwd-group id="kwd-group-aa8013b917b0497a8cdbaf4aad924fff">
        <title>Keywords</title>
        <kwd>LST</kwd>
        <kwd>Landsat 8</kwd>
        <kwd>Level 1</kwd>
        <kwd>Band 10</kwd>
        <kwd>band 4</kwd>
        <kwd>band 5</kwd>
      </kwd-group>
      <funding-group>
        <funding-statement>Nil</funding-statement>
      </funding-group>
    </article-meta>
  </front>
  <body>
    <sec>
      <title id="t-630ca65bbe12">1 Introduction</title>
      <p id="t-422721eb4f6f">A comprehensive understanding of land surface temperature is crucial for monitoring environmental changes and informed urban planning. Chikmagalur Taluk, nestled within the Western Ghats, presents a unique and dynamic landscape undergoing rapid transformations due to urbanization and climatic variations. This report embarks on an in-depth analysis of LST fluctuations in Chikmagalur Taluk, utilizing Landsat 8 Level 1 data to unravel the thermal dynamics that influence the region's ecosystems and urban development <xref id="x-000950a901fe" rid="R238433831159379" ref-type="bibr">1</xref>.</p>
    </sec>
    <sec>
      <title id="title-b227d024987a40258c67cb58b2838fe5">2 Study Area</title>
      <p id="paragraph-9dee2d9e242e4e06b8b2f4e219b6d7fe"> This particular study area is located in Chikmangalur District of Karnataka, which represents a diverse and ecologically diverse region characterized by varying topograph , land cover and land uses. This encompasses an area of 1,073 Sq.Km. Chikmangalur Taluk is located amidst the Western Ghats. ﻿This taluk﻿ experiences a subtropical climate influenced by its elevation and proximity to Western Ghats, with varying temperature across different seasons and elevations. Given the landscape, land cover and land use, and various socio-economic activities, ﻿﻿this taluk provides an ideal </p>
      <fig id="f-9bc8af401588" orientation="portrait" fig-type="graphic" position="anchor">
        <label>Figure 0 </label>
        <graphic id="g-46ecb834dc0a" xlink:href="https://typeset-prod-media-server.s3.amazonaws.com/article_uploads/d67f4f0a-60b2-4d91-913b-938a23985cf0/image/1b9ed68b-b1da-4edd-99f6-ac09afed5010-upicture1.jpg"/>
      </fig>
      <p id="p-72579a7c8a8f"/>
      <sec>
        <title id="title-a0ed215fc1b04b16bfced317cb911b9f">2.1 Data</title>
        <p id="paragraph-811d5f5798704497a4027b2810787840">The foundation of this research is a robust dataset comprising a spectrum of geospatial information:</p>
        <p id="paragraph-5fa0d9e74eaf4094823615ad0b384c6f">Satellite Imagery and MTL File:</p>
        <p id="paragraph-8069f7a4dddd40c2a57259b029ecb10d">Our primary data source is the Landsat 8 Level 1 satellite imagery, offering multispectral data with exceptional spatial resolution. This dataset was meticulously collected over multiple time periods, enabling an in-depth exploration of temporal LST changes. Accompanying this imagery, the Metadata (MTL) file plays a pivotal role. It contains critical information about sensor calibration, sun angles, and other metadata essential for accurate radiometric and atmospheric correction. This metadata file ensures the scientific rigor and accuracy of the study, allowing us to transform raw digital numbers into radiance values and, subsequently, into brightness temperature. The MTL file and satellite imagery together constitute the bedrock of our geospatial analysis <xref id="xref-15a40af2b3da43e7b2bd642292f3ba1e" rid="R238433831159376" ref-type="bibr">2</xref>.</p>
      </sec>
    </sec>
    <sec>
      <title id="title-4ff9e3fa843143e5aa535d9d934c4f31">3 Methodology</title>
      <p id="paragraph-5d8f10daaa14473ea767acb1266a5ca3">Our methodology comprises a series of meticulously orchestrated steps, each designed to ensure the precision and reliability of our findings. </p>
      <fig id="figure-a6edc2379a004b3992a9ad863dbbf63a" orientation="portrait" fig-type="graphic" position="anchor">
        <label>Figure 1 </label>
        <caption id="caption-0863b10dd7a247c791a3ade8cb2deb11">
          <title id="title-f09504cf219b44d79c8a3af5cfe6da23">
            <bold id="s-5d3d3b4145ca">Methodology Chart</bold>
          </title>
        </caption>
        <graphic id="graphic-4893d54e5a70483b890ee27540a5d9f3" xlink:href="https://s3-us-west-2.amazonaws.com/typeset-prod-media-server/df6aa85c-71fc-4568-b65a-ddaa4735fed9image2.jpeg"/>
      </fig>
      <p id="paragraph-03fd65757aab42c3b14c1d6faf00caf4">The Landsat 8 imagery we collected underwent rigorous preprocessing, which included radiometric calibration, atmospheric correction, and geometric rectification. This preparatory phase guaranteed the dataset's readiness for subsequent analytical processes <xref id="xref-8e0cdf606ddf42f292c00790b1d73e8e" rid="R238433831159374" ref-type="bibr">3</xref>.</p>
      <sec>
        <title id="title-383ad64e96e04b4f82cd50952eef1c63">3.1 TOA Radiance Calculation</title>
        <p id="paragraph-0e7580febd1a48fe856629908deca59b">Following preprocessing, we calculated the Top-of-Atmosphere (TOA) radiance. This conversion endowed the dataset with radiometric precision, rendering it physically meaningful and amenable to temperature estimation.</p>
        <p id="paragraph-358b600bc9f54de4b3c1f946686cd47b">
          <bold id="strong-6fca193199aa4cf1ad9f36373108def3">TOA radiance, Lλ = ML x Qcal + AL-Qi (1)</bold>
        </p>
        <p id="paragraph-98072e0b49f24cbeac0273e47c23b22f">• ML = Specific multiplicative scaling factor of each band. Value obtained from    the MTL metadata file under the name of "RADIANCE_MULT_BAND_X". </p>
        <p id="paragraph-9b0207f1208f474fadc104b14d11751e">• Qcal = Is the band or the cut of it.</p>
        <p id="paragraph-40734b2f3a5443fe8a6d1018304f2154">• AL = Value included in the MTL metadata "Radiance_Add_Band_X", where X    corresponds to the number of the band</p>
      </sec>
      <sec>
        <title id="title-6c76cc79d1b24d85ad0f25e11e175398">3.2 Conversion to degrees (Brightness Temperature</title>
        <p id="paragraph-c75b05862b5c406dbbf7ece4a5bc6d48">The radiometric data was further transformed into degrees Kelvin (K), a vital conversion for accurate land surface temperature estimation. This transformation primed the dataset for precise temperature calculations.</p>
        <p id="paragraph-c004960ccd1c4ffb9fcdecb8adedd8f5">
          <bold id="strong-ddda51dfdedb4cafaf1f6420d08300c9"> </bold>
          <bold id="strong-bc755a4f1185457097034ee6bb6c4c6f">BT = [K2/ In (K1/ Lλ + 1)] - 273.15 (2)</bold>
        </p>
        <p id="paragraph-4fffebf1f80e4066980a672cc7d8cdfc">K1 and K2 = Conversion constants, included in the metadata (K1_CONSTANT_BAND_x and K2_CONSTANT_BAND_x) apply to each band, 10 and 11.</p>
      </sec>
      <sec>
        <title id="title-a1adb32f49e54d739723d946563c17c6">3.3 NDVI Calculation</title>
        <p id="paragraph-8182658c9fb045f186309b57f4e04eaa">We computed the Normalized Difference Vegetation Index (NDVI), a pivotal step in assessing regional vegetation health and land use patterns. The NDVI values, ranging from -1 to 1, played a critical role in evaluating green cover within the study area.</p>
        <p id="paragraph-bc3db5eda1bd47029af235a59f8da7cf">
          <bold id="strong-1bd316e2b03b49d4a7c09c8772dc8c58">NDVI = (Near infrared – Red)/ (Near infrared + Red) (3)</bold>
        </p>
      </sec>
      <sec>
        <title id="title-31bd9228d62d4883aa1024f48795527b">3.4 Proportion of Vegetation (PV</title>
        <p id="paragraph-ba71bdff13ed4932b049233e9044b0a2">Building upon the derived NDVI values, we assessed the proportion of vegetation (PV), quantifying the extent of vegetative cover across the study area. This metric provided valuable insights for our land surface temperature assessment.</p>
        <p id="paragraph-2de9bcd848564e4ba3cc9ecf68d7dc68">
          <bold id="strong-91875290f2f244b098609573fb347c50">Pv = Square ((NDVI - NDVImin) / (NDVImax - NDVImin)) (4)</bold>
        </p>
      </sec>
      <sec>
        <title id="title-c1dda98c36044fafaaa1987681f5b1fb">3.5 Land Surface Emissivity</title>
        <p id="paragraph-75f356bf37e44e929c741959f106e731">The determination of land surface emissivity values was a fundamental step, as these values significantly influenced temperature estimations. Emissivity values were based on surface properties and land cover types, ensuring the temperature estimates accurately reflected the region's unique characteristics. </p>
        <p id="paragraph-53eac51812024495aeb2c9a961025a0b">
          <bold id="strong-d4c1ad96d1b34802948951cd8ca4d831">e = m Pv + n </bold>
          <bold id="strong-bca35608f21549c599e8c17d843ae7cb">(5)</bold>
        </p>
        <p id="paragraph-6506f2d8254d458380466f9be171f1f5">• m = value of emissivity of vegetation, in this case 0.004 was used </p>
        <p id="paragraph-5c059c2468414175b99f3416ec93c3d5">• Pv = corresponds to the percentage of vegetation </p>
        <p id="paragraph-04def79687614643883c24fc5ab081c5">• n = Soil emissivity value, in this case 0.986 was used </p>
        <p id="paragraph-32ae3569aea24739af2ed7d4a330af98">Therefore, <bold id="strong-9274bab2ba994dda8ef8f44822a6c68a"> LSE = 0.004 * Pv + 0.986 (7)</bold></p>
      </sec>
      <sec>
        <title id="title-c2da5d02f4674a869df94c8b88c98616">3.6 Land Surface Temperature (LST estimation</title>
        <p id="paragraph-695d22920c9448dcb3b95461b11a98a7">The pinnacle of our methodology was the estimation of land surface temperature (LST). This process involved the application of the Radiative Transfer Equation (RTE), which integrated brightness temperature, emissivity values, and atmospheric properties. The final output was a spatial distribution of LST values across Chikmagalur Taluk <xref id="xref-d2f8fb2ae45a41e3b80eab8cac8ca086" rid="R238433831159375" ref-type="bibr">4</xref>.</p>
        <p id="paragraph-0c546fdd9e664708b7028190bd97d0b5">
          <bold id="strong-7d2c104f24f14f2ca66f702137c4d364">LST = BT / 1 + w (BT / p) * Ln (ε)  (8)</bold>
        </p>
        <p id="paragraph-d0c600fe6ff54c2c9bb8d56b3db6c92a">• BT = Brightness temperature</p>
        <p id="paragraph-3e9826cc4ec04a1f85ee3cf4f8e57b5c">• w = Length of the emitted radiation (band 10 or 11 as the case may be) </p>
        <p id="paragraph-9bf605c8108a434282fbbd08c3cb752d">• p = Constant value obtained by the formula h * c / s that when substituting the values is 1.438 * 10^ -34Js and results in 14,380 </p>
        <p id="paragraph-6708c725930547dfb1443b5464eaf355">• ε is LSE</p>
        <fig id="f-d70d6914cb6d" orientation="portrait" fig-type="graphic" position="anchor">
          <label>Figure 2 </label>
          <caption id="c-9e17353aabd3">
            <title id="t-9870448ec8df">
              <bold id="s-6f0868cebed7">LST of 2015 Chikamangalur Taluk</bold>
            </title>
          </caption>
          <graphic id="g-8f2ec4a9cdab" xlink:href="https://typeset-prod-media-server.s3.amazonaws.com/article_uploads/d67f4f0a-60b2-4d91-913b-938a23985cf0/image/40a34369-73a7-4884-9e96-c7573db1aeb5-upicture1.jpg"/>
        </fig>
        <p id="p-73e27c52f1ef"/>
        <fig id="f-94f47a8fb23d" orientation="portrait" fig-type="graphic" position="anchor">
          <label>Figure 3 </label>
          <caption id="c-d153ba85de55">
            <title id="t-4ceb80fb7bbd">
              <bold id="s-af166b23f9cc">LST of 2023 Chikamangalur Taluk</bold>
            </title>
          </caption>
          <graphic id="g-6ac43b450f16" xlink:href="https://typeset-prod-media-server.s3.amazonaws.com/article_uploads/d67f4f0a-60b2-4d91-913b-938a23985cf0/image/c7df20f0-f5a5-4f58-96c8-425d68fa3c03-upicture2.jpg"/>
        </fig>
        <p id="p-ae9143e66a06"/>
      </sec>
    </sec>
    <sec>
      <title id="title-ce7636bda5a84786957bfa02db6c8e5d">4 Results and conclusions</title>
      <p id="paragraph-573635fbb414403e8a0aa033b6aa6186">Our meticulous geospatial analysis revealed a diverse pattern of LST changes in Chikmagalur Taluk. Notably, the eastern part of the Taluk exhibited localized decreases in temperature, while the hill regions in the west and central areas experienced significant temperature increases.</p>
      <p id="paragraph-68551dc90fe44fb7bbdf8fdb2aa0a7f6">These temperature increases, particularly in the hill regions, corresponded with recent environmental incidents in Chikmagalur. Reports of deforestation, increased land clearing for agriculture, urban infrastructure projects, and a surge in tourism have raised concerns. These incidents may be contributing to localized urban heat island effects and the observed temperature changes. While our analysis doesn't establish causation, it underscores the need for further investigations to comprehensively understand these temperature variations and their potential implications for the region <xref id="xref-6595b3ddf7cb48288e269ca0229a2fdd" rid="R238433831159377" ref-type="bibr">5</xref>.</p>
      <fig id="f-8f5d79fe3f13" orientation="portrait" fig-type="graphic" position="anchor">
        <label>Figure 4 </label>
        <caption id="c-474c4739aed9">
          <title id="t-cf98f864ae5e">
            <bold id="s-fcd9d20ec170">Change Detection of LST</bold>
          </title>
        </caption>
        <graphic id="g-93c7a71add46" xlink:href="https://typeset-prod-media-server.s3.amazonaws.com/article_uploads/d67f4f0a-60b2-4d91-913b-938a23985cf0/image/7592ee87-d925-4def-a7be-5e915682fadf-upicture3.jpg"/>
      </fig>
      <p id="p-9983379ac2d0"/>
    </sec>
  </body>
  <back>
    <ref-list>
      <title>References</title>
      <ref id="R238433831159379">
        <element-citation publication-type="journal">
          <person-group person-group-type="author">
            <name>
              <surname>Sobrino</surname>
              <given-names>José A.</given-names>
            </name>
            <name>
              <surname>Jiménez-Muñoz</surname>
              <given-names>Juan C.</given-names>
            </name>
            <name>
              <surname>Paolini</surname>
              <given-names>Leonardo</given-names>
            </name>
            <collab/>
          </person-group>
          <article-title>Land surface temperature retrieval from LANDSAT TM 5</article-title>
          <source>Remote Sensing of Environment</source>
          <year>2004</year>
          <volume>90</volume>
          <issue>4</issue>
          <fpage>434</fpage>
          <lpage>440</lpage>
          <issn>0034-4257</issn>
          <publisher-name>Elsevier BV</publisher-name>
          <uri>https://dx.doi.org/10.1016/j.rse.2004.02.003</uri>
        </element-citation>
      </ref>
      <ref id="R238433831159376">
        <element-citation publication-type="journal">
          <person-group person-group-type="author">
            <name>
              <surname/>
              <given-names>T. V Ramachandra</given-names>
            </name>
            <name>
              <surname>Aithal</surname>
              <given-names>B H</given-names>
            </name>
            <name>
              <surname>Sanna</surname>
              <given-names>D D</given-names>
            </name>
            <collab/>
          </person-group>
          <article-title>Land Surface Temperature Analysis in an Urbanising Landscape through Multi-Resolution Data</article-title>
          <source>Journal of Space Science &amp; Technology</source>
          <year>2012</year>
          <volume>1</volume>
          <issue>1</issue>
          <fpage>1</fpage>
          <lpage>10</lpage>
          <uri>https://wgbis.ces.iisc.ac.in/energy/paper/stmj_lsta/land_surface_temp_analysys.pdf</uri>
        </element-citation>
      </ref>
      <ref id="R238433831159374">
        <element-citation publication-type="inproceedings">
          <person-group person-group-type="author">
            <name>
              <surname>Akhoondzadeh</surname>
              <given-names>M</given-names>
            </name>
            <name>
              <surname>Saradjian</surname>
              <given-names>M R</given-names>
            </name>
            <collab/>
          </person-group>
          <person-group person-group-type="editor">
            <etal/>
          </person-group>
          <article-title>COMPARISON OF LAND SURFACE TEMPERATURE MAPPING USING MODIS AND ASTER IMAGES IN SEMI-ARID AREA</article-title>
          <year>2008</year>
          <publisher-name>The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences</publisher-name>
          <publisher-loc>Beijing</publisher-loc>
          <uri>https://www.isprs.org/proceedings/xxxvii/congress/8_pdf/9_wg-viii-9/03a.pdf</uri>
        </element-citation>
      </ref>
      <ref id="R238433831159375">
        <element-citation publication-type="journal">
          <person-group person-group-type="author">
            <name>
              <surname>Garouani</surname>
              <given-names>Manal El</given-names>
            </name>
            <name>
              <surname>Amyay</surname>
              <given-names>Mhamed</given-names>
            </name>
            <name>
              <surname>Lahrach</surname>
              <given-names>Abderrahim</given-names>
            </name>
            <name>
              <surname>Oulidi</surname>
              <given-names>Hassane Jarar</given-names>
            </name>
            <collab/>
          </person-group>
          <article-title>Land Surface Temperature in Response to Land Use/Cover Change Based on Remote Sensing Data and GIS Techniques: Application to Saïss Plain, Morocco</article-title>
          <source>Journal of Ecological Engineering</source>
          <year>2021</year>
          <volume>22</volume>
          <issue>7</issue>
          <fpage>100</fpage>
          <lpage>112</lpage>
          <issn>2299-8993</issn>
          <publisher-name>Wydawnictwo Naukowe Gabriel Borowski (WNGB)</publisher-name>
          <uri>https://dx.doi.org/10.12911/22998993/139065</uri>
        </element-citation>
      </ref>
      <ref id="R238433831159377">
        <element-citation publication-type="journal">
          <person-group person-group-type="author">
            <name>
              <surname>Twumasi</surname>
              <given-names>Yaw A.</given-names>
            </name>
            <name>
              <surname>Merem</surname>
              <given-names>Edmund C.</given-names>
            </name>
            <name>
              <surname>Namwamba</surname>
              <given-names>John B.</given-names>
            </name>
            <name>
              <surname>Mwakimi</surname>
              <given-names>Olipa S.</given-names>
            </name>
            <name>
              <surname>Ayala-Silva</surname>
              <given-names>Tomas</given-names>
            </name>
            <name>
              <surname>Frimpong</surname>
              <given-names>Diana B.</given-names>
            </name>
            <name>
              <surname>Ning</surname>
              <given-names>Zhu H.</given-names>
            </name>
            <name>
              <surname>Asare-Ansah</surname>
              <given-names>Abena B.</given-names>
            </name>
            <name>
              <surname>Annan</surname>
              <given-names>Jacob B.</given-names>
            </name>
            <name>
              <surname>Oppong</surname>
              <given-names>Judith</given-names>
            </name>
            <name>
              <surname>Loh</surname>
              <given-names>Priscilla M.</given-names>
            </name>
            <name>
              <surname>Owusu</surname>
              <given-names>Faustina</given-names>
            </name>
            <name>
              <surname>Jeruto</surname>
              <given-names>Valentine</given-names>
            </name>
            <name>
              <surname>Petja</surname>
              <given-names>Brilliant M.</given-names>
            </name>
            <name>
              <surname>Okwemba</surname>
              <given-names>Ronald</given-names>
            </name>
            <name>
              <surname>McClendon-Peralta</surname>
              <given-names>Joyce</given-names>
            </name>
            <name>
              <surname>Akinrinwoye</surname>
              <given-names>Caroline O.</given-names>
            </name>
            <name>
              <surname>Mosby</surname>
              <given-names>Hermeshia J.</given-names>
            </name>
            <collab/>
          </person-group>
          <article-title>Estimation of Land Surface Temperature from Landsat-8 OLI Thermal Infrared Satellite Data. A Comparative Analysis of Two Cities in Ghana</article-title>
          <source>Advances in Remote Sensing</source>
          <year>2021</year>
          <volume>10</volume>
          <issue>04</issue>
          <fpage>131</fpage>
          <lpage>149</lpage>
          <issn>2169-267X, 2169-2688</issn>
          <publisher-name>Scientific Research Publishing, Inc.</publisher-name>
          <uri>https://dx.doi.org/10.4236/ars.2021.104009</uri>
        </element-citation>
      </ref>
    </ref-list>
  </back>
</article>
