<!DOCTYPE article PUBLIC "-//NLM//DTD JATS (Z39.96) Journal Archiving and Interchange DTD v1.0 20120330//EN" "JATS-archivearticle1.dtd">
<article xmlns:xlink="http://www.w3.org/1999/xlink">
  <front>
    <journal-meta />
    <article-meta>
      <title-group>
        <article-title>Algorithmic Fairness in Geo-intelligence Workflows through Causality</article-title>
      </title-group>
      <contrib-group>
        <contrib contrib-type="author">
          <string-name>Brian K. Masinde</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Caroline M. Gevaert</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Michael H. Nagenborg</string-name>
          <xref ref-type="aff" rid="aff0">0</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Marc van den Homberg</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <contrib contrib-type="author">
          <string-name>Jaap A. Zevenbergen</string-name>
          <xref ref-type="aff" rid="aff1">1</xref>
        </contrib>
        <aff id="aff0">
          <label>0</label>
          <institution>Department of Philosophy, University of Twente</institution>
          ,
          <addr-line>Enschede</addr-line>
          ,
          <country country="NL">The Netherlands</country>
        </aff>
        <aff id="aff1">
          <label>1</label>
          <institution>Faculty of Geo-information Science and Earth Observation, University of Twente</institution>
          ,
          <addr-line>Enschede</addr-line>
          ,
          <country country="NL">The Netherlands</country>
        </aff>
      </contrib-group>
      <abstract>
        <p>In this paper, we investigate how causality (causal inference) can be used to detect bias and ensure fairness in geo-intelligence workflows . We investigate the usefulness of such a causality-based approach in the context of an early warning system that predicts building damage at municipality levels in The Philippines. We use directed acyclic graphs to reason about the causal relationships in the model case study and quantify the relationships using structural equation modelling. Mediation analysis is also used to validate the causal relationships between variables. We find cases of confounder bias and Simpsons paradox that could potentially bias the damage predictions. However we note that the objective and outcome variable in the early warning system needs to be defined in a manner that allows for more nuanced investigation on fairness (i.e., from damage assessment to impact assessment).</p>
      </abstract>
      <kwd-group>
        <kwd>eol&gt;Biases</kwd>
        <kwd>algorithmic fairness</kwd>
        <kwd>geo-intelligence</kwd>
        <kwd>disaster early warning systems</kwd>
      </kwd-group>
    </article-meta>
  </front>
  <body>
    <sec id="sec-1">
      <title>1. Introduction</title>
      <p>fairness and transparency as data biases and technical limitations can afect the reliability and
trustworthiness of geo-intelligence workflows.</p>
      <p>
        In this paper, we investigate how causality (causal inference) can be used to promote fairness
by detecting bias in the data and design of geo-intelligence workflows. As a case/study we
consider an early warning system that predicts building damage by tropical cyclones (TC) at
municipality levels in The Philippines. Fairness is a pertinent issue in disaster
response/management because of the limited resources and the lives and livelihoods at stake. We explore how
causality can be used for algorithmic fairness since non-causal (associational) models have been
shown to replicate, reinforce and propagate biases in observational data. Furthermore, some
biases and statistical anomalies cannot be addressed using association/correlation language
used in non-causal models (e.g., confounder bias [Figure 3] and the Simpson’s paradox) [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ].
      </p>
      <p>
        We follow Friedman and Nissenbaum’s [
        <xref ref-type="bibr" rid="ref5">5</xref>
        ] definition of a bias in computer systems. By their
definition, a bias is a systematic, and unfair discrimination of individuals or groups of individuals.
Friedman and Nissenbaum develop a framework of understanding bias in computer systems
as either preexisting bias (e.g., data generated from biased society or institutions), technical
bias (technical limitations of the algorithms) and or emergent bias (i.e., biases changing real
world use such as changing societal values). In the context of our case study, the data does not
contain the quintessential sensitive variables (e.g., gender, race/ethnicity). However, the case
study uses building typologies as predictor variables which are often considered an indicator
of socio-economic well-being. Therefore a biased damage assessment in this context would
consequently cause unfair distribution of aid resources to municipalities with a high number of
building typologies vulnerable to tropical cyclones.
      </p>
      <p>
        Although there is literature on using causality for quantifying disaster risk (e.g., [
        <xref ref-type="bibr" rid="ref6 ref7">6, 7</xref>
        ]), to
