- Research
- Open access
- Published:
Social inequalities and COVID-19 mortality between neighborhoods of Bariloche city, Argentina
International Journal for Equity in Health volume 22, Article number: 198 (2023)
Abstract
Background
The COVID-19 pandemic has shown how intraurban inequalities are likely to reinforce health and social inequalities. Studies at small area level help to visualize social inequialities hidden in large areas as cities or regions.
Aim
To describe the spatial patterning of COVID-19 death rates in neighborhoods of the medium-sized city of Bariloche, Argentina, and to explore its relationship with the socioeconomic characteristics of neighborhoods.
Methods
We conducted an ecological study in Bariloche, Argentina. The outcome was counts of COVID-19 deaths between June 2020 and May 2022 obtained from the surveillance system and georeferenced to neighborhoods. We estimated crude- and age-adjusted death rates by neighborhood using a Bayesian approach through a Poisson regression that accounts for spatial-autocorrelation via Conditional Autoregressive (CAR) structure. We also analyzed associations of age-adjusted death rates with area-level socioeconomic indicators.
Results
Median COVID-19 death rate across neighborhoods was 17.9 (10th/90th percentile of 6.3/35.2) per 10,000 inhabitants. We found lower age-adjusted rates in the city core and western part of the city. The age-adjusted death rate in the most deprived areas was almost double than in the least deprived areas, with an education-related relative index of inequality (RII) of 2.14 (95% CI 1.55 to 2.96).
Conclusion
We found spatial heterogeneity and intraurban variability in age-adjusted COVID-19 death rates, with a clear social gradient, and a higher burden in already deprived areas. This highlights the importance of studying inequalities in health outcomes across small areas to inform placed-based interventions.
Introduction
By mid-2022 the Latin American region, which houses only 8.4% of the world population [1], had experienced 26.6% of the global COVID-19 deaths [2], making it one of the most affected regions by the pandemic [3]. Cities were especially vulnerable to viral spread, especially among residents living in segregated neighborhoods [1, 4, 5] due in part to working conditions, overcrowding, and poor access to quality health services. Moreover, co-existing health conditions linked to social disadvantage also magnified the severity of COVID-19 [1, 4].
This was noteworthy in Latin America, where almost 80% of the population lives in cities [6] and deep-rooted inequalities affect cities in this region. Countries as Argentina, that went through a long and strict lockdown, as a mean of avoiding healthcare services collapse [7], showed a greater impact on slums dwellers or other lower income groups [8]. These spatial differences could be explained by contextual characteristics where social distance and isolation were difficult to maintain over time. Aside from material deprivation, other explanations include labour informality, dependency on public transport, and limited possibilities to work from home [9, 10].
Studies assessing the impact of social inequities in COVID-19 incidence, hospitalization, and mortality in small areas within cities have been documented in high-income countries [11,12,13,14,15,16], but they are not common in middle-income countries, partly due to challenges in obtaining COVID-19 data at small levels of disagregation such as neighborhoods. Nonetheless, these studies have demonstrated a strong association between socioeconomic vulnerability and the degree of COVID-19 severity and death [1, 8, 17,18,19,20].
Researchers have widely documented that the risk of developing severe complications and dying from COVID-19 increases sharply with age [21]. However, the magnitude differs across countries: younger people are more likely to die from COVID-19 in middle-income countries compared to high-income countries [21, 22]. Some explanations of these differences may be that the young population in middle-income countries has a higher prevalence of pre-existing diseases, a higher proportion of informal employment, and higher overcrowding, the latter disproportionately affecting the youngest and also related to a more rapid spread of the disease [21].
Research on social inequalities and the differential impact of the pandemic in Latin American cities are more focused in comparing relatively large areas, which may hide meaningful differences across neighborhoods. The development of statistical modeling techniques, such as Bayesian modeling, has facilitated the investigation of the relationships between small areas and helped obtain reliable estimates, particularly when working with small numbers in small areas such as neighborhoods [23,24,25].
Studies focused on small areas could help policymakers to better understand how these health events are distributed and clustered within cities and target effective interventions at microsocial levels to improve health. Therefore, our aim was to analyze spatial patterning of COVID-19 death rates in neighborhoods in Bariloche city and describe its relationship with socioeconomic characteristics.
Methods
Study setting
We conducted this study as part of the SALud URBana en America Latina (SALURBAL) project [26]. The study area was Bariloche, a city of 134,000 residents in 2018, located in the southern region of Argentina, surrounded by a national park with large green areas, part of which is along the southern margins of the Nahuel Huapi Lake. Known worldwide as a tourist attraction, it is home to nearly half of its population lives in neighborhoods with poor and critical living conditions and high income inequalities [27].
Design and data sources
We performed an ecological study of 159 census tracts from Bariloche, known locally as radios censales. Census tracts, henceforth neighborhoods, are the smallest geostatistical census units defined by the national statistics bureau (INDEC) composed by about 300 households. We used confirmed deaths of COVID-19 that were systematically reported to the National Health Surveillance System (SNVS for its acronym in Spanish). These information was complemented by local authorities through the systematic review of other registries, such as civil registration and hospital records [28]. Each record relates to a single death containing information such as age and place of residence. We obtained data from the 2010 census to proxy socioeconomic characteristics of neighborhoods, compiled a by the SALURBAL project.
Outcome
The outcome was counts of COVID-19 deaths among residents of Bariloche city, spanning the period between June 2020 (first COVID-19 death in Bariloche) and May 2022. Deaths were georeferenced using the street address shapefile from OpenStreetMap in QGIS 3.18® and geocoded to the corresponding census tract using the shapefile of polygons provided by the INDEC [29]. A total of 1.9% of COVID-19 deaths (10 out of 529) could not be georeferenced and were excluded from the analysis. Confidentiality has been protected according to the Protection of Personal Data Law, No. 25,326, and the use of Sensitive Data and the Statistical Confidenciality Law, No. 17.622. Individual information was only used for georeferenciation; all data for analyses was de-identified and aggregated as counts of deaths by age-groups per neighborhood.
We obtained the 2010 census population for each of the neighborhoods by age from REDATAM (acronym of REtrieval of DATa for small Areas by Microcomputer), a software developed by CELADE used in Latin America for processing census data [30]. We also obtained the projected population estimates for 2018 (details in supplementary material, Appendix 1), the latest year for which the data was available, by age at the city level (not by neighborhood) from Dirección de Estadística y Censos for the Río Negro province [31]. Applying the age-and neighborhood-specific population proportion from the census 2010 population to the projected 2018 population, we approximated the age-and neighborhood-specific population for 2018, as done in prior studies examining small-area inequalities in mortality in other cities of Argentina [32, 33]. Because our study window spans 2 years, we doubled the estimated 2018 population by age and neighborhood as the population denominator in calculating mortality rates.
Given the role of age in determining mortality, we estimated both crude and age-adjusted death rates due to COVID-19 cumulatively through the end of the study period. We performed direct age-adjustment using the WHO world standard population 2000–2025 [34].
Exposures
The key exposure variable is the proportion of households in the neighborhood with unmet basic needs, a measure of poverty used in Latin America based on census data [35]. This variable is defined as the percentage of the households with at least one of the following indicators: (a) households with critical overcrowding (3 or more people per room) (b) inconvenient housing units (defined as tenancy piece, precarious households or other type) (c) housing with inadequate services (defined as lack of toilet) (d) households with school-aged (6 to 12 years) children not attending school, and (e) households with high economic dependence (households which have 4 or more individuals per working member and with the household head who has not completed the third grade of primary school) [36].
In addition, we also included the following measures of area-level socioeconomic status (SES): completed high school (% of the population aged 25 years and above with completed high school), piped water inside the dwelling (% of households with piped water inside the dwelling), overcrowding (% of households with more than 3 person per room), school attendance (% of population aged 15 to 17 years attending school), sewage connection (% of households connected to a sewage system of any type), and unemployment rate (% of the population 15 years or older who are unemployed and are actively seeking employment).
These variables were selected since they have been associated with other outcomes in small areas in Latin American cities in previous studies [27, 32, 33, 37, 38], are good measures of area-level socio-economic status (SES), are directly linked to COVID-19 transmission [18], and are commonly available from census data at small area level. They were all retrieved from the 2010 census [30].
Statistical analysis
We aggregated the deaths in each neighborhood into 4 age groups (0–39 years, 40–59 years, 60–79 years, 80 years and above) and modeled the age-specific mortality rates via a Poisson model with the number of deaths by age group as the outcome and population by age group as the offset. We used a Bayesian approach, accounting for spatial dependency, which enables the borrowing of information from neighboring areas and is particularly useful in obtaining reliable estimates when the number of events (sample sizes) is small [39].
Specifically, we modeled the number of COVID-19 deaths in age group \(a\) for small area \(i\), denoted \({y}_{ia}\), as being Poisson-distributed with mean \({n}_{ia}{\uplambda }_{ia}\), where \({n}_{ia}\) and \({\uplambda }_{ia}\) denote the population size and COVID-19 mortality rate, respectively. To model \({\uplambda }_{ia}\), we assume:
where \({\beta }_{0a}\) are the age-specific intercepts corresponding to the overall log mortality rate in each of the four age groups, \({z}_{i}\) is spatially structured random effect (clustering) modeled by conditional autoregressive (CAR) distribution and \({\phi }_{ia}\) is an unstructured, zero-mean random effect (heterogeneity). We assigned vague, uninformative priors with large variances for each of the parameters of interest, thereby allowing inference on these parameters to be driven primarily by the data.
The models were run using the R2WinBUGS package in R [40] and 4,000 iterations’ worth of samples were obtained for inference after first discarding pre-convergence samples as burn-in and thinning samples to reduce autocorrelation (Software code available in supplementary material, Appendix 2). We then used the posterior samples of age-specific mortality rates derived from the model to calculate crude and age-adjusted rates. To obtain the crude rates for each neighborhood, we aggregated the posterior samples of age-specific mortality rates based on the age distribution within the neighborhood. For age-adjusted rates, we used the population weights from the WHO world 2000–2025 age distribution. As a result, we obtained 4000 estimates of crude- and age-adjusted mortality rates for each neighborhood. We present the median of these estimates in a choropleth map of the crude and age-adjusted mortality rates.
To study the association between age-adjusted mortality rates and socio-economic characteristics of the neighborhoods, we fit two sets of simple linear regression models (a) with the logarithm of the age-adjusted mortality rate as the outcome and (b) mortality rate per 10,000 as the outcome. We used seven predictor variables that proxy the socio-economic context of neighborhoods, one at a time, for easier interpretation. We parameterized explanatory variables by creating the deciles of each of them and re-coding the value of the exposure variable on a 0–1 scale, such that values within the 1st decile are recoded to 0, those in the 2nd decile are recoded to 1/9, and so on. The estimated slope from regression analysis with this parameterization of explanatory variables and mortality rate as the outcome variable is the slope index of inequality (SII), a measure of the absolute inequality in the mortality rate between areas with high values of explanatory variables (9th decile) and areas with low values of explanatory variables (1st decile), while the model with log mortality rate as the outcome variable results in the relative index of inequality (RII), a relative inequality (ratio) measure. We ran the regression analyses for each of the 4000 samples and pooled the regression coefficients and standard errors with Rubin’s formula [41], using mitlml R package.
Bivariate choropleth maps were created to visualize associations between the outcome (age-adjusted COVID-19 mortality rates) and each socio-economic characteristic.
Results
From June 2020 to May 2022, a total of 529 COVID-19 deaths were registered in Bariloche. The median number of deaths and population per neighborhood were 3 and 789, respectively. A median of 90.0% (10th percentile(p10)- 90th percentile (p90): 78.4–100) of the population aged 15 to 17 years was attending school, but a lower proportion (median 55.9%; p10-p90 25.6–82.9) of individuals over 25 years completed high school. A median of 4.8% (p10-p90: 0.6–22.7) of the population had unmet basic needs and 94.3% (p10-p90: 58.6–100) of the households were connected to a sewage system revealing considerable heterogeneity between neighborhoods (Table 1).
Figure 1 displays maps of modeled crude (top panel) and age-adjusted (bottom panel) COVID-19 death rates, obtained from Bayesian spatial modeling. Crude death rates across neighborhoods ranged from 5.2 to 41.8 deaths per 10,000 inhabitants, and age-adjusted rates from 8.5 to 31.3 per 10,000 inhabitants (Fig. 1). Death rates were not randomly distributed showing a spatial pattern with clear clustering within the city. The crude death rates were higher in the core area of the city extending to the immediate western neighborhoods. In contrast, the distribution of age-adjusted mortality rates was a mirror image of the crude rates, with higher rates in the southeastern neighborhoods, and lower rates in the core area and westward in the neighborhoods along the shoreline.
We found that all SES variables examined were associated with age-adjusted COVID-19 mortality rates: neighborhoods with lower SES had higher COVID-19 death rates than those with higher SES. Figure 2 shows the correlation between each socioeconomic measure and the age-adjusted mortality rate. As expected, education, sewage, and water inside the dwellings were negatively associated with the outcome (higher values in each exposure, which imply higher SES, were associated with lower rates). On the other hand, overcrowding, unemployment, and unmet basic needs were positively associated (higher values in each of these exposures, implying lower SES, were associated with higher rates).
The spatial patterning of the different SES measures presented a similar pattern to the spatial distribution of age-adjusted COVID-19 mortality (eFigure 1). The bivariate maps in Fig. 3 simultaneously show the spatial distribution of COVID-19 age-adjusted mortality rates and the two selected SES characteristics (unmet basic needs and overcrowding). The rest of the exposures are available in supplementary material (eFigure 2). We were able to identify two areas of spatial concordance in the city. One area in the southeastern neighborhoods with high mortality and lower SES, and an opposite one, in the core area and western neighborhoods, with lower mortality and higher SES.
Table 2 shows the magnitude of these associations. Education, sewage and water indicators were reverse-coded to facilitate comparison with other indicators. In general, we found that people living in the most deprived areas presented higher mortality rates than those living in the least deprived areas. All variables, except for sewage connection, exhibited a statistically significant relationship with the age-adjusted mortality rate. The SII (difference in mortality rate between the 9th and 1st decile of each variable) ranged from a low of 4.83 (95% CI -0.87 to 10.54) for sewage connection to a high of 13.39 (95% CI 7.17 to 19.60) for completed high school education. The RII (ratio of mortality rates between the 9th and 1st decile of each variable) ranged from a low of 1.31 (95% CI 0.95 to 1.78) for sewage connection to a high of 2.14 (95% CI 1.55 to 2.96) for completed high school education.
Discussion
Summary of the main results
We used georeferenced death records to characterize urban inequalities in COVID-19 mortality in a medium-sized Latin American city. We found important spatial patterning of COVID-19 mortality and also found that this spatial patterning was closely linked to social conditions. People living in the most deprived areas had between 1.3- and 2-times higher COVID-19 death rates than those living in the least deprived areas, depending on the SES characteristic studied. Studies of spatial variation in health outcomes across small areas can shed light on the presence and magnitudes of health inequities [42].
Spatial distribution
We found that COVID-19 mortality rates were not randomly distributed across the city. Spatial modeling provides us with a means to detect and quantify patterns and investigate the relationship between potential exposures with the studied outcome [43]. In addition to identifying strong spatial patterning, we also documented a strong association of spatial differences with social differences. Hence, space becomes a privileged analytical category due to its potential to mediate social dimensions and health, since it is where social relations produce and reproduce [44].
We were able to both identify the spatial patterning of COVID-19 mortality and quantify its relationship with social factors, identifying a clear social gradient: the more deprived the area of residence, the higher the mortality rates. Different measures of SES were included because each one captures a different dimension and may reflect different pathways linking social conditions to COVID-19 although they are all interrelated. Specifically, educational level and unemployment showed the strongest association with COVID-19 mortality rates. Lower educational level and higher unemployment are related to poorer working conditions which may favor virus transmission by impacting the ability to stay home for fear of being fired, the lack of paid sick leave, work in the informal economy, or work in jobs that do not facilitate working from home or social distancing [13, 45]. Unmet basic needs and overcrowding showed a slightly weaker but still substantial association with COVID-19 deaths. These variables proxy precarious living conditions, which make social distancing physically impossible since multiple persons must share the same room. Poor quality housing is associated with higher risks of mortality from COVID-19 and these are more likely to be located in deprived areas and inhabited by people with lower incomes [45].
The consistency with different SES measures linked to COVID-19 deaths lends strength to the idea that social inequality is a key driver of COVID-19 mortality. Each indicator is correlated with many other population attributes and may act through many interrelated mechanisms [46]. Pre-existing social and health conditions have a substantial impact on mortality rates from COVID-19 [4, 45]. People from more disadvantaged social classes have a higher prevalence of chronic conditions (hypertension, cardiovascular diseases, diabetes, etc.) [37, 47] and are therefore more vulnerable to severe COVID-19. Deprivation-related inequalities in the mortality rate from COVID-19 showed a pattern that was similar to that observed for mortality from other causes across Latin American cities and across larger areas within cities [38, 48]. Moreover, as in this study, area-level education for small areas was also the indicator with the strongest association with life expectancy across neighborhoods in CABA and Córdoba, two other Argentinian cities [32, 33].
An interesting finding of our study is the impact that age-adjustment had on the spatial patterning of risks of death. Before age-adjustment, death rates were higher in the least deprived areas of the city, but this changed dramatically after age adjustment. The risk of dying from COVID-19 is strongly related to age, with higher mortality rates at older ages. Consequently, the rates in the two areas may appear different, but this may be due to different age distributions rather than to a difference in the underlying risk of COVID-19. Age-adjustment is necessary to accurately describe differences in risks that are unrelated to the age distribution such as those conferred by social conditions. Social conditions may also interact with age. For example, younger individuals in poorer neighborhoods may be at a substantially higher risk of being exposed to COVID-19 infection than younger individuals in less deprived neighborhoods because of the type of jobs they have access to (informal economy, inability to isolate). In addition, because they may be more likely to reside in multigenerational and overcrowded households, they may be more likely to transmit infection to older relatives thus magnifying transmission and severity with implications for mortality rates [21].
Strengths and limitations
This study has several strengths. First, the availability of geo-referenceable individual-level data on COVID-19 deaths over two years, with less than 2% missing georeferencing. To our knowledge, this is the first study examining COVID-19 data at the neighborhood level in a mid-sized Argentinian city. Second, we were able to conduct direct age-standardization of the outcome, an analysis that is often lacking in COVID-19 research due to the lack of age-specific mortality data. Age-standardization allowed us to uncover a substantially different spatial pattern. The COVID-19 pandemic has highlighted the importance of descriptive epidemiology in responding to serious public health crises, and the challenges in characterizing basic epidemiological metrics [49, 50]. Rigorous descriptive analysis may provide valuable insights into understanding the social patterning of health fundamental to discuss how social organization shapes the distribution of health and disease [46]. Moreover, descriptive studies at small areas reveal disparities in indicators that are obscured when using larger units. This allows for a more precise description of the social reality of an area [44], and the potential identification of vulnerable groups for targeted interventions.
Our study has several limitations. First, small numbers can result in imprecise (statistically unstable) estimates for accurate comparisons [43]. We addressed this limitation by employing a spatial Bayesian framework that smooths mortality rates by borrowing information from age groups and neighboring areas to improve the precision of local estimates, resulting in the smoothing of extreme values [23, 39]. Another limitation was the impossibility of stratifying the analysis by gender due to the low number of deaths, even though gender is an important dimension that would have presented differential risks [51, 52]. Finally, the social characteristics of neighborhoods rely on latest available census (year 2010) and therefore did not align with the period of the pandemic. Unfortunately, the closest census in time (2022) does not yet have available data at the neighborhood level. Oiur analysis therefore relies on the assumption that the relative ranking of social characteristics by neighborhood has remained stable over time. This assumption was studied in previous work for other countries from SALURBAL project, and found that socioeconomic rankings did not change much across time [48].
Policy recommendations
Studies comparing small areas such as neighborhoods can inform placed-based interventions as well as broader policies aimed at addressing structural drivers of health inequities. Identifying areas where inequalities are combined can target resourses to improve the situation of those at higher risk. Systematic geocoding of all deaths would allow cities to describe and track inequities, and to identify possible pathways for developing effective policy actions within and outside the healthcare system, allowing contextualized problems to be identified and quantified. We should use the accumulated knowledge gained throughout the pandemic to improve the health of populations, particularly those living in the most disadvantaged areas.
Conclusions
The COVID-19 pandemic has not only exposed the underlying inequalities in society but has magnified them [45]. The most vulnerable populations, such as those living in urban slums and informal settlements, were the worst affected by the immediate and long-term effects of the pandemic [4]. Deep-rooted inequalities in cities, manifested in overcrowded households and unsafe working conditions, made it almost impossible to comply with the lockdowns and stay-at-home measures.
The spatial differences that we observed are driven by the embodiment of historical and social context [34]. Spatial emplacement and segregation are linked to social status, ethnocultural background, and the interrelated effects of inequality and racism, as it happens in Bariloche and other cities. The associations found in the spatial distribution of mortality from an acute illness, such as COVID-19, are possibly related to a long-term historical process, given that the social processes leading to health differences based on social position, race or inmigration involve multiple pathways and long casual chains [46].
Availability of data and materials
This data cannot be shared openly. The data use agreement does not allow us to share the vital registration data openly. For the rest of the data, we are preparing an open data platform that is planned to be launched in late 2023. The SALURBAL project welcomes queries from anyone interested in learning more about its dataset and potential access to data. To learn more about SALURBAL’s dataset, visit https://drexel.edu/lac/ or contact the project at salurbal@drexel.edu.
References
CEPAL. Mortalidad por COVID-19 y las desigualdades por nivel socioeconómico y por territorio. CEPAL. 2021. Available from: https://www.cepal.org/es/enfoques/mortalidad-covid-19-desigualdades-nivel-socioeconomico-territorio. [Cited 2022 Nov 8].
WHO Coronavirus (COVID-19) Dashboard. Available from: https://covid19.who.int. [Cited 2022 Nov 9].
Bilal U. Has the pandemic exacerbated urban inequality in Latin America? Ciudades Sostenibles. 2022. Available from: https://blogs.iadb.org/ciudades-sostenibles/en/has-the-pandemic-exacerbated-urban-inequality-in-latin-america/. [Cited 2022 Nov 11].
Boza-Kiss B, Pachauri S, Zimm C. Deprivations and Inequities in Cities Viewed Through a Pandemic Lens. Frontiers in Sustainable Cities. 2021;3. Available from: https://www.frontiersin.org/article/10.3389/frsc.2021.645914. [Cited 2022 Jun 7].
United Nations. Policy Brief: COVID-19 in an Urban World - World | ReliefWeb. 2020. Available from: https://reliefweb.int/report/world/policy-brief-covid-19-urban-world. [Cited 2022 Jun 7].
World Bank Group. The World Bank Data. 2022. Urban population (% of total population) - Latin America & Caribbean. Available from: https://data.worldbank.org/indicator/SP.URB.TOTL.IN.ZS?locations=XJ. [Cited 2022 Sept 19].
Ise A, Villalba S, Clementi L, Carrizo S. Extra long Argentinian lockdown: revising the energy regime. Global Transitions. 2021;1(3):43–54.
Macchia A, Ferrante D, Battistella G, Mariani J, de Quirós FGB. COVID-19 among the inhabitants of the slums in the city of Buenos Aires: a population-based study. BMJ Open. 2021;11(1):e044592.
Nwosu CO, Kollamparambil U, Oyenubi A. Socio-economic inequalities in ability to work from home during the coronavirus pandemic. Econ Labour Relat Rev. 2022;33(2):290–307.
Chowkwanyun M, Reed AL. Racial health disparities and Covid-19 - caution and context. N Engl J Med. 2020;383(3):201–3.
Bilal U, Tabb LP, Barber S, Diez Roux AV. Spatial inequities in COVID-19 testing, positivity, confirmed cases, and mortality in 3 U.S. Cities : an ecological study. Ann Intern Med. 2021;174(7):936–44.
Glodeanu A, Gullón P, Bilal U. Social inequalities in mobility during and following the COVID-19 associated lockdown of the Madrid metropolitan area in Spain. Health Place. 2021;1(70):102580.
Marí-Dell’Olmo M, Gotsens M, Pasarín MI, Rodríguez-Sanz M, Artazcoz L, Garcia de Olalla P, et al. Socioeconomic inequalities in COVID-19 in a European Urban area: two waves, two patterns. Int J Environ Res Public Health. 2021;18(3):1256.
Bilal U, Jemmott JB, Schnake-Mahl A, Murphy K, Momplaisir F. Racial/ethnic and neighbourhood social vulnerability disparities in COVID-19 testing positivity, hospitalization, and in-hospital mortality in a large hospital system in Pennsylvania: A prospective study of electronic health records. Lancet Regional Health - Americas. 2022;1(10):100220.
Bilal U, Mullachery PH, Schnake-Mahl A, Rollins H, McCulley E, Kolker J, et al. Heterogeneity in Spatial Inequities in COVID-19 Vaccination Across 16 Large US Cities. Am J Epidemiol. 2022;191(9):1546–56.
Wachtler B, Michalski N, Nowossadeck E, Diercke M, Wahrendorf M, Santos-Hövener C, et al. Socioeconomic inequalities and COVID-19 – A review of the current international literature. J Health Monit. 2020;5(Suppl 7):3–17.
Bermudi PMM, Lorenz C, de Aguiar BS, Failla MA, Barrozo LV, Chiaravalloti-Neto F. Spatiotemporal ecological study of COVID-19 mortality in the city of São Paulo, Brazil: Shifting of the high mortality risk from areas with the best to those with the worst socio-economic conditions. Travel Med Infect Dis. 2021;1(39):101945.
Bilal U, Alfaro T, Vives A. COVID-19 and the worsening of health inequities in Santiago, Chile. Int J Epidemiol. 2021. https://doi.org/10.1093/ije/dyab007.
Arceo-Gomez EO, Campos-Vazquez RM, Esquivel G, Alcaraz E, Martinez LA, Lopez NG. The income gradient in COVID-19 mortality and hospitalisation: An observational study with social security administrative records in Mexico. Lancet Regional Health – Americas. 2022;6. Available from: https://www.thelancet.com/journals/lanam/article/PIIS2667-193X(21)00111-3/fulltext. [Cited 2022 Nov 9].
Marconi AM, Castillo Salgado C, Sarrouf EB, Zamora RJ, Irurzun AM, Islam N. Socioeconomic inequities of COVID-19 mortality in vulnerable Comunas of the City of Buenos Aires. Sci Rep. 2023;13(1):13642.
Chauvin JP, Fowler A, Herrera L. N. The Younger Age Profile of COVID-19 Deaths in Developing Countries. Inter-American Development Bank; 2020. Available from: https://publications.iadb.org/en/node/29388. [Cited 2022 Jul 1].
Demombynes G. COVID-19 Age-Mortality Curves Are Flatter in Developing Countries. Washington, DC: World Bank; 2020. Available from: https://openknowledge.worldbank.org/handle/10986/34028. [Cited 2022 Jul 1].
Moraga P. Small area disease risk estimation and visualization using R. R J. 2018;10(1):495–506.
Kong AY, Zhang X. The use of small area estimates in place-based health research. Am J Public Health. 2020;110(6):829–32.
Quick H, Waller LA, Casper M. Multivariate spatiotemporal modeling of age-specific stroke mortality. Ann Appl Stat. 2017;11(4):2165–77.
Quistberg DA, Diez Roux AV, Bilal U, Moore K, Ortigoza A, Rodriguez DA, et al. Building a Data Platform for Cross-Country Urban Health Studies: the SALURBAL Study. J Urban Health. 2019;96(2):311–37.
Perner MS, Alazraqui M. Bariloche: contrastes ocultos en un paraíso alpino. Desigualdades socio-sanitarias en el espacio urbano. Revista Chilena de Salud Pública. 2021;25(2):139–52.
Rearte A, Moisés MS, Rueda DV, Laurora MA, Marucco AF, Pennini VA, et al. Exceso de mortalidad por todas las causas en el contexto de la pandemia del COVID-19 en Argentina, 2020. Revista Argentina de Salud Pública. 2021;13. Available from: https://rasp.msal.gov.ar/index.php/rasp/article/view/672. [Cited 2022 Jun 17].
INDEC. Unidades Geoestadísticas - Cartografía y códigos geográficos del Sistema Estadístico Nacional. 2023. Available from: https://www.indec.gob.ar/indec/web/Institucional-Indec-Codgeo. [Cited 2023 Aug 16].
INDEC, REDATAM. CEPAL/CELADE. Base de datos - REDATAM. 2015. Available from: https://redatam.indec.gob.ar/argbin/RpWebEngine.exe/PortalAction?BASE=CPV2010A. [Cited 2022 May 16].
INDEC RN. Proyecciones y Estimaciones. 2013. Proyecciones y Estimaciones. Available from: https://estadisticaycensos.rionegro.gov.ar/?contID=55285. [Cited 2022 May 16].
López SR, Tumas N, Bilal U, Moore KA, Acharya B, Quick H, et al. Intraurban socioeconomic inequalities in life expectancy: a population-based cross-sectional analysis in the city of Córdoba, Argentina (2015–2018). BMJ Open. 2022;12(9):e061277.
Trotta A, Bilal U, Acharya B, Harrison Q, Moore K, Perner MS, et al. Spatial inequities in life expectancy in small areas of Buenos Aires, Argentina 2015–2017. J Urban Health. 2023;100:577–90.
Ahmad OB, Boschi Pinto C, Lopez A, Murray C, Lozano R, Inoue M. Age Standardization of Rates: A New WHO Standard. 2001;(31). Available from: https://www.researchgate.net/publication/238744905_Ahmad_OB_Boschi-Pinto_C_Lopez_AD_et_al_2000_Age_Standardization_of_Rates_A_New_WHO_Standard_GPE_Discussion_Paper_Series_No_31_World_Health
Feres J, Mancero X. El método de las necesidades básicas insatisfechas (NBI) y sus aplicaciones en América Latina. 2001. p. 7.
INDEC. Necesidades Básicas Insatisfechas (Version ampliada con datos departamentales). DINREP; 2014. Available from: http://www2.mecon.gov.ar/hacienda/dinrep/Informes/archivos/NBIAmpliado.pdf.
Perner MS, Alazraqui M, Amorim LDAF. Social inequalities between neighborhoods and cardiovascular disease: a multilevel analysis in a Latin American city. Cien Saude Colet. 2022. Available from: http://www.cienciaesaudecoletiva.com.br/en/articles/desigualdades-sociais-entre-vizinhancas-e-doencas-cardiovasculares-analise-multinivel-em-uma-cidade-latinoamericana/18262. [Cited 2022 May 5].
Bilal U, Hessel P, Perez-Ferrer C, Michael YL, Alfaro T, Tenorio-Mucha J, et al. Life expectancy and mortality in 363 cities of Latin America. Nat Med. 2021;27(3):463–70.
Blangiardo M, Cameletti M. Spatial and Spatio-temporal Bayesian Models with R-INLA. United Kingdom: Wiley; 2015.
Sturtz S, Ligges U, Gelman A. R2WinBUGS : A Package for Running WinBUGS from R. J Stat Soft. 2005;12(3). Available from: http://www.jstatsoft.org/v12/i03/. [Cited 2022 Nov 16].
Rubin DB. Multiple Imputation for Nonresponse in Surveys | Wiley. 2004. 258 p. Available from: https://www.wiley.com/en-us/Multiple+Imputation+for+Nonresponse+in+Surveys-p-9780471655749. [Cited 2022 Nov 10].
Diez Roux AV, Mair C. Neighborhoods and health. Ann N Y Acad Sci. 2010;1186:125–45.
Waller L, Gotway C. Applied Spatial Statistics for Public Health Data. New Jersey: John Wiley & Sons, Ltd; 2004.
Aquino R, Gouveia N, Texeira MG, Costa M da C, Barreto ML. Estudos Ecológicos Desenho de Dados Agregados. In: Epidemiologia Saúde Fundamentos Metódos Aplicações. Rio de Janeiro: Guanabara Koogan; 2012. p. 175–85.
Marmot M, Allen J, Goldblatt P, Herd JME. Build Back Fairer: The COVID-19 Marmot Review. The Pandemic, Socioeconomic and Health Inequalities in England. Insitutte of Health Equity; 2020. Available from: https://www.health.org.uk/publications/build-back-fairer-the-covid-19-marmot-review. [Cited 2022 Jun 9].
Diez Roux AV. Social epidemiology: past, present, and future. Annu Rev Public Health. 2022;43:79–98.
Tumas N, Rodríguez López S, Bilal U, Ortigoza AF, Diez Roux AV. Urban social determinants of non-communicable diseases risk factors in Argentina. Health Place. 2021;29:102611.
Bilal U, Alazraqui M, Caiaffa WT, Lopez-Olmedo N, Martinez-Folgar K, Miranda JJ, et al. Inequalities in life expectancy in six large Latin American cities from the SALURBAL study: an ecological analysis. Lancet Planetary Health. 2019;3(12):e503–10.
Fox MP, Murray EJ, Lesko CR, Sealy-Jefferson S. On the need to revitalize descriptive epidemiology. Am J Epidemiol. 2022;191(7):1174–9.
Diez Roux AV. Population Health in the Time of COVID-19: Confirmations and Revelations. Milbank Q. 2020;98(3):629–40.
Mi J, Zhong W, Huang C, Zhang W, Tan L, Ding L. Gender, age and comorbidities as the main prognostic factors in patients with COVID-19 pneumonia. Am J Transl Res. 2020;12(10):6537–48.
Doerre A, Doblhammer G. The influence of gender on COVID-19 infections and mortality in Germany: Insights from age- and gender-specific modeling of contact rates, infections, and deaths in the early phase of the pandemic. PLoS ONE. 2022;17(5):e0268119.
Acknowledgments
We want to acknowledge the healthcare workers from Hospital Zonal Bariloche. SALURBAL acknowledges the contributions of many different agencies in generating, processing, facilitating access to data, or assisting with other aspects of the project.
Funding
This research was funded in part by the Wellcome Trust [205177/Z/16/Z]. UB was supported by the Office of the Director of the National Institutes of Health under award number DP5OD26429. For the purpose of Open Access, the author has applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.
Author information
Authors and Affiliations
Contributions
M.S.P. wrote the main manuscript text. M.A. and A.T. contributed to the study conception, data analysis, interpretation, writing, and critical review. B.A. contributed with the methodological approach, data analysis, prepared the figures and helped with the interpretation. H.Q. contributed with the methodological approach and interpretation. U.B. and A.D.R. contributed to the study design, interpretation, writing, and critical review. RB contributed with the interpretation and critical review. All authors approved the final version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Data used for this study was aggregated to ensure data confidentiality required in Argentina and the data agreement policy of the SALURBAL project.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 2: eFigure 1.
Spatial distribution of socio-economic variables in Bariloche, Argentina. Note: Data from 2010 census. Legends show the quartiles of the variables (in percentages). Definition of each of these variables are provided in the exposure subsection of the methods section in the main manuscript. In purple scale maps, higher percentages indicates a better socioeconomic characteristics. In orange scale maps, higher percentages indicates a worse socioeconomic characteristics.
Additional file 3: eFigure 2.
Spatial distribution of COVID-19 age-adjusted mortality rates and socioeconomic variables in Bariloche, Argentina. A first set of bivariate maps with the spatial distribution of COVID-19 age-adjusted mortality rates with the two measures of education (school attendance and completed high school), piped water, or sewage. These maps are in a purple/light blue/green range, where darker the light blue, better the SES (higher education, piped water, or sewage), and darker the green, higher the mortality rates. In the bivariate map, the extremes of high concordance areas are dark-green and light blue. Dark green represents the most unequal areas (high/low in the legend), with high COVID-19 age-adjusted mortality rates and more deprived areas: low percentage of education, piped water, and sewage. Light-blue areas (low/high in the legend), the other extreme, low COVID-19 age-adjusted mortality rates and less deprived areas: high percentage of education, piped water, and sewage. In the botton, a map, with the spatial distribution of COVID-19 age-adjusted mortality rates with unemployment, same colours and interpretation as Figure 3. These maps is in a range of brown/light blue/orange colors, where the darker the light blue the worse the SES (higher the unemployment), and the darker the orange, the higher the mortality rate. The bivariate pallet shows areas of high concordance in brown and in light orange. Brown (high/high in the legend) represent the more unequal areas, high COVID-19 age-adjusted mortality rates and more deprived areas (high unemployment); light orange in the other extreme (low/low in the legend), low COVID-19 age-adjusted mortality rates and less deprived areas (low unemployment).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Perner, M.S., Trotta, A., Bilal, U. et al. Social inequalities and COVID-19 mortality between neighborhoods of Bariloche city, Argentina. Int J Equity Health 22, 198 (2023). https://doi.org/10.1186/s12939-023-02019-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12939-023-02019-w