Skip to main content


Mapping child maltreatment risk: a 12-year spatio-temporal analysis of neighborhood influences

Article metrics



‘Place’ matters in understanding prevalence variations and inequalities in child maltreatment risk. However, most studies examining ecological variations in child maltreatment risk fail to take into account the implications of the spatial and temporal dimensions of neighborhoods. In this study, we conduct a high-resolution small-area study to analyze the influence of neighborhood characteristics on the spatio-temporal epidemiology of child maltreatment risk.


We conducted a 12-year (2004–2015) small-area Bayesian spatio-temporal epidemiological study with all families with child maltreatment protection measures in the city of Valencia, Spain. As neighborhood units, we used 552 census block groups. Cases were geocoded using the family address. Neighborhood-level characteristics analyzed included three indicators of neighborhood disadvantage—neighborhood economic status, neighborhood education level, and levels of policing activity—, immigrant concentration, and residential instability. Bayesian spatio-temporal modelling and disease mapping methods were used to provide area-specific risk estimations.


Results from a spatio-temporal autoregressive model showed that neighborhoods with low levels of economic and educational status, with high levels of policing activity, and high immigrant concentration had higher levels of substantiated child maltreatment risk. Disease mapping methods were used to analyze areas of excess risk. Results showed chronic spatial patterns of high child maltreatment risk during the years analyzed, as well as stability over time in areas of low risk. Areas with increased or decreased child maltreatment risk over the years were also observed.


A spatio-temporal epidemiological approach to study the geographical patterns, trends over time, and the contextual determinants of child maltreatment risk can provide a useful method to inform policy and action. This method can offer a more accurate description of the problem, and help to inform more localized prevention and intervention strategies. This new approach can also contribute to an improved epidemiological surveillance system to detect ecological variations in risk, and to assess the effectiveness of the initiatives to reduce this risk.


Child maltreatment is a major social, public health, and human rights problem, with severe, far-reaching, and long-lasting consequences. Its impact on victims’ physical, mental, and reproductive health, behavioral problems, or education attainment; its role in the intergenerational transmission of violence; and its elevated costs for the criminal, health, and social welfare systems, poses a high burden on society [1,2,3,4]. Child maltreatment is a global phenomenon, and its prevalence in high-income countries remains high and is considered a leading cause of health inequality and social injustice [4,5,6,7]. From a public health perspective, child maltreatment is, however, a preventable problem as potentially risk-modifying factors can be identified and targeted in preventive interventions.

Child maltreatment research has typically focused on individual or family risk factors, but ‘place’ also matters in understanding prevalence variations and inequalities in child maltreatment risk [8,9,10,11,12]. A growing body of research is increasingly recognizing the importance of the context in which families live, linking neighborhood characteristics and processes—such as poverty, disorder and crime, immigrant concentration, social impoverishment, or diminished social control—to child maltreatment [13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]. However, studying neighborhood influences on child maltreatment presents important challenges, and there are still some shortcomings in the available literature, which we address in the present study.

Child maltreatment and neighborhood risk factors are not equally distributed spatially, and it is to some extent surprising that the use of spatial analysis techniques and disease mapping methods to analyze geographical patterns of child maltreatment risk, and whether these patterns are associated with neighborhood-level explanatory variables, has been almost non-existing. In this regard, and despite the substantial body of research showing an association between neighborhood characteristics and child maltreatment, most studies examining neighborhood influences fail to take into account the implications of the spatial dynamics of neighborhoods. For example, important issues such as the similarity and influence between neighboring areas, are not appropriately addressed in most studies analyzing neighborhood influences on child maltreatment [10, 11, 25]. Neighborhood risk factors that have been associated with child maltreatment are usually clustered in space, and therefore a spatial analytical approach is particularly appropriate for the study of their influences on the spatial variations of child maltreatment. Furthermore, when analyzing neighborhood influences on ecological variations in child maltreatment risk, the temporal dimension must also be taken into account. Neighborhoods characteristics may change over time and, therefore, their influence on child maltreatment risk can also change over time [10]. In this regard, adding the temporal dimension is a key in identifying and tracking trends in child maltreatment risk—for example, stable high or low risk areas, or areas of increasing or decreasing risk over time. Finally, most studies exploring neighborhood influences on child maltreatment do not provide small-area-specific risk estimations, which limits their relevance to inform more localized intervention and prevention strategies. A new generation of ecological studies needs to take into account the spatial and temporal dimensions to better understand small-area variations in child maltreatment risk.

Bayesian spatio-temporal modeling provides an adequate methodological framework to overcome the above limitations in studying neighborhood influences on child maltreatment risk variations. This analytical approach, allows to incorporate geographical and temporal information to provide more reliable area-specific risk estimates than other non-Bayesian methods, by addressing important methodological issues such as modeling small area counts, spatial auto-correlation, or overdispersion, that can bias estimates if ignored [33,34,35,36,37,38,39]. Although Bayesian spatio-temporal disease mapping is common in other public health and epidemiological areas, this approach is relatively new in the area of family violence [40,41,42,43,44]. So far, only a handful of studies conducted in the US have addressed the spatial and temporal dimensions to study neighborhood influences on child maltreatment with appropriate methodologies [40, 42, 44]. However, these studies are usually low-resolution ones, using larger geographical areas, such as counties, zip codes, or census track, which somehow limits their potential to inform more localized interventions. On the other hand, high-resolution studies using small-area geographical units offer a finer neighborhood characterization, and provide specific risk estimations for small areas, which increases their potential to inform highly localized policies targeting high-risk areas [36, 38, 39].

