Land potential assessment and trend-analysis using 2000-2021 FAPAR monthly time-series at 250 m spatial resolution

The article presents results of using remote sensing images and machine learning to map and assess land potential based on time-series of potential Fraction of Absorbed Photosynthetically Active Radiation (FAPAR) composites. Land potential here refers to the potential vegetation productivity in the hypothetical absence of short-term anthropogenic influence, such as intensive agriculture and urbanization. Knowledge on this ecological land potential could support the assessment of levels of land degradation as well as restoration potentials. Monthly aggregated FAPAR time-series of three percentiles (0.05, 0.50 and 0.95 probability) at 250 m spatial resolution were derived from the 8-day GLASS FAPAR V6 product for 2000-2021 and used to determine long-term trends in FAPAR, as well as to model potential FAPAR in the absence of human pressure. CCa 3 million training points sampled from 12,500 locations across the globe were overlaid with 68 bio-physical variables representing climate, terrain, landform, and vegetation cover, as well as several variables representing human pressure including: population count, cropland intensity, nightlights and a human footprint index. The training points were used in an ensemble machine learning model that stacks three base learners (extremely randomized trees, gradient descended trees and artificial neural network) using a linear regressor as meta-learner. The potential FAPAR was then projected by removing the impact of urbanization and intensive agriculture in the covariate layers. The results of strict cross-validation show that the global distribution of FAPAR can be explained with an R2 of 0.89, with the most important covariates being growing season length, forest cover indicator and annual precipitation. From this model, a global map of potential monthly FAPAR for the recent year (2021) was produced, and used to predict gaps in actual vs. potential FAPAR. The produced global maps of actual vs. potential FAPAR and long-term trends were each spatially matched with stable and transitional land cover classes. The assessment showed large negative FAPAR gaps (actual lower than potential) for classes: urban, needle-leave deciduous trees, and flooded shrub or herbaceous cover, while strong negative FAPAR trends were found for classes: urban, sparse vegetation and rainfed cropland. On the other hand, classes: irrigated or post-flooded cropland, tree cover mixed leaf type, and broad-leave deciduous showed largely positive trends. The framework allows land managers to assess potential land degradation from two aspects: as an actual declining trend in observed FAPAR and as a difference between actual and potential vegetation FAPAR.

Evaluating low-cost monitoring designs for PM2.5 exposure assessment with a spatiotemporal modeling approach

Determining the most feasible and cost-effective approaches to improving PM2.5 exposure assessment with low-cost monitors (LCMs) can considerably enhance the quality of its epidemiological inferences. We investigated features of fixed-site LCM designs that most impact PM2.5 exposure estimates to be used in long-term epidemiological inference for the Adult Changes in Thought Air Pollution (ACT-AP) study. We used ACT-AP collected and calibrated LCM PM2.5 measurements at the two-week level from April 2017 to September 2020 (N of monitors [measurements] = 82 [502]). We also acquired reference-grade PM2.5 measurements from January 2010 to September 2020 (N = 78 [6186]). We used a spatiotemporal modeling approach to predict PM2.5 exposures with either all LCM measurements or varying subsets with reduced temporal or spatial coverage. We evaluated the models based on a combination of cross-validation and external validation at locations of LCMs included in the models (N = 82), and also based on an independent external validation with a set of LCMs not used for the modeling (N = 30). We found that the model's performance declined substantially when LCM measurements were entirely excluded (spatiotemporal validation R2 [RMSE] = 0.69 [1.2 μg/m3]) compared to the model with all LCM measurements (0.84 [0.9 μg/m3]). Temporally, using the farthest apart measurements (i.e., the first and last) from each LCM resulted in the closest model's performance (0.79 [1.0 μg/m3]) to the model with all LCM data. The models with only the first or last measurement had decreased performance (0.77 [1.1 μg/m3]). Spatially, the model's performance decreased linearly to 0.74 (1.1 μg/m3) when only 10% of LCMs were included. Our analysis also showed that LCMs located in densely populated, road-proximate areas improved the model more than those placed in moderately populated, road-distant areas.

A comprehensive review of the development of land use regression approaches for modeling spatiotemporal variations of ambient air pollution: A perspective from 2011 to 2023

