Skip to main content

Spatiotemporal association of low birth weight with Cs-137 deposition at the prefecture level in Japan after the Fukushima nuclear power plant accidents: an analytical-ecologic epidemiological study

A Letter to the Editor to this article was published on 25 November 2020

A Letter to the Editor to this article was published on 25 November 2020

Abstract

Background

Perinatal mortality increased in contaminated prefectures after the Fukushima Daichi Nuclear Power Plant (FDNPP) accidents in Japan in 2011. Elevated counts of surgeries for cryptorchidism and congenital heart malformations were observed throughout Japan from 2012 onward. The thyroid cancer detection rate (2011 to 2016) was associated with the dose-rate at the municipality level in the Fukushima prefecture. Since the birth weight is a simple and objective indicator for gestational development and pregnancy outcome, the question arises whether the annual birth weight distribution was distorted in a dose-rate-dependent manner across Japan after Fukushima.

Methods

The Japanese Ministry of Health, Labour, and Welfare provides prefecture-specific annual counts for 26.158 million live births from 1995 to 2018, of which 2.366 million births (9.04%) with weights < 2500 g. Prefecture-specific spatiotemporal trends of the low birth weight proportions were analyzed. Logistic regression allowing for level-shifts from 2012 onward was employed to test whether those level-shifts were proportional to the prefecture-specific dose-rates derived from Cs-137 deposition in the 47 Japanese prefectures.

Results

The overall trend of the low birth weight prevalence (LBWp) in Japan discloses a jump in 2012 with a jump odds ratio (OR) 1.020, 95%-confidence interval (1.003,1.037), p-value 0.0246. A logistic regression of LBWp on the additional dose-rate after the FDNPP accidents adjusted for prefecture-specific spatiotemporal base-line trends yields an OR per μSv/h of 1.098 (1.058, 1.139), p-value < 0.0001. Further adjusting the logistic regression for the annual population size and physician density of the prefectures, as well as for the counts of the dead, the missing, and the evacuees due to earthquake and tsunami (as surrogate measures for medical infrastructure and stress) yields an OR per μSv/h of 1.109 (1.032, 1.191), p-value 0.0046.

Conclusions

This study shows increased low birth weight prevalence related to the Cs-137 deposition and the corresponding additional dose-rate in Japan from 2012 onward. Previous evidence suggesting compromised gestational development and pregnancy outcome under elevated environmental ionizing radiation exposure is corroborated.

Peer Review reports

Introduction

Low birth weight (LBW) is defined as having a birth weight of < 2500 g. LBW is an objective and reliable indicator used as a comprehensive demographic reporting measure of fetal development and pregnancy outcome [1,2,3,4]. The OECD provides an international comparison of the LBW prevalence (LBWp), which shows that Japan is among the 5 countries with the highest LBWp in the range of 9 to 10% [5]. Environmental pollutants are consistently linked to untoward pregnancy outcome and reductions in birth weight [6,7,8,9,10,11,12,13]. LBW has been suggested as an indicator of genetic detriment caused by mutation in humans exposed to ionizing radiation [14]. Analyses of birth weight and duration of pregnancies in relation to maternal age, parity, and infant survival indicated that non-survivors were significantly lighter at birth than survivors [15]. LBW is closely linked to fetal and perinatal mortality and morbidity [16]. It has been reported to be associated with disorders in perinatal periods, in childhood, and in adulthood [17, 18]. Studies in Great Britain showed that people who had low birth weight were at increased risk of coronary heart disease and the disorders related to it [19]. Animal and human studies have shown that the LBW proportions increase with toxic exposure and with radiation exposure [20,21,22]. Smoking increases the LBWp in a dose-dependent manner [23], possibly due to elevated radionuclides in tobacco [24]. Females subject to pelvic radiotherapy experience an increased risk of pre-term delivery and LBW among their offspring [25]. Treatment of female childhood cancer patients may entail restricted fetal growth and pre-term births [26]. LBW was reported after dental radiography during pregnancy [27]. A cohort study in China identified multiple risk factors of LBW including radiation exposure of fathers [28]. A natural experiment in Taiwan revealed that prenatal exposure to a continuous low-dose radiation reduced the gestational length and increased the LBW proportion [29]. In Belarus, increased LBW prevalence was reported from the highly Chernobyl-contaminated regions Gomel and Mogilev [30]. In the Ukraine, detrimental radiation-dose dependent outcomes in neonates were observed [31]. Temporarily elevated LBWp was seen in Sweden after Chernobyl [32]. Detrimental effects in humans are supported by a recent animal study: In wild Japanese monkeys (Macaca fuscata), body weight growth rate and proportional head size were significantly lower in fetuses conceived after the Fukushima disaster [33].