Based on a social disorganization theoretical framework [43, 45,46,47,48] in this study we analyze the influence of a set of neighborhood characteristics on local patterns of substantiated child maltreatment over a 12-year period at the small-area level. Neighborhood-level characteristics analyzed include three indicators of neighborhood concentrated disadvantage—neighborhood economic status, neighborhood education level, and levels of policing activity—, immigrant concentration, and residential instability. As far as we are aware, the present study is the first to conduct a high-resolution small-area study on the spatio-temporal epidemiology of child maltreatment risk using Bayesian spatio-temporal modeling and disease mapping methods. Our study is conducted in a European city which also provides a ground for comparison to US cities where most of the research on neighborhood influences on child maltreatment has been conducted.



The study was conducted in the city of Valencia, the third largest city of Spain with a population of 736,580. As neighborhood proxies, we used the 552 census block groups into which the city is divided. This was the minimum administrative unit available, with an average of 1334 residents, and ranging from 630 to 2845. To capture temporal trends, data for 12 years were used, from 2004 to 2015. Data for cases of substantiated child maltreatment were provided by Valencia Child Protection Services. Covariates for each census block group were provided by the Valencia Statistics Office and the Valencia Police Department.

Outcome variable

Number of families with child maltreatment protection measures. Data of all families with child protection measures from 2004 to 2015 were collected—no computerized and systematized data were available before 2004. The data we used correspond to all official cases of substantiated child maltreatment in the city of Valencia during the period of the study, to which a child protection measure is associated after the maltreatment and the risk for the child is established. Child maltreatment refer to any type of child maltreatment, including physical, psychological, or sexual abuse, as well as neglect. Perpetrators were parents or legal tutors–child abuse by other parties, such as non-related adults, are dealt by other agencies and are not considered in this study. The data provided by the Child Protection Services did not allow to distinguish between child maltreatment types or perpetrators. Protection measures are issued by the Child Protection Services for all substantiated cases of child maltreatment, and may include a range of measures such as home visiting, family support programs, family or residential foster care, or adoption. Data for this study were 1799 families with child maltreatment protection measures. To avoid data dependency, cases were ‘unique’ families, meaning that a family was only included the first time they received a child maltreatment protection measure. Data were geocoded using the family address, and we counted the cases in each census block group for each of the 12 years in the study period.


Neighborhood concentrated disadvantage. Economic status: Neighborhood economic status was measured using the average cadastral property value. This value is set by the City Hall and is used to establish local taxes. Education level: The average education level of neighborhood residents in each census block group was measured on a 4-point scale, where 1 = less than primary education, 2 = primary education, 3 = secondary education, 4 = college education. Policing activity: An indirect measure of public disorder and crime was measured through police officers’ assessments of their policing activity. Senior police officers provided an index of policing activity that included interventions in drug-related crime, public disorder such as drunkenness and fights, vandalism, homeless people and truancy. The index was structured as a 5-item scale, and each item ranged from 0—very low level of interventions—to 4—very high level of interventions. This measure of policing activity has been associated in previous studies with other types of family violence such as intimate partner violence [43], as well as with a number of neighborhood-level characteristics, such as neighborhood disorder, low socioeconomic status, and high immigrant concentration [49, 50]. Cronbach’s alpha was .74.

Immigrant concentration: Using census data, immigrant concentration was measured as the percentage of immigrant population in each census block group.

Residential instability: An index of residential mobility, based on census data, was used as the proportion—rate per 1000 inhabitants—of the population who had moved into or out of each census block group during the previous year—for example, residential instability for 2015 captures all population movements occurred in 2014.

Statistical analysis

The outcome variable was the number of families with child maltreatment protection measures, corresponding to all substantiated cases of child maltreatment in the city during the period of the study. We use, therefore, a conditionally independent Poisson distribution based on the count of families in each census block group in the 12 years of the study:

$$y_{it} |\eta_{it} \sim\;Po\left( {E_{it} \exp (\eta_{it} )} \right), \quad i = 1, \ldots ,552 \quad t = 1, \ldots ,12$$

where \(E_{it}\) is a fixed quantity that accounts for the expected number of families with child maltreatment protection measures, in proportion to the total number of families, in census block group \(i\) in year t; \(\eta_{it}\) is the log-relative risk for every area and year.

We used different models for \(\eta_{it}\) with an increasing level of complexity, from a Poisson regression model to a spatio-temporal autoregressive model. First, model 1 only included all covariates—that is, economic status, education level, policing activity, residential instability, and immigrant concentration.

$$\eta_{it} = \mu + X_{it} \beta$$

where \(\mu\) is the intercept, \(X_{it}\) is the vector of covariates, and \(\beta\) is a vector of regression coefficients.

Model 2 was specified as a spatial model; it included all covariates and added unstructured and structured random effects. The unstructured random effect accounted for spatial heterogeneity or overdispersion, while the structured random effect accounted for the spatial effect:

$$\eta_{it} = \mu + X_{it} \beta + \phi_{i} + \theta_{i}$$

where \(\phi_{i}\) represents the spatially structured term, and \(\theta_{i}\) the spatially unstructured term.

Model 3 included the previous terms and incorporated an unstructured temporal effect:

$$\eta_{it} = \mu + X_{it} \beta + \phi_{i} + \theta_{i} + \alpha_{t}$$

where \(\alpha_{t}\) accounts for the temporal heterogeneity. This model, however, does not account for past cases of child maltreatment—that is, temporal dependency.

Finally, model 4 included a spatio-temporal effect. To this end, we followed an autoregressive approach [51], combining autoregressive time series and spatial modeling. We defined a spatio-temporal structure in which the relative risks are both spatially and temporally dependent.

$$\eta_{i1} = \mu + X_{it} \beta + \alpha_{1} + \left( {1 - \rho^{2} } \right)^{ - 1/2} \cdot\left( {\phi_{i1} + \theta_{i1} } \right)$$
$$\eta_{it} = \mu + X_{it} \beta + \alpha_{t} + \rho \cdot\left( {\eta_{{i\left( {t - 1} \right)}} - \mu - \alpha_{t - 1} } \right) + \phi_{it} + \theta_{it}$$

The first equation defines the log-relative risk for the first year observed (2004) and the second equation defines the log-relative risk for the following years. In both, \(\alpha_{t}\) is the mean deviation of the risk in the year t, \(\rho\) represents the temporal correlation between years, and \(\phi_{it}\) and \(\theta_{it}\) refer to structured and unstructured spatial random effects, respectively.

Models were specified following a Bayesian approach. Therefore, we assigned appropriate prior distributions for all parameters. We assigned vague Gaussian distributions for the fixed effects \(\beta\); \(\mu\) was specified as an improper uniform distribution; unstructured effects were modeled as a normal distribution \(N\left( {0,\sigma^{2} } \right)\) in the different models (\(\theta\) and \(\alpha\)). Structured effects (\(\phi\)) were specified by a conditional spatial autoregressive (CAR) model [52] defined as follows:

$$\phi_{i} |\phi_{ - i} \sim\;N\left( {\frac{1}{{n_{i} }}\mathop \sum \limits_{j\sim\;i} \phi_{j} ,\frac{{\sigma_{\phi }^{2} }}{{n_{i} }}} \right)$$

where \(n_{i}\) is the number of neighboring areas of each census block group i, \(\phi_{ - i}\) represents the values of the \(\phi\) vector except the component i, \(\sigma_{\phi }\) is the standard deviation parameter, and \(j\sim\;i\) indicates all units j that are neighboring areas of census block group i. Finally, and following the structure of the hierarchical Bayesian models, hyperparameters \(\sigma\) were specified by uniform distributions \(U\left( {0,1} \right)\) in the models.

Bayesian estimations were performed using Markov Chain Monte Carlo simulation techniques with the software R and the WinBUGS package. 100,000 iterations were generated, discarding the first 10,000 as a burn-in period. Models were compared by the Deviance Information Criterion (DIC) [53]. This measure of fit assumes that models with smaller DIC should be preferred to models with larger DIC. Following this criterion, the model with smaller DIC was chosen. Differences in DIC between 5 and 10 are considered substantial; whereas differences of more than 10 units clearly indicate that the model with the higher DIC should be ruled out.

To ensure robustness of the results, we checked convergence with the convergence diagnosis \(\hat{R}\) [54] which was near to 1.0 for all parameters. A sensitivity analysis was also performed on prior distributions of hyperparameters, with consistent results.


Table 1 summarizes the descriptive statistics of the variables in the study. The neighborhood-level economic status, based on the average cadastral property value, had a mean of 26,320 € [standard deviation (SD) = 13,046], with wide variation across neighborhoods, ranging from 7943€ to 98,560€. The average neighborhood education level, corresponded to secondary education (Mean = 3; SD = .33). Policing activity had a mean of 7.16 (SD = 3.99), again with wide variations across neighborhoods ranging from 0 to 19. The mean of neighborhood residential instability was 200 (SD = 65.96) meaning that, in average, 200 people moved into or out of a census block group in a specific year. Neighborhood immigrant concentration had a mean of 13.3%, ranging from just 1 to 51%, which means that in some neighborhoods over half of the population were immigrant. Finally, the outcome variable, families with child protection measures, ranged from 0 cases of substantiated child maltreatment in some neighborhoods to a maximum of 7 families per census block group in a single year with child protection measures after child maltreatment was substantiated.

Table 1 Variables (mean, standard deviation, minimum and maximum values) at the census block group and year level

After conducting the four Bayesian Poisson regression models, we analyzed the DIC values (Table 2). Model 1, which included only the covariates, showed the worst fit (DIC = 8517.9). Once we introduced both structured and unstructured spatial random effects in model 2, the DIC decreased significantly (DIC = 8164.9). In model 3, which included an unstructured temporal effect, the DIC slightly increased to 8166.8. Finally, model 4, an autoregressive model, despite being the most complex had the lowest DIC value (8126.1), 38 units lower than model 2, and was therefore chosen as the final model. The sign of the covariate estimations (positive or negative) remains invariant in the different models, ensuring the stability of the effects.