the best of our knowledge there are none so far on assessing biases and fairness using causality
in this context. Though our case study is on aggregated data and does not contain the typical
problematic sensitive variables (e.g., gender, ethnicity), it is still important to ensure that there
are no biases in the early warning system and that it is fair. In this paper we present preliminary
results on detecting biases through causality.
      </p>
    </sec>
    <sec id="sec-2">
      <title>2. The Case Study and Data</title>
      <p>
        The Philippines is prone to tropical cyclones (characterized by high speed winds and heavy
rainfall) which cause loss of lives and infrastructure damage. While there are accurate models
that predict the characteristics of tropical cyclones (e.g., wind speed and rainfall) it is still
a challenge to predict and quantify the impact it would have on people. Because of
socioeconomic inequalities that influence choice of building materials, often this impact is measured
by assessing building damage. Because of these challenges (high frequency of tropical cyclone
events and cascading efects), geo-intelligence workflows are increasingly leveraged to guide
Impact-based Forecasting (IbF) that in turn inform disaster anticipatory actions such as
Forecastbased Financing (FbF) [
        <xref ref-type="bibr" rid="ref8">8</xref>
        ]. 510, an initiative of The Netherlands Red Cross developed an early
warning system for The Philippines. Their model and data serves as our case study.
      </p>
      <p>The observation data is from previous 39 tropical cyclone incidents in The Philippines. The
data is aggregated at the municipality level (1486 municipalities), therefore the variables are an
aggregate representation of the hazard efects in each municipality. Hazard variables include
wind speed and amount of rainfall and physical vulnerability variables are captured by house
typology numbers (grouped by roof type and wall type). Geographical variables are also included
(e.g., coastal length). The resulting early warning system under study is a regression model that
predicts the damage percentage in municipalities identifying the most impacted (see Table 1 for
metrics).</p>
    </sec>
    <sec id="sec-3">
      <title>3. Methods: Causality</title>
      <p>
        Causal inference is a methodology that aims at quantifying cause and efect from observational
data. Causal inference has found applications in algorithmic fairness [e.g., 9, 10, 11] and in
explainable AI [
        <xref ref-type="bibr" rid="ref12">12</xref>
        ]. There are diferent methods on causal inference such as Bayesian Networks
and Structural Equation Models (SEMs). BNs and SEMs are popular because of the ease of
representing causal relationships using graphs, specifically directed acylic graphs (DAGs) which
do not have feedback loops [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ]. In this paper we use SEMs (Equation 1) because they can also
be used in mediation analysis to estimate and validate causal chain graphs (Figure 2).
      </p>
      <p>= (, ),  = 1, . . . ,</p>
      <p>Where  is a function model (e.g., linear regression models) with  as parent nodes of 
and  a random error.</p>
      <p>
        Mediation analysis aims to quantify both the direct efects and indirect efects of an exposure
variable to an outcome. This is done in two steps, one regressing on the outcome variable on the
exposure variable (Equation 2) and two regressing the outcome on both exposure and mediator
variable (Equation 3). The diference of the exposure coeficients in the two model gives the
efects of the exposure (E) interfered by mediator variable (M) (Equation 4) [
        <xref ref-type="bibr" rid="ref13">13</xref>
        ].
      </p>
      <p>=  0 +   ×  + 
 =  0 +  1 ×  +  2 ×  +</p>
      <p>=   −  1</p>
      <p>
        To detect potential biases through causality in the early warning system we first define a DAG
informed by literature (e.g., [
        <xref ref-type="bibr" rid="ref14">14</xref>
        ]) of how tropical cyclones cause damage. This is however limited
by the available variables in the data. We use correlation analysis to identify relationships,
mediation analysis to confirm chain relationships (Figure 2), and linear regression models to
identify cases of Simpsons paradox. Importantly, the DAG help identify potential confounder
biases. Secondly, we conduct a path analysis, quantifying the efects of variables on the outcome.
We use SEMs with the functions specified as linear regression models to quantify the direct
efects among the variables. These steps are a precursor to conducting further causal fairness
checks, for example path-specific counterfactual fairness as illustrated by Chiappa [
        <xref ref-type="bibr" rid="ref15">15</xref>
        ].
(1)
(2)
(3)
(4)
      </p>
    </sec>
    <sec id="sec-4">
      <title>4. Preliminary results</title>
      <p>Figure 1 shows our causal representation of the tropical cyclones early warning system. It is
informed by literature of how tropical cyclones cause infrastructure damage [e.g., 14] and this
damage varies because of the building typologies1. With many variables to consider, we focus
on the ones that have a strong correlation with the outcome variable and among each other.
Figure 5 shows the correlation between the variables. Furthermore, departure from the initial
associational ML design to a causal model leads to removing variables that do not have a causal
link (or interpretation) to the other variables. For example, the number of households registered
for social protection can be associated with higher damage levels but there’s no causal link/path
between this variable and the outcome.</p>
      <p>The DAG is based on the causal reasoning that natural hazards have a direct causal efect on