Land use regression (LUR) models are widely used in epidemiological and environmental studies to estimate humans' exposure to air pollution within urban areas. However, the early models, developed using linear regressions and data from fixed monitoring stations and passive sampling, were primarily designed to model traditional and criteria air pollutants and had limitations in capturing high-resolution spatiotemporal variations of air pollution. Over the past decade, there has been a notable development of multi-source observations from low-cost monitors, mobile monitoring, and satellites, in conjunction with the integration of advanced statistical methods and spatially and temporally dynamic predictors, which have facilitated significant expansion and advancement of LUR approaches. This paper reviews and synthesizes the recent advances in LUR approaches from the perspectives of the changes in air quality data acquisition, novel predictor variables, advances in model-developing approaches, improvements in validation methods, model transferability, and modeling software as reported in 155 LUR studies published between 2011 and 2023. We demonstrate that these developments have enabled LUR models to be developed for larger study areas and encompass a wider range of criteria and unregulated air pollutants. LUR models in the conventional spatial structure have been complemented by more complex spatiotemporal structures. Compared with linear models, advanced statistical methods yield better predictions when handling data with complex relationships and interactions. Finally, this study explores new developments, identifies potential pathways for further breakthroughs in LUR methodologies, and proposes future research directions. In this context, LUR approaches have the potential to make a significant contribution to future efforts to model the patterns of long- and short-term exposure of urban populations to air pollution.

Improved daily estimates of relative humidity at high resolution across Germany: A random forest approach

The lack of readily available methods for estimating high-resolution near-surface relative humidity (RH) and the incapability of weather stations to fully capture the spatiotemporal variability can lead to exposure misclassification in studies of environmental epidemiology. We therefore aimed to predict German-wide 1 × 1 km daily mean RH during 2000-2021. RH observations, longitude and latitude, modelled air temperature, precipitation and wind speed as well as remote sensing information on topographic elevation, vegetation, and the true color band composite were incorporated in a Random Forest (RF) model, in addition to date for capturing the temporal variations of the response-explanatory variables relationship. The model achieved high accuracy (R2 = 0.83) and low errors (Root Mean Square Error (RMSE) of 5.07%, Mean Absolute Percentage Error (MAPE) of 5.19% and Mean Percentage Error (MPE) of - 0.53%), calculated via ten-fold cross-validation. A comparison of our RH predictions with measurements from a dense monitoring network in the city of Augsburg, South Germany confirmed the good performance (R2 ≥ 0.86, RMSE ≤ 5.45%, MAPE ≤ 5.59%, MPE ≤ 3.11%). The model displayed high German-wide RH (22y-average of 79.00%) and high spatial variability across the country, exceeding 12% on yearly averages. Our findings indicate that the proposed RF model is suitable for estimating RH for a whole country in high-resolution and provide a reliable RH dataset for epidemiological analyses and other environmental research purposes.

High-resolution spatiotemporal modeling of daily near-surface air temperature in Germany over the period 2000-2020

The commonly used weather stations cannot fully capture the spatiotemporal variability of near-surface air temperature (T_air), leading to exposure misclassification and biased health effect estimates. We aimed to improve the spatiotemporal coverage of T_air data in Germany by using multi-stage modeling to estimate daily 1 × 1 km minimum (T_min), mean (T_mean), maximum (T_max) T_air and diurnal T_air range during 2000-2020. We used weather station T_air observations, satellite-based land surface temperature (LST), elevation, vegetation and various land use predictors. In the first stage, we built a linear mixed model with daily random intercepts and slopes for LST adjusted for several spatial predictors to estimate T_air from cells with both T_air and LST available. In the second stage, we used this model to predict T_air for cells with only LST available. In the third stage, we regressed the second stage predictions against interpolated T_air values to obtain T_air countrywide. All models achieved high accuracy (0.91 ≤ R² ≤ 0.98) and low errors (1.03 °C ≤ Root Mean Square Error (RMSE) ≤ 2.02 °C). Validation with external data confirmed the good performance, locally, i.e., in Augsburg for all models (0.74 ≤ R² ≤ 0.99, 0.87 °C ≤ RMSE ≤ 2.05 °C) and countrywide, for the T_mean model (0.71 ≤ R² ≤ 0.99, 0.79 °C ≤ RMSE ≤ 1.19 °C). Annual T_mean averages ranged from 8.56 °C to 10.42 °C with the years beyond 2016 being constantly hotter than the 21-year average. The spatial variability within Germany exceeded 15 °C annually on average following patterns including mountains, rivers and urbanization. Using a case study, we showed that modeling leads to broader T_air variability representation for exposure assessment of participants in health cohorts. Our results indicate the proposed models as suitable for estimating nationwide T_air at high resolution. Our product is critical for temperature-based epidemiological studies and is also available for other research purposes.