Table 2 Results of different spatial and spatio-temporal regression Bayesian models for child maltreatment risk. Posterior mean, standard deviation (SD) and the 95% credible interval (CI) of all parameters

For model relevance, we consider posterior probability distributions of the regression coefficients (β) of being different from zero. The probability of being positive was higher than 90% for policing activity and immigrant concentration. The probability of being negative was higher than 90% for economic status and education level. These results indicate that the risk of substantiated child maltreatment was particularly high in neighborhoods with low economic status and education level, and with high levels of policing activity and immigrant concentration. Residential instability, however, did not show a relevant association with substantiated child maltreatment.

To know the relative influence of the four neighborhood variables that were relevant for the final model, our results can be interpreted in terms of odds ratios expressed as \(\exp (\upbeta \cdot\Delta {\text{X}})\). For example a 10,000€ increase in economic status decreases the relative risk of substantiated child maltreatment by 17%; a 0.1 increase in education level decreases the relative risk by 15%; a 5 unit increase in policing activity increases the relative risk by 17%; finally, a 10% increase in immigrant concentration increases the relative risk by 9%.

Bayesian spatio-temporal modeling allows area-specific risks of substantiated child maltreatment to be mapped and differences analyzed over the years. Figure 1 shows the relative risk for each year of the study. These maps show areas with higher (> 1) or lower (< 1) than average risk. In some areas, the relative risk is more than twice the average, which reflects very high-risk levels of substantiated child maltreatment.

Fig. 1

Maps of relative risks of child maltreatment by census block group in each year of study, Valencia, Spain, 2004–2015

These maps also reveal common patterns over the years, showing areas with higher levels of relative risk at the periphery, especially in the eastern part of the city. The parameter \(\rho\) (\(\rho\) = .90) indicated a high temporal correlation between a particular year and the previous one.

Autoregressive models can be used to represent temporal paths of relative risk in different census block groups, thus identifying areas with stable risks and areas with changes in risk over time. Figure 2 shows the most stable areas. Results showed both chronic spatial patterns of high risk of substantiated child maltreatment during the years analyzed, as well as stability over time in areas of low risk. The most stable low risks are located in the city center, while peripheral areas present more stable high risks. Figure 3 shows areas with increased or decreased substantiated child maltreatment risk over the years.

Fig. 2

Temporal paths of relative risk in areas with stable high risk (above), and stable low risk (below). Relative risk values greater than 1 indicate higher risk than the city average. Relative risk values lower than 1 indicate lower risk than the city average

Fig. 3

Temporal paths of relative risk in areas with increasing and decreasing child maltreatment risk, respectively


This study showed that the neighborhoods where the families live matter in understanding spatial and temporal variations in child maltreatment risk [8,9,10,11,12]. Our results showed that neighborhoods with low levels of economic and educational status, with high levels of policing activity, and high immigrant concentration had higher risk of substantiated child maltreatment. This study illustrates how the unequal distribution of neighborhood risk factors in an urban area is linked to the unequal spatial distribution of substantiated child maltreatment risk in the city. These results support the view that ‘place’ matters in relation to child maltreatment, as these ecological variations reflect important inequalities in substantiated child maltreatment risk [10,11,12]. Among the explanatory mechanisms that have been proposed to explain this link are: social impoverishment—that is, lack of trust and support networks in the neighborhood—, and diminished social control; social isolation from mainstream values regarding what is acceptable in parent–child relationships; and high levels of parental stress [8, 10,11,12, 19, 27].

Previous studies have linked these neighborhoods characteristics—that is, neighborhood socioeconomic disadvantage, disorder and crime, and immigrant concentration—indicative of social disorganization to child maltreatment rates [10, 11, 27, 29]. However, the present study is, as far as we are aware, the first to analyze this link using a high-resolution Bayesian spatio-temporal approach to study the influence of neighborhood characteristics on small-area variations in substantiated child maltreatment risk. Using this approach, we were able to analyze small-area variations in substantiated child maltreatment risk over a 12-year period, which provided the possibility to identify and track risk trends. Our results, in addition, show that ‘time’ is also a key element to analyze variations in child maltreatment risk. Studies that do not take into account the temporal dimension may misinterpret risk estimations due to data aggregation—for example increasing or decreasing risks in certain neighborhoods would not be detected. In our study, introducing a spatio-temporal autoregressive structure clearly improved the model fit, and provided more reliable area-specific risk estimates. Thus, our results showed chronic spatial patterns of high risk of substantiated child maltreatment during the years analyzed. We found areas where risks were over five-times higher than the city average, and this high risk was stable over the years. Results suggest that in certain city areas substantiated child maltreatment risk can become ‘endemic’. On the other hand, risks in other areas were between two and four times lower than the city average, and these lower risks were also stable over the years. In this regard, low-risk city areas could be considered to be providing a more protective social environment for the well-being of children, as opposed to the very high-risk ones—for example, differences in risk between some of the low risk areas and some of the high-risk ones are nearly tenfold. Finally, spatio-temporal analyses allowed us to identify areas where substantiated child maltreatment risk increased or decreased over the years analyzed.

