Introduction. The best known scaling relation for earthquakes is the Gutenberg-Richter law (G-R).1) It states that the (cumulative) number of earthquakes with magnitude greater than M occurring in a specified area and time is given by N(> M) ~ 10 -bM [1] where b is a constant. Numerous publications have examined the spatial and temporal variations of the bvalue, but it is currently considered 2) that it is generally a constant varying only slightly from region to region being in the range 0.8 ≤ b ≤ 1.2. It has been recognized that G-R belongs to a broad range of natural phenomena that exhibit fractal scaling.2),3) Such a scaling reflects that the number of earthquakes (occurring in a specified area and time) with rupture areas greater than A is given by 2)The main aim of the present paper could be described as an attempt towards understanding the origin of the aforementioned constancy of the b-or D-value. Nowadays, it is generally accepted that main shocks can be considered in the frame of critical phenomena 5) (e.g. In the frame of a single fault representation, the stick-slip mechanism is usually considered, in which the friction properties play the prominent role; we then additionally take into account that earthquakes occur in large areas characterized by a diversity of fault sizes and depths. Along these lines, the role of the geometry of the fault profiles, for example, in earthquake dynamics has been Abstract: It is explained, from first Principles, why in the Gutenberg-Richter law (stating that the cumulative number of earthquakes N(> M) with magnitude greater than M is given by N(> M) ~ 10 -bM ) the so called b-value is usually found to be around unity varying only slightly from region to region. The explanation is achieved just by applying the analysis in the natural time domain, without using any adjustable parameter.