Estimating intervention effects on infectious disease control: The effect of community mobility reduction on Coronavirus spread

Understanding the effects of interventions, such as restrictions on community and large group gatherings, is critical to controlling the spread of COVID-19. Susceptible-Infectious-Recovered (SIR) models are traditionally used to forecast the infection rates but do not provide insights into the causal effects of interventions. We propose a spatiotemporal model that estimates the causal effect of changes in community mobility (intervention) on infection rates. Using an approximation to the SIR model and incorporating spatiotemporal dependence, the proposed model estimates a direct and indirect (spillover) effect of intervention. Under an interference and treatment ignorability assumption, this model is able to estimate causal intervention effects, and additionally allows for spatial interference between locations. Reductions in community mobility were measured by cell phone movement data. The results suggest that the reductions in mobility decrease Coronavirus cases 4 to 7 weeks after the intervention.

Image segmentation and separation of spectrally similar dyes in fluorescence microscopy by dynamic mode decomposition of photobleaching kinetics

Background

Image segmentation in fluorescence microscopy is often based on spectral separation of fluorescent probes (color-based segmentation) or on significant intensity differences in individual image regions (intensity-based segmentation). These approaches fail, if dye fluorescence shows large spectral overlap with other employed probes or with strong cellular autofluorescence.

Results

Here, a novel model-free approach is presented which determines bleaching characteristics based on dynamic mode decomposition (DMD) and uses the inferred photobleaching kinetics to distinguish different probes or dye molecules from autofluorescence. DMD is a data-driven computational method for detecting and quantifying dynamic events in complex spatiotemporal data. Here, DMD is first used on synthetic image data and thereafter used to determine photobleaching characteristics of a fluorescent sterol probe, dehydroergosterol (DHE), compared to that of cellular autofluorescence in the nematode Caenorhabditis elegans. It is shown that decomposition of those dynamic modes allows for separating probe from autofluorescence without invoking a particular model for the bleaching process. In a second application, DMD of dye-specific photobleaching is used to separate two green-fluorescent dyes, an NBD-tagged sphingolipid and Alexa488-transferrin, thereby assigning them to different cellular compartments.

Conclusions

Data-based decomposition of dynamic modes can be employed to analyze spatially varying photobleaching of fluorescent probes in cells and tissues for spatial and temporal image segmentation, discrimination of probe from autofluorescence and image denoising. The new method should find wide application in analysis of dynamic fluorescence imaging data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04881-x.

Forest tree species distribution for Europe 2000-2020: mapping potential and realized distributions using spatiotemporal machine learning

This article describes a data-driven framework based on spatiotemporal machine learning to produce distribution maps for 16 tree species (Abies alba Mill., Castanea sativa Mill., Corylus avellana L., Fagus sylvatica L., Olea europaea L., Picea abies L. H. Karst., Pinus halepensis Mill., Pinus nigra J. F. Arnold, Pinus pinea L., Pinus sylvestris L., Prunus avium L., Quercus cerris L., Quercus ilex L., Quercus robur L., Quercus suber L. and Salix caprea L.) at high spatial resolution (30 m). Tree occurrence data for a total of three million of points was used to train different algorithms: random forest, gradient-boosted trees, generalized linear models, k-nearest neighbors, CART and an artificial neural network. A stack of 305 coarse and high resolution covariates representing spectral reflectance, different biophysical conditions and biotic competition was used as predictors for realized distributions, while potential distribution was modelled with environmental predictors only. Logloss and computing time were used to select the three best algorithms to tune and train an ensemble model based on stacking with a logistic regressor as a meta-learner. An ensemble model was trained for each species: probability and model uncertainty maps of realized distribution were produced for each species using a time window of 4 years for a total of six distribution maps per species, while for potential distributions only one map per species was produced. Results of spatial cross validation show that the ensemble model consistently outperformed or performed as good as the best individual model in both potential and realized distribution tasks, with potential distribution models achieving higher predictive performances (TSS = 0.898, R2 logloss = 0.857) than realized distribution ones on average (TSS = 0.874, R2 logloss = 0.839). Ensemble models for Q. suber achieved the best performances in both potential (TSS = 0.968, R2 logloss = 0.952) and realized (TSS = 0.959, R2 logloss = 0.949) distribution, while P. sylvestris (TSS = 0.731, 0.785, R2 logloss = 0.585, 0.670, respectively, for potential and realized distribution) and P. nigra (TSS = 0.658, 0.686, R2 logloss = 0.623, 0.664) achieved the worst. Importance of predictor variables differed across species and models, with the green band for summer and the Normalized Difference Vegetation Index (NDVI) for fall for realized distribution and the diffuse irradiation and precipitation of the driest quarter (BIO17) being the most frequent and important for potential distribution. On average, fine-resolution models outperformed coarse resolution models (250 m) for realized distribution (TSS = +6.5%, R2 logloss = +7.5%). The framework shows how combining continuous and consistent Earth Observation time series data with state of the art machine learning can be used to derive dynamic distribution maps. The produced predictions can be used to quantify temporal trends of potential forest degradation and species composition change.