Since several radiation inducible genetic effects, which can be associated with general radiation exposure were observed in all of Japan after the Fukushima nuclear power plant accidents [34,35,36,37,38], an increase in the LBW proportion in whole Japan was also conceivable. This argument is supported by the specific observation of increased thyroid cancer in children and adolescents in the Fukushima prefecture, which can be related to I-131 contamination in a dose-dependent manner [39,40,41,42]. Among the investigations after the Fukushima nuclear accidents, there are reports that LBW in humans is increasing and reports that deny the increase. In the following, we shortly address two reports that are questionnaire-based surveys with a response rate in the 50% range and one survey of a small number of births in one clinic in Fukushima [43,44,45]. Questionnaire-based studies are prone to selection bias and studies with small populations (mostly in clinical settings) may likely entail type-2 errors [46]. A questionnaire-based pregnancy and birth survey was conducted by the Radiation Medical Science Center for the Fukushima Health Management Survey [43]. In this study, an increase of the LBW proportion is documented in the combined two most eastern regions Iwaki and Soso compared to the five central and western regions of the Fukushima prefecture: OR 1.163, p-value 0.0723. This observation is further supported by an increase of the stillbirth proportion in Soso and Iwaki with OR 1.923, p-value 0.1321. Since this study [43] had a participation rate of below 60%, it is likely that significant effects would be obtained with lager populations considered during longer periods. Maternal and perinatal data (2008 to 2015) were retrospectively collected for singleton live births at a hospital located 23 km from the Fukushima nuclear power plant [44]. In 1101 births, LBWp was compared pre- and post-disaster. There was no increased LBWp in any year from 2011 onward. However, with 4 years before/after the accident, i.e., 140 births per year, which means about 10 LBW-births per year, it was unlikely to receive a meaningful result, i.e., there is a large type-2 error probability in this study [44]. A more recent investigation considered 12,804 maternal outcomes during 2011–2014 in the Fukushima Prefecture [45]. However, this study neither analyzed perinatal outcomes with distance from the nuclear accident nor chronological factors. Therefore, it is unclear whether increases of LBW are due to a temporary cause of the earthquake/tsunami or due to radiation exposure. These surveys cover the Fukushima Prefecture only incompletely over short periods. In the Miyagi Prefecture, the overall rate of LBW infants was reported to be 8.7%, which tended to be lower than LBWp in 2012 of 9.3 and 9.8% in 2013 [47]. In conclusion, the trends of the LBW prevalence in most of the Japanese prefectures have not yet been scrutinized, although detailed spatiotemporal data is publicly available.

In the present study, we analyzed data of the Japanese governments ‘Demographical Survey’, which accounts for all live births and all LBW children registered in Japan excluding births to parents living abroad. Therefore, not only Fukushima Prefecture but also the whole country with differently contaminated prefectures [48] was targeted, and statistical accuracy is guaranteed by using official practically complete long-term data from 1995 to 2018, i.e., 16 years (1995 to 2010) before and 7 years (2012 to 2018) after the nuclear power plant accidents in Fukushima in March 2011. Since radioactive contamination was much more comprehensively measured and documented in the prefecture Fukushima [49] compared to the rest of Japan [48], we put both measurement regimes in perspective. We propose a rescaling of the overall Japanese contamination measurements, and we study the possible association of the LBW prevalence with radioactivity at the prefecture level.

Methods

Vital statistics and psycho-social stressors

The Japanese Statistics Bureau publishes demographical information compiled by the Ministry of Health, Labor, and Welfare. Statistics include the annual numbers of live births and the annual counts of children with a low birth weight of < 2500 g (LBW), see Table 1 or the internet platform Vital Statistics of Japan. We investigated the spatiotemporal distribution of 26.158 million live births, of which 2.366 million births (9.04%) with weights < 2500 g, across 47 Japanese prefectures from 1995 to 2018. A considerable proportion of the Japanese population was physically and psychosocially affected to a significant degree by the Great East Japan Earthquake and subsequent tsunami [51]. Therefore, the counts of earthquake related deaths, the counts of the dead and missing after earthquake and tsunami, as well as the numbers of evacuees to and within any prefecture were obtained from official sources [52] and served as additional explicit ecological confounding variables in our logistic regression models. Since medical supply including medical information may also impact the general health-behavior and thus the prevalence of LBW, annual physician density by prefecture in Japan was deployed as a surrogate confounder variable in the spatiotemporal logistic regression models, see Table 2.

Table 1 Annual live births, live births with low birth weight (LBW: birth weight < 2500 g), and LBW proportions (LBWp) in Japan stratified by exposure status of prefectures; see Table 3 and Fig. 2; Table 1 excludes 3774 live births (including 197 LBWs, 5.2%) to Japanese parents in foreign countries; see https://www.mhlw.go.jp/english/database/db-hw/vs01.html
Table 2 Mean annual population (1000) for the Japanese prefectures 1995 to 2018, dead and missing after earthquake and tsunami 2011, earth-quake related deaths, evacuated persons within or to the prefectures (sources: National Police Agency March 8, 2019; Reconstruction Agency September 30, 2018), mean annual live births, mean annual LBW and LBWp, jump OR in LBWp from 2012 onward, and 95%-CI for jump OR in LBWp from 2012 to 2018; see https://stats-japan.com/t/kiji/10343

Cs-137 deposition

Yasunari et al. published average prefecture-specific Cs-137 deposition after the Fukushima nuclear power plant accidents for the 47 prefectures of Japan [48]; see Table 3. Yasunari et al. assumed that half of the total cesium deposited was Cs-134. Therefore, it is easily possible to re-scale all calculations in this paper to total cesium in place of Cs-137. The Yasunari et al. data understate the true Cs-137 deposition, which underestimation may be supported by the following aspects:

  • Yasunari et al.’s data based on measurements restricted to March 20th to April 19th, 2011.

  • Yasunari et al. report a value of 24.7 kBq/m2 Cs-137 for Fukushima prefecture (see Table 3), whereas the UNSCEAR data set 2013/2014 documents a mean value of 153.957 kBq/m2 Cs-137, which amounts to a factor 6 underestimation of the deposition in Fukushima by Yasunari et al. [48], see the Excel-file provided by UNSCEAR as referenced in Table 3 [49].

  • It is implausible that Fukushima prefecture would be less contaminated than Miyagi and Ibaraki prefectures, see Table 3.

  • The Yasunari et al. data decay with r-3.27 at distance r form the FDNPP, see Fig. 1, whereas a theoretical decay law of r-1.42 is expected according to UNSCEAR [50], and as empirically confirmed for the Fukushima prefecture [40].

Table 3 Distances of the centers of the Japanese prefectures’ area polygons from the FDNPP, Cs-137 deposition in the Japanese prefectures after the Fukushima nuclear power plant accidents as of March 2011 according to [48], rescaled Cs-137 deposition according to [48,49,50], and dose-rate [μSv/h] derived from the rescaled deposition; see https://www.unscear.org/docs/publications/2013/UNSCEAR_2013_Annex-A_Attach_C-2.xls and http://www.pnas.org/content/108/49/19530.full
Fig. 1
figure 1

