EconBase
← Back to paper

Social and Genetic Ties Drive Skewed Cross-Border Media Coverage of Disasters

Extracted main text — title through conclusion, appendix excluded. This is what our citation measures are computed over, published so the extraction can be checked by eye.

80,626 characters · 15 sections · 64 citation commands

Rendered from LaTeX for readability, not typeset faithfully. Citation keys are highlighted; maths is left as source; figures, tables and equation environments are summarised rather than reproduced; unrecognised commands are greyed out so nothing is silently dropped. Email addresses are removed.

Social and Genetic Ties Drive Skewed Cross-Border Media Coverage of Disasters

abstractClimate change is increasing the frequency and severity of natural disasters worldwide. Media coverage of these events may be vital to generate empathy and mobilize global populations to address the common threat posed by climate change. Using a dataset of 466 news sources from 123 countries, covering 135 million news articles since 2016, we apply an event study framework to measure cross-border media activity following natural disasters. Our results shows that while media attention rises after disasters, it is heavily skewed towards certain events, notably earthquakes, accidents, and wildfires. In contrast, climatologically salient events such as floods, droughts, or extreme temperatures receive less coverage. This cross-border disaster reporting is strongly related to the number of deaths associated with the event, especially when the affected populations share strong social ties or genetic similarities with those in the reporting country. Achieving more balanced media coverage across different types of natural disasters may be essential to counteract skewed perceptions. Further, fostering closer social connections between countries may enhance empathy and mobilize the resources necessary to confront the global threat of climate change.

Introduction

Climate change is a pressing global risk that affects all of humanity, both directly and indirectly, often through complex and unforeseen mechanisms IPCC2018, Stern2007Economics. As climate change intensifies, natural disasters become more frequent and severe, posing significant threats to societies worldwide Field2012Managing, Coumou2012Decade, IPCC2021Climate. Yet, traditional media systems frequently appear “blind” to climate-induced disasters that are less visually dramatic, such as floods, storms, and extreme temperatures, focusing instead on events like earthquakes and wildfires Boykoff2009WeSpeakForTheTrees, Doulton2009TenYears.

Accurate and comprehensive media coverage is essential for raising public awareness of climate change and encouraging global cooperation Boykoff2007Climate, Painter2013Climate, Schafer2015Climate. Media reporting helps shape perceptions of climate-related risks and can influence public support for mitigation and adaptation policies Pralle2009AgendaSetting, McCombs2004SettingAgenda, Gavin2009Climate. However, the media often devotes disproportionate attention to certain events or regions, potentially skewing public perceptions Soroka2019BadNews, Zhuravskaya2020PoliticalEffects, Prat2011MediaInfluence.

Global action on climate change requires coordination and cooperation across borders Hulme2009WhyWeDisagree, Beck2009WorldRisk. Media coverage not only shapes domestic opinions but also influences global policy debates and social norms, making transnational perspectives essential for fully understanding the interplay between climate issues and public discourse schafer2018media, entman2009projections. Empathy—encouraged by comprehensive media reporting that clearly links climate-related disasters and their human consequences to their underlying causes—can play a key role in building global solidarity Batson2011Altruism, Hoffman2001Empathy, deWaal2009AgeOfEmpathy. Empathy is considered a pinnacle of human moral development and helps motivate prosocial behavior, moral responsibility, and cooperative actions Hoffman2001Empathy, Tomasello2019BecomingHuman, Decety2014ComplexRelation, Goleman1995Emotional. Media reporting that potentially enhances empathy towards those affected by climate change-induced disasters could be important for fostering global cooperation.

Despite the increasing prevalence of climate change-induced natural disasters, and the importance of empathy and global cooperation, there is limited research on cross-border media reporting of these events and potential biases therein. Most empirical work has focused on individual countries, particularly the United States, providing a deep understanding of the political economy of media, economics, and politics Gentzkow2006MediaBias, DellaVigna2007FoxNews, Bursztyn2020Misinformation. However, transnational perspectives are essential to address global challenges like climate change.

Despite the increasing prevalence of climate change-induced disasters and the importance of empathy in promoting global cooperation, relatively few studies have examined cross-border media reporting on these events. Most existing work focuses on single countries, often the United States, offering insights into the political economy of media and politics Gentzkow2006MediaBias, DellaVigna2007FoxNews, Bursztyn2020Misinformation. Yet, addressing global challenges like climate change demands a broader, transnational perspective.

In this study, we analyze whether media sources linked to one country report more on another country following a natural disaster. We use a unique open-source intelligence data resource, the Europe Media Monitor (EMM) Atkinson2009NearRealTime, Steinberger2009JRC,jacquet2020covid, developed by the European Commission’s Joint Research Centre. The EMM compiles data from billions of publicly available online sources, employing natural language processing algorithms to extract key information—such as article content and event locations—from this vast media universe.

Combining the EMM dataset with the Emergency Events Database (EM-DAT) database of global natural disasters EMDAT2021, we construct a dyadic panel capturing cross-country media reporting before and after each of 3,769 disasters. Our dataset includes 466 sources from 123 countries, amounting to over 135 million articles indexed since 2016.

We estimate whether media outlets in one country systematically increase coverage of another country following a disaster, using a co-occurrence approach of article counts mentioning affected countries. With these estimates, we study how coverage varies by disaster type and by the severity, measured in fatalities, of the events. We then explore to what extent ties between countries correlate with the intensity and patterns of coverage. This analysis quantifies the extent of how existing biases in newsworthiness and global connectedness skew the portrayal of natural disasters.

Our findings suggest that media coverage is heavily skewed towards events that are visually salient, such as earthquakes and wildfires. In contrast, climatologically significant disasters like floods, droughts, storms, and extreme temperatures receive far less attention. This pattern appears to be remarkably consistent across countries, suggesting that it is a general feature of global media reporting. This aligns with findings that visually compelling imagery is crucial in eliciting public engagement and shaping perceptions of climate issues o2014climate

The number of deaths caused by a natural disaster appears to be one of the strongest predictors of the average increase in media coverage. The media in most countries exhibit a noticeable death gradient, where the degree of reporting increases with the number of fatalities. However, the influence of death count on media coverage varies significantly between countries, indicating that this relationship is not uniform.

When decomposing what drives differences in reporting patterns between countries, we find that historical measures of country connectedness, on their own, do not adequately explain these variations. Instead, a combination of natural disaster characteristics and measures of country-to-country connections plays a key role. Specifically, the explanatory power of the number of disaster deaths is significantly enhanced when considering the strength of social ties between countries. Genetic similarity also appears to be influential, possibly reflecting deeply rooted in-group preferences shaped by evolutionary and sociocultural factors simpson2015beyond, choi2007coevolution. Such parochial altruism and the moral significance of in-group membership may make certain disasters more ‘relatable’ and thus more newsworthy. This suggests that human connections, whether cultural, social, or historical, play an important role in fostering empathy and shaping how the media respond to disasters in other countries.

This study contributes to a deeper understanding of how skewed media reporting may shape perceptions of climate-related risks and highlights the need for more balanced media reporting to foster global awareness, empathy, and cooperation in addressing climate change. By shedding light on the systematic patterns and biases in media reporting, we hope to inform more effective strategies for addressing the global challenge of climate change.

Results

Media coverage following a natural disaster