Spatiotemporal modeling of land subsidence using a geographically weighted deep learning method based on PS-InSAR

The demand for water resources during urbanization forces the continuous exploitation of groundwater, resulting in dramatic piezometric drawdown and inducing regional land subsidence (LS). This has greatly threatened sustainable development in the long run. LS modeling helps understanding the factors responsible for the ongoing loss of land elevation and hence enhances the development of prevention strategies. Data-driven LS models perform well with fewer variables and faster convergence than physically-based hydrogeological models. However, the former models often cannot simultaneously reflect the temporal nonlinearity and spatial correlation (SC) characteristics of LS under complex variables. We proposed a LS spatiotemporal model which considers both nonlinear and spatial correlations between LS and groundwater level change of exploited aquifers. It is based on deep learning method and LS time series detected by permanent scatterer-interferometric synthetic aperture radar (PS-InSAR). The LS time series and hydrogeological properties are constructed as a spatiotemporal dataset for model training. The spatiotemporal LS model, geographically weighted long short-term memory (GW-LSTM), is constructed by integrating SC with LSTM. This latter is a deep recurrent neural network approach incorporating sequential data. The model is validated by a case study in the Beijing plain. The results show that the accuracy of the proposed model can be greatly improved considering the spatial correlation between LS and influencing factors. Furthermore, the comparison between the LSTM and GW-LSTM models reveals that groundwater level variation is not a unique causation of LS in the study area. The developed model deals with the spatiotemporal characteristics of LS under multiple variables and can be used to predict LS under different scenarios of groundwater level variations for the purpose of monitoring and providing evidence to support the prevention of future LS.

Spatial temporal distribution of COVID-19 risk during the early phase of the pandemic in Malawi

Background

COVID-19 has been one of the greatest challenges the world has faced since the second world war. This study aimed at investigating the distribution of COVID-19 in both space and time in Malawi.

Methods

The study used publicly available data of COVID-19 cases for the period from 2 April 2020 to 28 October 2020. Semiparametric spatial temporal models were fitted to the number of monthly confirmed cases as an outcome data, with time and district as independent variables, where district was the spatial unit, while accounting for sociodemographic factors.

Results

The study found significant effects of location and time, with the two interacting. The spatial distribution of COVID-19 risk showed major cities being at greater risk than rural areas. Over time, the COVID-19 risk was increasing then decreasing in most districts with the rural districts being consistently at lower risk. High proportion of elderly people was positively associated with COVID-19 risk (β = 1.272, 95% CI [0.171, 2.370]) than low proportion of elderly people. There was negative association between poverty incidence and COVID-19 risk (β = −0.100, 95% CI [−0.136, −0.065]).

Conclusion

Future or present strategies to limit the spread of COVID-19 should target major cities and the focus should be on time periods that had shown high risk. Furthermore, the focus should be on elderly and rich people.

Comparison of Machine Learning and Land Use Regression for fine scale spatiotemporal estimation of ambient air pollution: Modeling ozone concentrations across the contiguous United States

BACKGROUND

Spatial linear Land-Use Regression (LUR) is commonly used for long-term modeling of air pollution in support of exposure and epidemiological assessments. Machine Learning (ML) methods in conjunction with spatiotemporal modeling can provide more flexible exposure-relevant metrics and have been studied using different model structures. There is however a lack of comparisons of methods available within these two modeling frameworks, that can guide model/algorithm selection in air quality epidemiology.

OBJECTIVE

