We numerically study the two-dimensional, area preserving, web map. When the map is governed by ergodic behavior, it is, as expected, correctly described by Boltzmann-Gibbs statistics, based on the additive entropic functional $S_{BG}\left[p\left(x\right)\right]=-k\int dxp\left(x\right)lnp\left(x\right)$. In contrast, possible ergodicity breakdown and transitory sticky dynamical behavior drag the map into the realm of generalized $q$ statistics, based on the nonadditive entropic functional $S_q\left[p\left(x\right)\right]=k\frac{1-\int dx{\left[p\left(x\right)\right]}^q}{q-1}$ ($q\in R;S_1=S_{BG}$).

We statistically describe the system (probability distribution of the sum of successive iterates, sensitivity to the initial condition, and entropy production per unit time) for typical values of the parameter that controls the ergodicity of the map. For small (large) values of the external parameter $K$, we observe $q$-Gaussian distributions with $q=1.935\cdots$ (Gaussian distributions), like for the standard map. In contrast, for intermediate values of $K$, we observe a different scenario, due to the fractal structure of the trajectories embedded in the chaotic sea. Long-standing non-Gaussian distributions are characterized in terms of the kurtosis and the box-counting dimension of chaotic sea.

G. Ruiz, U. Tirnakli, E.P. Borges, C. Tsallis, Statistical characterization of discrete conservative systems: The web map, Physical Review E 96 (2017) 042158