Detecting and tracking these inequalities in substantiated child maltreatment risk is important for better-informed prevention and intervention strategies. One of the advantages of the high-resolution approach used in the present study is that it provides small-area-specific risk estimates pointing to areas of stable excess risk, as well as to areas of increasing or decreasing substantiated child maltreatment risk. Identifying and tracking these risk trends at the small-area level can provide a more useful method to inform policy and action, as compared to other low-resolution approaches—for example, counties or large census tracks. Interventions targeting areas of excess risk of child maltreatment—for example, increasing or chronic high-risk areas—can have an important preventive potential. Neighborhood-level interventions, as opposed to a person-centered approach, can reach a larger number of families, providing a more cost-effective strategy, as not only the individuals are the subject of the intervention, but also the context where these families live [10, 55, 56]. The spatio-temporal epidemiological approach used in this study not only provides a powerful method to map and track variations in risk, but it can also contribute to a surveillance system to assess the effectiveness of intervention and prevention strategies by monitoring changes in risk over time across different city areas [6, 12].

Limitations to our study include the type of data used, namely, only officially reported and substantiated cases of child maltreatment under the supervision of Child Protection Services. Child maltreatment is, however, still largely underreported and underestimated, as many cases do not come to the attention of these services, or are unsubstantiated after reporting [4, 44]. On the other hand, it is important to note that families living in high risk neighborhoods may be more visible to authorities and therefore can be more susceptible to be reported and substantiated, as these residential areas may lead to a higher surveillance by social welfare or law enforcement agencies, as compared to other residential areas [10, 11, 42, 44]. A related issue is the potential problem of neighborhood selection bias, whereby families with higher risk of being investigated by child protection services, either choose or are forced to live in these high-risk neighborhoods [10]. The question here is to which extent the higher risk of substantiated child maltreatment is the result of the influence of neighborhood-level factors or the self-selection of families with certain characteristics in certain neighborhoods. Although this line of criticism tends to give more preeminence to individual-based explanations, than to other higher-order explanations such as neighborhood mechanisms, a substantial body of research supports the link between neighborhood concentrated disadvantage and the spatial inequality in a wide variety of outcomes, including violence, crime, or health. This body of literature led Sampson to conclude that “spatially inscribed social differences, …, constitute a family of neighborhood effects that are pervasive, strong, cross-cutting, and paradoxically stable even as they are changing in manifest form” [55, p. 6].

Regarding the covariates used in our study, the measurement of policing activity could be biased by self-report. However, this measure has adequate psychometric characteristics [49, 50] and previous studies have shown that police perceptions capture their valuable experience, are correlated with police records, and can provide important information to identify high crime areas [50, 57,58,59]. Nevertheless, clearly, future studies would benefit from using more objective crime reports. Another limitation in our study is that a number of potentially relevant neighborhood-level measures were not available for a 12-year spatio-temporal analysis. For example, socioeconomic measures such as people living below the poverty line, rates of unemployment, or income, or other covariates tapping potentially relevant neighborhood processes—such as collective efficacy, social networks, neighborhood disorder, neighborhood social norms, or alcohol outlets—, were not available [29,30,31, 60,61,62,63,64,65,66]. The modifiable areal unit problem is also a potential limitation. We are confident, however, that the high-resolution approach in this study, using the smallest geographical administrative unit available—census block groups are walkable areas with a small number of city blocks—, was particularly adequate to capture neighborhood influences on small-area variations in risk [7].

Regarding the implications for policy and action, the neighborhood conditions linked to substantiated child maltreatment in the present study are risk factors that cluster in space and, therefore, can be targeted in more focused preventive interventions. Although potentially modifiable, some of these factors, such as poor housing, or high levels of crime, or high levels of immigrant concentration are clearly difficult to change. However, other type of focused neighborhood-level interventions can address indirectly these factors. For example, urban planning and environmental approaches such as urban redevelopment and revitalization—for instance, providing new infrastructures, greening vacant lots, improving access to services, or increasing community programs—, have been shown to improve the quality of life of residents, and reduce crime, drug use, and violence in disadvantage communities [67,68,69,70].

Finally, this study was conducted in a medium-sized European city. Although similar neighborhood effects in relation to other type of offenses have been observed in European and American cities despite their differences in culture and organization [43, 55, 71, 72], future cross-national research with a similar approach would help to strengthen and generalize our results.