damage of buildings. The hazard components - wind speed and total rainfall2 - have a strong
correlation (see Figure 5H) and statistical significance to the outcome (damage). Since every
tropical cyclone event is diferent and there might be other characteristics not captured in the
data, we consider the event itself as a parent node to wind speed and total rainfall. The event
and municipalities are parent nodes of distance which causes the municipalities to experience
diferent levels of wind speed and rainfall. For example, Table 6 shows that wind speed is indeed
a significant mediator between distance and damage. Tables 4 and 3 show that the significance
of rainfall as a predictor of damage changes depending on whether or not we account for the
specific tropical cyclone events in the model. This is therefore a case of the Simpsons paradox.
Furthermore, here we consider distance to express other hazard efects on damage and therefore
it becomes a confounder variable. The hazard components may not be entirely representative of
tropical cyclones since it does not also consider storm surges as a side efect of the wind speed.</p>
      <p>
        For simplicity, we do not add paths between hazard nodes and build typologies nodes because
this introduces a time component. That is over time communities adjust to impact tropical
cyclones rebuilding damaged buildings with stronger materials. Mediation analysis however
confirms that build typologies can be mediators of efect of hazard on damage levels in our data
(see Table 5). Furthermore, the relationship between build typologies could be bidirectional but
doing so in our model introduces cycles and structural causal models do not account for cyclic
relationships [
        <xref ref-type="bibr" rid="ref4">4</xref>
        ].
      </p>
    </sec>
    <sec id="sec-5">
      <title>5. Discussion</title>
      <p>In this causal re-interpretation of the efect of tropical cyclones on damage levels in
municipalities, based on our data, we find the problems of confounder bias (distance confounds both
hazard components and the outcome damage levels). In addition, we find a case of Simpsons
paradox with the direct efect of rainfall on damage. Following the correct causal path ways to
establish efect of hazard on damage plays an important role in ensuring fairness. Path analysis
is a precursor to carrying out further fairness tests (e.g., counterfactual analysis). In this case
the weaker build typologies (e.g., light roof, light walls) are sometimes proxies for vulnerable</p>
      <sec id="sec-5-1">
        <title>1Stronger building materials are more resilient to natural hazards 2We exclude rainfall measurements at 6 hours and 24 hours because of high correlation to total rainfall, see Figure 5</title>
        <p>socio-economic groups of people and hence can be considered as sensitive variables. Because
of the identified complexities (i.e., confounder bias and Simpsons paradox) further analysis is
required to quantify the efect of these biases on the accuracy of damage prediction and how
they would afect fair distribution of aid resources. It is however apparent that the outcome
variable and the objectives also need to be designed or specified in a fair manner that ensures
fairness. In this context we note that predicting impact instead of damage would allow for a
deeper engagement on fairness.</p>
      </sec>
    </sec>
    <sec id="sec-6">
      <title>6. Acknowledgments</title>
      <p>Authors received funding from the Netherlands Organization for Scientific Research
(NWOMVI) and UNICEF (grant number MVI.19.007). Authors also benefit from collaborations with
510, an Initiative of the Netherlands Red Cross.</p>
      <p>Authors also acknowledge Anna Manchens for comments on mediation reasoning.</p>
    </sec>
    <sec id="sec-7">
      <title>A. Basic DAGs</title>
      <sec id="sec-7-1">
        <title>Chain graph conditional independence.</title>
        <p>⊥ |
(5)</p>
      </sec>
    </sec>
    <sec id="sec-8">
      <title>C. Tables</title>
      <p>Mean Absolute Error (MAE)
Root Mean Square Error (RMSE)</p>
      <p>Description
Percentage of damage buildings observed in the previous typhoons
Hazard - wind speed (/)
Hazard - total rainfall () during the tropical cyclone
Geography - Distance from the epicentre of the tropical cyclone
Percentage of strong roof and strong wall build types
Percentage of strong roof and salvage wall build types
Percentage of strong roof and light wall build types
Percentage of salvage roof and strong wall build types
Percentage of salvage roof and salvage wall build types
Percentage of salvage roof and light wall build type
Percentage of light roof and strong wall build type
Percentage of light roof and salvage wall build type</p>
      <p>Percentage of light roof and light wall build types
