In Kang (2015) we calculated the acceleration of cosmic-ray electrons at weak spherical shocks that are expected to form in the cluster outskirts, and estimated the diffuse synchrotron radiation emitted by those electrons. There we demonstrated that, at decelerating spherical shocks, the volume integrated spectra of both electrons and radiation deviate significantly from the test-particle power-laws predicted for constant planar shocks, because the shock compression ratio and the flux of inject electrons decrease in time. In this study, we consider spherical blast waves propagating through a constant density core surrounded by an isothermal halo with ρ ∝ r ⁻ⁿ in order to explore how the deceleration of the shock affects the radio emission from accelerated electrons. The surface brightness profile and the volumeintegrated radio spectrum of the model shocks are calculated by assuming a ribbon-like shock surface on a spherical shell and the associated downstream region of relativistic electrons. If the postshock magnetic field strength is about 0.7 or 7 µ G, at the shock age of ∼ 50 Myr, the volume-integrated radio spectrum steepens gradually with the spectral index from α _inj to α _inj + 0.5 over 0.1–10 GHz, where α _inj is the injection index at the shock position expected from the diffusive shock acceleration theory. Such gradual steepening could explain the curved radio spectrum of the radio relic in cluster A2266, which was interpreted as a broken power-law by Trasatti et al. (2015) , if the relic shock is young enough so that the break frequency is around 1 GHz. Cosmological hydrodynamic simulations have shown that shock waves may form due to supersonic flows in the baryonic intracluster medium (ICM) during the formation of the large scale structure in the Universe (e.g., Ryu et al. 2003 ; Kang et al. 2007 ; Vazza et al. 2009 ; Skillman et al. 2011 ). The time evolution of such shocks in numerical simulations (e.g., Vazza et al. 2012 ) indicates that shock surfaces behave like spherical bubbles blowing out from the cluster center during major episodes of mergers or infalls from adjacent filaments. Shock surfaces could last for the order of the dynamical time scale of clusters, t _dyn ∼ 1 Gyr. This implies that some cosmological shock waves are associated with merger induced outflows, and hence shock bubbles are likely to expand into the cluster outskirts where the density profile decreases with radial distance from the cluster center. Observational evidence for shock bubbles can be found in the so-called “radio relics” detected in the outskirts of galaxy clusters, which are interpreted as diffuse synchrotron emitting structures containing relativistic electrons accelerated at weak ICM shocks ( M_s ∼ 2 − 4) (e.g., van Weeren et al. 2010 , 2011 ; Nuza et al. 2012 ; Feretti et al. 2012 ; Brüggen et al. 2012 ; Brunetti & Jones 2014 ). Some radio relics, for instance, the “Sausage relic” in cluster CIZA J2242.8+5301 and the “Toothbrush” relic in cluster 1RXS J0603.3+4214, have thin arc-like shapes of ∼ 50 kpc in width and ∼ 1 − 2 Mpc in length ( van Weeren et al. 2010 , 2012 ). They could be represented as a ribbon-like structure on a spherical shell and the associated downstream volume of radio-emitting electrons, projected onto the sky plane (e.g., van Weeren et al. 2010 ; Kang et al. 2012 ). As shown in Figure 1 , the viewing depth of a relic shock structure can be parameterized by the extension angle 2 ψ , while the physical width of the postshock volume of radiating electrons is mainly determined by the advection length, ∆ l _adv ≈ u ₂ · t _age , for low energy electrons or the cooling length, ∆ l _cool ( γ_e ) ≈ u ₂ · t _rad ( γ_e ) for high energy electrons. Here u ₂ is the flow speed behind the shock and t _rad ( γ_e ) is the radiative cooling time scale for electrons with Lorentz factor, γ_e . It remains largely unknown how such a ribbon-like structure can be formed by the spherical outflows in galaxy clusters. We note, however, a recent study by Shimwell et al. (2015) who suggested that the uniform arc-like shape of some radio relics may trace the underlying region of pre-existing seed electrons remaining from an old radio lobe. According to diffusive shock acceleration (DSA) theory, cosmic-ray (CR) particles can be generated via Fermi 1st order process at collisionless shocks ( Bell 1978 ; Drury 1983 ; Malkov & Drury 2001 ). In the testparticle regime, the CR energy spectrum at the position of the shock has a power-law form, N ( E ) ∝ E ^−s , where s = ( σ + 2)/( σ − 1) and σ = ρ ₂ / ρ ₁ is the shock compression ratio. Hereafter, we use the sub- scripts ‘1’ and ‘2’ to denote the conditions upstream and downstream of shock, respectively. Then the synchrotron radiation spectrum due to these CR electrons has a power-law form of j_ν ( x_s ) ∝ ν − ^αinj , where is the injection index at the shock. Moreover, the volumeintegrated synchrotron spectrum downstream of a planar shock becomes a simple power-law of J_ν ∝ ν ^−A_ν with A_ν = α _inj +0.5 above a break frequency, since electrons cool via synchrotron and inverse-Compton (iC) losses behind the shock (e.g., Kang 2011 ). Such spectral characteristics are commonly used to infer the shock Mach number of observed radio relics (e.g., van Weeren et al. 2010 ; Stroe et al. 2014 ) In Kang (2015) (Paper I), we calculated the electron acceleration at spherical shocks similar to Sedov-Taylor blast waves with M_s ∼ 2.5 − 4.5, which expand into a hot uniform ICM. We found that the electron energy spectrum at the shock location reaches a steady state defined by the instantaneous shock parameters. Hence the spatially resolved, synchrotron radiation spectra at the shock can be described properly by the test-particle DSA predictions for steady planar shocks. However, the volume integrated spectra of both electrons and radiation evolve differently from those of planar shocks and exhibit some nonlinear signatures, depending on the time-dependence of the shock parameters. For instance, the shock compression ratio σ and the injection flux of CR electrons decrease, as the shock expands and slows down, resulting in some curvatures in both electron and radiation spectra. Magnetic fields play key roles in DSA at collisionless shocks and control the synchrotron cooling and emission by relativistic electrons. The observed magnetic field strength is found to decrease from ∼ 1 − 10 µ G in the core region to ∼ 0.1 − 1 µ G in the periphery of clusters ( Feretti et al. 2012 ). On the other hand, it is well established that magnetic fields can be ampli- fied via resonant and non-resonant instabilities induced by CR protons streaming upstream of strong shocks ( Bell 1978 ; Lucek & Bell 2000 ; Bell 2004 ). Recently Caprioli & Spitkovsky (2014) have shown that the magnetic field amplification (MFA) factor scales with the Alfvénic Mach number, M_A , and the CR proton acceleration efficiency as < δB / B > ² ∼ 3 M_A Here δB is the turbulent magnetic fields perpendicular to the mean background magnetic fields, ρ ₁ is the upstream gas density, and P _cr,2 is the downstream CR pressure. For typical cluster shocks with 2 M_s 5 and 10 M_A 25 ( Ryu et al. 2003 ), the MFA factor due to the streaming stabilities is expected be rather small but not negligible, < δB / B > ² ∼ 0.3 − 3. However, it has not yet been fully understood how magnetic fields may be amplified both upstream and downstream of a weak shock in high beta ICM plasmas with β_p = P_g / P_B ∼ 100. Therefore, in Paper I, we considered several models with decaying postshock magnetic fields, in which the downstream magnetic field, B_d ( r ), decreases behind the shock with a scale height of 100-150 kpc. Varying B_d ( r ) profiles were found to have little impact on the spatial distribution of the electron energy spectrum, because the baseline cooling rate of electrons is set by iC scattering off the cosmic background photons. Further, electrons in a broad range of γ_e contribute to j_ν ( r ), which smooths out any spatial variations. We note, however, that the B ² dependence of the synchrotron emissivity ( j_ν ∝ N_eB ² ) can become significant in some cases. Moreover, any nonlinear features due to the spatial variations of N_e ( r , γ_e ) and B_d ( r ) are mostly averaged out, leaving only subtle signatures in the volume integrated spectrum, J_ν . For the case with a constant background density (e.g., MF1-3 model in Paper 1), the shock speed decreases approximately as u_s ∝ t ^−3/5 . In fact, it has not been examined, through cosmological hydrodynamic simulations, whether these shock bubbles would accelerate or decelerate as they propagate through the cluster outskirts. In this study, we have performed additional DSA simulations, in which the initial Sedov blast wave travels through a background medium with ρ_u ∝ r ⁻ⁿ , where r is the radial distance from the cluster center and n = 2 − 4. This effectively mimics a blast wave that expands into a constant-density core surrounded by an isothermal halo with a decreasing pressure pro- file. In these new runs, the spherical shock decelerates much more slowly than u_s ∝ t ^−3/5 , and so the nonlinear effects due to the deceleration of the shock speed are expected to be reduced compared to the uniform density models considered in Paper I. We have also calculated the projected surface brightness, I_ν , which depends on the three dimensional structure of the shock surface and the viewing direction. Because of the curvature in the structure of the model shock, synchrotron emissions from downstream electrons with different ages contribute to the surface brightness along a given line-of-sight (i.e., projection effects). So the observed spatial profile of I_ν ( R ) is calculated by assuming the geometrical configuration described in Figure 1 . We note that this model with a ribbon-like shock surface gave rise to radio flux pro- files, S_ν ( R ), that were consistent with those of the the Sausage relic and the double relics in ZwCl0008.8+5215 (e.g., van Weeren et al. 2010 , 2011 ; Kang et al. 2012 ). In paper I, we demonstrated that the spectral index of the volume-integrated spectral index increases gradually from A_ν = α _inj to A_ν = α _inj + 0.5 over a broad frequency range, ∼ (0.1 − 10) ν _br , where the break frequency is ν _br ∼ 0.4 GHz at the shock age of about 50 Myr for postshock magnetic fields B ₂ ∼ 0.7 or 7 µ G (see Equation 9). Here, we explore whether such a transition can explain the broken power-law spectra observed in the radio relic in A2256 ( Trasatti et al. 2015 ). In the next section we describe the numerical calculations. The DSA simulation results of blast wave models with different postshock magnetic field profiles and with different background density profiles will be discussed in Section 3. A brief summary will be given in Section 4. In order to calculate DSA of CR electrons at spherical shocks, we have solved the time-dependent diffusion-convection equation for the pitch-angle-averaged phase space distribution function for CR electrons, f_e ( r, p, t ) = g_e ( r, p, t ) p ⁻⁴ , in the one-dimensional (1D) spherically symmetric geometry: where u ( r, t ) is the flow velocity, y = ln( p/m_ec ), m_e is the electron mass, c is the speed of light, and D ( r, p ) is the spatial diffusion coefficient ( Skilling 1975 ). We adopt a Bohm-like diffusion coefficient with a weaker non-relativistic momentum dependence The cooling term b ( p ) = − dp/dt accounts for electron synchrotron and iC losses, and the cooling time scale for electrons is given as (e.g., Kang 2011 ). The ‘effective’ magnetic field strength takes account for radiative losses due to both synchrotron and iC processes, where B _rad = 3.24 µ G(1+ z ) ² corresponds to the cosmic background radiation at redshift z ( Schlickeiser 2002 ). In this study, we set z = 0.2 as a reference epoch and so B _rad = 4.67 µ G. Here we adopt a simple phenomenological injection model, in which the particles with p p _inj ≈ 5.34 m_pu_s ( t )/ σ are allowed to cross the shock front and get injected into the CR populations ( Kang et al. 2002 ). For spherically expanding shocks considered here, the injection momentum decreases in time as the shock decelerates. Since we are concerned with synchrotron emission from electrons alone, so long as we do not care about the exact units of flux, we can set the electron to proton ratio K_e/p = 1 in the simulations. More detailed discussion on the pre-acceleration and injection of electrons can be found in Paper I. At weak shocks in the test-particle limit, the CR feedback becomes negligible and the background flow, u ( r, t ), is governed by the standard gasdynamic conservation equations in 1D spherical coordinates ( Kang & Jones 2006 ). The test-particle version of CRASH (Cosmic-Ray Amr SHock) code in a comoving spherical grid was used. For the initial shock structure, we adopt a Sedov-Taylor similarity solution propagating in a uniform ICM with the following parameters: the ICM density, n _H,1 = 10 ⁻³ cm ⁻³ , the ICM temperature, T ₁ = 5 × 10 ⁷ K, the initial shock radius, r_s,i = 0.78Mpc, and the initial shock speed, u_s,i = 4.5 × 10 ³ km s ⁻¹ with the sonic Mach number, M_s,i = 4.3 at the onset of the simulations. The shock parameters change in time as the spherical shock expands out, depending on the upstream conditions in the cluster outskirts. The synchrotron emissivity, j_ν ( r ), at each shell is calculated (in units of erg cm ⁻³ s ⁻¹ Hz ⁻¹ str ⁻¹ ), using the electron distribution function, f_e ( r, p, t ), and the magnetic field profile, B ( r, t ), from the DSA simulations. Then the radio intensity or surface brightness, I_ν (erg cm ⁻² s ⁻¹ Hz ⁻¹ str ⁻¹ ), is calculated by integrating j_ν along the path length, h , as shown in Figure 1 : Here R is the distance behind the projected shock edge in the plane of the sky and d = r_s − r is the distance of a shell behind the shock, where ( r_s − d ) ² = ( r_s − R ) ² + h ² . The extension angle ψ = 10 ^◦ is assumed in this study. The volume integrated emissivity, J_ν = ∫ j_ν ( r ) dV , is calculated by integrating j_ν over the downstream volume specified with r_s and ψ , as shown in Figure 1 . The spectral indices of j_ν ( d ), J_ν , and I_ν ( R ) are defined as follows: estimated between ν_i and ν _i+1 . We chose the following four frequencies at the source, ν ₁ = 240 MHz, ν ₂ = 600 MHz, ν ₃ = 1.4 GHz, and ν ₄ = 3.0 GHz in Figures 2 and 4 . Then the redshifted frequency for objects at a redshift z is ν _obs , i = ν_i /(1 + z ). As in the previous study of Kang et al. (2012) , we adopt the postshock magnetic field strength B ₂ ∼ 7 µ G in order to model the observed width of the Sausage relic (∆ l ∼ 55 kpc) as the cooling length of γ_e ∼ 10 ⁴ electrons (∼ u ₂ t _rad ), while the preshock magnetic field strength is chosen to be B ₁ ∼ 2 − 3 µ G. Since B ₁ cannot be constrained directly from observations, it is adjusted so that B ₂ becomes about 7 µ G after considering MFA or compression of the perpendicular components of magnetic fields across the shock. However, we will show in the next section that the observed width can be modeled with a weaker magnetic field of B ₂ ∼ 0.7 µ G (see the discussion below Equation 9). We consider several models whose characteristics are summarized in Table 1 . As described in Paper I, we adopt the following downstream magnetic field profile, B_d ( r ), in MF models: For these models, the upstream density and temperature are assumed to be uniform. In Figure 2 , we compare these models at the shock age of 47 Myr. The evolution of the shock, i.e., r_s ( t ) ∝ t ^2/5 and u_s ( t ) ∝ t ^−3/5 , is identical in these three models. Synchrotron emissivity scales with the electron energy spectrum and the magnetic field strength as j_ν ∝ N_e ( γ_e ) B ² . Of course, the evolution of N_e ( γ_e ) in each MF model also depends on the assumed pro- file of B_d ( r ) through DSA and synchrotron cooling. In MF1 model where B_d ( r ) = B ₂ , the effects of shock deceleration are dramatically revealed in the spatial distribution of j_ν ( d ) (at the low frequency of 240MHz). It increases downstream because of higher shock compression and higher electron injection flux at the earlier epoch. On the other hand, at the higher frequency of 1.4 GHz, j_ν ( d ) decreases behind the shock due to the fast synchrotron/iC cooling of high energy electrons. The downstream increase of j_ν ( d ) at low frequencies is softened in the models with decaying postshock magnetic fields, compared to MF1 model, because of B ² dependence of the synchrotron emissivity. In MF3 model, in which the downstream magnetic field decreases with the gas density behind the shock, the spatial distributions of j_ν ( d ) at all three frequencies decrease downstream. Thus the downstream distribution of j_ν ( d ) depends on the shock speed evolution, u_s ( t ), and the postshock B_d ( r ) as well as the chosen frequency. In Paper I, we indicated that signatures imprinted on synchrotron emission, j_ν ( d ), and its volume integrated spectrum, J_ν , due to different postshock magnetic field profiles could be rather subtle to detect. That expectation arose because we examined the behavior of log j_ν ( r ) in paper I (see Figure 6 there), while we plot νj_ν ( d ) in Figure 2 . As can be seen in the lower two rows in Figure 2 , the surface brightness profile, I_ν ( R ), is affected by projection effects as well as the evolution of u_s ( t ) and the spatial variation of B_d ( r ). For instance, the gradual increase of I_ν just behind the shock up to R _inf ≈ 15 kpc is due to the increase of the path length, and its inflection point, R _inf = r_s (1 − cos ψ ), depends on the value of ψ . Beyond the inflection point, the path length decreases but j_ν ( d ) may increase or decrease depending on u_s ( t ) and B_d ( r ), resulting in a range of spatial profiles of I_ν ( R ). Figure 2 also demonstrates that the spectral indices, α _νi−νi+1 ( d ) and B _νi−νi+1 ( R ), at all three frequencies decrease behind the shock and do not show significant variations among the different MF models other than faster steepening of both indices for more rapidly decaying B_d ( r ) profiles. In summary, at Sedov-Taylor type spherical shocks decelerating with u_s ∝ t ^−3/5 , the energy spectrum, N_e ( r , γ_e ), of low energy electrons increases downstream behind the shock. As a result, the spatial distributions of the radio emissivity, j_ν ( d ), and the surface brightness, I_ν ( R ), at low radio frequencies (< 1 GHz) could depend significantly on the postshock magnetic field profile. At high radio frequencies, such dependence becomes relatively weaker, because the width of the postshock spatial distribution of high energy electrons is much narrower. So, the magnetic field profile far downstream has less influence on synchrotron emission. On the other hand, the spectral indices, α_ν and B_ν , are relatively insensitive to these variations. In BD models, we assume that the initial blast wave propagates into an isothermal halo with a different density profile for r > r_s,i : The upstream temperature is T ₁ = 5×10 ⁷ K, and B ₁ = 2 µ G and B ₂ = 7 µ G for these models. In fact MF1 and BD1 models are identical. We also consider BD2b and BD3b models, in which the downstream magnetic field profile is the same as MF2 , i.e., In the so-call beta model for isothermal ICMs, the gas distribution can be represented by ρ ( r ) ∝ r ^−3β in the outskirts of galaxy clusters ( Sarazin 1988 ). So BD2 model corresponds to the beta model with β ≈ 2/3, which is consistent with typical X-ray brightness profile of observed X-ray clusters. Recall that in our simulations the spherical blast wave propagating through a uniform ICM is adopted for the initial conditions. So we are effectively considering a spherical blast wave that propagates first in a constant-density core (i.e., ρ = ρ ₀ for r < r_s,i ) and then expands into a surrounding isothermal halo with ρ ∝ r ⁻ⁿ (for r > r_s,i ). With the different background density (or gas pressure) profile, the shock speed evolves differently as it expands outward. Figure 3 shows how the shock radius, r_s ( t ), shock speed, u_s ( t ), the sonic Mach number, M_s ( t ), and the DSA spectral index at the shock position, α _inj , vary in time in BD models. As expected, the shock decelerates much more slowly if the background pressure decreases outward as in BD2 and BD3 models. In BD3 model, for instance, the shock speed decreases less than 10% during 100 Myr. Hence these BD models allow us to explore the dependence of radio spectral properties for a range of the time evolution of the shock. Figure 4 compares the surface brightness profile, I_ν ( R ), at ν ₁ = 240 MHz (black solid lines), ν ₂ = 600 MHz (red dotted), and ν ₃ = 1.4 GHz (blue dashed) in models with different ρ_u and B_d . In BD1 model, the shock parameters are: r_s = 0.83Mpc, u_s = 4.1 × 10 ³ km s ⁻¹ and M_s = 3.9 at t ₁ = 13 Myr; r_s = 0.91Mpc, u_s = 3.6 × 10 ³ km s ⁻¹ and M_s = 3.4 at t ₂ = 33 Myr; r_s = 1.0Mpc, u_s = 3.1 × 10 ³ km s ⁻¹ and M_s = 2.9 at t ₃ = 67 Myr. The width of the radio structure at low frequencies, ∆ l _adv ∼ t _age u_s / σ , increases with the shock age, while the amplitude of the surface brightness decreases in time (from left to right in Figure 4 ). For high frequencies (> 1 GHz), however, the width asymptotes to the cooling length, ∆ l _cool ∼ t _rad ( γ_e ) u_s / σ . So the postshock magnetic field strength could be inferred from t _rad ( γ_e ) with high frequency observations. Since the synchrotron emission from monoenergetic electrons with γ_e in the postshock magnetic field of B ₂ peaks at ν _peak ≈ 0.3(3 eB ₂ /4 πm_ec ) ≈ 0.88 GHz( B ₂ /7 µ G)( γ_e /10 ⁴ ) ² , according to Equation (3), the cooling time scale of the electrons that emit mostly at ν _peak becomes where Here note that the quantity, ( B _e,2 /8.4 µ G) ⁻² ( B ₂ /7 µ G) ^1/2 , has an identical value for two values of B ₂ , e.g., 0.7 and 7 µ G. So the downstream cooling length of the relic shocks becomes degenerate for the two possible values of B ₂ . At low frequencies (< 1 GHz) the surface brightness is affected by the evolution of u_s ( t ) and the spatial variation of B_d ( r ) as well as the projection effects, resulting in a wide range of spatial profiles of I_ν ( R ). This suggests that, if the radio surface brightness can be spatially resolved at several radio frequencies over ∼ (0.1 − 10) GHz, we may extract the time evolution of u_s ( t) as well as the shock age from low frequency observations. So it would be useful to compare multifrequency radio observations with the intensity modeling that accounts for the shock evolution, magnetic field profile and projection effects. The comparison of I_ν ( R ) of MF1 and MF2 models in Figure 2 demonstrates that the decaying postshock magnetic field could lessen the shock deceleration signatures in the surface brightness profile. If we compare I_ν ( R ) of BD2 and BD2b models in Figure 4 , on the other hand, the difference in their profiles is much smaller than that between MF1 and MF2 models. This is because the the shock deceleration effects and the ensuing downstream increase of N_e ( r , γ_e ) are relatively milder in BD2 and BD2b models. In Figure 5 , we compare the volume integrated spectrum, J_ν , and its spectral index, A_ν , at three different shock ages, t _age = 13 Myr (black solid), 33 Myr (red dotted), and 67 Myr (blue dashed) for the same set of models shown in Figure 4 . For the test-particle power-law at steady planar shocks, the radio index A_ν is the same as α _inj for at the source. For ν > ν _br , A_ν is expected to steepen to α _inj + 0.5 due to synchrotron/iC cooling of electrons. Again note that the quantity, ( B _e,2 /8.4 µ G) ⁻⁴ ( B ₂ /7 µ G) can have the same value for two possible value of B ₂ , and so the break frequency at a given shock age becomes identical as well. The observed (redshifted) break frequency corresponds to ν _br,obs = ν _br /(1 + z ). If the shock age is t _age ∼ 67 Myr and B ₂ ∼ 0.7 or 7 µ G, for example, the observed break frequency becomes about 120 MHz for objects at z = 0.2. The smallest possible index is α _inj = 0.5 at strong shocks, so A_ν = α _inj + 0.5 ≥ 1.0 above the break frequency. As can be seen in Figure 5 , the transition from α _inj to α _inj + 0.5 occurs gradually over about two orders of magnitude in frequency range. If the shock is 30-70 Myr old for the shock parameters considered here, the volume-integrated radio spectrum is expected to steepen gradually from 100 MHz to 10 GHz, instead of a sharp broken power-law with the break frequency at ν _br,obs . Such gradual steepening may explain why the volume-integrated radio spectra of some observed relic shocks might be interpreted as a broken power-law with A_ν < 1.0 at low frequencies. For example, according to a recent observation of the relic in A2256, the observed spectral index is between 351 and 1369 MHz and increases to between 1369 and 10450 MHz ( Trasatti et al. 2015 ). Similar indices were obtained for our J_ν between 355 and 1413 MHz, and between 1413 MHz and 10 GHz, using the DSA simulation results. The best fits are obtained for BD2b and BD3b models. In BD2b model, at t _age = 33 Myr (red dotted line in Figure 5 ), these spectral indices are and In BD3b model, again at t _age = 33 Myr, they are and Thus we may explain the curvature in the integrated spectrum of the observed radio relics around 1GHz by the electron population radiatively cooling behind a cluster shock, if we assume that the age of the shock (or the electron acceleration duration) is about 30 Myr, which is relatively young compared to the dynamical time scales of typical clusters. In the case of BD3b model , the spectral index behaves very similarly to that of a plane shock case (see Figure 4 of Paper 1), since the shock speed is more or less constant in time. Departures from the predictions for the test-particle planar shock are the most severe in BD1 ( MF1 ) model, while it becomes relatively milder in BD2 and BD3 models with decreasing halo density profiles. As shown in Paper I, any variations in the spatial distributions of f_e ( r ) and B_d ( r ) are averaged in the volume integrated quantities such as J_ν . So signatures imprinted on the volume-integrated emission due to different B_d ( r ) (e.g., between BD2 and BD2b models) would be too subtle to detect. We have performed time-dependent DSA simulations for cosmic-ray (CR) electrons at decelerating spherical shocks with parameters relevant for weak cluster shocks: u_s ≈ (3.0 − 4.5) × 10 ³ km s ⁻¹ and M_s ≈ 3.0−4.3. Several models with different postshock magnetic field profiles ( MF1-3 ) and different upstream gas density profiles ( BD1-3 ) were considered as summarized in Table 1 . Using the synchrotron emissivity, j_ν ( r ), calculated from the CR electron energy spectra at these model shocks, the radio surface brightness pro- file, I_ν ( R ), and the volume integrated spectrum, J_ν , were estimated by assuming a ribbon-like shock structure described in Figure 1 . At low frequencies (< 1 GHz) the surface brightness is affected by the evolution of u_s ( t ) and the spatial variation of B_d ( r ) as well as projection effects, resulting in a wide range of spatial profiles of I_ν ( R ) (see Figures 2 and 4 ). At high frequencies (> 1 GHz), for a given geometrical structure of the shock, such dependences become relatively weaker, because the width of the postshock spatial distribution of high energy electrons is much narrower and so the magnetic field profile far downstream has less influence on synchrotron emission. For low frequency observations, the width of radio relics increases with the shock age as ∆ l _adv ∼ t _age u_s / σ , while it asymptotes to the cooling length, ∆ l _cool ∼ t _rad ( γ_e ) u_s / σ , for high frequencies. If the surface brightness can be spatially resolved at multi-frequency observations over ∼ (0.1 − 10) GHz, we may extract significant information about the time evolution of u_s ( t ), the shock age, t _age , and the postshock magnetic field strength, B_d ( r ), through the detail modeling of DSA and projection effects. For instance, one may infer two possible values of the postshock magnetic field strength from the observed width of radio relics, provided that the geometrical structure (e.g., ψ ) of the radio structure is known. The spectral index of I_ν ( R ), however, behaves rather similarly in all the models considered here. If the postshock magnetic field strength is about 0.7 or 7 µ G, at the shock age of ∼ 30 Myr, the volume integrated radio spectrum has a break frequency, ν _br ∼ 1 GHz, and steepens gradually with the spectral index from α _inj to α _inj + 0.5 over the frequency range of 0.1- 10 GHz (see Figure 5 ). Thus, we suggest that such a curved spectrum could explain the observed spectrum of the relic in cluster A2256 ( Trasatti et al. 2015 ).