Figure (ref)(a) displays the average increase in the volume of news articles attributed to a country after that country experiences a natural disaster. This figure focuses on the type of natural disaster and shows the average estimated increase in reporting within three days following the event. We find that, after a natural disaster occurs, media reporting by countries not directly affected rises unevenly, depending on the type of disaster. The increase is most pronounced for earthquakes, technological disasters (such as, for example transportation accidents), and wildfires. Specifically, the average increase in media reporting on a country affected by an earthquake, technological accident, wildfire, or storm is 7.4% (b = 7.38%, 95% CI [6.49,13.03], p < 0.001), 2.87% (b = 2.87%, 95% CI [2.54,4.93], p < 0.001), 2.08% (b = 2.08%, 95% CI [0.42,3.92], p < 0.05), and 1.33% (b =1.33, 95% CI [0.79,1.79], p <0.01), respectively. In contrast, we do not detect a robust and systematic increase in reporting when a country is affected by floods, droughts, or extreme temperatures.

To check the robustness of these findings, Appendix Figure (ref) explores a more flexible time window of -7 to +7 days relative to the start date of the natural disaster. This wider window accounts for potential anticipatory reporting (e.g., for floods) and confirms that the main pattern remains intact. Thus, focusing on the three-day window after the disaster start date provides a suitable estimate of pooled reporting effects.

To verify that the skewed reporting pattern by disaster type is consistent across reporting countries, Figure (ref)(b) plots the country-specific estimated increases for each disaster type (vertical axis) against the country-specific estimated increase for earthquakes (horizontal axis). Each point represents a reporting country, and the fitted regression line is shown without an intercept for visual ease. Although there is considerable heterogeneity, the overall pattern—more pronounced reporting on earthquakes, technological disasters, and wildfires—is systematic across countries.

Fatalities are an important driver of media reporting

Figure (ref) highlights another salient factor influencing the increase in media reporting after a natural disaster: the number of fatalities. Panel (a) of Figure (ref) presents the results of a non-linear regression, showing that the average increase in media coverage rises substantially when a natural disaster causes more than 20 deaths. This suggests threshold effects may influence editorial decisions on whether and how much to report on foreign disasters.

In Figure (ref)(b), we examine how this “death gradient” in media reporting varies across different reporting countries. The horizontal axis plots, for each reporting country, the average estimated increase in reporting following a natural disaster, regardless of the disaster’s characteristics. The vertical axis shows the country-specific death gradient, obtained from a linear regression estimated separately for each reporting country. This death gradient estimate is normalized by the country’s average reporting increase (shown on the horizontal axis), allowing us to relate the relative strength of the death gradient to a country’s baseline reporting level.

The results indicate that the positive impact of fatalities on media coverage is universal. However, countries that have lower baseline reporting increases show a more pronounced death gradient. In other words, the number of deaths drives reporting to a greater degree for countries that, on average, devote less coverage to foreign disasters. This finding suggests that, for media sources attached to these countries, the severity of a disaster (as measured by fatalities) is a particularly important cue for deciding whether to cover it.

Naturally, this result does not imply that the public necessarily pays attention to the detected statistical patterns in reporting. Rather, it indicates that in certain countries, the likelihood of reporting on a foreign disaster is heavily skewed by whether the disaster has caused a substantial number of immediate casualties.

Disaster characteristics and bilateral country connections jointly explain notable reporting differences

So far, our analysis has focused on how disaster characteristics—specifically type and fatalities—explain variations in reporting patterns across countries. We have found substantial heterogeneity: while disaster type and fatalities significantly influence media coverage, their impact is not uniform across different reporting countries.

This complexity suggests that additional factors, such as the nature of connections between countries, may interact with disaster attributes. Appendix Figure (ref) presents univariate relationships showing that various measures of country connectedness—including historical ties, cultural or linguistic similarities, values-based alignments, social connections, and genetic similarities—are generally positively associated with increased reporting on disasters in the other country. For example, country-pairs with stronger social ties or similar languages or religious compositions tend to report more on each other’s disasters. However, these associations are relatively small, indicating that connectedness measures alone do not account for much of the observed variation in reporting patterns.

To better understand these interactions and the underlying heterogeneity, we employ a random forest regression to decompose the variation in the estimated reporting increases $\hat{\beta}_{k(j),i}$. The random forest approach is well-suited for this task because it can capture non-linearities and interaction effects. In learning the random forest, we focus on the set of dyadic features that are consistently available—this choice excludes genetic dissimilarity data, which is only available for a smaller subset of country pairs.

Panel A of Figure (ref) shows how well the random forest models perform using three different sets of features: (1) the dyadic measures of country connectedness, (2) the characteristics of the natural disasters, and (3) a combination of both sets. When we use only dyadic connectedness measures, the best cross-validated random forest explains about 0.5% of the variation in reporting increases. Using only disaster characteristics—namely, disaster type, number of fatalities, and disaster duration—improves the explained variation to around 12%. Combining both sets of predictors raises the $R^2$ to about 19%. Although a large portion of the variation remains unexplained, this marked improvement suggests that country-connectedness measures interact strongly with disaster characteristics to shape cross-border reporting patterns. Appendix Figure (ref) confirms that this pattern is not an artifact of the chosen functional form for estimating $\hat{\beta}_{k(j),i}$, nor is it driven by the particular dependent variable specification—dyadic occurrence counts—used in the analysis.

Panel B of Figure (ref) compares variable importance measures for the disaster characteristics in models that include and exclude the dyadic connectedness features. Here, we report scaled permutation importance scores to ensure comparability. The results show that incorporating country connectedness increases the relative importance of the number of fatalities as a predictive factor. In other words, the explanatory power of fatalities in shaping media coverage appears to be enhanced when we account for how connected the countries are. This suggests that the perceived “newsworthiness” of a natural disaster, as proxied by fatalities, interacts with country-level linkages to influence media reporting decisions.

Panel C in Figure (ref) explores which specific country-connectedness features interact most strongly with the number of fatalities. To do this, we estimate a series of linear models, each introducing more stringent sets of fixed effects, and convert dyadic features into z-scores to facilitate comparisons. Two variables stand out: the degree of social connections between countries and their genetic similarity. Although genetic similarity data are available for fewer pairs of countries, its importance is still notable. By contrast, factors such as colonial ties, trade linkages, or religious or linguistic similarity of the populations appear to have much weaker effects on cross-border media reporting.

Overall, these findings indicate that social connections between the peoples of different countries are a primary factor enhancing media coverage of natural disasters abroad, especially as the number of fatalities increases.