Estimate
0.005
0.235
0.240
0.021
Average Causal Mediation Efects(ACME)
Average Direct Efects (ADE)
Total Efect
Prop. Mediated</p>
    </sec>
  </body>
  <back>
    <ref-list>
      <ref id="ref1">
        <mixed-citation>
          [1]
          <string-name>
            <given-names>S.</given-names>
            <surname>Boeke</surname>
          </string-name>
          , M. van den Homberg, A. Teklesadik,
          <string-name>
            <given-names>J.</given-names>
            <surname>Fabila</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D.</given-names>
            <surname>Riquet</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Alimardani</surname>
          </string-name>
          ,
          <article-title>Towards predicting rice loss due to typhoons in the philippines</article-title>
          ,
          <source>The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences</source>
          <volume>42</volume>
          (
          <year>2019</year>
          )
          <fpage>63</fpage>
          -
          <lpage>70</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref2">
        <mixed-citation>
          [2]
          <string-name>
            <surname>M. J. Van den Homberg</surname>
          </string-name>
          ,
          <string-name>
            <surname>C. M. Gevaert</surname>
            ,
            <given-names>Y. Georgiadou,</given-names>
          </string-name>
          <article-title>The changing face of accountability in humanitarianism: Using artificial intelligence for anticipatory action</article-title>
          ,
          <source>Politics and Governance</source>
          <volume>8</volume>
          (
          <year>2020</year>
          )
          <fpage>456</fpage>
          -
          <lpage>467</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref3">
        <mixed-citation>
          [3]
          <string-name>
            <given-names>C. M.</given-names>
            <surname>Gevaert</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M.</given-names>
            <surname>Carman</surname>
          </string-name>
          ,
          <string-name>
            <given-names>B.</given-names>
            <surname>Rosman</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Georgiadou</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Soden</surname>
          </string-name>
          ,
          <article-title>Fairness and accountability of ai in disaster risk management: Opportunities and challenges</article-title>
          ,
          <source>Patterns</source>
          <volume>2</volume>
          (
          <year>2021</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref4">
        <mixed-citation>
          [4]
          <string-name>
            <given-names>J.</given-names>
            <surname>Pearl</surname>
          </string-name>
          , Causality: Models,
          <string-name>
            <surname>Reasoning</surname>
          </string-name>
          , and Inference, 2nd ed., Cambridge Univeristy Press,
          <year>2009</year>
          .
        </mixed-citation>
      </ref>
      <ref id="ref5">
        <mixed-citation>
          [5]
          <string-name>
            <given-names>B.</given-names>
            <surname>Friedman</surname>
          </string-name>
          ,
          <string-name>
            <given-names>H.</given-names>
            <surname>Nissenbaum</surname>
          </string-name>
          ,
          <article-title>Bias in computer systems</article-title>
          ,
          <source>ACM Transactions on information systems (TOIS) 14</source>
          (
          <year>1996</year>
          )
          <fpage>330</fpage>
          -
          <lpage>347</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref6">
        <mixed-citation>
          [6]
          <string-name>
            <given-names>S.</given-names>
            <surname>Xu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Dimasaka</surname>
          </string-name>
          ,
          <string-name>
            <given-names>D. J.</given-names>
            <surname>Wald</surname>
          </string-name>
          ,
          <string-name>
            <given-names>H. Y.</given-names>
            <surname>Noh</surname>
          </string-name>
          ,
          <article-title>Seismic multi-hazard and impact estimation via causal inference from satellite imagery</article-title>
          ,
          <source>Nature Communications</source>
          <volume>13</volume>
          (
          <year>2022</year>
          )
          <fpage>7793</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref7">
        <mixed-citation>
          [7]
          <string-name>
            <given-names>H.</given-names>
            <surname>Burton</surname>
          </string-name>
          ,
          <article-title>Causal inference on observational data: Opportunities and challenges in earthquake engineering</article-title>
          ,
          <source>Earthquake Spectra</source>
          <volume>39</volume>
          (
          <year>2023</year>
          )
          <fpage>54</fpage>
          -
          <lpage>76</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref8">
        <mixed-citation>
          [8]
          <string-name>
            <given-names>A.</given-names>
            <surname>Teklesadik</surname>
          </string-name>
          , M. van den Homberg,
          <article-title>Forecasting impacts of tropical cyclones with machine learning: A case study in the philippines</article-title>
          , in: EGU General Assembly Conference Abstracts,
          <year>2022</year>
          , pp.
          <fpage>EGU22</fpage>
          -
          <lpage>12917</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref9">
        <mixed-citation>
          [9]
          <string-name>
            <given-names>X.</given-names>
            <surname>Wang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Y.</given-names>
            <surname>Zhang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Zhu</surname>
          </string-name>
          ,
          <article-title>A brief review on algorithmic fairness</article-title>
          ,
          <source>Management System Engineering</source>
          <volume>1</volume>
          (
          <year>2022</year>
          )
          <article-title>7</article-title>
          .
        </mixed-citation>
      </ref>
      <ref id="ref10">
        <mixed-citation>
          [10]
          <string-name>
            <given-names>J. R.</given-names>
            <surname>Loftus</surname>
          </string-name>
          ,
          <string-name>
            <given-names>C.</given-names>
            <surname>Russell</surname>
          </string-name>
          ,
          <string-name>
            <given-names>M. J.</given-names>
            <surname>Kusner</surname>
          </string-name>
          ,
          <string-name>
            <given-names>R.</given-names>
            <surname>Silva</surname>
          </string-name>
          ,
          <article-title>Causal reasoning for algorithmic fairness</article-title>
          , arXiv preprint arXiv:
          <year>1805</year>
          .
          <volume>05859</volume>
          (
          <year>2018</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref11">
        <mixed-citation>
          [11]
          <string-name>
            <given-names>C.</given-names>
            <surname>Su</surname>
          </string-name>
          ,
          <string-name>
            <given-names>G.</given-names>
            <surname>Yu</surname>
          </string-name>
          ,
          <string-name>
            <given-names>J.</given-names>
            <surname>Wang</surname>
          </string-name>
          ,
          <string-name>
            <given-names>Z.</given-names>
            <surname>Yan</surname>
          </string-name>
          ,
          <string-name>
            <given-names>L.</given-names>
            <surname>Cui</surname>
          </string-name>
          ,
          <article-title>A review of causality-based fairness machine learning (</article-title>
          <year>2022</year>
          ).
        </mixed-citation>
      </ref>
      <ref id="ref12">
        <mixed-citation>
          [12]
          <string-name>
            <surname>Y.-L. Chou</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          <string-name>
            <surname>Moreira</surname>
            ,
            <given-names>P.</given-names>
          </string-name>
          <string-name>
            <surname>Bruza</surname>
            ,
            <given-names>C.</given-names>
          </string-name>
          <string-name>
            <surname>Ouyang</surname>
            ,
            <given-names>J.</given-names>
          </string-name>
          <string-name>
            <surname>Jorge</surname>
          </string-name>
          ,
          <article-title>Counterfactuals and causability in explainable artificial intelligence: Theory, algorithms</article-title>
          , and applications,
          <source>Information Fusion</source>
          <volume>81</volume>
          (
          <year>2022</year>
          )
          <fpage>59</fpage>
          -
          <lpage>83</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref13">
        <mixed-citation>
          [13]
          <string-name>
            <given-names>S. J.</given-names>
            <surname>Jung</surname>
          </string-name>
          ,
          <article-title>Introduction to mediation analysis and examples of its application to real-world data</article-title>
          ,
          <source>Journal of Preventive Medicine and Public Health</source>
          <volume>54</volume>
          (
          <year>2021</year>
          )
          <fpage>166</fpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref14">
        <mixed-citation>
          [14]
          <string-name>
            <surname>U. M. K. Eidsvig</surname>
            ,
            <given-names>K.</given-names>
          </string-name>
          <string-name>
            <surname>Kristensen</surname>
            ,
            <given-names>B. V.</given-names>
          </string-name>
          <string-name>
            <surname>Vangelsten</surname>
          </string-name>
          ,
          <article-title>Assessing the risk posed by natural hazards to infrastructures</article-title>
          ,
          <source>Natural Hazards and Earth System Sciences</source>
          <volume>17</volume>
          (
          <year>2017</year>
          )
          <fpage>481</fpage>
          -
          <lpage>504</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref15">
        <mixed-citation>
          [15]
          <string-name>
            <given-names>S.</given-names>
            <surname>Chiappa</surname>
          </string-name>
          ,
          <article-title>Path-specific counterfactual fairness</article-title>
          ,
          <source>in: Proceedings of the AAAI conference on artificial intelligence</source>
          , volume
          <volume>33</volume>
          ,
          <year>2019</year>
          , pp.
          <fpage>7801</fpage>
          -
          <lpage>7808</lpage>
          .
        </mixed-citation>
      </ref>
      <ref id="ref16">
        <mixed-citation>
          <string-name>
            <surname>P-value &lt;</surname>
          </string-name>
          2e-
          <fpage>16</fpage>
          *** &lt;
          <fpage>2e</fpage>
          -
          <lpage>16</lpage>
          *** &lt;
          <fpage>2e</fpage>
          -
          <lpage>16</lpage>
          *** &lt;
          <fpage>2e</fpage>
          -
          <lpage>16</lpage>
          ***
        </mixed-citation>
      </ref>
    </ref-list>
  </back>
</article>