key: cord-0928518-yza7ljlp authors: Baser, Onur title: Population density index and its use for distribution of Covid-19: A case study using Turkish data date: 2020-10-13 journal: Health Policy DOI: 10.1016/j.healthpol.2020.10.003 sha: 639d3896b95533fba6289f131f438706af4dbb72 doc_id: 928518 cord_uid: yza7ljlp Since March 2020, many countries around the world have been experiencing a large outbreak of a novel coronavirus (2019-nCoV). Because there is a higher rate of contact between humans in cities with higher population weighted densities, Covid-19 spreads faster in these areas. In this study, we examined the relationship between population weighted density and the spread of Covid-19. Using data from Turkey, we calculated the elasticity of Covid-19 spread with respect to population weighted density to be 0.67 after controlling for other factors. In addition to the density, the proportion of people over 65, the per capita GDP, and the number of total health care workers in each city positively contributed to the case numbers, while education level and temperature had a negative effect. We suggested a policy measure on how to transfer health care workers from different areas to the areas with a possibility of wide spread. An ongoing outbreak of a novel coronavirus (2019-nCov) was identified only a few days after the World Health Organization (WHO) was alerted about a cluster of pneumonia of unknown aetiology in the city of Wuhan, China on 31 December 2019. [1] The outbreak appears to have started from a single zoonotic transmission event or multiple zoonotic transmission events at a wet market in Wuhan where game animals and meat were sold [2] and is quickly approaching 25 million confirmed cases worldwide. [3] . The start date for the epidemic, total cases, and fatalities were different for each country. The country with the highest number of cases as of August 28, 2020 is the United States, with more than 6 million confirmed cases, followed by Brazil and India. There were more than 3.5 million confirmed cases in Europe. In all of these regions, crowded cities are the epicenters for the disease. [4] What sets these cities apart from rural areas are their high population densities. Population weighted density can be described as a weighted average of density across the tracts, where tracts are not weighted by land area but by population. [5] Density is one of the most fundamental characteristics of an urban area. [6] However, raw population density, simply population divided by count, is not a good measure of the density at which the population lives. [ Angeles is actually denser than New York, but it is hemmed in by mountains, limiting how far the commuting zone can reach. However, according to population weighted density, an average New Yorker lives in a census tract with more than 12,400 people within a kilometer square. That is three times more than the density of Los Angeles County. [7] The population weighted density of several of European cities, such as Barcelona (24, 600) , Madrid (18, 600) , Valencia (17, 300) and Paris (13, 300) are much higher than New York City. The population weighted density of Rome (8, 900) , Berlin (8, 200) and London(8,000) is also relatively high. [8] Especially in pandemies such as coronavirus, where human contact is the main reason for spread, populationweighted densities are better measure than conventional densities, because the variation in density across the subareas matters more than the density in total area. Before New York City placed restrictions on its residents in order to combat disease spread, the number of cases in New York City was close to 20 times the number of cases in The aim for this paper was first to derive population weighted density for the cities in Turkey and the districts for its major three cities in 2020. Then using the data points in April 2020, we analyzed the relationship with the density and the spread of coronavirus in those cities controlling for cities' education level, wealth, health care J o u r n a l P r e -p r o o f force, temperature, and demographics. We suggest that, in conjunction with information about a cities' number of health care workers and fatality statistics in August 2020, population weighted density can also inform us on how to mobilize each region health care force from low spread risk areas to high spread risk areas. Let D be the density of the urban area, which is the total population, P divided by the total area A: = Let is the population and is the area of subareas, by definition P = ∑ and A = ∑ . Therefore the density for each areas is = . Population-weighted density is the mean of the subareas densities weighted by the population of the subareas: Ottensmann [9] showed that the difference between population-weighted density and conventional density is a simple function of the variance in density across the census subareas and conventional density. Craig [5] suggested the amount of differences will depend on the variation in density across the subareas. We would expect similar results for the areas that have been defined in such a ways that they do not include sparsely settled territory. In the USA, this measure has been part of the national statistics since 2010, but it has not been used yet in Turkey. By using population values from the Turkish Statistical Institution and area values from several websites that use Google Earth, we calculated population weighted density for each city in Turkey. [10] Istanbul, with a population of more than 15 million, was ranked as the first city according to population weighted density. On average, residents in Istanbul live with 16, 757 people around their 1 kilometer square. Istanbul population weighted density is lower than Barcelona, Madrid and Valencia. However, it is higher than Paris population density and almost double that of London. Figure 1 shows a map of Turkey with respect to the population weighted density of each of its regions. Izmir was the second dense population in Turkey according to population weighted index, although it is in third J o u r n a l P r e -p r o o f place when it comes to population and raw density. Ankara was the third dense population, although according to raw density it was as eighth. (Table A1 ) For the three major cities, we also calculated population density of the each district (Table A2, Table A3, Table A4 ) as well as their corresponding population weighted density maps. For our knowledge, this was the first time the population weighted densities are calculated for Turkey and these three cities. ( Figure A1 ) Turkey's Health Ministry has released only limited data on the spread of the virus and announced the number of Covid-19 cases in individual cities on two occasions, on April 1st and April 4th. [11] As of August 28, 2020, there were more than 265,000 cases in Turkey and Istanbul accounted for almost 60 percent of confirmed cases of Covid-19. Izmir and Ankara have been declared growing hot spots, as was predicted by population weighted density. For each city, we calculated the difference in the case numbers to proxy for the spread of the disease in each city. There was a strong positive relationship between population weighted density and and the spread of the disease. ( Figure 2 ) Correlation coefficient was calculated as 0.97 with p − value < 0.0001. To determine the relationship between these variables and our outcome variable, corona spread, we consider the following model: where j indexes m cities. All of our variables are at the city level. The spread, our outcomes variable, defined as log of differences in case numbers in a given city measured two different time period. However, note that estimating the model on individuals and clustering standard errors by city would yield the same coefficients and standard errors as estimating city means using analytic weights and standard errors robust to heteroskedasticity. where ̃ and ̃ obtained by multiplying each row of [n × (k + 1)] matrix X and row of [n × 1] matrix y by √ , being the number of individuals contributing to the average. For standard errors, the variance-covariance matrix can be calculated as The variables such as population weighted density, total health care workers, and city's per capita GDP in logarithmic form allowed us to measure elasticity. The other explanatory variables were used as level forms, thus providing semi-elasticity measures. On average, a person in Turkey lives with 3,868 people within 1 kilometer square. The average education level is around 7.5 years and 9.12% of the population is 65 years old and over. There is a slightly higher male population than female population and the per capita income is $9,745. There are about 1.9 doctors, 2.37 nurses and 2.22 other health care workers per 1000 people in Turkey. (Table 1 ) Pandemics spread through the movement of and interaction between infected people, and these interactions occur more frequently in places with high population weighted densities. Therefore, it has been assumed that during Figure 3 : The Expected Percentage increase in district cases in Istanbul relative to Istanbul average. pandemics such as Covid-19, density is associated with higher rates of transmission, infection, and mortality. [14, 15] After controlling for other factors, weighted regression yield that the elasticity of population weighted density with respect to the growth of coronaspread is calculated as 0.67. (Table 2) 1% increase in population-weighted-density increased the growth of the disease spread by .67%. For each district in Istanbul, we measured percentage changes in expected cases relative to the district with the average population density. We choose Besiktas district as a reference since population weighted density for district is approximately equal to the density of Istanbul as a whole. For example, Avcilar district's population weighted density is around 10% higher than Kadikoy, so the growth of Covid-19 cases for Avcilar will be 6.7% higher than Kadikoy. (Figure 3 ) As expected, there were significant associations between socioeconomic factors and the growth of the spread. As is consistent with the literature, cities with higher-than-average education levels had significantly lower Covid-19 infection rates. This trend can be explained by the facts that people with higher education levels have a better understanding of the virus and take shelter-in-place restrictions more seriously. They are also more likely to be able to work from home. [16] The proportion of the population aged 65 years and older was also positively associated with the growth rate in cases. Each percentage point increased the growth rate by 11%. This may be due to the fact that people aged 65 years and older have weaker immune systems than the rest of the population. [17] Since the higher health care workers in the city is related with the higher number of testing, we found strong positive correlation between total number of health care workers and the spread of the disease. Cities with 1% higher health care workers were associated with .84% higher disease spread. These findings are consistent with previous studies in the US that found statewide testing is the most significant predictor of the county infection In epidemiological terms, "flattening the curve"refers to the implementation of measures that slow the rate at which people are infected by the virus, thus lessening the burden on medical professionals and the health care system. The "curve" refers to the projected number of people who will come into contact with COVID-19 over a period of time. As more people contract the virus, the infection curve rises. If it rises too quickly, then the health care system risks becoming overloaded, which can lead to hospitals running out of the supplies they need to help infected people recover. How many people a given patient is likely to infect is defined by the reproductive number. Decreasing this number is the ultimate goal in fighting the pandemic. If it is less than one, then group of infected people would be generating For any health care system to perform well it depends on the availability of a sufficient number of skilled health care workers. Furthermore, it is crucial that these health workers be mobile, since urban areas are first to experience spikes in Covid-19 cases. For example, in the first months of the Covid-19 epidemic, health care workers flew from San Francisco to New York City, where there was a higher need for services. [26] If such mobility is available, we suggest that attention to population weighted density can direct how this transfer might be done. We graphed the cities in Turkey with low risk spread (proxied by low population density) with a high number of health care workers. (Figure 4) The big green cities like Yozgat, Sivas, Tokat are the cities with low population weighted density with relatively high number of health care workers, and the brown small cities like Istanbul, Batman, Kocaeli and Yalova are the ones with high population weighted density with relatively low number health care workers. Possible transfer of health care workers can be done from cities marked by big green bubbles to cities marked by small brown bubbles. Note that healthcare workers are one major resource in the health systems, but for treating patients with Covid it would not be enough to shift workers without providing/shifting other care resources such as beds, reanimation equipment. Population weighted density can also be used as a tool to evaluate the success of fighting with pandemic. We compiled the death data from 8 different cities in Turkey. These were the only 8 cities that make the death data available on a daily basis online. [27] We graphed the series of daily death since March 25, 2020 (two weeks after J o u r n a l P r e -p r o o f For eight cities we selected, we calculated the increase amount of death relative to average of last three years. We then divided this number with population weighted density to determine, controlling for the density, the proportion death per density measure. (Table 3 ) Consistent with previous research, cities with higher population densities tend to have lower death rates, possibly they enjoyed a higher level of development including better health care systems. The recent analysis found that after controlling for factors such as metropolitan size, education, race and age, doubling the activity density was associated with an 11.3 percent lower death rate. [18] The impact of population density on the rate of spread of emerging highly contagious infectious diseases has rarely been studied. The current Covid-19 pandemic allows us to investigate these relationships. Our study uses a regression model to study the impact of population weighted density on Covid-19 spread in Turkey while controlling for key compounding. We found that population weighted density is one of the most significant predictors of infection rates. However, counties with higher densities have significantly lower virus-related mortality rates, possibly due to superior health care systems. J o u r n a l P r e -p r o o f Most of the data is preliminary during the outbreak and the data that we can use in Turkey is limited. We were able to update our analysis for the death statistics by August 30th. However, the Turkish ministry of health has not published case numbers by cities after April 2, 2020. Because of this, the results in Table 2 have not been updated beyond April 2, 2020. However, our results were consistent with the previous research. [14, [16] [17] [18] [19] [20] [21] Additionally, we are unable to measure the number of patients with underlying conditions in each city, which would decrease the survival rate of each of these regions. We have seen death toll raise in a cities like Zonguldak where most of the population suffers from chronic lung diseases due to work conditions in the minings located in the city. We have used the death statistics to measure the success of dealing with the outbreak, but any death statistics in the midst of pandemic are tricky to pin down and must be considered preliminary. We have seen most of the countries are improving their death statistics, which they now acknowledge incomplete. In order to control and manage the outbreak, our analysis suggests that population weighted density can be a useful tool. According to previous research, there are many advantages of compact development. It is associated with with open space preservation, higher innovation and overall economic productivity, greater social capital, less likelihood of obesity and related chronic diseases and increased overall life expectancy. [28] [29] [30] [31] [32] . However, compact development can be big "enemy" in the coronavirus fight. [15] High density means that people in those areas live very differently from other people. Those who live or work in or near the city shop and commute differently: they are far more likely to walk or take public transit than the rest of the people. The disease spreads faster in the areas with high population weighted areas then elsewhere simply because there is so much human contact. At this stage, particular attention should be given to the prevention of spreading in the highest dense areas directed by population weighted density. The concept of "flattening the curve" ultimately assumes that the same number of people will contract the coronavirus whether or not the curve is steep or flattened. If the curve is flattened, however, there is less stress placed upon the health care system, which results in better health care access for those who are sick. Acknowledgments: None J o u r n a l P r e -p r o o f J o u r n a l P r e -p r o o f The role of superspreading in Middle East respiratory syndrome coronavirus (MERS-CoV) transmission A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster Novel Coronavirus (2019-nCoV) situation report 9 Averaging population density Density and creativity in us regions Density Comparing the densities of australian, european, canadian, and new Zealand cities On population-weighted density Turkish statistical instituion Turkish meteorolgy instituion Econometric analysis of cross section and panel data mit press Spatial distribution and trend analysis of current status of covid-19 in nepal and global future preventive perspectives Density is new york city big "enemy" in the coronavirus fight The social determinants of health and pandemic h1n1 2009 influenza severity Clinical features of covid-19 in elderly patients: A comparison with young and middle-aged patients Does density aggravate the covid-19 pandemic? early findings and lessons for planners High temperature and high humidity reduce the transmission of covid-19 Rich at risk: socio-economic drivers of covid-19 pandemic spread Gender differences in patients with covid-19: Focus on severity and mortality Calculating virus spread The real time world statistics Emerged ha and na mutants of the pandemic influenza h1n1 viruses with increasing epidemiological significance in taipei and kaohsiung, taiwan Economic effects of the 1918 influenza pandemic These san francisco doctors flew to new york to fight the coronavirus and they have a warning for us E-devlet turkey Compactness versus sprawl: A review of recent evidence from the united states The effectiveness of urban containment regimes in reducing exurban sprawl Evidence of the impacts of urban sprawl on social capital Relationship between urban sprawl and physical activity, obesity, and morbidity-update and refinement Associations between urban sprawl and life expectancy in the united states