Deposition of Cs-137 in the 47 Japanese prefectures according to Yasunari et al. [48] by the prefectures’ distances from the Fukushima Daichi Nuclear Power Plant (FDNPP); gray circles: original Yasunari et al. data; black dots: deposition for Fukushima corrected and remaining depositions rescaled to a decay of r−2 with distance r, see Table 3

Since strong underestimation of radiation exposure would exaggerate any dose-specific radiation risk estimates, we suggest and propagate a correction and a rescaling of the Yasunari et al. deposition data for all of Japan despite the disadvantages of these data listed above. The rationale behind this is that the Yasunari et al. data, while restricted to a narrow time frame nevertheless reflect a valid mutual relative exposure status amongst the prefectures. To this end, we firstly increased the original deposition value 24,718.4 Bq/m2 of the Fukushima prefecture by a factor of 4.3 to 106,867.0 Bq/m2 based on the MEXT/UNSCEAR data [49] excluding 20 locations in the immediate vicinity of FDNPP with more than 2.0E+ 6 Bq/m2 with less likely importance for public exposure. Secondly, we rescaled the deposition data with the original decay-rate r-3.27 to a decay of r-2.00, which is a compromise between the theoretical decay r-1.42 by UNSCEAR [50] and the Yasunari decay r-3.27. The rescaling details and results are depicted in Fig. 1 and listed in Table 3. Figure 2 shows a geographic region value plot for the rescaled Cs-137 deposition in the Japanese prefectures.

Fig. 2
figure 2

Geographic region value plot of the decadic logarithm for the rescaled Cs-137 deposition in 47 Japanese prefectures after the Fukushima nuclear power plant accidents as of March 2011 [48], see Table 3; indication of the positions of the earthquake epi-center, the FDNPP, and a 300 km geo-circle around FDNPP; for the prefecture codes see Table 2 or Table 3

Dose-rate (μSv/h) derived from Cs-137 deposition

The United Nations Scientific Committee on the Effects of Atomic Radiation (UNSCEAR) published Cs-137 deposition and corresponding dose-rate readings at 2180 locations in the Fukushima prefecture or close to the borders of the Fukushima prefecture [49]. These data were provided by the Government of Japan as described in the report titled ‘Summarized version of the results of the research on distribution of radioactive substances discharged by the accident at TEPCO’s Fukushima Daiichi NPP’. The Japan Atomic Energy Authority (JAEA) conducted the survey with cooperation of universities and research institutes. The Ministry of Education, Culture, Sports, Science, and Technology in Japan (MEXT) was responsible for the measurements and their validity. UNSCEAR reviewed and published the dataset [49]. The single dose-rate readings resulting from all relevant deposited radionuclides including Cs-134 range from 0.040 μSv/h to 54.800 μSv/h, with mean 1.259 μSv/h and median 0.40 μSv/h. The single C-137 measurements range from 590 Bq/m2 to 15,450,928 Bq/m2, with mean 153,957 Bq/m2 and median 39,714 Bq/m2. A 2nd degree regression of the dose-rate on the Cs-137 deposition allows the translation of fallout to dose-rate. The functional details of this association are presented in Fig. 3 and the resulting dose-rates are listed in Table 3. Unfortunately, fine-resolution contamination data as available for the Fukushima prefecture, does not exist for all of Japan.

Fig. 3
figure 3

Association of the dose-rate [μSv/h] at 1 m height with the Cs-137 deposition [kBq/m2] in and near the Fukushima prefecture for 2180 positive deposition measurements and 2175 positive dose-rate readings; see [49] and Fig. 5A in [40] https://www.unscear.org/docs/publications/2013/UNSCEAR_2013_Annex-A_Attach_C-2.xls

Statistical methods

A powerful method to assess data that contains spatial as well as temporal information is spatiotemporal (logistic) regression [38, 53,54,55,56,57]. The basic idea is to adjust a regression model for region-specific trend functions and to allow for local or global level-shifts at certain points in time, or, preferably, to allow for the heights of local drops or jumps to be proportional to the contamination level of the regional strata. The advantage of this spatiotemporal method is optimization of adjustment and minimization of confounding by considering partial trends of regional units. Values of the outcome variable (here LBW) in those partial trends are modeled and compared within the same regional stratum, as the target variable describing the characteristic of interest varies from year to year. Information on several regional units is combined in a global spatiotemporal model, giving rise to tests of local or global change-points in time as well as spatial trends in the outcome variable with regionally determined contamination or exposure. In the present study, we considered the 47 prefectures of Japan with their documented individual annual LBW proportions (1995 to 2018), their original [48] and rescaled Cs-137 depositions, and the associated dose-rates compiled in Table 3. We tested whether possible changes of LBWp from 2012 onward were dependent on Cs-137 fallout at the prefecture level. We further adjusted the spatiotemporal logistic regression model for the annual population size of the prefectures (1995 to 2015), for the number of physicians per 1000 population, as well as for the counts of the dead, the missing, and the evacuees due to the earthquake and the tsunami as surrogate measures for the available prefecture-specific medical infrastructure and the stress associated with this triple catastrophe and its aftermath, see Table 2. For data processing, statistical analyses, and results display, we used Microsoft Excel 2016 (Office 365), R 3.5.1 (Version 2017-10-04), Wolfram MATHEMATICA 11.3, and mostly SAS/STAT software 9.4 (SAS Institute Inc.: SAS/STAT User’s Guide, Version 9.4, Cary NC: SAS Institute Inc.,© 2002–2012).

Results

Figure 4A shows the annual marginal LBW distribution for all of Japan 1995 to 2018, excluding the births to Japanese parents in foreign countries (number of births abroad: 3774, LBWp = 0.052). See the columns ‘Total’ in Table 1 for the corresponding absolute counts and the LBW proportions (LBWp). As a first step, we fit to this overall LBWp a smooth 4th degree polynomial allowing for a change-point in 2012 after the Fukushima nuclear power plant accidents. This approach discloses a significant jump in 2012 with a jump odds ratio (OR) 1.020, 95%-CI (1.003, 1.037), p-value 0.0246, see Fig. 4A. In Japan, people consider that the reportedly healthy lean structure of women led to reduced weight gain during pregnancy and contributed to the rise in LBW prior to 2007, because low BMI and poor weight gain are risk factors for LBW. However, the increase in LBWp reached its overall maximum in 2007 (see Fig. 4A), and subsequently LBWp slightly fell or remained nearly constant due to a change in health-awareness and behavior [58].