Our 12-year study showed that there are important spatial inequalities in substantiated child maltreatment risk across the city areas and over the years. Our study illustrates that a spatio-temporal epidemiological approach to study the geographical patterns, trends over time, and the contextual determinants of substantiated child maltreatment risk, can provide a useful method that can be of help to better inform policy and action. This methodological approach—that uses data that can be routinely collected—, can offer a more accurate description of the problem, and help to design new local prevention and intervention strategies. A spatio-temporal approach can also contribute to an improved epidemiological surveillance system to detect ecological variations in risk, and to assess the effectiveness of the initiatives to reduce this risk.


  1. 1.

    Gilbert R, Kemp A, Thoburn J, Sidebotham P, Radford L, Glaser D, Macmillan HL. Recognising and responding to child maltreatment. Lancet. 2009;373:167–80.

  2. 2.

    Reading R, Bissell S, Goldhagen J, Harwin J, Masson J, Moynihan S, Parton N, Pais MS, Thoburn J, Webb E. Promotion of children’s rights and prevention of child maltreatment. Lancet. 2009;373:332–43.

  3. 3.

    Fang X, Brown DS, Florence CS, Mercy JA. The economic burden of child maltreatment in the United States and implications for prevention. Child Abuse Negl. 2012;36:156–65.

  4. 4.

    World Health Organization. European report on preventing child maltreatment. Copenhagen: Regional Office for Europe, World Health Organization; 2013.

  5. 5.

    Finkelhor D, Turner HA, Shattuck A, Hamby SL. Violence, crime, and abuse exposure in a national sample of children and youth: an update. JAMA Pediatr. 2013;167:614–21.

  6. 6.

    Finkelhor D, Turner HA, Shattuck A, Hamby SL. Prevalence of childhood exposure to violence, crime, and abuse: results from the National Survey of Children’s Exposure to Violence. JAMA Pediatr. 2015;169:746–54.

  7. 7.

    Stoltenborgh M, Bakermans-Kranenburg MJ, Alink LR, IJzendoorn MH. The prevalence of child maltreatment across the globe: review of a series of meta-analyses. Child Abuse Rev. 2015;24:37–50.

  8. 8.

    Garbarino J, Sherman D. High-risk neighbourhoods and high-risk families: the human ecology of child maltreatment. Child Dev. 1980;51:188–98.

  9. 9.

    Zuravin SJ. The ecology of child abuse and neglect: review of the literature and presentation of data. Violence Vict. 1989;4:101–20.

  10. 10.

    Freisthler B, Merritt DH, LaScala EA. Understanding the ecology of child maltreatment: a review of the literature and directions for future research. Child Maltreat. 2006;11:263–80.

  11. 11.

    Coulton CJ, Crampton DS, Irwin M, Spilsbury JC, Korbin JE. How neighborhoods influence child maltreatment: a review of the literature and alternative pathways. Child Abuse Negl. 2007;31:1117–42.

  12. 12.

    Petersen AC, Joseph J, Feit M. New directions in child abuse and neglect research. Washington, DC: National Academics Press; 2014.

  13. 13.

    Coulton C, Korbin J, Su M, Chow J. Community level factors and child maltreatment rates. Child Dev. 1995;66:1262–76.

  14. 14.

    Korbin JE, Coulton CJ, Chard S, Platt-Houston C, Su M. Impoverishment and child maltreatment in African American and European American neighborhoods. Dev Psychopathol. 1998;10:215–33.

  15. 15.

    Drake B, Pandey S. Understanding the relationship between neighborhood poverty and specific types of child maltreatment. Child Abuse Negl. 1996;20:1003–18.

  16. 16.

    Coulton CJ, Korbin JE, Su M. Neighborhoods and child maltreatment: a multilevel study. Child Abuse Negl. 1999;23:1019–40.

  17. 17.

    Ernst JS. Community-level factors and child maltreatment in a suburban county. Soc Work Res. 2001;25:133–42.

  18. 18.

    Gracia E, Musitu G. Social isolation from communities and child maltreatment: a cross-cultural comparison. Child Abuse Negl. 2003;27:153–68.

  19. 19.

    Molnar BE, Buka SL, Brennan RT, Holton JK, Earls F. A multilevel study of neighborhoods and parent-to-child physical aggression: results from the Project on Human Development in Chicago Neighborhoods. Child Maltreat. 2003;8:84–97.

  20. 20.

    Paulsen DJ. No safe place: Assessing spatial patterns of child maltreatment victimization. J Aggress Maltreat Trauma. 2003;8:63–85.

  21. 21.

    Freisthler B. A spatial analysis of social disorganization, alcohol access, and rates of child maltreatment in neighborhoods. Child Youth Serv Rev. 2004;26:803–19.

  22. 22.

    Gracia E, Herrero J. Perceived neighborhood social disorder and residents’ attitudes toward reporting child physical abuse. Child Abuse Negl. 2006;30:357–65.

  23. 23.

    Freisthler B, Gruenewald PJ, Ring L, LaScala EA. An ecological assessment of the population and environmental correlates of childhood accident, assault and child abuse injuries. Alcohol Clin Exp Res. 2008;32:1969–75.

  24. 24.

    Freisthler B, Bruce E, Needell B. Understanding the geospatial relationship of neighborhood characteristics and rates of maltreatment for black, hispanic, and white children. Soc Work. 2007;52:7–16.

  25. 25.

    Klein S. The availability of neighborhood early care and education resources and the maltreatment of young children. Child Maltreat. 2011;16:300–11.

  26. 26.

    Maguire-Jack K. Multilevel investigation into the community context of child maltreatment. J Aggress Maltreat Trauma. 2014;23:229–48.

  27. 27.

    Freisthler B, Maguire-Jack K. Understanding the interplay between neighborhood structural factors, social processes, and alcohol outlets on child physical abuse. Child Maltreat. 2015;20:268–77.

  28. 28.

    Kim B, Maguire-Jack K. Community interaction and child maltreatment. Child Abuse Negl. 2015;41:146–57.

  29. 29.

    Molnar BE, Goerge RM, Gilsanz P, Gilsanz P, Hill A, Subramanian SV, Holton JK, Duncan DT, Beatriz ED, Beardslee WR. Neighborhood-level social processes and substantiated cases of child maltreatment. Child Abuse Negl. 2016;51:41–53.

  30. 30.

    Maguire-Jack K, Showalter K. The protective effect of neighborhood social cohesion in child abuse and neglect. Child Abuse Negl. 2016;52:29–37.

  31. 31.

    Fujiwara T, Yamaoka Y, Kawachi I. Neighborhood social capital and infant physical abuse: a population-based study in Japan. Int J Ment Health Syst. 2016;10:13.

  32. 32.

    Maguire-Jack K, Font SA. Intersections of individual and neighborhood disadvantage: implications for child. Child Youth Serv Rev. 2016;72:44–51.

  33. 33.

    Bernardinelli L, Clayton D, Pascutto C, Montomoli C, Ghislandi M, Songini M. Bayesian analysis of space-time variation in disease risk. Stat Med. 1995;14:2433–43.

  34. 34.

    Waller LA, Gotway CA. Applied spatial statistics for public health data. Hoboken: John Wiley and Sons; 2004.

  35. 35.

    Haining R, Law J, Griffith D. Modelling small area counts in the presence of overdispersion and spatial autocorrelation. Comput Stat Data Anal. 2009;53:2923–37.

  36. 36.

    Jonker MF, Congdon PD, van Lenthe FJ, Donkers B, Burdorf A, Mackenbach JP. Small-area health comparisons using health-adjusted life expectancies: a Bayesian random-effects approach. Health Place. 2013;23:70–8.

  37. 37.

    Congdon P. Applied Bayesian modelling. West Sussex: Wiley; 2014.

  38. 38.

    Law J, Quick M, Chan P. Bayesian spatio-temporal modeling for analysing local patterns of crime over time at the small-area level. J Quant Criminol. 2014;30:57–78.

  39. 39.

    Quick M, Law J, Luan H. The influence of on-premise and off-premise alcohol outlets on reported violent crime in the region of Waterloo, Ontario: applying Bayesian spatial modeling to inform land use planning and policy. Appl Spat Anal Policy. 2017;10:435–54.

  40. 40.

    Freisthler B, Weiss RE. Using Bayesian space-time models to understand the substance use environment and risk for being referred to child protective services. Subst Use Misuse. 2008;43:239–51.

  41. 41.

    Cunradi CB, Mair C, Ponicki W, Remer L. Alcohol outlets, neighborhood characteristics, and intimate partner violence: ecological analysis of a California city. J Urban Health. 2011;88:191–200.

  42. 42.

    Freisthler B, Kepple NJ, Holmes MR. The geography of drug market activities and child maltreatment. Child Maltreat. 2012;17:144–52.

  43. 43.

    Gracia E, López-Quílez A, Marco M, Lladosa S, Lila M. The spatial epidemiology of intimate partner violence: do neighborhoods matter? Am J Epidemiol. 2015;182:58–66.

  44. 44.

    Thurston H, Freisthler B, Bell J, Tancredi D, Romano PS, Miyamoto S, Joseph JG. The temporal-spatial distribution of seriously maltreated children. Spat Spatiotemporal Epidemiol. 2017;20:1–8.

  45. 45.

    Shaw CR, McKay HD. Juvenile delinquency and urban areas. Chicago: University of Chicago Press; 1942.

  46. 46.

    Bursik RJ, Grasmick HG. Neighborhood and crime. New York: Lexington Books; 1993.

  47. 47.

    Sampson RJ, Groves WB. Community structure and crime: testing social disorganization theory. AJS. 1989;94:774–802.

  48. 48.

    Sampson RJ, Morenoff JD, Gannon-Rowley T. Assessing, “neighborhood effects”: social processes and new directions in research. Annu Rev Sociol. 2002;28:443–78.

  49. 49.

    Marco M, Gracia E, Martín-Fernández M, López-Quílez A. Validation of a Google Street View-based neighborhood disorder observational scale. J Urban Health. 2017;94:190–8.

  50. 50.

    Marco M, Gracia E, López-Quílez A. Linking neighborhood characteristics and drug-related police interventions: a Bayesian spatial analysis. ISPRS Int J Geo-Inf. 2017;6:65–77.

  51. 51.

    Martínez-Beneito MA, López-Quilez A, Botella-Rocamora P. An autoregressive approach to spatio-temporal disease mapping. Stat Med. 2008;27:2874–89.

  52. 52.

    Besag J, York J, Mollié A. Bayesian image restoration, with two applications in spatial statistics. Ann Inst Stat Math. 1991;43:1–20.

  53. 53.

    Spiegelhalter DJ, Best NG, Carlin BP, Van der Linde A. Bayesian measures of model complexity and fit. J R Stat Soc Ser B Stat Methodol. 2002;64:583–639.

  54. 54.

    Gelman A, Carlin JB, Stern HS, Rubin DB. Bayesian data analysis. 2nd ed. Boca Raton: CRC Press; 1990.

  55. 55.

    Sampson RJ. Great American city: Chicago and the enduring neighborhood effect. Chicago: University of Chicago Press; 2012.

  56. 56.

    Díez-Roux AV, Mair C. Neighborhoods and health. Ann NY Acad Sci. 2010;1186:125–45.

  57. 57.

    Ratcliffe JH, McCullagh MJ. Chasing ghosts? Police perception of high crime areas. Br J Criminol. 2001;41:330–41.

  58. 58.

    Haining R, Law J. Combining police perceptions with police records of serious crime areas: a modelling approach. J R Stat Soc Ser A Stat Soc. 2007;170:1019–34.

  59. 59.

    Rengert GF, Pelfrey WV. Cognitive mapping of the city center: comparative perceptions of dangerous places. In: Weisburd D, McEwen JT, editors. Crime mapping and crime prevention. Monsey: Criminal Justice Press; 1997. p. 193–218.

  60. 60.

    Gracia E, López-Quílez A, Marco M, Lladosa S, Lila M. Exploring neighborhood influences on small-area variation in intimate partner violence risk: a Bayesian random-effect modeling approach. Int J Environ Res Public Health. 2014;11:866–82.

  61. 61.

    Sampson RJ, Raudenbush SW, Earls F. Neighborhoods and violent crime: a multilevel study of collective efficacy. Science. 1997;277:918–24.

  62. 62.

    Freisthler B, Midanik LT, Gruenewald PJ. Alcohol outlets and child physical abuse and neglect: applying routine activities theory to the study of child maltreatment. J Stud Alcohol. 2004;65:586–92.

  63. 63.

    Freisthler B, Needell B, Gruenewald PJ. Is the physical availability of alcohol and illicit drugs related to neighborhood rates of child maltreatment? Child Abuse Negl. 2005;29:1049–60.

  64. 64.

    Freisthler B, Gruenewald PJ, Remer LG, Lery B, Needell B. Exploring the spatial dynamics of alcohol outlets and child protective services referrals, substantiations, and foster care entries. Child Maltreat. 2007;12:114–24.

  65. 65.

    Morton CM, Simmel C, Peterson NA. Neighborhood alcohol outlet density and rates of child abuse and neglect: moderating effects of access to substance abuse services. Child Abuse Negl. 2014;38:952–61.

  66. 66.

    Daley D, Bachmann M, Bachmann BA, Pedigo C, Bui M, Coffman J. Risk terrain modeling predicts child maltreatment. Child Abuse Negl. 2016;62:29–38.

  67. 67.

    Kelaher M, Warr DJ, Tacticos T. Evaluating health impacts: results from the neighbourhood renewal stratgey in Victoria, Australia. Health Place. 2010;16:861–7.

  68. 68.

    Garvin EC, Cannuscio CC, Branas CC. Greening vacant lots to reduce violent crime: a randomised controlled trial. Inj Prev. 2013;19:198–203.

  69. 69.

    Mehdipanah R, Malmusi D, Muntaner C, Borrell C. An evaluation of an urban renewal program and its effects on neighborhood resident’s overall wellbeing using concept mapping. Health Place. 2013;23:9–17.

  70. 70.

    Linton SL, Jennings JM, Latkin CA, Kirk GD, Mehta SH. The association between neighborhood residential rehabilitation and injection drug use in Baltimore, Maryland, 2000–2011. Health Place. 2014;28:142–9.

  71. 71.

    Le Galès P, Zagrodzki M. Cities are back in town: The US/Europe comparison. Paris: Cahier Européen, Centre d’Etudes Européennes de Sciences Po; 2006.

  72. 72.

    Summers AA, Cheshire PC, Senn L. Urban Change in the United States and Western Europe: comparative analysis and policy. Washington, DC: The Urban Institute Press; 1999.