\paragraph{Visualising these patterns} To make our findings and their implications more salient, we develop an interactive visualization that illustrates cross-country differences in how media report on natural disasters.\footnote{This visualization, along with other plots that facilitate data exploration, is available at \url{https://www.trfetzer.com/megeo-dyadic-news/}.} Specifically, we use the final random forest model to generate out-of-sample predictions, $\widehat{\beta}_{i,j}$, for a suite of hypothetical yet common natural disasters occurring in selected locations. Each hypothetical disaster is assigned varying numbers of fatalities to capture different severity levels.

This procedure enables us to systematically plot, for a given reporting country, how the estimated reporting increase $\widehat{\beta}_{i,j}$ changes with the disaster’s death toll. Figure (ref) illustrates this approach by focusing on two reporting countries—Germany and India—and three hypothetical disaster locations: Bangladesh, Italy, and Mexico. In this example, each scenario is a storm with a varying number of fatalities.

The fitted values allow us to make like-for-like comparisons across countries and scenarios. For instance, from the perspective of the German media, a storm causing 100 fatalities in Bangladesh elicits roughly the same relative reporting increase as a storm with about 43 fatalities in Italy or 62 in Mexico. Conversely, for Indian media, a storm with 45 fatalities in Bangladesh yields the same approximate coverage as a storm with 100 fatalities in Mexico or Italy. These simulations highlight a core observation: although the absolute severity of a disaster matters, factors like cultural proximity, social connectedness, and genetic similarity significantly modulate how foreign media perceive and prioritize an event.

By conducting such counterfactual comparisons, we show that two disasters with the same death toll can generate markedly different levels of coverage, depending on the countries involved and the specific audiences. Hence, what may appear “close” or “relevant” to one set of media outlets—and, by extension, their audiences—can be strikingly different for another. This insight aligns with our broader finding that media coverage is not solely driven by the raw number of casualties but also by deeper social, cultural, and evolutionary ties between populations.

Discussion

Our study provides evidence that cross-border media reporting after natural disasters is heavily skewed toward certain event types—particularly earthquakes and wildfires—while underreporting other climate-relevant disasters such as floods, droughts, storms, and extreme temperatures. This imbalance in coverage may have far-reaching implications for public understanding of climate change, the growth of global empathy, and the development of international cooperation to address this urgent challenge.

Climate negotiations are often fraught with tensions, exacerbated by continued efforts to cast doubt on the scientific consensus. Media reporting may play a crucial role in fostering public awareness of both, the scientific consensus as well as the consequences of climate change. Rising geopolitical tensions, which risk fragmenting the global information sphere, may provide avenues for countries less concerned about the consequences of climate change to suppress reporting on adverse climate impacts, such as natural disasters. This suppression of information risks undermining popular support for climate action at a critical time.

This concern is heightened by the key role that elite messaging and media narratives play in influencing public perceptions of climate change. For example, merkley2020party highlights how partisan media can foster skepticism. Social norms and elite endorsements remain critical for garnering public support for mitigation policies. rinscheid2020public stresses that descriptive norms and elite cues can shape climate policy approval, emphasizing the need for careful, strategic communication. In an era defined by fragmented media and geopolitical discord, effectively using these tools is vital to counter skepticism and build momentum for global climate cooperation.

Climate change presents risks not just through direct impacts, but also via complex, cascading effects IPCC2018, Beck2009WorldRisk. The increasing prevalence of climate-related disasters underscores the urgency of global action Field2012Managing, Coumou2012Decade. However, traditional media often appears “blind” to many climate-induced natural disasters, especially those that lack dramatic visuals or immediate effects Boykoff2009WeSpeakForTheTrees, Doulton2009TenYears. This underreporting may limit public appreciation of the full, multifaceted nature of climate risks Schmidt2013MediaAttention, Painter2013Climate, Moser2010Communicating.\footnote{How climate risks and responsibilities are communicated can strongly influence both public engagement and policy outcomes Dryzek2011Oxford, stevenson2012discursive. Empirical work on climate change communication suggests that more balanced and inclusive reporting can foster more robust public debate and enhance global solidarity moser2016reflections.}

Our findings also show that increases in reporting are strongly skewed toward disasters that cause significant fatalities, implying potential threshold effects in drawing media attention. This relationship, however, varies considerably across countries. Research on the link between media coverage and public awareness, especially regarding global challenges like climate change, is extensive Pralle2009AgendaSetting, McCombs2004SettingAgenda. Media bias—where certain events are overemphasized—has a substantial influence on public understanding of critical issues Zhuravskaya2020PoliticalEffects, Prat2011MediaInfluence, Lacy2001EffectsOfContent. Such editorial decisions align with long-established principles of newsworthiness, where sudden, dramatic events are more likely to receive extensive foreign coverage galtung1965structure. Similar patterns have been documented more broadly in journalism studies, where certain event attributes systematically increase the likelihood of extensive reporting wahl2009handbook. The adage “if it bleeds, it leads” captures the media’s tendency to focus on dramatic events Soroka2019BadNews, Harcup2001NewsValues, which may distort perceptions of climate change by emphasizing sudden catastrophes over gradual, equally consequential changes.

A significant factor shaping these reporting biases is the interaction between disaster characteristics and country connectedness. Social ties, as measured by social connectedness, are particularly influential: media outlets devote more coverage to disasters in countries that share closer social bonds, even if the number of fatalities is the same. Genetic similarity between populations plays another central role. We can not distinguish whether this is due to genetic similarity being associated with stronger social ties and the associated improved information flows or whether it may capture a structural factor that may increase relatedness.

Empathy is crucial for motivating prosocial behavior and global cooperation, and is considered central to human moral and cognitive development Hoffman2001Empathy, Tomasello2019BecomingHuman, Decety2014ComplexRelation.\footnote{Theories of empathy, such as Martin Hoffman's "empathic arousal" Hoffman2001Empathy, suggest that empathy is triggered by the ability to emotionally and cognitively relate to the suffering of others.} Media coverage can evoke empathy by bringing distant suffering closer, both emotionally and cognitively Kearney2013CrossCulturalEmpathy, Decety2010Empathy. However, when coverage disproportionately focuses on dramatic disasters, it may trigger selective empathy\footnote{Empathy alone may not guarantee constructive collective action, particularly when it is narrowly focused or selective. Research suggests that expanding the circle of moral concern and understanding group identity dynamics is essential for fostering broader empathy brewer2007importance, while cautioning against uncritically assuming empathy will always yield beneficial outcomes bloom2017against}, overshadowing victims of less visually striking events like floods or droughts. Similarly, skewed coverage toward disasters with higher death tolls or in countries deemed more relatable may distort the global picture. Such biases in empathy distribution could affect public engagement with climate issues and shape policy preferences Batson2011Altruism, Goleman1995Emotional.

Media-induced empathy—or its absence— may have wide-reaching implications for global climate efforts. Natural disasters are tangible manifestations of climate change, and accurate, balanced reporting on these events can foster empathy and solidarity deWaal2009AgeOfEmpathy. If media outlets focus narrowly on certain types of disasters, they risk skewing public perceptions of climate change’s risks and reducing empathy for underreported events, potentially weakening the impetus for collective action. Since empathy drives prosocial behavior and collective action, correcting these biases in coverage is essential for building the global cooperation needed to address climate change Hulme2009WhyWeDisagree, Beck2009WorldRisk.

Empathy in response to media coverage is influenced not only by the type of disaster but also by the cultural, social, and ancestral proximity between countries. Studies in political economy show how ancestral links and transnational relationships shape economic and political ties Burchardi2019Migrant, Bahar2018Migration. Thus, media reporting on disasters may serve as a conduit through which “soft power” is either amplified or diminished Nye2004SoftPower. Our data indicate that simple dyadic measures of connectedness alone do not fully explain these patterns. Instead, the interaction between disaster characteristics and connectedness, especially social ties and genetic similarity, significantly influences cross-border media coverage. These findings suggest that media reporting may disproportionately highlight disasters that are likely to generate empathy and a sense of relatedness, reinforcing unequal attention to certain regions or disaster types Petrova2011NewspapersPoliticalCompetition, Fan2013Ownership.

Our study has several limitations. The EMM dataset may not uniformly cover all media sources, and EM-DAT may have incomplete records for less severe disasters or those in underrepresented regions Jones2022Evaluating, Below2009DisasterDatabase. Further research is needed to clarify the causal mechanisms behind the associations we observe.

Future investigations should identify root causes of media bias in disaster coverage and its direct effects on public perceptions, empathy, and policy preferences. While this study provides a statistical inference approach to detect average patterns across different disasters, experimental or quasi-experimental studies could directly test how exposure to various disaster types affects climate awareness and empathy. Additionally, examining the influence of social media and other information channels could yield a more comprehensive understanding of how information spreads in the digital era Allcott2017SocialMediaFakeNews, Vosoughi2018SpreadOfMisinformation, gargfetzer2024. For instance, the observed gradient in reporting tied to social connections may be driven by the ease of information flow rather than culturally conditioned empathy.

Despite these uncertainties, our findings highlight the need for more balanced and comprehensive media coverage of climate-related disasters to enhance global awareness, empathy, and cooperation. Media organizations have a responsibility to report on a wide range of events, emphasizing the links between climate change and natural disasters. Policymakers and educators can also strengthen public understanding of climate risks, working to counteract reporting biases Boykoff2011WhoSpeaksForClimate, Moser2007MoreBadNews. By promoting accurate, inclusive coverage, it may be possible to deepen global solidarity and inspire the collective action required to mitigate the escalating risks of climate change.

Methods

Primary Data Sources

We draw on two primary data sources. First, we use the Europe Media Monitor (EMM), an open-source intelligence platform developed by the Joint Research Centre of the European Commission Atkinson2009NearRealTime, Steinberger2009JRC. The EMM continuously scrapes publicly available online sources—covering text and images—from thousands of news outlets in up to 70 languages, accumulating billions of items since the early 2000s. Employing web scraping and natural language processing (NLP) algorithms, it extracts relevant information such as content and geographic references. The EMM operates under a legal mandate from the European Commission, making it a robust and comprehensive dataset for analyzing cross-country media coverage.

Second, we incorporate data from the Emergency Events Database (EM-DAT) EMDAT2021, maintained by the Centre for Research on the Epidemiology of Disasters (CRED). EM-DAT records global natural disasters, providing details such as type, location, duration, and number of fatalities. Widely used in disaster risk research and policy Below2009DisasterDatabase, EM-DAT offers a reliable reference for identifying and classifying disaster events since 2016.

\paragraph{Dataset Construction} The EMM capacity has grown rapidly with the expansion of the internet from its inception with a large number of sources that have been (semi) consistently indexed and made accessible through an OpenSearch interface. A prime challenge with media observatory work leveraging primary sources across countries is the lack of data consistency across content providers over longer duration. That is, news sources may enter- and exit a news index, rendering longitudinal comparisons as highly difficult. The EMM capacity, due to its scale and breadth, made it possible to define a consistent set of sources. For the period from 2016 to 2023 we have identified a sample of 466 sources that are attached to 123 countries that appear consistently indexed. We proceeded in the source selection using the following filtering criteria. A source is considered if:

itemize• The source must have been consistently indexed by EMM since 2016, producing near-daily articles. This ensures stable long-term coverage. • The source must be a primary content provider rather than a content aggregator or platform (e.g., not Google News or MSN News).

Since the EMM has extensive coverage in Europe, including local and regional outlets, we further restrict the sample to achieve a more balanced country-level representation. For each country, we include at most five sources. To determine the top sources for each country, we use web traffic and search rank data from Semrush. \footnote{For example, for the US the listing is provided on \url{https://www.semrush.com/trending-websites/us/newspapers}. We retrieve the associated rank and traffic data for each country.} For countries where fewer than five sources have available traffic data, we rank remaining sources by their total indexed articles between 2016 and 2023 and select the top ones.

Applying these filters results in a sample of 466 sources associated with 123 countries. Each source provides a consistent time series of articles from 2016 to 2023, yielding 134 million articles and allowing cross-country comparisons unaffected by source turnover.

\paragraph{Geographic attribution} The geographic aboutness of a news source is based on the country of attribution of a source. That is, the country $i$ that is attached to a source $s$ is static. The country that a specific news item may cover or relate to is retrieved from the EMM's in house geographic parsing solution. The EMM in-house approach leverages both a combination of country-specific gazetters as well as named entity recognition to identify names of places, organisations or people in the original news articles. The EMM geographic attribution pipeline has a range of disambiguation techniques in place to handle ambiguous place names. Each indexed article is tagged with a geographic aboutness or more if one has been detected. For the purpose of geographic attribution of an article to a country $j$ the country that is attached to the most frequently occurring locations is considered as the articles geographic aboutness.

\paragraph{Language handling} Our sample includes 466 sources publishing in 44 languages. Figure (ref) tabulates the number of sources by language: about 22% of sources are in English, and the top 10 languages cover 72% of all sources. EMM is explicitly multilingual with a skew towards European languages. The indexing focuses on the title, the first paragraph of a text as this would traditionally contain the most relevant and important information around a natural disaster event. For languages that are not natively supported, the EMM machine translates the content into English.

\paragraph{Dyadic count measure} Having mapped each source $s$ to a country $i(s)$ and attributed a geographic location $j$ to each article, we construct a measure of country-pair-level article counts as a time series. To achieve this, we use OpenSearch queries to generate a structured dataset where, for each source $s$, each date $t$, and each possible destination country $j$, we count the number of articles indexed by that source that reference country $j$ on that specific date. This approach yields a fully dyadic dataset, capturing all source-destination-country combinations over time.

\paragraph{Natural disaster keyword based measure} In addition to raw dyadic article counts, we refine our dependent variable using a keyword-based filter. This method relies on the EMM’s categorization engine, which employs keyword indexing developed by media analysts at the Joint Research Centre. These filters focus specifically on natural disaster-related content. While our primary analysis uses raw country-pair-specific dyadic counts to capture general attention, the keyword-based measure allows us to isolate and validate patterns specifically tied to natural disasters.

This approach may introduce idiosyncrasies or measurement errors. However, our empirical design addresses potential concerns by centering the analysis on robust, statistically detectable reporting increases. Details of the statistical adjustments and validation steps are described below.

Empirical Design of Event Study

Our analysis employs an event study design centered around each natural disaster event. Let $E_{k(j),t}$ denote a disaster $k$ occurring in country $j$ at time $t$. Given that the start and end dates of natural disasters may be imprecise, we account for this with flexible time windows in our model.

The dependent variable in our analysis dataset $y_{s(i)jt}$ provides for a dyadic daily count dataset measuring the number of times a news source $s$ attached to country $i(s)$ mentions or refers to news covering any other country $j$ on day $t$ either across all articles or, in natural disaster keyword filtered articles. In our case this dyadic dataset has dimensions $466 $ news sources $\times$ 123 countries $\times 8$ years $\times 365$ days = 167,368,560.

Since the focus of this paper is specifically on spillovers, we remove the observations $y_{s(i)jt}$ where $s(i) = j$. That is, we exclude own country pairs to ensure our media reporting differential is capturing the differential increase in media reporting in other countries, not the countries where the natural disaster is happening.

For each event $E_{k(j)t}$, we define a 7-day window around the disaster’s start and end dates. This accommodates events of varying durations and accounts for both lead and lag effects. For example, if a disaster begins in country $j$ on day $t$ and ends $n$ days later, the dataset for this event, $y_{s(i)jt_k}$, includes all dates from $t_k - 7$ to $t_k + n + 7$. We then estimate the following model:

equation[equation omitted — 223 chars of source]

This event study design employs a two-way fixed effects estimation that “centers” the data, removing idiosyncratic differences in reporting levels specific to each source-destination pair, around each event, while also controlling for common time shocks. More concretely, the fixed effects $\alpha_{s(i)j}$ capture baseline differences in reporting patterns from a news source $s(i)$ to a particular destination country $j$. This is important because some media outlets in country $i$ may consistently report more on country $j$ regardless of any disaster event, potentially due to linguistic familiarity, cultural ties, or historical interest. By including $\alpha_{s(i)j}$, we eliminate these persistent differences in baseline reporting, ensuring that any identified increase in coverage can be attributed to the occurrence of a natural disaster rather than long-standing reporting biases. Additionally, this centering also helps correct for possible misattributions in geographic references, which can arise due to place name ambiguities or incomplete disambiguation.

The day-level fixed effects $\eta_t$ capture common shocks or global news cycles that might influence reporting uniformly across all sources and destinations. For example, a major geopolitical event could affect worldwide news patterns on a given day. By including $\eta_t$, we control for these time-specific fluctuations, ensuring that our measured effects are not driven by day-to-day variance in global news focus.

The coefficient of interest, $\beta_{k(j),i}$, is identified through an interaction that isolates when, where, and for how long a natural disaster influences foreign reporting. Specifically, $\mathbbm{1}(j=d)$ flags the country $j$ as the one experiencing the disaster, and $\mathbbm{1}(t \in [t_k, t_k+\tau])$ indicates whether the current day $t$ falls within $\tau$ days of the disaster’s start date $t_k$. Multiplying these indicators ensures that $\beta_{k(j),i}$ measures the average change in reporting on country $j$ by media from country $i$ strictly in the wake of that specific disaster $k$. For our main analysis, we focus on $\tau=3$ days, but we also consider alternative windows, such as $\tau=7$ or $\tau=1$, to examine how quickly reporting spikes and then decays. Conceptually, $\beta_{k(j),i}$ represents the statistically detectable attention increase that occurs within a short time frame around the disaster event, relative to a baseline period (typically the two weeks surrounding the disaster window). This helps distinguish between normal reporting fluctuations and those specifically triggered by the occurrence of a natural disaster.

It is worth noting that we estimate a complete distribution of $\beta_{k(j),i}$ coefficients across all natural disasters $k$ and for each reporting country $i$. Given the size of our dataset—up to $3,769 \times 123 = 463,587$ potential estimates—this approach yields a rich portrayal of how different disasters and reporting countries interact. When we apply a $\log(1+y)$ transformation to the dependent variable, each estimated coefficient represents the relative increase in reporting on a particular $(i,j)$ dyad after a given disaster $k$. Additionally, we employ a shrinkage procedure, setting estimates to zero if their associated t-statistic is below 1.65. Standard errors are clustered two-way by source and destination to address potential correlation in the error structure. These refined estimates $\hat{\beta}_{k(j),i}$ are then used to explore patterns and heterogeneities in media reporting increases, examining how factors such as disaster type, fatalities, and country connectedness shape the observed reporting responses.

\paragraph{Functional form} In analyzing the dependent variable—the number of articles covering a given destination country—we consider three alternative transformations to accommodate zero-inflation and overdispersion. First, we use counts in their raw form, but the sparsity and overdispersion of the dyadic data present challenges. Second, we apply a $log(1+y)$ transformation, which moderates the impact of extreme values and handles zeros naturally by adding one before taking logs. Although adding one is arbitrary, it typically yields interpretable results. Third, we consider the inverse hyperbolic sine (IHS) transformation, which, like the log transformation, manages skewness and accommodates zeros without arbitrary shifts. In practice, we find that the choice of transformation does not qualitatively affect our main findings. For clarity and ease of interpretation, we therefore focus on the $log(1+y)$ transformation.

\paragraph{Characterizing $\hat{\beta}_{k(j),i}$} The core output of our estimation is the set of coefficients $\hat{\beta}_{k(j),i}$, each capturing the increment in media coverage from country $i$ to country $j$ following a disaster $k$. These coefficients reflect the excess attention that arises beyond baseline reporting patterns. Because the regression includes source-destination fixed effects and day fixed effects, the resulting $\hat{\beta}_{k(j),i}$ values can be interpreted as the part of reporting that is specifically attributable to the natural disaster shock, net of persistent reporting tendencies and global time trends. In other words, they represent statistically detectable boosts in coverage that can be linked to the occurrence of a particular disaster in a given country.

\paragraph{Relaxation of the sharp event timing assumption} While many disasters can be viewed as having a relatively clear starting point, others may have a more gradual onset or prolonged duration. To address this, we relax the assumption of a sharp event start. Instead of focusing solely on a window defined by $\tau$ days around the start date, we consider a more flexible model:

equation[equation omitted — 223 chars of source]

where $\kappa \in [-7,7]$ indexes days relative to the reported start date of the disaster. This approach allows us to trace how reporting evolves daily before, during, and after the event rather than simply aggregating over a fixed window. In this specification, we pool the coefficients $\beta_{k(j),\kappa}$ across reporting countries to simplify inference and visualize the average time profile of media attention. For comparison, we scale these flexible-time estimates so that they are directly comparable to the disaster-type-specific averages presented in Figure (ref), allowing us to see whether the identified patterns are robust to the timing assumptions.

Studying structure in the reporting increase estimates

Once we have obtained the full set of $\hat{\beta}_{k(j),i}$ estimates, we conduct a series of analyses to understand their underlying structure. We explore three main dimensions of heterogeneity: (1) differences across types of natural disasters, (2) variations linked to the severity of the disaster, particularly the number of fatalities, and (3) the role of country connectedness—encompassing geographic proximity, historical ties, cultural similarities, social networks, genetic relatedness, and economic relationships—in shaping the magnitude and distribution of reporting increases. Through these exercises, we aim to identify regularities in how foreign media respond to disasters and the conditions under which attention is most pronounced.

Disaster characteristics

In examining how disaster characteristics affect reporting patterns, we focus on two key attributes: the type of disaster and the number of associated fatalities.

\paragraph{Disaster type} To analyze how different disaster types influence reporting, we first compute the average estimated increase in reporting, $\hat{\beta}_{k(j),i}$, for each disaster type. Panel A of Figure (ref) displays these averages by disaster type. In Panel B we present reporting-country-specific results, comparing each disaster type’s average increase in reporting to that observed for earthquakes. We then fit a regression without an intercept to these relative increases. This approach demonstrates that the skew toward reporting on certain disaster types, notably earthquakes, is a broad, system-wide feature of the media landscapes represented in our data.

\paragraph{Number of fatalities} We also explore the extent to which the number of fatalities that are associated with an event may affect the estimated reporting increase. To do so, we compute a binned scatter plot dividing the empirical distribution of casualties across the $3,769 $ natural disasters into 50 bins and compute the average of the estimated reporting increase $\hat{\beta}_{k,i}$ for the natural disasters that fall into this bin. This is presented in Figure (ref) Panel A. The pattern is robust to controlling for natural disaster type.

To further explore differences in how countries’ media respond to disaster fatalities, we estimate the following regression model:

equation[equation omitted — 101 chars of source]

where $\nu_i$ represents the country-specific “death gradient” in reporting. We also estimate a version of the model without the fatalities measure, isolating $\xi_i$, which captures the average reporting increase for each reporting country, irrespective of disaster characteristics. Interpreting $\xi_i$ as a baseline scaling factor, we can relate the death gradient $\nu_i$ to this baseline. Figure (ref) Panel B plots the ratio $\frac{\hat{\nu}_i}{\hat{\xi}_i}$ against $\hat{\xi}_i$, revealing how sensitive each country’s media is, in relative terms, to the number of fatalities. This analysis allows us to assess the importance of mortality in shaping reporting increases across different national media environments.

Degree of country-connectedness

Next, we examine whether the increase in media reporting by country $i$, from a source $s(i)$, following a disaster $k(j)$ in country $j$ is influenced by how connected countries $i$ and $j$ are along various dimensions. We denote $w_{ij}$ as a measure of connectedness between countries $i$ and $j$. We consider a broad range of such measures, detailed below.

\paragraph{Measures of country connectedness} To quantify $w_{ij}$, we consider a suite of measures capturing geographic proximity, historical and cultural ties, linguistic and religious similarities, social and genetic connectedness, preference and attitudinal overlaps, and economic linkages (notably trade). Many of these measures are drawn from the CEPII Gravity Database CEPIIgravity2022, while others come from specialized data sources such as spolaore2024barriers, Bailey2018FacebookSCI, and Stojkoski2024Measuring.

\paragraph{Geography} Geographic connectedness is represented by the great-circle distance between the most populous cities of $i$ and $j$, providing a quantitative measure of proximity. We also include contiguity, a binary indicator for whether two countries share a land border, based on political boundaries as of 2020. Geographic proximity can enhance information flows and mutual interest, potentially increasing cross-border reporting.

\paragraph{Historical and cultural connectedness} We then consider a range of dyadic features that capture shared histories, languages, and cultural traits:

itemize• Same country A binary measure of whether two countries have been part of the same political entity in the past (e.g., former colonies or unified states). The underlying data comes from CEPII data. • Religious similarity, based on indices from spolaore2024barriers, which quantify overlap in the religious composition of the two countries. The underlying raw data stems from the World Religion Database (WRD). The measure accounts for hierarchical relationships between religions that have shared origins among Christianity, Islam, and Judaism and is weighted by the population. The measure intends to capture the expected overlap in religious affiliation between two randomly selected individuals from country $i$ and $j$ respectively. • Linguistic similarity, based on indices from spolaore2024barriers, which measures the degree to which the languages spoken in $i$ and $j$ are related (e.g., through shared official languages or linguistic roots). The underlying data stems from Ethnologue dataset which creates language family trees. The measure captures the "Normalized Tree Distance" calculating the shortest path connecting two languages along the tree, scaled relative to the root. The score ranges from 0 to 1 and is further adjusted to weight for the population in a country that speaks each language. • Cultural similarity, based on data from spolaore2024barriers, that measures similarity of cultures as an index constructed from the World Values Survey (WVS) Inglehart2014WVS, Welzel2013Freedom. The set of questions that are used for country-pair similarity calculation varies and is explained in more detail in their paper. We focus on the 2014 measure of cultural similarity related to the WVS from then. Countries appear more similar if the population of respondents in the WVS share a similar response profile. • Colony A binary indicator of whether the two countries have ever had a colonial relationship with one another.

Social and genetic connectedness reflects interpersonal and biological relationships.

Social Connectedness Index (SCI) from Bailey2018FacebookSCI quantifies the strength of interpersonal ties between two countries measued using data from Facebook friendships. Formally: \[ \text{SCI}_{ij} = \frac{\text{FB Friends}_{ij}}{\text{Users}_i \times \text{Users}_j}, \] where $\text{FB Friends}_{ij}$ is the total number of Facebook friendships between users in $i$ and $j$, and $\text{Users}_i$ and $\text{Users}_j$ are the total users in each country. We adjust this measure to focus on foreign ties by excluding same-country connections: \[ \text{Foreign Connection Share}_{i,j} = \frac{\text{SCI}_{i,j}}{\sum_{k \neq i} \text{SCI}_{i,k}}. \]

Genectic similarity We include measures of genetic similarity from Pemberton2013, which capture the extent of shared ancestry between the populations of $i$ and $j$. We also consider genetic similarity relative to the USA, which benchmarks the similarity of $i$ and $j$ to a common reference point, helping contextualize patterns of genetic relatedness.

To assess preferences and attitudinal similarities, we use data from the Global Preferences Survey (GPS) Falk2018GPS and the Global Climate Change Survey (GCCS) andre2024climate. The GPS provides country-level measures of six key economic preferences: patience, risk-taking, positive reciprocity, negative reciprocity, altruism, and trust.\footnote{Patience reflects time preference, measured through hypothetical choices between immediate and delayed rewards and self-assessment of willingness to delay gratification. Risk-taking captures willingness to take risks, assessed via hypothetical lottery choices and self-assessment. Positive reciprocity measures the tendency to respond kindly to friendly actions, based on scenarios involving returning favors. Negative reciprocity assesses the propensity to retaliate against unfair treatment, including costly punishment. Altruism is derived from hypothetical donation decisions and self-assessed willingness to give without expecting anything in return. Trust is based on self-reported belief in the benevolence of others' intentions.} To measure similarity or dissimilarity we compute a pairwise Euclidean distances as:

\[ d_{ij} = \sqrt{\sum_{k=1}^{n} (x_{ik} - x_{jk})^2}, \]

where $x_{ik}$ and $x_{jk}$ are standardized values of preference $k$ for countries $i$ and $j$, and $n$ is the number of preferences.

We proceed in a similar fashion to measure differences in climate change beliefs and policy preferences. The GCCS offers data on climate change beliefs and preferences from nearly 130,000 individuals in 125 countries. Key variables include willingness to pay to fight global warming, beliefs about others' willingness to contribute, social norms regarding climate action, and views on government responsibility.\footnote{Specifically, variables include: (1) Willingness to Pay (WTP) 1%—a binary indicator if respondents are willing to contribute 1% of their household income monthly to fight global warming; (2) Willingness to Pay Positive Amount—a binary indicator if respondents are willing to contribute a positive amount less than 1% if not willing to contribute 1%; (3) Willingness to Pay Index—a composite measure coding 1 for those willing to contribute at least 1%, 0.5 for those willing to contribute a positive amount less than 1%, and 0 for those not willing to contribute; (4) Belief about Others' WTP—respondents' estimates of how many out of 100 people in their country are willing to contribute at least 1%; (5) Social Norms—a binary indicator if respondents believe people in their country should try to fight global warming; (6) Government Responsibility—a binary indicator if respondents think the national government should do more to fight global warming.} We aggregate responses to the country level and compute pairwise Euclidean distances between countries based on these variables.

Economic connectedness is captured through trade relationships. We compute measures of goods trade by calculating the average share of a country's exports that go to a destination and the average market share of imports into a country that come from a source. For digital services trade, we use estimates from Stojkoski2024Measuring, who provide bilateral trade flows in digital products by combining corporate revenue data, consumption patterns, machine learning, and optimal transport methods.\footnote{Stojkoski2024Measuring estimate bilateral trade in digital products by collecting corporate revenue data from large digital firms and combining it with country consumption patterns. They use machine learning to augment missing data and optimal transport methods to allocate revenues to consumption patterns, allowing them to estimate trade flows in digitally ordered and delivered products such as video games, digital advertising, and digital intermediation services.}

In leveraging these measures, we aim to capture different aspects of connectedness between countries $i$ and $j$. Geographic and ethnographic measures account for physical proximity and shared history or language, facilitating media coverage through ease of communication and mutual interest. Cultural and preference similarities reflect shared values and attitudes, potentially influencing empathy or interest in news from each other. Social connectedness indicates the strength of interpersonal relationships across borders, which may drive media reporting through shared social networks. Economic connectedness through trade reflects the interdependence of countries' economies, which might affect media attention due to economic interests. By including these diverse measures in our analysis, we comprehensively capture the factors that may influence the extent to which media in country $i$ reports on natural disasters occurring in country $j$.

\paragraph{Country connectedness and differential media reporting} We conduct two sets of analyses to explore how country connectedness measures $w_{ij}$ relate to variations in media reporting increases. In the first set, we examine the sign of the relationship between $w_{ij}$ and the average reporting increase $\hat{\beta}_{k(j),i}$, without conditioning on any specific disaster characteristics. This helps us understand whether certain forms of connectedness are generally associated with higher or lower media coverage, although it may not fully explain the observed heterogeneity due to potential interactions. For clarity and comparability, we standardize the features $w_{ij}$ by converting them into z-scores.

In the second set of analyses, we focus on the interaction between $w_{ij}$ and the number of disaster-related deaths. This allows us to determine whether certain types of country connectedness modulate how strongly media coverage responds to disaster severity.

To study the sign of $w_{ij}$ on reporting, we estimate:

equation[equation omitted — 130 chars of source]

To examine the interaction between $w_{ij}$ and disaster fatalities, we estimate:

equation[equation omitted — 171 chars of source]

Here, we estimate six different specifications of these models, each adding progressively more control variables. These controls are captured by $x_{ijk}$ and include:

itemize$x_{ijk} = \emptyset$ in this case the coefficient $\nu$ captures the average difference in the reporting increase $ \hat{\beta}_{k(j),i} $ that we can attribute to the feature $w_{ij}$ unconditionally. • $x_{ijk} = \alpha_j$ in this case we remove potential level differences in the reporting elasticity on countries $j$ that may be attributable, e.g. to the countries $j$ unique global position. In Fetzermediamultiplier2024 we noted for example that there are notable level differences in reporting for specific countries that may also imply differential elasticities $ \hat{\beta}_{k(j),i} $ due to more attention being devoted to disasters in country $j$, on average. • $x_{ijk} = \gamma_i$ this removes potential level differences in the reporting increase in country $i$. This may capture potential distorting effects that the specific media source sample that we have available for country $i$ may have on the reporting elasticity. • $x_{ijk} =\alpha_j + \gamma_i$ this removes both the potential level differences in the estimated average reporting increase associated with the country of a news source $i$ or the country in which a disaster occurs $j$$x_{ijk} = \delta_k$ removes a disaster-specific level shifter that may affect the increase in the reporting increase, on average, for all sources in a common fashion. This can be thought of as the most general way of controlling for disaster-level characteristics. • $s_{ijk} = \alpha_j + \gamma_i + \delta_k$ combines all the previous additional control variables

The results of these analyses are presented in Appendix Figure (ref) for the sign of $w_{ij}$, and in Figure (ref) for the interaction with disaster deaths. In these figures, we plot the main coefficient $\hat{\nu}$ from equation (ref) and focus on the interaction term $\hat{\gamma}$ from specification (ref). The features have been turned into z-scores to ensure that the point estimates are more comparable. The stability of the point estimates across the specifications is a sign of the robustness of the interactive relationship to different variation. Note that in all cases, the effect of disaster deaths themselves is accounted for.

Each row in these figures displays six points, corresponding to the six specifications with increasingly stringent controls. By examining the consistency of point estimates across different sets of controls, we identify features of $w_{ij}$ that robustly correlate with media reporting increases even after accounting for multiple confounding factors. Emphasizing consistently positive or negative estimates helps us identify systematic relationships between country connectedness and cross-border media coverage of natural disasters.

Fitting a random forest

Finally, we fit random forests to the estimated $\beta_{k(j),i}$ values to assess how well the included features explain observed heterogeneity in reporting increases. This approach helps us determine the extent to which linear or non-linear interactions between dyadic connectedness measures and disaster characteristics can account for variation in media coverage.\footnote{Further, it is instructive to understand how much of the variation in the estimated $\hat{\beta}_{k(j),i}$ can be captured with the 15 combined features that we use to study the heterogeneity.}

We use the ranger package in R Wright2017Ranger, training 1,000 trees with a minimum node size of 30 and at least 10 instances per terminal node. We apply scaled permutation importance to evaluate the contribution of each predictor. We fit the random forest to the estimated $\hat{\beta}_{k(j),i}$ and carry out the analysis using three sets of features separately. We run three models:

enumerate• Dyadic Variables Only: This model uses only the dyadic country-pair measures of connectedness as predictors listed above. We exclude variables that are not available for most country pairs. This affects the genetic similarity measures and the survey-based measures from the Global Preferences Survey Falk2018GPS as well as the Global Climate Change Survey (GCCS) Andre2023GCCS. This leaves us with 10 variables. • Disaster Characteristics Only: This model uses only disaster-specific characteristics as predictors. These include the disaster type (categorical variable), the logarithm of the number of deaths, and the disaster duration in days. • Combined Model: This model includes both the dyadic country-pair measures and the disaster characteristics as predictors, totaling 15 variables.

Panel A of Figure (ref) plots the $R^2$ of each model. Appendix Figure (ref) provides the same plot for models where the $\beta_{k(j),i}$ have been trained learned with different transformations of the dependent variable or using the more refined natural disaster-relevant article counts. We note that the pattern remains the same throughout.

Panel B of Figure (ref) plots the variable importance for the main disaster variables comparing the full model that includes the dyadic characteristics with the best model that only includes the three disaster characteristics. As indicated we use permutation importance involves randomly permuting the values of a predictor variable and measuring the decrease in model performance (increase in mean squared error) resulting from the permutation. This method is less biased towards variables with many unique values or higher cardinality, making it suitable for our analysis.

We note that the increase in goodness of fit of the model that includes the dyadic characteristics is mostly driven by a sharp increase in the variable importance of the total fatalities measure, suggesting that the dyadic characteristics strongly inform the heterogeneity in reporting between countries on natural disasters in a way that is interactively shaped by the number of fatalities.

Panel C of (ref) explores which dyadic characteristic appears to most notably influence the reporting differences in the interaction with the number of disaster fatalities using the linear regression model outlined in specification (ref). In order to allow for comparisons across the different features we converted the different dyadic measures into z-scores before interacting them with the number of fatalities associated with a given natural disaster $k$. This allows for more direct comparisons. As highlighted in the discussion, the analysis suggests strongly that it is particular social- and genetic connectedness between countries that appear to shape the natural disaster deaths reporting gradient that we detected.

Visualization in an interactive tool

To deepen our understanding of the estimated results, we also build interactive visualizations that illustrate several core findings. We focus on two main components.

\paragraph{Visualisation of model estimates $\beta_{ji(k)}$ } We visualize the model observations, in particular the role of fatalities in shaping reporting, by simulating a large set of disasters common to a region with different numbers of fatalities. Using these hypothetical disasters, we construct fitted values from the random forest model, namely \(\widetilde{\beta}_{ji(k)}\), for each simulated natural disaster. We then compute the average \(\bar{\beta} = \widetilde{\beta}_{ji(k)}\) for each country-pair, grouping outcomes by bins of total deaths and absorbing the type of disaster dimension. To explore the empirical distribution, we plot these values with reference to a reporting country’s percentile statistics. Specifically, let \(P_{z}(\beta_{ij} \mid j)\) indicate the \(z\)-th percentile of the empirical distribution conditional on a reporting country \(j\).

We describe two conceptual viewpoints. First, the country-of-disaster view centers on how reporting countries \(j\) vary in their attention to a specific disaster-affected country \(i\). The data in this view are normalized within the focal disaster-affected country \(i\), across different reporting countries \(j\). This normalization adjusts for differences in absolute attention levels; for instance, some disasters may attract more global coverage due to size, media focus, or geopolitical relevance. Formally, the normalization is: \[ \text{Normalized}_{i} \, \widehat{\beta_{ij}} = \frac{2 \times \bigl(\beta_{ij} - P_{5}(\beta_{ij} \,\vert\, j)\bigr)}{P_{95}(\beta_{ij} \,\vert\, j) - P_{5}(\beta_{ij} \,\vert\, j)} - 1, \] where \(P_{5}\) and \(P_{95}\) are the 5th and 95th percentiles, respectively.

Second, the country-of-reporting view shifts the focus to how a single reporting country \(j\) distributes its coverage across multiple disaster-affected countries \(i\). Here, we normalize across the different countries \(i\) that are hit by disasters, allowing us to see if the same reporting country devotes disproportionate attention to certain regions or disaster types. The normalization in this view is: \[ \text{Normalized}_{j} \, \widehat{\beta_{ij}} = \frac{2 \times \bigl(\beta_{ij} - P_{5}(\beta_{ij} \,\vert\, i)\bigr)}{P_{95}(\beta_{ij} \,\vert\, i) - P_{5}(\beta_{ij} \,\vert\, i)} - 1. \] This two-tiered perspective allows us to examine whether biases or preferences in coverage emerge in different dimensions (reporting country vs.\ disaster-affected country). For instance, a reporting country might systematically prioritize certain regions based on geopolitical, cultural, or economic ties. Conversely, a disaster-affected country might draw universal attention from most sources because of high visibility or severity.

We also restrict 90% of the \(\widehat{\beta}_{ij}\) observations to lie within the \([-1,1]\) range by trimming at the 5th and 95th percentiles. This avoids excessive sensitivity to extreme outliers and centers the analysis on the main mass of the distribution. In practical terms, this procedure reveals how much coverage is “typical” or “excessive” for a given scenario.

\paragraph{Path of random forests} In addition to estimating a single fully saturated random forest, we also run a series of \(2^{10}-1\) random forests that progressively add more dyadic features alongside the baseline disaster characteristics. This best-subset selection approach traces out the path of potential random forest models, each with a distinct combination of connectedness variables. Ultimately, each specification produces an \(R^2\) measure reflecting how much additional variation in \(\hat{\beta}_{k(j),i}\) the included predictors can explain. Appendix Figure (ref) shows the full distribution of these goodness-of-fit values, highlighting the role of specific features such as social connectedness, which often considerably boosts the model’s predictive accuracy.

\paragraph{Counterfactual estimates} Finally, we aim to make the analysis more tangible by conducting counterfactual experiments that gauge how “close” or “distant” certain disasters may appear to foreign media. Our approach builds on the observation that social and family ties between countries can alter perceptions of disaster severity. Suppose a particular country \(i\) experiences a natural disaster of a specific type with a given number of fatalities. From the perspective of a different country \(j\), one might ask how many fatalities a similar disaster in yet another country would require to elicit the same relative increase in media coverage. This question helps evaluate the relative proximity or connectedness between different countries’ populations.

To implement this exercise, we leverage the random forest model that incorporates both disaster attributes and dyadic features. We first simulate a set of common disaster types and country locations. Next, we vary the death toll \(c\) (e.g., from 10 to 300) for each hypothetical event and compute the fitted values \(\widehat{\beta}_{ij}(c)\). This allows us to obtain a matrix of predicted coverage levels across countries and severity levels, enabling us to systematically evaluate patterns of reporting. Figure (ref) illustrates this concept by plotting the predicted reporting increases for storms in Bangladesh, Mexico, and Italy, respectively.

The results highlight the notion of equivalent attention. For example, we ask: how many fatalities in one country would generate the same level of media attention as a disaster with 100 fatalities in another country? From the perspective of German media, a storm with 100 fatalities in Bangladesh elicits a similar reporting increase to a storm with approximately 43 fatalities in Italy or 62 in Mexico. Conversely, for Indian media, 45 fatalities in Bangladesh yield a level of attention comparable to that generated by 100 fatalities in either Mexico or Italy. These counterfactual scenarios make salient the often skewed nature of reporting, where equivalent attention depends not only on the raw severity of the disaster but also on the social, cultural, and geographical ties between the countries involved.

By systematically exploring these scenarios, we show that coverage patterns are highly interactive. They reflect both the absolute severity of the disaster and the relational connectedness between reporting and affected countries. Equivalent attention thus underscores how shared histories, cultural affinities, or even economic ties can modulate media responses to ostensibly similar events, highlighting non-linear effects and biases inherent in cross-border reporting.

Figures and Tables

landscape\begin{figure}[H] \caption{Increases in media reporting following a natural disaster are strongly driven by the disaster type} \captionsetup{singlelinecheck=off,font=scriptsize} \caption*{Note: Figure displays the estimated increase in dyadic media coverage across 3,769 natural disasters that occurred since 2016 by natural disaster type. We construct an event-level panel dataset that captures dyadic news coverage across the 466 covering 123 countries. For each natural disaster, ocurring in a country $i$ we construct a 7 day window around the disaster start and the disaster end date. Around this event dataset we attach dyadic daily dataset that measures the number of media articles of a news source $s$ attached to country $j(s)$ that mentions or refers to news covering any other country. Throughout we remove the distorting effects of any other time-varying shock and the average level differences in the propensity of a news source to cover any country. Panel A documents the average increase in news articles from sources $j(s)$ within 3 days after a disaster hit country $i$, relative to its news coverage of all other countries. Panel B documents the average increase in news articles from sources $j(s)$ at different points in time relative to three days before the start of the natural disaster event to allow for anticipation effects.} \end{figure}
landscape\begin{figure}[H] \caption{Media reporting increase strongly informed by the number of deaths} \captionsetup{singlelinecheck=off,font=scriptsize} \caption*{Note: Figure displays the estimated increase in dyadic media coverage across 3,769 natural disasters that occurred since 2016 by natural disaster type. We construct an event-level panel dataset that captures dyadic news coverage across the 466 sources in 123 countries. For each natural disaster occurring in a country $i$, we construct a 7-day window around the disaster start and end dates. For this event dataset, we attach a dyadic daily dataset measuring the number of media articles of a news source $s$ in country $j(s)$ that mentions or refers to news covering any other country. Throughout, we remove the distorting effects of any other time-varying shock and the average level differences in the propensity of a news source to cover any country. Panel A highlights that the estimated increase in news reporting, across disasters, is strongly and non-linearly increasing in the number of deaths. Panel B highlights that this death gradient varies across countries. In Panel B, the x-axis captures the average country-specific reporting increase, while the y-axis is the ratio of the estimated country-specific death gradient divided by the value of the x-axis to see the two in relations. } \end{figure}
landscape\begin{figure}[H] \caption{Measures of country-connectedness capture systematic variation in media reporting between countries when interacted with natural disaster characteristics} \captionsetup{singlelinecheck=off,font=scriptsize} \caption*{Note: Figure displays analysis of the estimated increase in dyadic media coverage across 3,769 natural disasters that occurred since 2016. The left panel presents the results from training a random forest that aims to decompose the variation in the estimated increase in media reporting $\hat{\beta}_{k(j),i}$ of a natural disaster $k$ occurring in country $j$ through media sources associated with country $i$. To explain the variation in this estimated reporting increase we consider either random forest that is grown using just measures of country-level connectedness labeled dyadic characteristics only, or, we consider three natural disaster characteristics (its duration, its type and the number of deaths) associated. The third set combines both of these. The left panel presents the goodness of fit of the random forests that are grown with this set of features, while the right panel displays the variable importance of each measure. The variable importance is scaled to the maximal permutation variable importance that is measured within each model. Appendix Figure (ref) highlights the sign of the relationship between a dyadic characteristic and the extent of media coverage.} \end{figure}
figure[figure omitted — 1,251 chars of source]