Fig. 4
figure 4

Low birth weight (LBW) proportion in Japan 1995 to 2018; 4th degree polynomial logistic regression trends allowing for jumps from 2012 onward; A: Japan; B: Japan excluding 10 exposed prefectures; C: 5 moderately exposed prefectures; D: 5 highly exposed prefectures, see Table 1 for the absolute counts, the relative frequencies, and the ISO codes of the prefectures

Considering the possibility that the jump height in 2012 may be associated with the Cs-137 fallout in the prefectures, we analyze and depict in a second step the behavior of LBWp in the 3 strata of prefectures according to Table 1: in the 37 least contaminated prefectures, in the 5 moderately contaminated prefectures (Yamagata, Saitama, Tokyo, Kanagawa, Chiba), and in the 5 most contaminated prefectures (Fukushima, Miyagi, Ibaraki, Tochigi, Iwate). Figures 4B-D display the result: the higher the fallout in the prefectures, the higher the jumps in LBWp from 2012 to 2018. The excess LBW counts and their respective 95%-confidence intervals corresponding to the jump odds ratios (ORs) in the LBWp trends in Fig. 4A-D are (A) 11,561 (1470, 21,793), (B) 5659 (− 1165, 12,585), (C) 3458 (778, 6172), and (D) 2484 (584, 4448), respectively.

Table 2 and Fig. 5 generalize and visualize the effects seen in Fig. 4 by listing and plotting the prefecture-specific level-shifts in 2012 in the LBWp trends against the average dose-rate in the prefectures. The combination of the 37 least contaminated prefectures in one group avoids an overly scattered picture for these regions in the left part of Fig. 5. The leftmost data point in Fig. 5 represents this group of the 37 low or only slightly contaminated prefectures. In Fig. 5, a variance weighted straight line regression of the individual jump odds ratios against the dose-rates discloses a significant linear relationship (R2 = 0.82) with slope 0.11 per μSv/h and p-value < 0.0001.

Fig. 5
figure 5

Odds ratios for the jumps in the low birth weight proportion trends (LBWp) from 2012 onward by prefecture-specific dose-rates derived from the rescaled Cs-137 deposition in the Japanese prefectures from March 20th to April 19th 2011; restricted linear regression yields trend p-value < 0.0001; the left data point summarizes and represents 37 slightly radiologically impacted prefectures, the 10 data points from the right represent the 10 prefectures with high to moderate pollution, see Table 1

A more direct approach is logistic regression of LBWp on the additional dose-rate after Fukushima adjusted for prefecture-specific spatiotemporal base-line trends [57]. This yields an OR per μSv/h of 1.098 (1.058, 1.139), p-value < 0.0001. By additionally adjusting this spatiotemporal logistic regression for the counts of the earthquake-related deaths, the dead and the missing after the tsunami, and the counts of evacuees within and to the prefectures as surrogate measures for the disaster-related stress from 2011 onward (see Table 2), as well as additionally taking into account the prefecture-specific population size and physician density per 1000 population as surrogate measures of general stress and available medical infrastructure, we obtain a somewhat larger, however less precise adjusted OR per μSv/h of 1.109 (1.032, 1.191), p-value 0.0046. The decreased precision resulting from this additional adjustment may be explained by variance inflation in the spatiotemporal logistic regression model due to partly correlated surrogate confounder measures. In summary, the increase of the background dose-rate by 1 μSv/h elevates the prevalence odds of low birth weight babies by approximately 10%. Note, 1 μSv/h translates to a dose of 8.8 mSv/year. Importantly, without the suggested rescaling of the Yasunari et al. exposure data, the dose-rate specific effect would be 50% in place of 10%, and this would likely be an overestimation of the radiation effect due to an obvious underestimation of the overall Cs-137 deposition and the associated dose-rate across Japan by Yasunari et al. [48].

To more directly assess and display the relative impacts of the earthquake and the tsunami versus the effects of the Cs-137 deposition and the associated dose-rate on LBW, we compared the three contaminated prefectures Fukushima, Iwate, and Miyagi, where the dead and missing persons due to earthquake and tsunami were numerous (n = 18,359), to the somewhat weaker contaminated Ibaraki, Tochigi, and Yamagata where only relatively few immediate deaths occurred and few persons were missing (n = 31), see Table 2. Figure 6A and B show that for the less versus strongly earthquake and tsunami impacted groups of prefectures, the LBWp jump heights are similar with largely overlapping 95%-CIs and with similar p-values. Therefore, the long-term increasing LBWp is essentially independent of the direct or protracted impact of earthquake and tsunami.

Fig. 6
figure 6

Low birth weight (LBW) proportion and parsimonious constant trends allowing for jumps in 2012 in 6 Japanese Fukushima exposed prefectures (2005 to 2018) stratified by tsunami impact; A: low tsunami impact in Ibaraki, Tochigi, and Yamagata; B: high tsunami impact in Fukushima, Iwate, and Miyagi

Discussion