Download references

Authors’ contributions

EG conceptualized the study, drafted the initial manuscript, and approved the final manuscript as submitted. AL-Q designed the analytic strategy and supervised the statistical analysis, drafted the initial manuscript, and approved the final manuscript as submitted. MM conducted statistical analysis, drafted the initial manuscript, and approved the final manuscript as submitted. ML designed and coordinated the data collection, and critically reviewed the manuscript. All authors read and approved the final manuscript.


We would like to thank the Child Protective Services, the Statistics Office, and the Police Department of the city of Valencia, for their support and assistance in collecting the data for this study.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets generated and analyzed during the current study are not publicly available due to personal data protection, but are available from the corresponding author on reasonable request.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.


This study was supported by the Spanish Ministry of Economy and Competitiveness [PSI2014-54561-P, and MTM2016-77501-P] and the Spanish Ministry of Education, Culture and Sports [FPU2013/00164].

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Correspondence to Enrique Gracia.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Gracia, E., López-Quílez, A., Marco, M. et al. Mapping child maltreatment risk: a 12-year spatio-temporal analysis of neighborhood influences. Int J Health Geogr 16, 38 (2017) doi:10.1186/s12942-017-0111-y

Download citation


  • Child maltreatment
  • Neighborhood influences
  • Bayesian spatio-temporal modeling
  • Disease mapping
  • Spatial inequality
  • Small-area study
  • Area-specific risk estimation