The present study compares thirteen algorithms for spatial/spatiotemporal modeling applied for daily maxima of 8-hour running averages of ambient ozone concentrations at spatial resolutions corresponding to census tracts, to support estimation of annual ozone design values across the contiguous US. These algorithms were selected from nine representative categories and trained using predictors that included chemistry-transport model predictions, meteorological factors, land use and land cover, and stationary and mobile emissions.

METHODS

To obtain the best predictive performance, model structures were optimized through a repeated coarse/fine grid search with expert knowledge. Six target-oriented validation strategies were used to prevent overfitting and avoid over-optimistic model evaluation results. In order to take full advantage of the power of different algorithms, we introduced tuning sample weights in spatiotemporal modeling to ensure predictive accuracy of peak concentrations, that is crucial for exposure assessments. In spatial modeling, four interpretation and visualization tools were introduced to explain predictions from different algorithms.

RESULTS

Nonlinear ML methods achieved higher prediction accuracy than linear LUR, and the improvements were more significant for spatiotemporal modeling (nearly 10%-40% decrease of predicted RMSE). By tuning the sample weights, spatiotemporal models can predict concentrations used to calculate ozone design values that are comparable or even better than spatial models (nearly 30% decrease of cross-validated RMSE). We visualized the underlying nonlinear relationships, heterogeneous associations and complex interactions from the two best performing ML algorithms, i.e., Random Forest and Extreme Gradient Boosting, and found that the complex patterns were relatively less significant with respect to model accuracy for spatial modeling.

CONCLUSION

Machine Learning can provide estimates that are actually more interpretable and practical than linear regression to improve accuracy in modeling human exposures. A careful design of hyperparameter tuning and flexible data splitting and validations is crucial to obtain reliable and stable results. Desirable/successful nonlinear models are expected to capture similar nonlinear patterns and interactions using different ML algorithms.

Spatiotemporal trends of PM2.5 concentrations in central China from 2003 to 2018 based on MAIAC-derived high-resolution data

Long-term PM_2.5 levels with high precision at fine spatiotemporal resolution are essential for quantitatively understanding the health risk of exposure to ambient fine particulate matter (PM_2.5) and making effective air pollution control policies. The emerging statistically derived PM_2.5 estimations from satellite remote sensing observations of aerosol optical depth (AOD) data are an effective alternative to reconstruct global, long-term, high spatiotemporal resolution PM_2.5 information. However, studies on PM_2.5 estimation and its application to exposure and health-related studies are limited in China due to the lack of historical in-situ measurements before 2013. In this study, we explored the long-term trends of PM_2.5 exposure in central China, a hotspot that has recently been experiencing severe particulate pollution, at the local scale. We first developed a spatiotemporal model incorporating periodical characteristics within the data to estimate daily concentrations of historical PM_2.5 at a fine scale of 1 km for 2003-2018. The linear effects of predictors including AOD, meteorological and land-use parameters and the non-linear interaction between AOD and meteorological parameters were considered in the modeling process. The most recently released high-resolution satellite aerosol product, Multi-Angle Implementation of Atmospheric Correction (MAIAC) was used to help to represent the fine-scale particle gradients. Our daily estimates correlated well with in-situ observations (cross-validation R² = 0.59), achieving precision comparable to previous statistical models. Through linking with gridded demographic data, the population-weighted PM_2.5 average during 2003 to 2018 was found to be high (62.23 μg/m³ for the whole domain) with obvious spatial variations and seasonality. An inverse U pattern was seen in the time series, with two inflection points around 2008 and 2015. Our model provides reliable particulate information with high spatial resolution and long-term temporal coverage, which can inform local-scale PM_2.5-related epidemiological studies and health-risk assessments for central China.

Phenology-adjusted dynamic curve number for improved hydrologic modeling