This study strengthens the evidence provided by previous investigations [20, 26,27,28,29,30, 32] that elevated exposure to ionizing radiation increases the prevalence of low birth weight children. The proportion of low birth weight babies in Japan (LBW < 2500 g) was increasing continuously from 1995 to a peak value in 2007, see Fig. 4A. Control of maternal weight gain in reproductive-age women in the general Japanese population appeared to be a major factor involved in the increase in LBW babies before 2007 [4]. Since 2007, the conventional practice of ‘suppressing weight gain during pregnancy to within 10 kg’ has been changed to not reducing the weight gain too much. This may be the reason why most of the prefectures have stopped increasing LBW since peaking around 2007. The revision may spread rapidly in prefectures with large cities and high physician density, and, therefore, the degree of spread may vary between the prefectures. This situation suggests the adjustment of the spatiotemporal LBWp logistic regression models not only for its estimable and known spatiotemporal base-line trend parameters and possible determinants, such as Cs-137 deposition and earthquake-related stress measures, but also for the annual population size and the physician density of the prefectures. In Table 2 we compiled these potential ecological confounders. However, the adjustment for the additional confounders (population counts, physician density, and the triple-disaster-related stress indicators from 2012 onward) does not change our effect estimate of approximately 10% per 1 μSv/h. The reason for this may be that all major prefecture-specific information is already captured by the prefecture-specific spatiotemporal base-line trends accounted for in the adopted spatiotemporal regression approach [57].

Our investigation disclosed a positive association between the Cs-137 deposition across Japan after Fukushima with the prevalence of low birth weight (< 2500 g). Therefore, previously reported epidemiological health detriment after Fukushima [34,35,36,37, 39,40,41, 59] can be generalized and corroborated. Nevertheless, the question whether ionizing radiation exposure of young people in reproductive age or peri-conceptional and embryonic radiation exposure impair fetal and post-natal development of offspring remains a controversial issue. There are articles in favor of and against this hypothesis, e.g., [60, 61]. A problem with statistically negative studies in the clinical setting is sample size - typically in the range of a few thousand [61] or less. Small sample sizes generally entail low statistical power implying large type-2 error probabilities. It may be rather improbable to detect relevant changes in low birth weight proportions, say in the order of 10%, with population sizes ranging in the thousands only. For example, a two-sided one-sample binomial test for testing a hypothetically increased LBWp of 0.11 against a typical null-LBWp of 0.10 (i.e. 10% increase) requires a sample size of 7248 to achieve a statistical power of at least 80%. For more realistic two-sample scenarios involving additional independent LBW-determinants entailing enhanced biological variability, the required sample sizes for obtaining meaningful results would be even larger. Therefore, it is of no surprise that no unequivocal evidence has been obtained yet. For example, in a study mentioned in the introduction [31], there was no statistically significant effect (p-value > 0.1) of fetal dose on birth weight in 2582 in-utero-exposed individuals from northern Ukraine for whom estimates of fetal thyroid I-131 dose were available. Because of this relatively small population size (n = 2582), the statistical power for detecting a relevant 10% increase in LBW prevalence, which prevalence is itself in the range of 10%, achieves only 40%. In contrast, our study with effective sample sizes in the order of 2,000,000 live births in the 5 moderately contaminated prefectures and 500,000 live births in the 5 highly contaminated prefectures from 2012 to 2018 (see Table 1) yields statistical powers of over 80% for detecting 2 and 4% increases in the LBW prevalence, after Fukushima, respectively.

Under the headline “Radiation-induced mutation rates in man”, UNSCEAR [14] emphasized already in the year 1958 “All the results obtained are subject to an inevitable sampling error which necessitates the collection of a very large amount of data. A number of quantitative characters, such as birth weight, size and various anthropometric measurements, as well as statistical data, such as neo-natal mortality, have been suggested and examined. Unfortunately, the precise genetic component in these variables is not known; on the contrary, they are known to be dependent upon factors which are economic (standard of living), demographic (age of parents, order of birth, etc.) and sociological (medical care).” The sample size issue addressed in this statement may be resolved when instead of at most thousands of births in clinical settings many millions of births in ecological epidemiological studies can be considered: After the nuclear accidents of Chernobyl and Fukushima, the populations of large regions or even whole countries have been exposed to additional ionizing radiation significantly elevating the existing background radiation by, e.g., 10% or above [54, 56, 62, 63]. Moreover, in a large-scale ecological design, as the one presented here, the socio-demographic and environmental determinants of the low birth weight prevalence can be considered similar in and comparable between the regional units (prefectures). The differences within and between the regional trends from 2012 onward can be assessed by spatiotemporal logistic regression adjusted for appropriately chosen base-line trend parameters and further LBW determinants and ecological confounders [57].

Conclusions

This study shows increased low birth weight prevalence across Japan related to the prefecture-specific dose-rate derived from Cs-137 deposition after Fukushima. One (1.0) μSv/h (equivalent to 8.8 mSv/year) increases the odds of observing low birth weight events by approximately 10%. Therefore, previous investigations suggesting compromised gestational development and impaired pregnancy outcome under elevated ionizing radiation levels have been corroborated by the present study. These findings, in the overall view, call for intensifying bio-physical research in exposure mechanisms and exposure pathways of natural or artificial ionizing radiation. Biological, epidemiological, and medical research should aim at clarifying the genetic and the carcinogenic consequences of enhanced radiation in the environment or in the workplace. Radiation-induced genetic effects may occur without immediately obvious link to spectacular incidents or accidents [63, 64]. Therefore, the legislator, the nuclear industry, and the nuclear and radio-pharmaceutical medicine must impose and exert even greater care when processing, employing, and disposing radioactive materials.

Availability of data and materials

The employed data has exclusively been published previously and/or it is contained in the Tables and in the Figures included in this paper.

Abbreviations

95%-CI or (. , .):

95%-confidence interval

μSv/h:

Micro-Sievert per hour

Cs-134:

Cesium-134

Cs-137:

Cesium-137

Bq:

Becquerel

FDNPP:

Fukushima Daiichi Nuclear Power Plant

FHMS:

Fukushima Health Management Survey

kBq/m2 :

Kilo-Becquerel per square meter

LBW:

Low birth weight

LBWp :

Low birth weight prevalence (or proportion)

MEXT:

Ministry of Education, Culture, Sports, Science, and Technology in Japan

mSv:

Milli-Sievert

OECD:

Organization for Economic Co-operation and Development

OR:

Odds Ratio

SAS:

Statistical Analysis System, software produced by SAS Institute Inc.

TEPCO:

Tokyo Electric Power Company

UNSCEAR:

United Nations Scientific Committee on the Effects of Atomic Radiation

WHO:

World Health Organization

References

  1. WHO: World Health Organization (WHO). Physical status: the use and interpretation of anthropometry. WHO technical report series no. 854. Geneva: WHO; 1995.

    Google Scholar 

  2. WHO. South East Asia regional NEONATAL-PERINATAL DATABASE, World Health Organization (South-East Asia region), working definitions; 2009. https://www.newbornwhocc.org/pdf/database.pdf, Accessed 15 November, 2019.

    Google Scholar 

  3. Hughes MM, Black RE, Katz J. 2500-g low birth weight cutoff: history and implications for future research and policy. Matern Child Health J. 2017;21(2):283–9.

    Google Scholar 

  4. Honami Y, Noriko K, Tetsuji Y. Current trends in low birth weight infants in Japan. J Natl Inst Public Health. 2014;63(1):2–16.

    Google Scholar 

  5. OECD. Health statistics; European Community health indicators (ECHI), OECD health at a glance; 2019. accessed through the HEIDI data tool; https://www.oecd.org/els/family/CO_1_3_Low_birth_weight.pdf. Accessed May 20, 2020.

    Google Scholar 

  6. Bobak M. Outdoor air pollution, low birth weight, and prematurity. Environ Health Perspect. 2000;108(2):173–6.

    CAS  Google Scholar 

  7. Dugandzic R, Dodds L, Stieb D, Smith-Doiron M. The association between low level exposures to ambient air pollution and term low birth weight: a retrospective cohort study. Environ Health. 2006;5(1):3.

    Google Scholar 

  8. Grazuleviciene R, Nieuwenhuijsen MJ, Vencloviene J, Kostopoulou-Karadanelli M, Krasner SW, Danileviciute A, et al. Individual exposures to drinking water trihalomethanes, low birth weight and small for gestational age risk: a prospective Kaunas cohort study. Environ Health. 2011;10(1):32.

    CAS  Google Scholar 

  9. Kloog I, Melly SJ, Ridgway WL, Coull BA, Schwartz J. Using new satellite based exposure methods to study the association between pregnancy pm2.5 exposure, premature birth and birth weight in Massachusetts. Environ Health. 2012;11(1):40.

    Google Scholar 

  10. Ruckart PZ, Bove FJ, Maslia M. Evaluation of contaminated drinking water and preterm birth, small for gestational age, and birth weight at marine Corps Base camp Lejeune, North Carolina: a cross-sectional study. Environ Health. 2014;13(1):99.

    Google Scholar 

  11. Seo JH, Ha EH, Kim OJ, Kim BM, Park HS, Leem JH, et al. Kim YJ: [environmental health surveillance of low birth weight in Seoul using air monitoring and birth data]. J Prev Med Public Health. 2007;40(5):363–70.

    Google Scholar 

  12. Morello-Frosch R, Jesdale BM, Sadd JL, Pastor M. Ambient air pollution exposure and full-term birth weight in California. Environ Health. 2010;9(1):44.

    Google Scholar 

  13. Genowska A, Jamiołkowski J, Szafraniec K, Stepaniak U, Szpak A, Pająk A. Environmental and socio-economic determinants of infant mortality in Poland: an ecological study. Environ Health. 2015;14(1):61.

    Google Scholar 

  14. UNSCEAR. Report of the United Nations Scientific committee on the effects of atomic radiation (UNSCEAR), general Assembly official records: thirteenth session supplement no. 17 (a/3838); ANEX H:the genetic effects of radiation; 1958. p. 180.

    Google Scholar 

  15. Karn MN, Penrose LS. Birth weight and gestation time in relation to maternal age, parity and infant survival. Ann Eugenics. 1952;16:147–64.

    CAS  Google Scholar 

  16. Coutinho PR, Cecatti JG, Surita FG, Costa ML, Morais SS. Perinatal outcomes associated with low birth weight in a historical cohort. Reprod Health. 2011;8:18.

    Google Scholar 

  17. Kaijser M, Bonamy AK, Akre O, Cnattingius S, Granath F, Norman M, et al. Perinatal risk factors for ischemic heart disease: disentangling the roles of birth weight and preterm birth. Circulation. 2008;117(3):405–10.

    Google Scholar 

  18. Kaijser M, Bonamy AK, Akre O, Cnattingius S, Granath F, Norman M, et al. Perinatal risk factors for diabetes in later life. Diabetes. 2009;58(3):523–6.

    CAS  Google Scholar 

  19. Barker DJ. Early growth and cardiovascular disease. Arch Dis Child. 1999;80(4):305–7.

    CAS  Google Scholar 

  20. Yoshimoto Y, Schull WJ, Kato H, Neel JV. Mortality among the offspring (F1) of atomic bomb survivors, 1946-85. J Radiat Res (Tokyo). 1991;32(4):327–51.

    CAS  Google Scholar 

  21. Tang FR, Loke WK, Khoo BC. Low-dose or low-dose-rate ionizing radiation-induced bioeffects in animal models. J Radiat Res. 2017;58(2):165–82.

    CAS  Google Scholar 

  22. Dehmel S, Nathan P, Bartel S, El-Merhie N, Scherb H, Milger K, et al. Intrauterine smoke exposure deregulates lung function, pulmonary transcriptomes, and in particular insulin-like growth factor (IGF)-1 in a sex-specific manner. Sci Rep. 2018;8(1):7547.

    Google Scholar 

  23. Voigt M, Hermanussen M, Wittwer-Backofen U, Fusch C, Hesse V. Sex-specific differences in birth weight due to maternal smoking during pregnancy. Eur J Pediatr. 2006;165(11):757–61.

    CAS  Google Scholar 

  24. Papastefanou C. Radioactivity of tobacco leaves and radiation dose induced from smoking. Int J Environ Res Public Health. 2009;6(2):558–67.

    CAS  Google Scholar 

  25. Wo JY, Viswanathan AN. Impact of radiotherapy on fertility, pregnancy, and neonatal outcomes in female cancer patients. Int J Radiat Oncol Biol Phys. 2009;73(5):1304–12.

    Google Scholar 

  26. Signorello LB, Cohen SS, Bosetti C, Stovall M, Kasper CE, Weathers RE, et al. Female survivors of childhood cancer: preterm birth and low birth weight among their children. J Natl Cancer Inst. 2006;98(20):1453–61.

    Google Scholar 

  27. Hujoel PP, Bollen AM, Noonan CJ, del Aguila MA. Antepartum dental radiography and infant low birth weight. JAMA. 2004;291(16):1987–93.

    CAS  Google Scholar 

  28. Chen S, Yang Y, Qv Y, Zou Y, Zhu H, Gong F, et al. Paternal exposure to medical-related radiation associated with low birthweight infants: a large population-based, retrospective cohort study in rural China. Medicine. 2018;97(2):e9565.

    Google Scholar 

  29. Tsou M, Liu J, Hammitt JK, et al. The effect of prenatal exposure to radiation on birth outcomes: exploiting a natural experiment in Taiwan. JER. 2019. https://0-doi-org.brum.beds.ac.uk/10.1007/s42973-019-00016-9.

  30. Petrova A, Gnedko T, Maistrova I, Zafranskaya M, Dainiak N. Morbidity in a large cohort study of children born to mothers exposed to radiation from Chernobyl. Stem Cells. 1997;15(Suppl 2):141–50.

    Google Scholar 

  31. Hatch M, Little MP, Brenner AV, Cahoon EK, Tereshchenko V, Chaikovska L, et al. Neonatal outcomes following exposure in utero to fallout from Chernobyl. Eur J Epidemiol. 2017;32(12):1075–88.

    CAS  Google Scholar 

  32. Ericson A, Kallen B. Pregnancy outcome in Sweden after the Chernobyl accident. Environ Res. 1994;67(2):149–59.

    CAS  Google Scholar 

  33. S-i H, Tsuchiya M, Ochiai K, Nakiri S, Nakanishi S, Ishii N, et al. Small head size and delayed body weight growth in wild Japanese monkey fetuses after the Fukushima Daiichi nuclear disaster. Sci Rep. 2017;7(1):3528.

    Google Scholar 

  34. Scherb H, Mori K, Hayashi K. Increases in perinatal mortality in prefectures contaminated by the Fukushima nuclear power plant accident in Japan: a spatially stratified longitudinal study. Medicine (Baltimore). 2016;95(38):e4958.

    Google Scholar 

  35. Scherb H, Mori K, Hayashi K. Comment on ‘perinatal mortality after the Fukushima accident’. J Radiol Prot. 2019;39(2):647–9.

    Google Scholar 

  36. Murase K, Murase J, Mishima A. Nationwide increase in complex congenital heart diseases after the Fukushima nuclear accident. J Am Heart Assoc. 2019;8(6):e009486.

    Google Scholar 

  37. Murase K, Murase J, Machidori K, Mizuno K, Hayashi Y, Kohri K. Nationwide increase in cryptorchidism after the Fukushima nuclear accident. Urology. 2018;118:65–70.

    Google Scholar 

  38. Korblein A, Kuchenhoff H. Perinatal mortality after the Fukushima accident: a spatiotemporal analysis. J Radiol Prot. 2019;39(4):1021–30.

    Google Scholar 

  39. Tsuda T, Tokinobu A, Yamamoto E, Suzuki E. Thyroid cancer detection by ultrasound among residents ages 18 years and younger in Fukushima, Japan: 2011 to 2014. Epidemiology. 2016;27(3):316–22.

    Google Scholar 

  40. Yamamoto H, Hayashi K, Scherb H. Association between the detection rate of thyroid cancer and the external radiation dose-rate after the nuclear power plant accidents in Fukushima, Japan. Medicine (Baltimore). 2019;98(37):e17165.

    CAS  Google Scholar 

  41. Kato T. Re: associations between childhood thyroid cancer and external radiation dose after the Fukushima Daiichi nuclear power plant accident. Epidemiology. 2019;30(2):e9–e11.

    Google Scholar 

  42. Toki H, Wada T, Manabe Y, Hirota S, Higuchi T, Tanihata I, et al. Relationship between environmental radiation and radioactivity and childhood thyroid cancer found in Fukushima health management survey. Sci Rep. 2020;10(1):4074.

    CAS  Google Scholar 

  43. Fujimori K, Kyozuka H, Yasuda S, Goto A, Yasumura S, Ota M, et al. Pregnancy and birth survey after the great East Japan earthquake and Fukushima Daiichi nuclear power plant accident in Fukushima prefecture. Fukushima J Med Sci. 2014;60(1):75–81.

    Google Scholar 

  44. Leppold C, Nomura S, Sawano T, Ozaki A, Tsubokura M, Hill S, et al. Birth outcomes after the Fukushima Daiichi nuclear power plant disaster: a long-term retrospective study. Int J Environ Res Public Health. 2017;14(5):542.

    Google Scholar 

  45. Kyozuka H, Fujimori K, Hosoya M, Yasumura S, Yokoyama T, Sato A, et al. The Japan environment and Children's study (JECS) in Fukushima prefecture: pregnancy outcome after the great East Japan earthquake. Tohoku J Exp Med. 2018;246(1):27–33.

    Google Scholar 

  46. Choi BC, Pak AW. A catalog of biases in questionnaires. Prev Chronic Dis. 2005;2(1):A13.

    Google Scholar 

  47. Sugawara J, Iwama N, Hoshiai T, Tokunaga H, Nishigori H, Metoki H, et al. Regional birth outcomes after the 2011 great East Japan earthquake and tsunami in Miyagi prefecture. Prehosp Disaster Med. 2018;33(2):215–9.

    Google Scholar 

  48. Yasunari TJ, Stohl A, Hayano RS, Burkhart JF, Eckhardt S, Yasunari T. Cesium-137 deposition and contamination of Japanese soils due to the Fukushima nuclear accident. Proc Natl Acad Sci. 2011;108(49):19530.

    CAS  Google Scholar 

  49. UNSCEAR. Report 2013, Volume I, United Nations Scientific Committee on the Effects of Atomic Radiation, REPORT TO THE GENERAL ASSEMBLY, SCIENTIFIC ANNEX A: Levels and effects of radiation exposure due to the nuclear accident after the 2011 great east-Japan earthquake and tsunami, http://www.unscear.org/docs/reports/2013/13-85418_Report_2013_Annex_A.pdf and https://www.unscear.org/docs/publications/2013/UNSCEAR_2013_Annex-A_Attach_C-2.xls. Accessed May 22, 2019.

  50. UNSCEAR. Report 2017. United Nations Scientific Committee on the Effects of Atomic Radiation, Report to the General Assembly, SCIENTIFIC ANNEXES A and B, https://www.unscear.org/unscear/en/publications/2017.html. Accessed November 5, 2019.

  51. Harada N, Shigemura J, Tanichi M, Kawaida K, Takahashi S, Yasukata F. Mental health and psychological impacts from the 2011 great East Japan earthquake disaster: a systematic literature review. Disaster Mil Med. 2015;1:17.

    Google Scholar 

  52. Hasegawa A, Ohira T, Maeda M, Yasumura S, Tanigawa K. Emergency responses and health consequences after the Fukushima accident; evacuation and relocation. Clin Oncol (R Coll Radiol). 2016;28(4):237–44.

    CAS  Google Scholar 

  53. Turner SL, Karahalios A, Forbes AB, Taljaard M, Grimshaw JM, Cheng AC, et al. Design characteristics and statistical methods used in interrupted time series studies evaluating public health interventions: protocol for a review. BMJ Open. 2019;9(1):e024096.

    Google Scholar 

  54. Scherb H, Weigelt E, Brüske-Hohlfeld I. European stillbirth proportions before and after the Chernobyl accident. Int J Epidemiol. 1999;28(5):932–40.

    CAS  Google Scholar 

  55. Scherb H, Weigelt E, Brüske-Hohlfeld I. Regression analysis of time trends in perinatal mortality in Germany, 1980-1993. Environ Health Perspect. 2000;108(2):159–65.

    CAS  Google Scholar 

  56. Scherb H, Weigelt E. Congenital malformation and stillbirth in Germany and Europe before and after the Chernobyl nuclear power plant accident. Environ Sci Pollut Res Special Issue. 2003;1:117–25.

    Google Scholar 

  57. Scherb H, Voigt K. Analytical ecological epidemiology: exposure-response relations in spatially stratified time series. Environmetrics. 2009;20(6):596–606.

    Google Scholar 

  58. Yorifuji T, Naruse H, Kashima S, Murakoshi T, Kato T, Inoue S, et al. Trends of preterm birth and low birth weight in Japan: a one hospital-based study. BMC Pregnancy Childbirth. 2012;12(1):162.

    Google Scholar 

  59. Kato T: Area Dose Response of Prevalent Childhood Thyroid Cancers after the Fukushima Nuclear Power Plant Accident. In., vol. Clinical Oncology and Research: Area Dose Response of Prevalent Childhood Thyroid Cancers: Science Repository, DOI: https://0-doi-org.brum.beds.ac.uk/10.31487/j.COR.2019.06.16, Accessed 6 January, 2020.

  60. Goldberg MS, Mayo NE, Levy AR, Scott SC, Poitras B. Adverse reproductive outcomes among women exposed to low levels of ionizing radiation from diagnostic radiography for adolescent idiopathic scoliosis. Epidemiology. 1998;9(3):271–8.

    CAS  Google Scholar 

  61. Mortazavi SM, Shirazi KR, Mortazavi G. The study of the effects of ionizing and non-ionizing radiations on birth weight of newborns to exposed mothers. J Nat Sci Biol Med. 2013;4(1):213–7.

    CAS  Google Scholar 

  62. Cardis E, Krewski D, Boniol M, Drozdovitch V, Darby SC, Gilbert ES, et al. Estimates of the cancer burden in Europe from radioactive fallout from the Chernobyl accident. Int J Cancer. 2006;119(6):1224–35.

    CAS  Google Scholar 

  63. Scherb H, Kusmierz R, Voigt K. Increased sex ratio in Russia and Cuba after Chernobyl: a radiological hypothesis. Environ Health. 2013;12:63.

    Google Scholar 

  64. Schmitz-Feuerhake I, Busby C, Pflugbeil S. Genetic radiation risks: a neglected topic in the low dose debate. Environ Health Toxicol. 2016;31:e2016001.

    Google Scholar 

Download references

Acknowledgements

We are most grateful to the reviewers for many valuable detailed suggestions improving our initial draft and for indicating pertinent additional references.

Funding

The authors declare that they have no funding for this study.

Author information

Authors and Affiliations

Authors

Contributions

KH conceived of the study, developed the study design, organized the data, and provided an initial manuscript. HS elaborated the statistical methodology, processed and analyzed the data, designed and produced the Figures, gathered additional literature. Both authors edited and created the final paper. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Hagen Scherb.

Ethics declarations

Ethics approval and consent to participate

Not applicable. Ethics approval and consent to participate are not required and not necessary, since only publicly available data and previously published information is being used.

Consent for publication

Not applicable. Only anonymous data is being used.

Competing interests

The authors declare that they have no conflicts of interest.

Additional information

Publisher’s Note

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

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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Scherb, H., Hayashi, K. Spatiotemporal association of low birth weight with Cs-137 deposition at the prefecture level in Japan after the Fukushima nuclear power plant accidents: an analytical-ecologic epidemiological study. Environ Health 19, 82 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s12940-020-00630-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s12940-020-00630-w

Keywords