We discuss the use of fractal analysis of geophysical big data as a precursory method for the estimation of current seismic hazard. We begin by analyzing GPS time series of land surface displacements recorded in various regions of the planet exhibiting different degrees of seismic activity. We show how the quantitative measure of chaoticity based on a modified detrended fluctuation analysis (DFA) method can be used to distinguish critical and noncritical regimes of evolution of the Earth’s crust. Having established the reference values for the measure of chaoticity, we then choose California as the case study region and apply modified DFA to over 1200 time series for determining local areas currently exposed to higher seismic hazard. We comprehensively analyze the results obtained and compare them with those produced by conventional risk assessment methods. We demonstrate consistency of our results and yet discover that the suggested approach allows us to quantify two qualitatively different mechanisms involved in the growth of seismic hazards. We propose a hypothesis which relates these mechanisms to the known mechanisms of earthquake generation that evince a growth of collective behavior in the Earth’s crust on the eve of an earthquake. In order to verify the hypothesis, we analyze the same GPS data by an independent statistical method, canonical coherence analysis (CCA), which permits us to quantify the degree of synchronized behavior of constituents of the Earth’s crust. The CCA method confirms our previously obtained estimates, thus supporting fractal analysis of GPS time series as a suitable precursory method for current seismic hazard assessment.