The Soil Conservation Service Curve Number (SCS-CN, or CN) is a widely used method to estimate runoff from rainfall events. It has been adapted to many parts of the world with different land uses, land cover types, and climatic conditions and successfully applied to situations ranging from simple runoff calculations and land use change assessment to comprehensive hydrologic/water quality simulations. However, the CN method lacks the ability to incorporate seasonal variations in vegetated surface conditions, and unnoticed landuse/landcover (LULC) change that shape infiltration and storm runoff. Plant phenology is a main determinant of changes in hydrologic processes and water balances across seasons through its influence on surface roughness and evapotranspiration. This study used regression analysis to develop a dynamic CN (CN_NDVI) based on seasonal variations in the remotely-sensed Normalized Difference Vegetation Index (NDVI) to monitor intra-annual plant phenological development. A time series of 16-day MODIS NDVI (MOD13Q1 Collection 5) images were used to monitor vegetation development and provide NDVI data necessary for CN_NDVI model calibration and validation. Twelve years of rainfall and runoff data (2001-2012) from four small watersheds located in the Konza Prairie Biological Station, Kansas were used to develop, calibrate, and validate the method. Results showed CN_NDVI performed significantly better in predicting runoff with calibrated CN_NDVI runoff increasing by approximately 0.74 for every unit increase in observed runoff compared to 0.46 for SCS-CN runoff and was more highly correlated to observed runoff (r = 0.78 vs. r = 0.38). In addition, CN_NDVI runoff had better NSE (0.53) and PBIAS (4.22) compared to the SCS-CN runoff (-0.87 and -94.86 respectively). In the validated model, CN_NDVI runoff increased by approximately 0.96 for every unit of observed runoff, while SCS-CN runoff increased by 0.49. Validated runoff was also better correlated to observed runoff than SCS-CN runoff (r = 0.52 vs. r = 0.33). These findings suggest that the CN_NDVI can yield improved estimates of surface runoff from precipitation events, leading to more informed water and land management decisions.

Spatial infectious disease epidemiology: on the cusp

Infectious diseases continue to pose a significant public health burden despite the great progress achieved in their prevention and control over the last few decades. Our ability to disentangle the factors and mechanisms driving their propagation in space and time has dramatically advanced in recent years. The current era is rich in mathematical and computational tools and detailed geospatial information, including sociodemographic, geographic, and environmental data, which are essential to elucidate key drivers of infectious disease transmission from epidemiological and genetic data. Indeed, this paradigm shift was driven by dramatic advances in complex systems approaches along with substantial improvements in data availability and computational power. The burgeoning output of infectious disease spatial modeling suggests that we are close to a fully integrated approach for early epidemic detection and intervention. This special collection in BMC Medicine aims to bring together a broad range of quantitative investigations that improve our understanding of the spatiotemporal transmission dynamics of infectious diseases in order to mitigate their impact on the human population.

Applying spatiotemporal models to monitoring data to quantify fish population responses to the Deepwater Horizon oil spill in the Gulf of Mexico

Quantifying the impacts of disturbances such as oil spills on marine species can be challenging. Natural environmental variability, human responses to the disturbance (e.g., fisheries closures), the complex life histories of the species being monitored, and limited pre-spill data can make detection of effects of oil spills difficult. Using long-term monitoring data from the state of Louisiana (USA), we applied novel spatiotemporal approaches to identify anomalies in species occurrence and catch rates. We included covariates (salinity, temperature, turbidity) to help isolate unusual events. While some species showed evidence of unlikely temporal anomalies in occurrence or catch rates, we found that the majority of the observed anomalies were also before the Deepwater Horizon event. Several species-gear combinations suggested upticks in the spatial variability immediately following the spill, but most species indicated no trend. Across species-gear combinations, there was no clear evidence for synchronous or asynchronous responses in occurrence or catch rates across sites following the spill. Our results are in general agreement to other analyses of monitoring data that detected small impacts, but in contrast to recent results from ecological modeling that showed much larger effects of the oil spill on fish and shellfish.

Satellite-based high-resolution PM2.5 estimation over the Beijing-Tianjin-Hebei region of China using an improved geographically and temporally weighted regression model

Ground fine particulate matter (PM2.5) concentrations at high spatial resolution are substantially required for determining the population exposure to PM2.5 over densely populated urban areas. However, most studies for China have generated PM2.5 estimations at a coarse resolution (≥10 km) due to the limitation of satellite aerosol optical depth (AOD) product in spatial resolution. In this study, the 3 km AOD data fused using the Moderate Resolution Imaging Spectroradiometer (MODIS) Collection 6 AOD products were employed to estimate the ground PM2.5 concentrations over the Beijing-Tianjin-Hebei (BTH) region of China from January 2013 to December 2015. An improved geographically and temporally weighted regression (iGTWR) model incorporating seasonal characteristics within the data was developed, which achieved comparable performance to the standard GTWR model for the days with paired PM_2.5- AOD samples (Cross-validation (CV) R² = 0.82) and showed better predictive power for the days without PM_2.5- AOD pairs (the R² increased from 0.24 to 0.46 in CV). Both iGTWR and GTWR (CV R² = 0.84) significantly outperformed the daily geographically weighted regression model (CV R² = 0.66). Also, the fused 3 km AODs improved data availability and presented more spatial gradients, thereby enhancing model performance compared with the MODIS original 3/10 km AOD product. As a result, ground PM2.5 concentrations at higher resolution were well represented, allowing, e.g., short-term pollution events and long-term PM2.5 trend to be identified, which, in turn, indicated that concerns about air pollution in the BTH region are justified despite its decreasing trend from 2013 to 2015.

Spatiotemporal analysis of brucellosis incidence in Iran from 2011 to 2014 using GIS

OBJECTIVE

To investigate the distribution and trends associated with brucellosis incidence rates in Iran from 2011 to 2014.

METHODS

The reported incidence rates of brucellosis for the years 2011-2014 were collected and entered into GIS 10.1. The Cochran-Armitage test for linear trends, choropleth maps, hot-spot analysis, and high-low clustering analysis were used to investigate patterns of the disease over the study period and by season, and to identify high-risk areas and any clustering of the disease. The significance level was set at p<0.05.

RESULTS

A total of 68493 cases of brucellosis were reported during the study period, giving an average brucellosis incidence rate for this period of 38.67/100000. In 2011, the highest rate of brucellosis was observed in Koohrang County of Chaharmahal-Bakhtiari Province, with 317/100 000. In the subsequent years, 2012-2014, Charuymaq County of East-Azerbaijan Province had incidence rates of 384, 534, and 583/100000, respectively. However, the incidence rate of the disease did not follow a linear trend (p<0.001). The maximum and minimum incidence rates of the disease occurred in mid-summer and mid-winter, respectively. The results of the hot-spot analysis showed that the distribution of the disease was highest in the mountainous areas of Iran, particularly along the Zagros mountain range and in most cities near the Zagros Mountains (p<0.01). In addition, the cluster analysis showed a clustering pattern in these high incidence areas (p<0.01).

CONCLUSIONS

There were significant differences in the geographic distribution of brucellosis, with the incidence rates being highest in most of the cities in the west and north-west of the country. The incidence of this disease also increased during the summer. It is important to take these patterns into account when allocating resources to combat this disease and to ensure that health programs and other interventions focus on the areas of greatest need.

Spatial Multiresolution Analysis of the Effect of PM2.5 on Birth Weights

Fine particulate matter (PM_2.5) measured at a given location is a mix of pollution generated locally and pollution traveling long distances in the atmosphere. Therefore, the identification of spatial scales associated with health effects can inform on pollution sources responsible for these effects, resulting in more targeted regulatory policy. Recently, prediction methods that yield high-resolution spatial estimates of PM_2.5 exposures allow one to evaluate such scale-specific associations. We propose a two-dimensional wavelet decomposition that alleviates restrictive assumptions required for standard wavelet decompositions. Using this method we decompose daily surfaces of PM_2.5 to identify which scales of pollution are most associated with adverse health outcomes. A key feature of the approach is that it can remove the purely temporal component of variability in PM_2.5 levels and calculate effect estimates derived solely from spatial contrasts. This eliminates the potential for unmeasured confounding of the exposure - outcome associations by temporal factors, such as season. We apply our method to a study of birth weights in Massachusetts, U.S.A from 2003-2008 and find that both local and urban sources of pollution are strongly negatively associated with birth weight. Results also suggest that failure to eliminate temporal confounding in previous analyses attenuated the overall effect estimate towards zero, with the effect estimate growing in magnitude once this source of variability is removed.

Spatiotemporal distribution modeling of PET tracer uptake in solid tumors

OBJECTIVE

Distribution of PET tracer uptake is elaborately modeled via a general equation used for solute transport modeling. This model can be used to incorporate various transport parameters of a solid tumor such as hydraulic conductivity of the microvessel wall, transvascular permeability as well as interstitial space parameters. This is especially significant because tracer delivery and drug delivery to solid tumors are determined by similar underlying tumor transport phenomena, and quantifying the former can enable enhanced prediction of the latter.

METHODS

We focused on the commonly utilized FDG PET tracer. First, based on a mathematical model of angiogenesis, the capillary network of a solid tumor and normal tissues around it were generated. The coupling mathematical method, which simultaneously solves for blood flow in the capillary network as well as fluid flow in the interstitium, is used to calculate pressure and velocity distributions. Subsequently, a comprehensive spatiotemporal distribution model (SDM) is applied to accurately model distribution of PET tracer uptake, specifically FDG in this work, within solid tumors.

RESULTS

The different transport mechanisms, namely convention and diffusion from vessel to tissue and in tissue, are elaborately calculated across the domain of interest and effect of each parameter on tracer distribution is investigated. The results show the convection terms to have negligible effect on tracer transport and the SDM can be solved after eliminating these terms.

CONCLUSION

The proposed framework of spatiotemporal modeling for PET tracers can be utilized to comprehensively assess the impact of various parameters on the spatiotemporal distribution of PET tracers.

Impact of meteorological factors on the spatiotemporal patterns of dengue fever incidence

Dengue fever is one of the most widespread vector-borne diseases and has caused more than 50 million infections annually over the world. For the purposes of disease prevention and climate change health impact assessment, it is crucial to understand the weather-disease associations for dengue fever. This study investigated the nonlinear delayed impact of meteorological conditions on the spatiotemporal variations of dengue fever in southern Taiwan during 1998-2011. We present a novel integration of a distributed lag nonlinear model and Markov random fields to assess the nonlinear lagged effects of weather variables on temporal dynamics of dengue fever and to account for the geographical heterogeneity. This study identified the most significant meteorological measures to dengue fever variations, i.e., weekly minimum temperature, and the weekly maximum 24-hour rainfall, by obtaining the relative risk (RR) with respect to disease counts and a continuous 20-week lagged time. Results show that RR increased as minimum temperature increased, especially for the lagged period 5-18 weeks, and also suggest that the time to high disease risks can be decreased. Once the occurrence of maximum 24-hour rainfall is >50 mm, an associated increased RR lasted for up to 15 weeks. A temporary one-month decrease in the RR of dengue fever is noted following the extreme rain. In addition, the elevated incidence risk is identified in highly populated areas. Our results highlight the high nonlinearity of temporal lagged effects and magnitudes of temperature and rainfall on dengue fever epidemics. The results can be a practical reference for the early warning of dengue fever.

Reduced-Rank Spatio-Temporal Modeling of Air Pollution Concentrations in the Multi-Ethnic Study of Atherosclerosis and Air Pollution

There is growing evidence in the epidemiologic literature of the relationship between air pollution and adverse health outcomes. Prediction of individual air pollution exposure in the Environmental Protection Agency (EPA) funded Multi-Ethnic Study of Atheroscelerosis and Air Pollution (MESA Air) study relies on a flexible spatio-temporal prediction model that integrates land-use regression with kriging to account for spatial dependence in pollutant concentrations. Temporal variability is captured using temporal trends estimated via modified singular value decomposition and temporally varying spatial residuals. This model utilizes monitoring data from existing regulatory networks and supplementary MESA Air monitoring data to predict concentrations for individual cohort members. In general, spatio-temporal models are limited in their efficacy for large data sets due to computational intractability. We develop reduced-rank versions of the MESA Air spatio-temporal model. To do so, we apply low-rank kriging to account for spatial variation in the mean process and discuss the limitations of this approach. As an alternative, we represent spatial variation using thin plate regression splines. We compare the performance of the outlined models using EPA and MESA Air monitoring data for predicting concentrations of oxides of nitrogen (NO _x )-a pollutant of primary interest in MESA Air-in the Los Angeles metropolitan area via cross-validated R². Our findings suggest that use of reduced-rank models can improve computational efficiency in certain cases. Low-rank kriging and thin plate regression splines were competitive across the formulations considered, although TPRS appeared to be more robust in some settings.

REGULARIZED IMAGE RECONSTRUCTION FOR PS MODEL-BASED CARDIOVASCULAR MRI

Real-time cardiovascular MRI is a useful and challenging dynamic imaging application. The partial separability (PS) model enables reconstruction of dynamic cardiac images from highly undersampled (k, t)-space data. However, the underlying PS model-based reconstruction problem is ill-conditioned, so regularization is often necessary to stabilize its solution. It has been shown that ℓ₁ regularization is useful for finding sparse solutions, and ℓ₂ regularization is widely used to incorporate anatomical constraints. An important practical question is which regularization scheme to use for PS model-based cardiovascular imaging. We address this problem by implementing both schemes and evaluating their performances in terms of reconstruction error, image artifacts, image noise, computation time, and performance characterizability. The ℓ₁-regularized results exhibit lower reconstruction error, artifact energy, and noise variance, while ℓ₂ regularization is much faster and produces predictable reconstruction results. This study indicates that the ℓ₁ scheme is preferable when image quality is the main concern.